Changeset 199

Timestamp:

2004-12-22T16:49:02+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE140 : Add the storage of mbathy in order to be 100% compatible with standard domwri.F90 output

File:

• trunk/NEMO/OPA_SRC/DOM/domwri_dimg.h90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DOM/domwri_dimg.h90

@@ -49 +49 @@
       INTEGER, DIMENSION(11:15)  :: irecv ! index of this array will be inumxxx
       REAL(wp) ::   zdate0
+      REAL(wp), DIMENSION(jpi,jpj) ::   &
+         zprt                       ! temporary array for bathymetry
+
       CHARACTER (len=21) ::      &
          clnam0 = 'mesh_mask',   &  ! filename (mesh and mask informations)
@@ -242 +245 @@
          WRITE(inum3,REC=irec) ff(:,:)
          irecv(inum3) = irecv(inum3) + jpnij 
-         
-            
+         !
+         ! mbathy
+         !       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)
+         irec = irecv(inum4) + (narea - 1 )
+         WRITE(inum4,REC=irec) zprt(:,:)
+         irecv(inum4) = irecv(inum4) + jpnij 
+
 #if defined key_s_coord
          !