Changeset 29 for branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_th_diff.f

Timestamp:

10/12/16 18:08:11 (8 years ago)

Author:

vancop

Message:

Add ice_sal_adv routine for GN13

File:

• branches/2016/dev_v3.20_2016_platelet/SOURCES/source_3.20/ice_th_diff.f (modified) (2 diffs)

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

@@ -722 +722 @@
       ! ice energy of melting
       CALL ice_th_enmelt(kideb, kiut, nlay_s, nlay_i) 
-      
+
+      !----------------------------------------------------------
+      ! Interfacial temperatures (for salinity advection scheme)
+      !----------------------------------------------------------
+
+      ! Uppermost ice temperature
+      IF ( ht_s_b(ji) .LE. 0. ) THEN 
+         t_i_int(ji,1) = t_su_b(ji)
+      ELSE
+         zk            = xkg + betak1*s_i_b(ji,1) 
+     &                 / MIN( -zeps , t_i_b(ji,1) - tpw ) 
+     &                 - betak2* ( t_i_b(ji,1) - tpw )
+         zk            = MAX( zk, zkimin , zkimin )
+         zdummy1       = zk      * ht_s_b(ji) * t_s_b(ji,1) + 
+     &                   ykn(ji) * deltaz_i_phy(1) * t_i_b(ji,1)
+         zdummy2       = zk      * ht_s_b(ji) +
+     &                   ykn(ji) * deltaz_i_phy(1)
+         t_i_int(ji,1) = zdummy1 / zdummy2
+      ENDIF
+
+      ! Inner ice temperatures
+      DO layer = 2, nlay_i
+         t_i_int(ji,layer) =   
+     &             t_i_b(ji,layer-1) + deltaz_i_phy(layer-1) / 2. * 
+     &             ( t_i_b(ji,layer) - t_i_b(ji,layer-1) ) /
+     &             ( z_i_phy(layer) - z_i_phy(layer-1) ) 
+      END DO
+
+      ! Basal ice temperature
+      t_i_int(ji,nlay_i+1) = t_bo_b(ji)
+
       IF ( ln_write ) THEN
          WRITE(numout,*) ' nconv : ', nconv
@@ -729 +759 @@
          WRITE(numout,*) ' t_s_b : ', ( t_s_b(ji,layer), 
      &                      layer = 1, nlay_s )
+         WRITE(numout,*) ' t_i_int : ', ( t_i_int(ji,layer),
+     &                      layer = 1, nlay_i+1 )
          WRITE(numout,*) ' t_i_b : ', ( t_i_b(ji,layer), 
      &                      layer = 1, nlay_i )
          WRITE(numout,*) ' t_bo_b : ', t_bo_b(ji)
+      
          WRITE(numout,*)
       ENDIF