Changeset 116 for trunk/NEMO/OPA_SRC

Timestamp:

2004-06-28T17:09:50+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE075 : Save the mbathy 2D array and the 3D scale factors e3t_ps and e3w_ps in partial steps in the mesh_zgr.nc file

File:

• trunk/NEMO/OPA_SRC/DOM/domwri.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/DOM/domwri.F90

@@ -2 +2 @@
    !!======================================================================
    !!                       ***  MODULE domwri  ***
-   !! Ocean initialization : write the ocean domain mesh and mask file(s)
+   !! Ocean initialization : write the ocean domain mesh ask file(s)
    !!======================================================================
 
@@ -88 +88 @@
       INTEGER  ::   itime           !  output from restini ???
       REAL(wp) ::   zdate0
+      REAL(wp), DIMENSION(jpi,jpj) ::   &
+         zprt                       ! temporary array for bathymetry 
+
       CHARACTER (len=21) ::      &
          clnam0 = 'mesh_mask',   &  ! filename (mesh and mask informations)
@@ -102 +105 @@
       CALL ymds2ju( 0, 1, 1, 0.e0, zdate0 )    ! calendar initialization
 
+!       note that mbathy has been modified in dommsk or in solver.
+!       it is the number of non-zero "w" levels in the water, and the minimum 
+!       value (on land) is 2. We define zprt as the number of "T" points in the ocean 
+!       at any location, and zero on land. 
+!
+      zprt = tmask(:,:,1)*(mbathy-1)
 
       SELECT CASE (nmsh)
@@ -179 +188 @@
          CALL restput( inum3, 'ff', jpi, jpj, 1, 0, ff )           !    ! coriolis factor
 
+         CALL restput( inum4, 'mbathy', jpi, jpj, 1, 0, zprt )
+
 # if defined key_s_coord
          !                                                         ! s-coordinate
@@ -197 +208 @@
          CALL restput( inum4, 'hdepw' , jpi, jpj, 1, 0, hdepw  ) 
 
-         CALL restput( inum4, 'e3tp'  , jpi, jpj, 1  , 0, e3tp   )  !    ! scale factors
-         CALL restput( inum4, 'e3wp'  , jpi, jpj, 1  , 0, e3wp   )
+         CALL restput( inum4, 'e3t_ps', jpi, jpj, jpk, 0, e3t_ps )  !    ! scale factors
          CALL restput( inum4, 'e3u_ps', jpi, jpj, jpk, 0, e3u_ps )
          CALL restput( inum4, 'e3v_ps', jpi, jpj, jpk, 0, e3v_ps )
+         CALL restput( inum4, 'e3w_ps', jpi, jpj, jpk, 0, e3w_ps )
 
          CALL restput( inum4, 'gdept' , 1, 1, jpk, 0, gdept )       !    ! reference z-coord.