Changeset 1565 for trunk/NEMO/LIM_SRC_3/limthd.F90

Timestamp:

2009-07-31T16:01:08+02:00 (15 years ago)

Author:

rblod

Message:

Use appropiate e3 for divergence computation and for sea-ice, see ticket #507

File:

• trunk/NEMO/LIM_SRC_3/limthd.F90 (modified) (2 diffs)

trunk/NEMO/LIM_SRC_3/limthd.F90

@@ -14 +14 @@
    USE phycst          ! physical constants
    USE dom_oce         ! ocean space and time domain variables
+   USE domvvl          ! Variable volume
    USE lbclnk
    USE in_out_manager  ! I/O manager
@@ -241 +242 @@
             ! Energy needed to bring ocean surface layer until its freezing
             ! qcmif, limflx
-            qcmif  (ji,jj) =  rau0 * rcp * fse3t(ji,jj,1) * ( t_bo(ji,jj) - (sst_m(ji,jj) + rt0) ) * ( 1. - zinda )
+            qcmif  (ji,jj) =  rau0 * rcp * fse3t_m(ji,jj,1)   &
+               &           * ( t_bo(ji,jj) - (sst_m(ji,jj) + rt0) ) * ( 1. - zinda )
 
             !  calculate oceanic heat flux (limthd_dh)