Changeset 11837 for NEMO/branches/2019/dev_r11756_SI3restart_XIOS/src/ICE/icerst.F90

Timestamp:

2019-10-30T15:10:36+01:00 (4 years ago)

Author:

andmirek

Message:

ticket #2323 read SI3 restart with XIOS

File:

• NEMO/branches/2019/dev_r11756_SI3restart_XIOS/src/ICE/icerst.F90 (modified) (6 diffs)

NEMO/branches/2019/dev_r11756_SI3restart_XIOS/src/ICE/icerst.F90

@@ -184 +184 @@
       ENDIF
 
+      lxios_sini = .FALSE.
       CALL iom_open ( TRIM(cn_icerst_indir)//'/'//cn_icerst_in, numrir, kdlev = jpl )
 
@@ -192 +193 @@
       IF( id0 > 0 ) THEN   ! == case of a normal restart == !
          !                 ! ------------------------------ !
-         
+         IF(.NOT.lxios_set) THEN
+            lrxios = lrxios.AND.lxios_sini
+         ELSE
+            IF(.NOT.lxios_sini) THEN
+               WRITE(ctmp1,*) 'E R R O R : XIOS can not be used to read SI3 restart from multiple files'
+               CALL ctl_stop( 'ice_rst_read:', ctmp1 )
+            ENDIF
+         ENDIF
+
+         lrsi3_nemo = .TRUE.                       ! SI3 restart written by NEMO
+
+         IF(lrxios) THEN
+           crixios_context = 'si3_rst'
+           CALL iom_init( crixios_context, TRIM(cn_icerst_indir)//'/'//cn_icerst_in, ld_tmppatch = .false. )
+           lxios_set = .TRUE.
+         ENDIF
+
          ! Time info
-         CALL iom_get( numrir, 'nn_fsbc', zfice )
-         CALL iom_get( numrir, 'kt_ice' , ziter )    
+         CALL iom_get( numrir, 'nn_fsbc', zfice, ldxios = lrxios )
+         CALL iom_get( numrir, 'kt_ice' , ziter, ldxios = lrxios )    
          IF(lwp) WRITE(numout,*) '   read ice restart file at time step    : ', ziter
          IF(lwp) WRITE(numout,*) '   in any case we force it to nit000 - 1 : ', nit000 - 1
@@ -210 +227 @@
 
          ! --- mandatory fields --- ! 
-         CALL iom_get( numrir, jpdom_autoglo, 'v_i'  , v_i   )
-         CALL iom_get( numrir, jpdom_autoglo, 'v_s'  , v_s   )
-         CALL iom_get( numrir, jpdom_autoglo, 'sv_i' , sv_i  )
-         CALL iom_get( numrir, jpdom_autoglo, 'a_i'  , a_i   )
-         CALL iom_get( numrir, jpdom_autoglo, 't_su' , t_su  )
-         CALL iom_get( numrir, jpdom_autoglo, 'u_ice', u_ice )
-         CALL iom_get( numrir, jpdom_autoglo, 'v_ice', v_ice )
+         CALL iom_get( numrir, jpdom_autoglo, 'v_i'  , v_i, ldxios = lrxios   )
+         CALL iom_get( numrir, jpdom_autoglo, 'v_s'  , v_s, ldxios = lrxios   )
+         CALL iom_get( numrir, jpdom_autoglo, 'sv_i' , sv_i, ldxios = lrxios  )
+         CALL iom_get( numrir, jpdom_autoglo, 'a_i'  , a_i, ldxios = lrxios   )
+         CALL iom_get( numrir, jpdom_autoglo, 't_su' , t_su, ldxios = lrxios  )
+         CALL iom_get( numrir, jpdom_autoglo, 'u_ice', u_ice, ldxios = lrxios )
+         CALL iom_get( numrir, jpdom_autoglo, 'v_ice', v_ice, ldxios = lrxios )
          ! 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
-            CALL iom_get( numrir, jpdom_autoglo, znam , z3d )
-            e_i(:,:,jk,:) = z3d(:,:,:)
-         END DO
+         IF(lrsi3_nemo) THEN
+!NEMO read inerface
+            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
+               CALL iom_get( numrir, jpdom_autoglo, znam , z3d )
+               e_i(:,:,jk,:) = z3d(:,:,:)
+            END DO
+         ELSE
+!XIOS read interface
+            CALL iom_get( numrir, jpdom_autoglo, "e_s" , e_s )
+            CALL iom_get( numrir, jpdom_autoglo, "e_i" , e_i ) 
+         ENDIF
          ! -- optional fields -- !
          ! ice age
          id1 = iom_varid( numrir, 'oa_i' , ldstop = .FALSE. )
          IF( id1 > 0 ) THEN                       ! fields exist
-            CALL iom_get( numrir, jpdom_autoglo, 'oa_i', oa_i )
+            CALL iom_get( numrir, jpdom_autoglo, 'oa_i', oa_i, ldxios = lrxios )
          ELSE                                     ! start from rest
             IF(lwp) WRITE(numout,*) '   ==>>   previous run without ice age output then set it to zero'
@@ -243 +267 @@
          id2 = iom_varid( numrir, 'a_ip' , ldstop = .FALSE. )
          IF( id2 > 0 ) THEN                       ! fields exist
-            CALL iom_get( numrir, jpdom_autoglo, 'a_ip' , a_ip )
-            CALL iom_get( numrir, jpdom_autoglo, 'v_ip' , v_ip )
+            CALL iom_get( numrir, jpdom_autoglo, 'a_ip' , a_ip, ldxios = lrxios )
+            CALL iom_get( numrir, jpdom_autoglo, 'v_ip' , v_ip, ldxios = lrxios )
          ELSE                                     ! start from rest
             IF(lwp) WRITE(numout,*) '   ==>>   previous run without melt ponds output then set it to zero'
@@ -255 +279 @@
             id4 = iom_varid( numrir, 't1_ice'  , ldstop = .FALSE. )
             IF( id3 > 0 .AND. id4 > 0 ) THEN         ! fields exist
-               CALL iom_get( numrir, jpdom_autoglo, 'cnd_ice', cnd_ice )
-               CALL iom_get( numrir, jpdom_autoglo, 't1_ice' , t1_ice  )
+               CALL iom_get( numrir, jpdom_autoglo, 'cnd_ice', cnd_ice, ldxios = lrxios )
+               CALL iom_get( numrir, jpdom_autoglo, 't1_ice' , t1_ice, ldxios = lrxios  )
             ELSE                                     ! start from rest
                IF(lwp) WRITE(numout,*) '   ==>>   previous run without conductivity output then set it to zero'
@@ -263 +287 @@
             ENDIF
          ENDIF
-
+         IF(lrxios) CALL iom_swap( TRIM(cxios_context) )
          CALL iom_delay_rst( 'READ', 'ICE', numrir )   ! read only ice delayed global communication variables