Changeset 8644 for branches/2017/dev_r8600_xios_write/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90

Timestamp:

2017-10-20T10:04:07+02:00 (7 years ago)

Author:

andmirek

Message:

ticket #1962 xios write functionality works

File:

• branches/2017/dev_r8600_xios_write/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90 (modified) (5 diffs)

branches/2017/dev_r8600_xios_write/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90

@@ -144 +144 @@
       CALL day( nit000 )
       !
+      IF( lwxios ) THEN
+! define variables in restart file when writing with XIOS
+          CALL set_rstw_var_active('kt')
+          CALL set_rstw_var_active('ndastp')
+          CALL set_rstw_var_active('adatrj')
+          CALL set_rstw_var_active('ntime')
+      ENDIF
+
    END SUBROUTINE day_init
 
@@ -341 +349 @@
                ndastp = NINT( zndastp )
                CALL iom_get( numror, 'adatrj', adatrj  )
-          CALL iom_get( numror, 'ntime', ktime )
-          nn_time0=INT(ktime)
+               CALL iom_get( numror, 'ntime', ktime )
+               nn_time0=INT(ktime)
                ! calculate start time in hours and minutes
-          zdayfrac=adatrj-INT(adatrj)
-          ksecs = NINT(zdayfrac*86400)        ! Nearest second to catch rounding errors in adatrj         
-          ihour = INT(ksecs/3600)
-          iminute = ksecs/60-ihour*60
-           
+               zdayfrac=adatrj-INT(adatrj)
+               ksecs = NINT(zdayfrac*86400)       ! Nearest second to catch rounding errors in adatrj        
+               ihour = INT(ksecs/3600)
+               iminute = ksecs/60-ihour*60
+                 
                ! Add to nn_time0
                nhour   =   nn_time0 / 100
                nminute = ( nn_time0 - nhour * 100 )
-          nminute=nminute+iminute
-          
-          IF( nminute >= 60 ) THEN
-             nminute=nminute-60
-        nhour=nhour+1
-          ENDIF
-          nhour=nhour+ihour
-          IF( nhour >= 24 ) THEN
-        nhour=nhour-24
-             adatrj=adatrj+1
-          ENDIF          
-          nn_time0 = nhour * 100 + nminute
-          adatrj = INT(adatrj)                    ! adatrj set to integer as nn_time0 updated          
+               nminute=nminute+iminute
+               
+               IF( nminute >= 60 ) THEN
+                  nminute=nminute-60
+                  nhour=nhour+1
+               ENDIF
+               nhour=nhour+ihour
+               IF( nhour >= 24 ) THEN
+                  nhour=nhour-24
+                  adatrj=adatrj+1
+               ENDIF      
+               nn_time0 = nhour * 100 + nminute
+               adatrj = INT(adatrj)                    ! adatrj set to integer as nn_time0 updated        
             ELSE
                ! parameters corresponding to nit000 - 1 (as we start the step loop with a call to day)
@@ -378 +386 @@
             ndastp = ndate0           ! ndate0 read in the namelist in dom_nam
             nhour   =   nn_time0 / 100
-       nminute = ( nn_time0 - nhour * 100 )
+            nminute = ( nn_time0 - nhour * 100 )
             IF( nhour*3600+nminute*60-ndt05 .lt. 0 )  ndastp=ndastp-1      ! Start hour is specified in the namelist (default 0)
             adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday
@@ -388 +396 @@
             WRITE(numout,*) '   date ndastp                                      : ', ndastp
             WRITE(numout,*) '   number of elapsed days since the begining of run : ', adatrj
-       WRITE(numout,*) '   nn_time0                                         : ',nn_time0
+            WRITE(numout,*) '   nn_time0                                         : ',nn_time0
             WRITE(numout,*)
          ENDIF
@@ -405 +413 @@
          CALL iom_rstput( kt, nitrst, numrow, 'adatrj' , adatrj, lxios = lwxios            )   ! number of elapsed days since
          !                                                                     ! the begining of the run [s]
-    CALL iom_rstput( kt, nitrst, numrow, 'ntime'  , REAL( nn_time0, wp), lxios = lwxios ) ! time
+         CALL iom_rstput( kt, nitrst, numrow, 'ntime'  , REAL( nn_time0, wp), lxios = lwxios ) ! time
          IF( lwxios ) CALL iom_swap(      cxios_context          )
       ENDIF