Changeset 1334 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90

Timestamp:

2009-03-03T15:07:48+01:00 (15 years ago)

Author:

smasson

Message:

complete work on time origin in outputs (ticket:335) + downward vertical axis (ticket:357)

File:

• trunk/NEMO/OPA_SRC/DIA/diaptr.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/DIA/diaptr.F90

@@ -418 +418 @@
 
       CHARACTER (len=40)       ::   clhstnam, clop                   ! temporary names
-      INTEGER                  ::   iline, it, ji                    !
+      INTEGER                  ::   iline, it, ji, itmod             !
       REAL(wp)                 ::   zsto, zout, zdt, zjulian   ! temporary scalars
       REAL(wp), DIMENSION(jpj) ::   zphi, zfoo
@@ -424 +424 @@
       
       ! define time axis
-      it = kt - nit000 + 1
+      it = kt
+      itmod = kt - nit000 + 1
 
       ! Initialization
@@ -481 +482 @@
          ! Horizontal grid : zphi()
          CALL histbeg(clhstnam, 1, zfoo, jpj, zphi,   &
-            1, 1, 1, jpj, 0, zjulian, zdt, nhoridz, numptr, domain_id=nidom )
+            1, 1, 1, jpj, nit000-1, zjulian, zdt, nhoridz, numptr, domain_id=nidom )
          ! Vertical grids : gdept_0, gdepw_0
          CALL histvert( numptr, "deptht", "Vertical T levels",   &
-            "m", jpk, gdept_0, ndepidzt )
+            "m", jpk, gdept_0, ndepidzt, "down" )
          CALL histvert( numptr, "depthw", "Vertical W levels",   &
-            "m", jpk, gdepw_0, ndepidzw )
+            "m", jpk, gdepw_0, ndepidzw, "down" )
          
          !  Zonal mean T and S
@@ -555 +556 @@
       ENDIF
 
-      IF( MOD( it, nf_ptr ) == 0 ) THEN
+      IF( MOD( itmod, nf_ptr ) == 0 ) THEN
 
          IF(lwp) THEN