Changeset 489

Timestamp:

2006-09-01T15:48:40+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_bugfix_052:CE+RB: correct time step parameter in histwrite call

File:

• trunk/NEMO/TOP_SRC/trcdit.F90 (modified) (14 diffs)

trunk/NEMO/TOP_SRC/trcdit.F90

@@ -146 +146 @@
 
       ! define time axis
-      it = kt - nit000 + 1
+      it = kt - nittrc000 + 1
 
 ! 1. Define NETCDF files and fields at beginning of first time step
@@ -152 +152 @@
 
       IF(ll_print)WRITE(numout,*)'trcdit_wr kt=',kt,' kindic ',kindic
-      IF(kt == nit000) THEN
+      IF(kt == nittrc000) THEN
 
 ! Compute julian date from starting date of the run
@@ -158 +158 @@
          CALL ymds2ju(nyear,nmonth,nday,0.0,zjulian)
          IF(lwp)WRITE(numout,*)' '  
-         IF(lwp)WRITE(numout,*)' Date 0 used :',nit000     &
+         IF(lwp)WRITE(numout,*)' Date 0 used :',nittrc000     &
        &     ,' YEAR ',nyear,' MONTH ',nmonth,' DAY ',nday   &
        &     ,'Julian day : ',zjulian    
@@ -177 +177 @@
 ! Vertical grid for tracer : gdept
          CALL histvert(nit5, 'deptht', 'Vertical T levels', &
-         &    'm', ipk, gdept, ndepit5)
+         &    'm', ipk, gdept_0, ndepit5)
 
 ! Index of ocean points in 3D and 2D (surface)
@@ -338 +338 @@
 
       ! define time axis
-      it = kt - nit000 + 1
+      it = kt - nittrc000 + 1
 
 ! Define the NETCDF files (one per tracer)
 !
       IF(ll_print)WRITE(numout,*)'trcdid kt=',kt,' kindic ',kindic
-      IF(kt == nit000) THEN
+      IF(kt == nittrc000) THEN
 
           DO jn=1,jptra
@@ -366 +366 @@
 ! Vertical grid for tracer trend - one per each tracer IF needed
                 CALL histvert(nit6(jn), 'deptht', 'Vertical T levels',  &
-                &    'm', ipk, gdept, ndepit6(jn)) 
+                &    'm', ipk, gdept_0, ndepit6(jn)) 
 
 
@@ -647 +647 @@
 
       ! define time axis
-      it = kt - nit000 + 1
+      it = kt - nittrc000 + 1
 
 ! 1. Define NETCDF files and fields at beginning of first time step
@@ -653 +653 @@
 
       IF(ll_print)WRITE(numout,*)'trcdii_wr kt=',kt,' kindic ',kindic
-      IF(kt == nit000) THEN
+      IF(kt == nittrc000) THEN
 
 ! Define the NETCDF files for additional arrays : 2D or 3D
@@ -671 +671 @@
 
           CALL histvert(nitd, 'deptht', 'Vertical T levels', &
-          &    'm', ipk, gdept, ndepitd)
+          &    'm', ipk, gdept_0, ndepitd)
 
 
@@ -729 +729 @@
           DO jn=1,jpdia2d
             cltra=ctrc2d(jn) ! short title for 2D diagnostic
-            CALL histwrite(nitd, cltra, kt, trc2d(:,:,jn), ndimt51    &
+            CALL histwrite(nitd, cltra, it, trc2d(:,:,jn), ndimt51    &
             &   ,ndext51)
           END DO
@@ -849 +849 @@
 
       ! define time axis
-      it = kt - nit000 + 1
+      it = kt - nittrc000 + 1
 
 ! 1. Define NETCDF files and fields at beginning of first time step
@@ -855 +855 @@
 
       IF(ll_print)WRITE(numout,*)'trcdib_wr kt=',kt,' kindic ',kindic
-      IF(kt == nit000) THEN
+      IF(kt == nittrc000) THEN
 
 ! Define the NETCDF files for biological trends
@@ -868 +868 @@
 ! Vertical grid for biological trends
           CALL histvert(nitb, 'deptht', 'Vertical T levels',  &
-          &    'm', ipk, gdept, ndepitb)
+          &    'm', ipk, gdept_0, ndepitb)
 
 ! Declare all the output fields as NETCDF variables
@@ -905 +905 @@
       DO jn=1,jpdiabio
          cltra=ctrbio(jn)  ! short title for biological diagnostic
-         CALL histwrite(nitb, cltra, kt, trbio(:,:,:,jn), ndimt50,ndext50)
+         CALL histwrite(nitb, cltra, it, trbio(:,:,:,jn), ndimt50,ndext50)
       END DO