Changeset 460 for trunk/NEMO/OPA_SRC/DIA/diawri.F90

Timestamp:

2006-05-10T19:10:25+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_051:RB: update diagnostics routines following the new coordinate definition (gdept -> gdept_0)

File:

• trunk/NEMO/OPA_SRC/DIA/diawri.F90 (modified) (8 diffs)

trunk/NEMO/OPA_SRC/DIA/diawri.F90

@@ -29 +29 @@
    USE flx_oce         ! sea-ice/ocean forcings variables
    USE diadimg         ! dimg direct access file format output
+   USE ioipsl
 
    IMPLICIT NONE
@@ -100 +101 @@
       !!   9.0  !  05-11  (V. Garnier) Surface pressure gradient organization
       !!----------------------------------------------------------------------
-      !! * Modules used
-      USE ioipsl
-
       !! * Arguments
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index
@@ -186 +184 @@
             &          0, zjulian, zdt, nh_T, nid_T, domain_id=nidom )
          CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept
-            &           "m", ipk, gdept, nz_T )
+            &           "m", ipk, gdept_0, nz_T )
          !                                                            ! Index of ocean points
          CALL wheneq( jpi*jpj*ipk, tmask, 1, 1., ndex_T , ndim_T  )      ! volume
@@ -199 +197 @@
             &          0, zjulian, zdt, nh_U, nid_U, domain_id=nidom )
          CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept
-            &           "m", ipk, gdept, nz_U )
+            &           "m", ipk, gdept_0, nz_U )
          !                                                            ! Index of ocean points
          CALL wheneq( jpi*jpj*ipk, umask, 1, 1., ndex_U , ndim_U  )      ! volume
@@ -212 +210 @@
             &          0, zjulian, zdt, nh_V, nid_V, domain_id=nidom )
          CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept
-            &          "m", ipk, gdept, nz_V )
+            &          "m", ipk, gdept_0, nz_V )
          !                                                            ! Index of ocean points
          CALL wheneq( jpi*jpj*ipk, vmask, 1, 1., ndex_V , ndim_V  )      ! volume
@@ -225 +223 @@
             &          0, zjulian, zdt, nh_W, nid_W, domain_id=nidom )
          CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw
-            &          "m", ipk, gdepw, nz_W )
+            &          "m", ipk, gdepw_0, nz_W )
  
 
@@ -539 +537 @@
       !!    "   !  05-11  (V. Garnier) Surface pressure gradient organization
       !!----------------------------------------------------------------------
-      !! * Modules used
-      USE ioipsl
-
       !! * Arguments
       CHARACTER (len=* ), INTENT( in ) ::   &
@@ -580 +575 @@
           1, jpi, 1, jpj, 0, zjulian, zdt, nh_i, id_i, domain_id=nidom )          ! Horizontal grid : glamt and gphit
       CALL histvert( id_i, "deptht", "Vertical T levels",   &    ! Vertical grid : gdept
-          "m", jpk, gdept, nz_i)
+          "m", jpk, gdept_0, nz_i)
 
       ! Declare all the output fields as NetCDF variables