Changeset 9367 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfosm.F90

Timestamp:

2018-02-28T17:23:20+01:00 (6 years ago)

Author:

mathiot

Message:

Add restart read/write via XIOS capability (#1953 and #1962 and twiki: 2017WP/Met_Office-1_Mirek_XIOSread). WARNING: need to upgrade XIOS to r1296 to compile

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfosm.F90 (modified) (4 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfosm.F90

@@ -1509 +1509 @@
      ghamv(:,:,:) = 0.
      !
+     IF( lwxios ) THEN
+        CALL iom_set_rstw_var_active('wn')
+        CALL iom_set_rstw_var_active('hbl')
+        CALL iom_set_rstw_var_active('hbli')
+     ENDIF
    END SUBROUTINE zdf_osm_init
 
@@ -1540 +1545 @@
         id1 = iom_varid( numror, 'wn'   , ldstop = .FALSE. )
         IF( id1 > 0 ) THEN                       ! 'wn' exists; read
-           CALL iom_get( numror, jpdom_autoglo, 'wn', wn )
+           CALL iom_get( numror, jpdom_autoglo, 'wn', wn, ldxios = lrxios )
            WRITE(numout,*) ' ===>>>> :  wn read from restart file'
         ELSE
@@ -1549 +1554 @@
         id2 = iom_varid( numror, 'hbli'   , ldstop = .FALSE. )
         IF( id1 > 0 .AND. id2 > 0) THEN                       ! 'hbl' exists; read and return
-           CALL iom_get( numror, jpdom_autoglo, 'hbl', hbl )
-           CALL iom_get( numror, jpdom_autoglo, 'hbli', hbli )
+           CALL iom_get( numror, jpdom_autoglo, 'hbl' , hbl , ldxios = lrxios )
+           CALL iom_get( numror, jpdom_autoglo, 'hbli', hbli, ldxios = lrxios  )
            WRITE(numout,*) ' ===>>>> :  hbl & hbli read from restart file'
            RETURN
@@ -1563 +1568 @@
      IF( TRIM(cdrw) == 'WRITE') THEN     !* Write hbli into the restart file, then return
         IF(lwp) WRITE(numout,*) '---- osm-rst ----'
-         CALL iom_rstput( kt, nitrst, numrow, 'wn'   , wn     )
-         CALL iom_rstput( kt, nitrst, numrow, 'hbl'   , hbl     )
-         CALL iom_rstput( kt, nitrst, numrow, 'hbli'   , hbli     )
+         CALL iom_rstput( kt, nitrst, numrow, 'wn'     , wn  , ldxios = lwxios )
+         CALL iom_rstput( kt, nitrst, numrow, 'hbl'    , hbl , ldxios = lwxios )
+         CALL iom_rstput( kt, nitrst, numrow, 'hbli'   , hbli, ldxios = lwxios )
         RETURN
      END IF