Changeset 81 for trunk

Timestamp:

09/17/13 12:51:04 (11 years ago)

Author:

smasson

Message:

force to write hdept and hdepw in the meshmask in any cases

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90

@@ -224 +224 @@
             CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r4 )
             CALL iom_rstput( 0, 0, inum4, 'gdepw', gdepw, ktype = jp_r4 )
-         ELSE                                                   !    ! 2D bottom depth
-            DO jj = 1,jpj   
-               DO ji = 1,jpi
-                  zprt(ji,jj) = gdept(ji,jj,mbkt(ji,jj)  ) * tmask(ji,jj,1)
-                  zprw(ji,jj) = gdepw(ji,jj,mbkt(ji,jj)+1) * tmask(ji,jj,1)
-               END DO
+         ENDIF
+         DO jj = 1,jpj   
+            DO ji = 1,jpi
+               zprt(ji,jj) = gdept(ji,jj,mbkt(ji,jj)  ) * tmask(ji,jj,1)
+               zprw(ji,jj) = gdepw(ji,jj,mbkt(ji,jj)+1) * tmask(ji,jj,1)
             END DO
-            CALL iom_rstput( 0, 0, inum4, 'hdept', zprt, ktype = jp_r4 )     
-            CALL iom_rstput( 0, 0, inum4, 'hdepw', zprw, ktype = jp_r4 ) 
-         ENDIF
+         END DO
+         CALL iom_rstput( 0, 0, inum4, 'hdept', zprt, ktype = jp_r4 )     
+         CALL iom_rstput( 0, 0, inum4, 'hdepw', zprw, ktype = jp_r4 ) 
          !
          CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0 )     !    ! reference z-coord.