Changeset 1929

Timestamp:

2010-06-15T09:41:10+02:00 (14 years ago)

Author:

smasson

Message:

close mesh files for all values of nn_msh, see ticket:655

File:

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

trunk/NEMO/OPA_SRC/DOM/domwri.F90

@@ -45 +45 @@
       !!      domhgr, domzgr, and dommsk. Note: the file contain depends on
       !!      the vertical coord. used (z-coord, partial steps, s-coord)
-      !!                    nmsh = 1  :   'mesh_mask.nc' file
+      !!            MOD(nmsh, 3) = 1  :   'mesh_mask.nc' file
       !!                         = 2  :   'mesh.nc' and mask.nc' files
-      !!                         = 3  :   'mesh_hgr.nc', 'mesh_zgr.nc' and
+      !!                         = 0  :   'mesh_hgr.nc', 'mesh_zgr.nc' and
       !!                                  'mask.nc' files
       !!      For huge size domain, use option 2 or 3 depending on your 
       !!      vertical coordinate.
+      !!
+      !!      if     nmsh <= 3: write full 3D arrays for e3[tuvw] and gdep[tuvw]
+      !!      if 3 < nmsh <= 6: write full 3D arrays for e3[tuvw] and 2D arrays 
+      !!                        corresponding to the depth of the bottom points hdep[tw]
+      !!      if 6 < nmsh <= 9: write 2D arrays corresponding to the depth and the
+      !!                        thickness of the bottom points hdep[tw] and e3[tw]_ps
       !!
       !! ** output file : 
@@ -241 +247 @@
       !                                     !        close the files 
       !                                     ! ============================
-      SELECT CASE ( nmsh )
+      SELECT CASE ( MOD(nmsh, 3) )
       CASE ( 1 )                
          CALL iom_close( inum0 )
@@ -247 +253 @@
          CALL iom_close( inum1 )
          CALL iom_close( inum2 )
-      CASE ( 3 )
+      CASE ( 0 )
          CALL iom_close( inum2 )
          CALL iom_close( inum3 )