Changeset 2068 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/IOM/restart.F90

Timestamp:

2010-09-06T17:56:51+02:00 (14 years ago)

Author:

mlelod

Message:

ticket: #663 ensuring restartability and conservation

File:

• branches/DEV_r1837_MLF/NEMO/OPA_SRC/IOM/restart.F90 (modified) (6 diffs)

branches/DEV_r1837_MLF/NEMO/OPA_SRC/IOM/restart.F90

@@ -26 +26 @@
    USE zdfmxl          ! mixed layer depth
    USE trdmld_oce      ! ocean active mixed layer tracers trends variables
+   USE domvvl          ! variable volume
 
    IMPLICIT NONE
@@ -38 +39 @@
 
    !! * Substitutions
+#  include "domzgr_substitute.h90"
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
@@ -123 +125 @@
       CALL iom_rstput( kt, nitrst, numrow, 'sshb'   , sshb    )
       IF( lk_vvl ) THEN
-         CALL iom_rstput( kt, nitrst, numrow, 'fse3t_b', fse3t_b )
+         CALL iom_rstput( kt, nitrst, numrow, 'fse3t_b', fse3t_b(:,:,:) )
       ENDIF
       !
@@ -157 +159 @@
       REAL(wp) ::   zrdt, zrdttra1
       INTEGER  ::   jlibalt = jprstlib
+      INTEGER  ::   jk                    ! dummy loop indices
       LOGICAL  ::   llok
       !!----------------------------------------------------------------------
@@ -195 +198 @@
       CALL iom_get( numror, jpdom_autoglo, 'sshb' , sshb  )
       IF( lk_vvl ) THEN
-         CALL iom_get( numror, jpdom_autoglo, 'fse3t_b', fse3t_b )
+         CALL iom_get( numror, jpdom_autoglo, 'fse3t_b', fse3t_b(:,:,:) )
       ENDIF
       !
@@ -224 +227 @@
          sshb (:,:)   = sshn (:,:)
          IF( lk_vvl ) THEN
-            fse3t_b(:,:,:) = fse3t_n(:,:,:)
+            DO jk = 1, jpk
+               fse3t_b(:,:,jk) = fse3t_n(:,:,jk)
+            ENDDO
          ENDIF
       ENDIF