Changeset 2084

Timestamp:

2010-09-10T16:29:12+02:00 (14 years ago)

Author:

cetlod

Message:

minor bug correction

File:

• branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/istate.F90 (modified) (3 diffs)

branches/DEV_r2006_merge_TRA_TRC/NEMO/OPA_SRC/istate.F90

@@ -75 +75 @@
       rhop (:,:,:) = 0.e0
       rn2  (:,:,:) = 0.e0 
+      ta   (:,:,:) = 0.e0    
+      sa   (:,:,:) = 0.e0
 
       IF( ln_rstart ) THEN                    ! Restart from a file
@@ -80 +82 @@
          neuler = 1                              ! Set time-step indicator at nit000 (leap-frog)
          CALL rst_read                           ! Read the restart file
+         CALL tra_swap                           ! swap 3D arrays (ta,sa)  in a 4D array
          CALL day_init                           ! model calendar (using both namelist and restart infos)
       ELSE
@@ -123 +126 @@
          ENDIF
 
-         ta(:,:,:) = 0.       ;     sa(:,:,:) = 0.
-         CALL tra_swap
+         CALL tra_swap                     ! swap 3D arrays (ta,sa)  in a 4D array
          CALL eos( tsb, rhd, rhop )        ! before potential and in situ densities
-         
          IF( ln_zps .AND. .NOT. lk_c1d )   &
                            CALL zps_hde( nit000, jpts, tsb, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient