Changeset 9416 for branches/2017/dev_merge_2017/NEMOGCM

Timestamp:

2018-03-22T12:42:04+01:00 (6 years ago)

Author:

clem

Message:

debug icetsni output by initializing t_si in iceistate.F90

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90 (modified) (3 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90

@@ -113 +113 @@
       ! init surface temperature
       DO jl = 1, jpl
-         t_su  (:,:,jl) = rt0 * tmask(:,:,1)
-         tn_ice(:,:,jl) = rt0 * tmask(:,:,1)
+         t_su   (:,:,jl) = rt0 * tmask(:,:,1)  ! temp at the surface
+         t_si   (:,:,jl) = rt0 * tmask(:,:,1)  ! temp at the ice-snow interface
+         cnd_ice(:,:,jl) = 0._wp               ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling)
       END DO
       !
@@ -355 +356 @@
          tn_ice (:,:,:) = t_su (:,:,:)
          t1_ice (:,:,:) = t_i (:,:,1,:)   ! initialisation of 1st layer temp for coupled simu
-         cnd_ice(:,:,:) = 0._wp           ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling)
 
          ! Melt pond volume and fraction
@@ -393 +393 @@
          tn_ice (:,:,:) = t_i (:,:,1,:)
          t1_ice (:,:,:) = t_i (:,:,1,:)   ! initialisation of 1st layer temp for coupled simu
-         cnd_ice(:,:,:) = 0._wp           ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling)
          
          a_ip(:,:,:)      = 0._wp