Changeset 14553 for NEMO/trunk/src/ABL/ablrst.F90

Timestamp:

2021-02-26T18:01:43+01:00 (3 years ago)

Author:

gsamson

Message:

merge ticket2628_r14502_abl_restart_xios branch into trunk; sette identical between r14502 and r14544; ticket #2628

File:

• NEMO/trunk/src/ABL/ablrst.F90 (modified) (6 diffs)

NEMO/trunk/src/ABL/ablrst.F90

@@ -45 +45 @@
       CHARACTER(len=50)   ::   clname   ! abl output restart file name
       CHARACTER(len=256)  ::   clpath   ! full path to abl output restart file 
+      CHARACTER(LEN=52)   ::   clpname  ! abl output restart file name including prefix for AGRIF
       !!----------------------------------------------------------------------
       !
@@ -74 +75 @@
             ENDIF
             !
-            CALL iom_open( TRIM(clpath)//TRIM(clname), numraw, ldwrt = .TRUE., kdlev = jpka, cdcomp = 'ABL' )
+            IF(.NOT.lwxios) THEN
+               CALL iom_open( TRIM(clpath)//TRIM(clname), numraw, ldwrt = .TRUE., kdlev = jpka, cdcomp = 'ABL' )
+            ELSE
+#if defined key_xios
+               cw_ablrst_cxt = "rstwa_"//TRIM(ADJUSTL(clkt))
+               IF( TRIM(Agrif_CFixed()) == '0' ) THEN
+                  clpname = clname
+               ELSE
+                  clpname = TRIM(Agrif_CFixed())//"_"//clname
+               ENDIF
+               numraw = iom_xios_setid(TRIM(clpath)//TRIM(clpname))
+               CALL iom_init( cw_ablrst_cxt, kdid = numraw, ld_closedef = .FALSE. )
+               CALL iom_swap( cxios_context )
+#else
+               CALL ctl_stop( 'Can not use XIOS in rst_opn' )
+#endif
+            ENDIF
             lrst_abl = .TRUE.
          ENDIF
@@ -103 +120 @@
 
       ! Write in numraw (if iter == nitrst)
-      ! ------------------ 
+      ! ------------------
       !                                                                        ! calendar control
       CALL iom_rstput( iter, nitrst, numraw, 'nn_fsbc', REAL( nn_fsbc, wp ) )      ! time-step 
       CALL iom_rstput( iter, nitrst, numraw, 'kt_abl' , REAL( iter   , wp ) )      ! date
-      CALL iom_delay_rst( 'WRITE', 'ABL', numraw )   ! save only abl delayed global communication variables
+
+      IF(.NOT.lwxios) CALL iom_delay_rst( 'WRITE', 'ABL', numraw )   ! save only abl delayed global communication variables
 
       ! Prognostic (after timestep + swap time indices = now timestep) variables
@@ -124 +142 @@
       ! ------------------
       IF( iter == nitrst ) THEN
-         CALL iom_close( numraw )
+         IF(.NOT.lwxios) THEN
+            CALL iom_close( numraw )
+         ELSE
+            CALL iom_context_finalize(      cw_ablrst_cxt          )
+            iom_file(numraw)%nfid       = 0
+            numraw = 0
+         ENDIF
          lrst_abl = .FALSE.
       ENDIF
@@ -146 +170 @@
       ENDIF
 
+      lxios_sini = .FALSE.
       CALL iom_open ( TRIM(cn_ablrst_indir)//'/'//cn_ablrst_in, numrar )
+
+      IF( lrxios) THEN
+          cr_ablrst_cxt = 'abl_rst'
+          IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for ABL'
+!         IF( TRIM(Agrif_CFixed()) == '0' ) THEN
+!            clpname = cn_ablrst_in
+!         ELSE
+!            clpname = TRIM(Agrif_CFixed())//"_"//cn_ablrst_in
+!         ENDIF
+          CALL iom_init( cr_ablrst_cxt, kdid = numrar, ld_closedef = .TRUE. )
+      ENDIF
 
       ! Time info
@@ -174 +210 @@
       CALL iom_get( numrar, jpdom_auto,'mxld_abl',mxld_abl(:,:,:           ) )
       CALL iom_get( numrar, jpdom_auto,    'pblh',    pblh(:,:             ) )
-      CALL iom_delay_rst( 'READ', 'ABL', numrar )   ! read only abl delayed global communication variables
+
+      IF(.NOT.lrxios) CALL iom_delay_rst( 'READ', 'ABL', numrar )   ! read only abl delayed global communication variables
 
    END SUBROUTINE abl_rst_read