Changeset 311 for trunk/NEMO/OPA_SRC/restart.F90

Timestamp:

2005-09-30T12:18:52+02:00 (19 years ago)

Author:

opalod

Message:

nemo_v1_update_017:RB: added extra outer halo (parameters jpr2di and jpr2dj) and the corresponding lbc_lnk_e for boundary conditions.It will be use for nsolv=4.

File:

• trunk/NEMO/OPA_SRC/restart.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/restart.F90

@@ -83 +83 @@
       REAL(wp), DIMENSION( 1) ::   zfice, zfblk   ! used only in case of ice & bulk
       REAL(wp), DIMENSION(10) ::   zinfo(10)
+      REAL(wp), DIMENSION(jpi,jpj) :: ztab 
       !!----------------------------------------------------------------------
       !!  OPA 9.0 , LOCEAN-IPSL (2005) 
@@ -159 +160 @@
          CALL restput( numwrs, 'hdivn'  , jpi, jpj, jpk, 0, hdivn   )
 
-         CALL restput( numwrs, 'gcx'    , jpi, jpj, 1  , 0, gcx     )   ! Read elliptic solver arrays
-         CALL restput( numwrs, 'gcxb'   , jpi, jpj, 1  , 0, gcxb    )
+         ztab(:,:) = gcx(1:jpi,1:jpj)
+         CALL restput( numwrs, 'gcx'    , jpi, jpj, 1  , 0, ztab    )   ! Read elliptic solver arrays
+         ztab(:,:) = gcxb(1:jpi,1:jpj)
+         CALL restput( numwrs, 'gcxb'   , jpi, jpj, 1  , 0, ztab    )
 # if defined key_dynspg_fsc
          CALL restput( numwrs, 'sshb'   , jpi, jpj, 1  , 0, sshb    )   ! free surface formulation (ssh)
@@ -243 +246 @@
       REAL(wp) ::   zdate0, zdt, zinfo(10)
       REAL(wp) ::   zdept(jpk), zlamt(jpi,jpj), zphit(jpi,jpj)
+      REAL(wp), DIMENSION(jpi,jpj) :: ztab 
 #   if defined key_ice_lim
       INTEGER  ::   ios1, ji, jj, jn
@@ -346 +350 @@
       CALL restget( inum, 'hdivn'  , jpi, jpj, jpk, 0, llog, hdivn   )
 
-      CALL restget( inum, 'gcxb'   , jpi, jpj, 1  , 0, llog, gcxb    )   ! Read elliptic solver arrays
-      CALL restget( inum, 'gcx'    , jpi, jpj, 1  , 0, llog, gcx     )
+      CALL restget( inum, 'gcxb'   , jpi, jpj, 1  , 0, llog, ztab    )   ! Read elliptic solver arrays
+      gcxb(1:jpi,1:jpj) = ztab(:,:) 
+      CALL restget( inum, 'gcx'    , jpi, jpj, 1  , 0, llog, ztab    )
+      gcx(1:jpi,1:jpj) = ztab(:,:) 
 # if defined key_dynspg_fsc
       CALL restget( inum, 'sshb'   , jpi, jpj, 1  , 0, llog, sshb    )   ! free surface formulation (ssh)