Changeset 8243 for branches/UKMO/test_moci_test_suite/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90

Timestamp:

2017-06-29T11:41:55+02:00 (7 years ago)

Author:

andmirek

Message:

#1914 working XIOS read, XIOS write and single processor read

File:

• branches/UKMO/test_moci_test_suite/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90 (modified) (5 diffs)

branches/UKMO/test_moci_test_suite/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90

@@ -19 +19 @@
    USE iom             ! IOM library
    USE lib_mpp         ! MPP library
+   USE iom_def, ONLY : lwxios
 
    IMPLICIT NONE
@@ -29 +30 @@
    LOGICAL, PUBLIC ::   ln_ref_apr   !: ref. pressure: global mean Patm (F) or a constant (F)
    REAL(wp)        ::   rn_pref      !  reference atmospheric pressure   [N/m2]
+   LOGICAL         ::   ln_apr_sio   ! single processor read flag
 
    REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   ssh_ib    ! Inverse barometer now    sea surface height   [m]
@@ -77 +79 @@
       IF( kt == nit000 ) THEN                   ! First call kt=nit000 !
          !                                      ! -------------------- !
+         ln_apr_sio = .FALSE.
          REWIND( numnam_ref )              ! Namelist namsbc_apr in reference namelist : File for atmospheric pressure forcing
          READ  ( numnam_ref, namsbc_apr, IOSTAT = ios, ERR = 901)
@@ -126 +129 @@
          IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields
          !
+         lspr = ln_apr_sio
          CALL fld_read( kt, nn_fsbc, sf_apr )               !* input Patm provided at kt + nn_fsbc/2
+         lspr = .false.
          !
          !                                                  !* update the reference atmospheric pressure (if necessary)
@@ -157 +162 @@
          IF(lwp) WRITE(numout,*) 'sbc_apr : ssh_ib written in ocean restart file at it= ', kt,' date= ', ndastp
          IF(lwp) WRITE(numout,*) '~~~~'
-         CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib )
+         IF( lwxios ) CALL iom_swap(      wxios_context          )
+         CALL iom_rstput( kt, nitrst, numrow, 'ssh_ibb' , ssh_ib, lxios = lwxios )
+         IF( lwxios ) CALL iom_swap(      cxios_context          )
       ENDIF
       !