Changeset 8933 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icerst.F90

Timestamp:

2017-12-07T12:14:23+01:00 (6 years ago)

Author:

clem

Message:

dev_CNRS_2017: raw commit of the Met-Office specifics for coupling with their Earth system model. It is not yet tested and debuging is expected

File:

• branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icerst.F90 (modified) (6 diffs)

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

@@ -19 +19 @@
    USE ice            ! sea-ice variables
    USE dom_oce        ! ocean domain
-   USE sbc_oce , ONLY : nn_fsbc
+   USE sbc_oce , ONLY : nn_fsbc, ln_meto_cpl
    USE icectl
    !
@@ -129 +129 @@
       CALL iom_rstput( iter, nitrst, numriw, 'a_i' , a_i  )
       CALL iom_rstput( iter, nitrst, numriw, 't_su', t_su )
-      !
       ! Melt ponds
       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, 'tempt_sl1' , z3d )
-      !
       ! Ice enthalpy (all ice layers)
       DO jk = 1, nlay_i 
@@ -146 +143 @@
          CALL iom_rstput( iter, nitrst, numriw, znam , z3d )
       END DO
-      !
       ! ice velocity
       CALL iom_rstput( iter, nitrst, numriw, 'u_ice', u_ice ) ! u_ice
       CALL iom_rstput( iter, nitrst, numriw, 'v_ice', v_ice ) ! v_ice
-      
+      ! fields needed for Jules coupling
+      IF (ln_meto_cpl) THEN
+         CALL iom_rstput( iter, nitrst, numriw, 'cnd_ice', cnd_ice )
+         CALL iom_rstput( iter, nitrst, numriw, 't1_ice' , t1_ice  )
+      ENDIF
+      !
+
       ! close restart file
       ! ------------------
@@ -207 +209 @@
       CALL iom_get( numrir, jpdom_autoglo, 'a_i' , a_i  )
       CALL iom_get( numrir, jpdom_autoglo, 't_su', t_su )
-      !
       ! Melt ponds
       id1 = iom_varid( numrir, 'a_ip_htc01' , ldstop = .FALSE. )
@@ -218 +219 @@
          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, 'tempt_sl1' , z3d )
       e_s(:,:,1,:) = z3d
-      !
       ! Ice enthalpy (all ice layers)
       DO jk = 1, nlay_i 
@@ -231 +230 @@
          e_i(:,:,jk,:) = z3d(:,:,:)
       END DO
-      !
       ! ice velocity
       CALL iom_get( numrir, jpdom_autoglo, 'u_ice', u_ice )
       CALL iom_get( numrir, jpdom_autoglo, 'v_ice', v_ice )
+      ! fields needed for Jules coupling
+      IF (ln_meto_cpl) THEN
+         CALL iom_get( numrir, jpdom_autoglo, 'cnd_ice', cnd_ice )
+         CALL iom_get( numrir, jpdom_autoglo, 't1_ice' , t1_ice  )
+      ENDIF
 
    END SUBROUTINE ice_rst_read