Changeset 1334 for trunk/NEMO/LIM_SRC_3/limwri.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/LIM_SRC_3/limwri.F90 (modified) (7 diffs)

trunk/NEMO/LIM_SRC_3/limwri.F90

@@ -112 +112 @@
          nice, nhorid, ndim, niter, ndepid
       INTEGER , SAVE ::      &
-         nicea, nhorida, nitera, ndimitd
+         nicea, nhorida, ndimitd
       INTEGER , DIMENSION( jpij ) , SAVE ::  &
          ndex51
@@ -141 +141 @@
          ENDIF
          zout     = nwrite * rdt_ice / nn_fsbc
-         niter    = 0
+         niter    = nit000 - 1
          zdept(1) = 0.
 
@@ -147 +147 @@
          zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment
          CALL dia_nam ( clhstnam, nwrite, 'icemod' )
-         CALL histbeg ( clhstnam, jpi, glamt, jpj, gphit, 1, jpi, 1, jpj, 0, zjulian, rdt_ice, nhorid, nice, domain_id=nidom )
-         CALL histvert( nice, "deptht", "Vertical T levels", "m", 1, zdept, ndepid)
+         CALL histbeg ( clhstnam, jpi, glamt, jpj, gphit, 1, jpi, 1, jpj, nit000-1, zjulian, rdt_ice,   &
+            &           nhorid, nice, domain_id=nidom )
+         CALL histvert( nice, "deptht", "Vertical T levels", "m", 1, zdept, ndepid, "down")
          CALL wheneq  ( jpij , tmask(:,:,1), 1, 1., ndex51, ndim)
 
@@ -169 +170 @@
          clop     = "ave(x)"
          zout     = nwrite * rdt_ice / nn_fsbc
-         nitera   = 0
          zdept(1) = 0.
 
@@ -175 +175 @@
          CALL histbeg ( clhstnama, jpi, glamt, jpj, gphit,         &
             1, jpi, 1, jpj,        & ! zoom
-            0, zjulian, rdt_ice,   & ! time
+            nit000 -1, zjulian, rdt_ice,   & ! time
             nhorida,               & ! ? linked with horizontal ...
             nicea , domain_id=nidom)                  ! file 
@@ -316 +316 @@
       END DO
 
-      IF ( ( nn_fsbc * niter + nit000 - 1 ) >= nitend .OR. kindic < 0 ) THEN
+      IF ( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN
          IF( lwp) WRITE(numout,*) ' Closing the icemod file '
          CALL histclo( nice )
@@ -375 +375 @@
          !     not yet implemented
 
-         IF ( ( nn_fsbc * niter + nit000 - 1 ) >= nitend .OR. kindic < 0 ) THEN
+         IF ( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN
             IF(lwp) WRITE(numout,*) ' Closing the icemod file '
             CALL histclo( nicea )