Changeset 34 for branches

Timestamp:

10/18/16 18:07:44 (8 years ago)

Author:

vancop

Message:

bug in ice_th_dh.f

File:

• branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_th_dh.f (modified) (1 diff)

branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_th_dh.f

@@ -328 +328 @@
      &                 - cpw*(tmelts-tpw) ) 
 
-         zE1 = - cpw * t_bo_b(ji) ! specific enthalpy of sea water
-         zE2 = q_i_b(ji,nlay_i+1) ! specific enthalpy of new sea ice
-
-         zdE = zE2 - zE1
-
-         dh_i_bott(ji)    =  - ddtb*(fc_bo_i(ji) + fbbqb(ji) )
-     &                    / zdE
+         zE1 = - cpw * ( t_bo_b(ji) - tpw )          ! specific enthalpy of sea water  <0
+         zE2 = - q_i_b(ji,nlay_i+1) / rhog ! specific enthalpy of new sea ice  <0
+
+         WRITE(numout,*) ' zE1 : ', zE1
+         WRITE(numout,*) ' zE2 : ', zE2
+
+         zdE = zE2 - zE1 ! <0
+
+         WRITE(numout,*) ' zdE : ', zdE
+
+         dh_i_bott(ji)    =  ddtb*(fc_bo_i(ji) + fbbqb(ji) )
+     &                    / ( zdE * rhog )
   
          ! salt flux due to initial salt entrapment (keep ?)