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

Timestamp:

2019-10-31T16:30:50+01:00 (4 years ago)

Author:

andmirek

Message:

Ticket #2323 write 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

@@ -51 +51 @@
       CHARACTER(len=50)   ::   clname   ! ice output restart file name
       CHARACTER(len=256)  ::   clpath   ! full path to ice output restart file 
+      CHARACTER(LEN=52)   ::   clpname   ! ocean output restart file name including prefix for AGRIF
       !!----------------------------------------------------------------------
       !
@@ -80 +81 @@
             ENDIF
             !
-            CALL iom_open( TRIM(clpath)//TRIM(clname), numriw, ldwrt = .TRUE., kdlev = jpl )
+            IF(.NOT.lwxios) THEN
+               CALL iom_open( TRIM(clpath)//TRIM(clname), numriw, ldwrt = .TRUE., kdlev = jpl )
+            ELSE
+#if defined key_iomput
+               cwixios_context = "si3_rstw_"//TRIM(ADJUSTL(clkt))
+               IF( TRIM(Agrif_CFixed()) == '0' ) THEN
+                  clpname = clname
+               ELSE
+                  clpname = TRIM(Agrif_CFixed())//"_"//clname
+               ENDIF
+               CALL iom_init( cwixios_context, TRIM(clpath)//TRIM(clpname), .false. )
+               CALL xios_update_calendar(nitrst)
+               CALL iom_swap(      cxios_context          )
+#else
+               clinfo = 'Can not use XIOS in rst_opn'
+               CALL ctl_stop(TRIM(clinfo))
+#endif
+            ENDIF
             lrst_ice = .TRUE.
          ENDIF
@@ -117 +135 @@
       ! ------------------ 
       !                                                                        ! calendar control
-      CALL iom_rstput( iter, nitrst, numriw, 'nn_fsbc', REAL( nn_fsbc, wp ) )      ! time-step 
-      CALL iom_rstput( iter, nitrst, numriw, 'kt_ice' , REAL( iter   , wp ) )      ! date
+      IF(lwxios) CALL iom_swap(      cwixios_context          )
+      CALL iom_rstput( iter, nitrst, numriw, 'nn_fsbc', REAL( nn_fsbc, wp ), ldxios = lwxios )      ! time-step 
+      CALL iom_rstput( iter, nitrst, numriw, 'kt_ice' , REAL( iter   , wp ), ldxios = lwxios )      ! date
       CALL iom_delay_rst( 'WRITE', 'ICE', numriw )   ! save only ice delayed global communication variables
 
       ! Prognostic variables
-      CALL iom_rstput( iter, nitrst, numriw, 'v_i'  , v_i   )
-      CALL iom_rstput( iter, nitrst, numriw, 'v_s'  , v_s   )
-      CALL iom_rstput( iter, nitrst, numriw, 'sv_i' , sv_i  )
-      CALL iom_rstput( iter, nitrst, numriw, 'a_i'  , a_i   )
-      CALL iom_rstput( iter, nitrst, numriw, 't_su' , t_su  )
-      CALL iom_rstput( iter, nitrst, numriw, 'u_ice', u_ice )
-      CALL iom_rstput( iter, nitrst, numriw, 'v_ice', v_ice )
-      CALL iom_rstput( iter, nitrst, numriw, 'oa_i' , oa_i  )
-      CALL iom_rstput( iter, nitrst, numriw, 'a_ip' , a_ip  )
-      CALL iom_rstput( iter, nitrst, numriw, 'v_ip' , v_ip  )
+      CALL iom_rstput( iter, nitrst, numriw, 'v_i'  , v_i, ldxios = lwxios   )
+      CALL iom_rstput( iter, nitrst, numriw, 'v_s'  , v_s, ldxios = lwxios   )
+      CALL iom_rstput( iter, nitrst, numriw, 'sv_i' , sv_i, ldxios = lwxios  )
+      CALL iom_rstput( iter, nitrst, numriw, 'a_i'  , a_i, ldxios = lwxios   )
+      CALL iom_rstput( iter, nitrst, numriw, 't_su' , t_su, ldxios = lwxios  )
+      CALL iom_rstput( iter, nitrst, numriw, 'u_ice', u_ice, ldxios = lwxios )
+      CALL iom_rstput( iter, nitrst, numriw, 'v_ice', v_ice, ldxios = lwxios )
+      CALL iom_rstput( iter, nitrst, numriw, 'oa_i' , oa_i, ldxios = lwxios  )
+      CALL iom_rstput( iter, nitrst, numriw, 'a_ip' , a_ip, ldxios = lwxios  )
+      CALL iom_rstput( iter, nitrst, numriw, 'v_ip' , v_ip, ldxios = lwxios  )
       ! Snow enthalpy
       DO jk = 1, nlay_s 
@@ -137 +156 @@
          znam = 'e_s'//'_l'//zchar1
          z3d(:,:,:) = e_s(:,:,jk,:)
-         CALL iom_rstput( iter, nitrst, numriw, znam , z3d )
+         CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios )
       END DO
       ! Ice enthalpy
@@ -144 +163 @@
          znam = 'e_i'//'_l'//zchar1
          z3d(:,:,:) = e_i(:,:,jk,:)
-         CALL iom_rstput( iter, nitrst, numriw, znam , z3d )
+         CALL iom_rstput( iter, nitrst, numriw, znam , z3d, ldxios = lwxios )
       END DO
       ! fields needed for Met Office (Jules) coupling
       IF( ln_cpl ) THEN
-         CALL iom_rstput( iter, nitrst, numriw, 'cnd_ice', cnd_ice )
-         CALL iom_rstput( iter, nitrst, numriw, 't1_ice' , t1_ice  )
-      ENDIF
+         CALL iom_rstput( iter, nitrst, numriw, 'cnd_ice', cnd_ice, ldxios = lwxios )
+         CALL iom_rstput( iter, nitrst, numriw, 't1_ice' , t1_ice, ldxios = lwxios  )
+      ENDIF
+      IF(lwxios) CALL iom_swap(      cxios_context          )
       !
 
@@ -156 +176 @@
       ! ------------------
       IF( iter == nitrst ) THEN
-         CALL iom_close( numriw )
+         IF(lwxios) THEN
+            CALL iom_context_finalize(      cwixios_context          )
+         ELSE
+            CALL iom_close( numriw )
+         ENDIF
          lrst_ice = .FALSE.
       ENDIF