Changeset 2145

Timestamp:

2010-10-04T15:36:14+02:00 (14 years ago)

Author:

rblod

Message:

Fixx sshn_b initialization, see ticket #723

File:

• trunk/NEMO/OPA_SRC/DYN/dynspg_ts.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/DYN/dynspg_ts.F90

@@ -595 +595 @@
             vn_b (:,:) = vn_b(:,:) * hvr(:,:)
          ENDIF
+         IF( iom_varid( numror, 'sshn_b', ldstop = .FALSE. ) > 0 ) THEN
+            CALL iom_get( numror, jpdom_autoglo, 'sshn_b' , sshn_b (:,:) )   ! filtered extrenal ssh
+         ELSE
+            sshn_b(:,:)=sshb(:,:)   ! if not in restart set previous time mean to current baroclinic before value   
+         ENDIF 
       ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN
-         CALL iom_rstput( kt, nitrst, numrow, 'un_b'  , un_b  (:,:) )   ! external velocity issued
-         CALL iom_rstput( kt, nitrst, numrow, 'vn_b'  , vn_b  (:,:) )   ! from barotropic loop
+         CALL iom_rstput( kt, nitrst, numrow, 'un_b'   , un_b  (:,:) )   ! external velocity and ssh
+         CALL iom_rstput( kt, nitrst, numrow, 'vn_b'   , vn_b  (:,:) )   ! issued from barotropic loop
+         CALL iom_rstput( kt, nitrst, numrow, 'sshn_b' , sshn_b(:,:) )   ! 
       ENDIF
       !