Changeset 389 for trunk/NEMO/OPA_SRC/restart_dimg.h90

Timestamp:

2006-03-09T18:22:04+01:00 (18 years ago)

Author:

opalod

Message:

RB:nemo_v1_update_038: first integration of Agrif :

• configuration parameters are just integer when agrif is used
• add call to agrif routines with key_agrif

File:

• trunk/NEMO/OPA_SRC/restart_dimg.h90 (modified) (10 diffs)

trunk/NEMO/OPA_SRC/restart_dimg.h90

@@ -62 +62 @@
        irecl8= jpi * jpj * wp
        WRITE(clres,'(a,i3.3)') 'restart.output.',narea
+#if defined key_AGRIF 
+       inum = Agrif_Get_Unit()      
+       If(.NOT. Agrif_root() ) clres = TRIM(Agrif_CFixed())//'_'//TRIM(clres)
+#endif 
        OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl8 )
 
@@ -166 +170 @@
 
        ! TKE arrays
+
 #if defined key_zdftke
          DO jk = 1, jpk
-            WRITE(inum,REC=irec) en(:,:,jk)   ;   irec = irec + 1 
+            WRITE(inum,REC=irec) en(:,:,jk) ; irec = irec + 1 
          END DO
 #endif
@@ -174 +179 @@
 #if defined key_ice_lim
           zfice(1) = FLOAT( nfice )                                      ! Louvain La Neuve Sea Ice Model
-          WRITE(inum,REC=irec) zfice(:)      ;   irec = irec + 1
-          WRITE(inum,REC=irec) sst_io(:,:)   ;   irec = irec + 1
-          WRITE(inum,REC=irec) sss_io(:,:)   ;   irec = irec + 1
-          WRITE(inum,REC=irec) u_io  (:,:)   ;   irec = irec + 1
-          WRITE(inum,REC=irec) v_io  (:,:)   ;   irec = irec + 1
+          WRITE(inum,REC=irec) zfice(:)     ; irec = irec + 1
+          WRITE(inum,REC=irec) sst_io(:,:)  ; irec = irec + 1
+          WRITE(inum,REC=irec) sss_io(:,:)  ; irec = irec + 1
+          WRITE(inum,REC=irec) u_io  (:,:)  ; irec = irec + 1
+          WRITE(inum,REC=irec) v_io  (:,:)  ; irec = irec + 1
 #    if defined key_coupled
-          WRITE(inum,REC=irec) alb_ice(:,:)  ;   irec = irec + 1
+          WRITE(inum,REC=irec) alb_ice(:,:)  ; irec = irec + 1
 #    endif
 #endif
 # if defined key_flx_bulk_monthly || defined key_flx_bulk_daily
           zfblk(1) = FLOAT( nfbulk )                                 ! Bulk
-          WRITE(inum,REC=irec) zfblk(:)   ;   irec = irec + 1
-          WRITE(inum,REC=irec) gsst(:,:)  ;   irec = irec + 1
+          WRITE(inum,REC=irec) zfblk(:)   ; irec = irec + 1
+          WRITE(inum,REC=irec) gsst(:,:)  ; irec = irec + 1
 # endif
 
@@ -229 +234 @@
     LOGICAL   :: lstop
 
-    REAL(wp), DIMENSION( 1) ::   zfice, zfblk   ! used only in case of ice & bulk
+      REAL(wp), DIMENSION( 1) ::   zfice, zfblk   ! used only in case of ice & bulk
     !!----------------------------------------------------------------------
 
@@ -268 +273 @@
     ! Open direct access file, with reclength for 2D wp fields
     WRITE(clres,'(a,i3.3)') 'restart.',narea
-
+#if defined key_AGRIF       
+    inum = Agrif_Get_Unit()
+    If(.NOT. Agrif_root() ) clres = TRIM(Agrif_CFixed())//'_'//TRIM(clres)
+#endif
     OPEN(inum,FILE=clres,FORM='UNFORMATTED', ACCESS='DIRECT', RECL=8 )
     READ(inum,REC=1)irecl8
@@ -278 +286 @@
     ! --------------
 
-    READ(inum,REC=1) irecl8, ino1, it1, isor1, ipcg1, itke1,   &
-       &             iice1, ibulk1, ios1, ios2, ios3, ios4,    &
-       &             idast1, adatrj0,  ipi,ipj,ipk,ipni,ipnj,ipnij,iarea
+
+    READ(inum,REC=1) irecl8, ino1, it1, isor1, ipcg1, itke1, &
+     &  iice1, ibulk1, ios1, ios2, ios3, ios4, &
+     &  idast1, adatrj0,  ipi,ipj,ipk,ipni,ipnj,ipnij,iarea
 
     ! Performs checks on the file
@@ -416 +425 @@
 
     ! TKE arrays
+
 #if defined key_zdftke
     IF ( itke1 == 1 ) THEN
@@ -433 +443 @@
     ! check if it was in the previous run
     IF ( ios1 == 1 ) THEN
-       READ(inum,REC=irec) zfice(:)      ;   irec = irec + 1
-       READ(inum,REC=irec) sst_io(:,:)   ;   irec = irec + 1
-       READ(inum,REC=irec) sss_io(:,:)   ;   irec = irec + 1
-       READ(inum,REC=irec) u_io  (:,:)   ;   irec = irec + 1
-       READ(inum,REC=irec) v_io  (:,:)   ;   irec = irec + 1
-# if defined key_coupled
-       READ(inum,REC=irec) alb_ice(:,:)   ;   irec = irec + 1
-# endif
+       READ(inum,REC=irec) zfice(:)    ; irec = irec + 1
+       READ(inum,REC=irec) sst_io(:,:) ; irec = irec + 1
+       READ(inum,REC=irec) sss_io(:,:) ; irec = irec + 1
+       READ(inum,REC=irec) u_io  (:,:) ; irec = irec + 1
+       READ(inum,REC=irec) v_io  (:,:) ; irec = irec + 1
+#  if defined key_coupled
+       READ(inum,REC=irec) alb_ice(:,:) ; irec = irec + 1
+#  endif
     ENDIF
     IF ( zfice(1) /= FLOAT(nfice) .OR. ios1 == 0 ) THEN
@@ -454 +464 @@
             END DO
          END DO
-# if defined key_coupled
+#    if defined key_coupled
          alb_ice(:,:) = 0.8 * tmask(:,:,1)
-# endif
+#    endif
     ENDIF
 #endif
@@ -489 +499 @@
     ENDIF
 
+
   END SUBROUTINE rst_read