Changeset 2678 for branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

Timestamp:

2011-03-09T16:02:46+01:00 (13 years ago)

Author:

rblod

Message:

Phasing branch dev_r2586_dynamic_mem with revision 2675 off the trunk

File:

• branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90 (modified) (2 diffs)

branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

@@ -79 +79 @@
       !!           profile of the ice/snow layers : z_i, z_s
       !!           total ice/snow thickness : ht_i_b, ht_s_b
+      !!
+      !! ** External : 
+      !!
+      !! ** References :
+      !!
+      !! ** History :
+      !!           (02-2003) Martin Vancoppenolle, Louvain-la-Neuve, Belgium
+      !!           (06-2005) Martin Vancoppenolle, 3d version
+      !!           (11-2006) Vectorized by Xavier Fettweis (UCL-ASTR)
+      !!           (04-2007) Energy conservation tested by M. Vancoppenolle
       !!------------------------------------------------------------------
       INTEGER , INTENT (in) ::  &
@@ -333 +343 @@
                END DO
             END DO
+         ENDIF
+
+         IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle (0.011/2=0.0055)
+            DO layer = 1, nlay_i-1
+               DO ji = kideb , kiut
+                  ztcond_i(ji,layer) = rcdic + 0.09*( s_i_b(ji,layer)   &
+                     + s_i_b(ji,layer+1) ) / MIN(-2.0*zeps,      &
+                     t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) - &
+                     0.0055* ( t_i_b(ji,layer) + t_i_b(ji,layer+1) - 2.0*rtt )  
+                  ztcond_i(ji,layer) = MAX(ztcond_i(ji,layer),zkimin)
+               END DO
+            END DO
+         ENDIF
+
+         IF ( thcon_i_swi .EQ. 0 ) THEN ! Untersteiner
             DO ji = kideb , kiut
                ztcond_i(ji,nlay_i)   = rcdic + zbeta*s_i_b(ji,nlay_i) / &