Changeset 9271 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icerst.F90

Timestamp:

2018-01-19T18:56:15+01:00 (6 years ago)

Author:

clem

Message:

first steps for having more than 1 snow layers in the ice (in theory). There is still icethd_dh.F90 routine to change

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icerst.F90 (modified) (2 diffs)

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

@@ -132 +132 @@
       CALL iom_rstput( iter, nitrst, numriw, 'a_ip', a_ip )
       CALL iom_rstput( iter, nitrst, numriw, 'v_ip', v_ip )
-!!gm dangerous !!!!!  ===>>>> better reading writing all snow layers !
-      ! Snow enthalpy (1st snow layer only)
-      z3d = e_s(:,:,1,:)
-      CALL iom_rstput( iter, nitrst, numriw, 'e_s_l01' , z3d )
-      ! Ice enthalpy (all ice layers)
+      ! Snow enthalpy
+      DO jk = 1, nlay_s 
+         WRITE(zchar1,'(I2.2)') jk
+         znam = 'e_s'//'_l'//zchar1
+         z3d(:,:,:) = e_s(:,:,jk,:)
+         CALL iom_rstput( iter, nitrst, numriw, znam , z3d )
+      END DO
+      ! Ice enthalpy
       DO jk = 1, nlay_i 
          WRITE(zchar1,'(I2.2)') jk
@@ -219 +222 @@
          v_ip(:,:,:) = 0._wp
       ENDIF
-!!gm dangerous !!!!!  ===>>>> better reading writing all snow layers !
-      ! Snow enthalpy (1st snow layer only)
-      CALL iom_get( numrir, jpdom_autoglo, 'e_s_l01' , z3d )
-      e_s(:,:,1,:) = z3d
-      ! Ice enthalpy (all ice layers)
-      DO jk = 1, nlay_i 
+      ! Snow enthalpy
+      DO jk = 1, nlay_s
+         WRITE(zchar1,'(I2.2)') jk
+         znam = 'e_s'//'_l'//zchar1
+         CALL iom_get( numrir, jpdom_autoglo, znam , z3d )
+         e_s(:,:,jk,:) = z3d(:,:,:)
+      END DO
+      ! Ice enthalpy
+      DO jk = 1, nlay_i
          WRITE(zchar1,'(I2.2)') jk
          znam = 'e_i'//'_l'//zchar1