Changeset 13334 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/OCE/DOM/domvvl.F90

Timestamp:

2020-07-22T16:20:32+02:00 (4 years ago)

Author:

jchanut

Message:

finish bypassing ocean/ice initialization with AGRIF, #2222, #2129

File:

• NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/OCE/DOM/domvvl.F90 (modified) (4 diffs)

NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/OCE/DOM/domvvl.F90

@@ -25 +25 @@
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE timing          ! Timing
+
+   USE agrif_oce       ! initial state interpolation
+   USE agrif_oce_interp 
 
    IMPLICIT NONE
@@ -803 +806 @@
       IF( TRIM(cdrw) == 'READ' ) THEN        ! Read/initialise 
          !                                   ! ===============
-         IF( ln_rstart ) THEN                   !* Read the restart file
-            CALL rst_read_open                  !  open the restart file if necessary
-            CALL iom_get( numror, jpdom_auto, 'sshn'   , ssh(:,:,Kmm), ldxios = lrxios    )
-            !
-            id1 = iom_varid( numror, 'e3t_b', ldstop = .FALSE. )
-            id2 = iom_varid( numror, 'e3t_n', ldstop = .FALSE. )
-            id3 = iom_varid( numror, 'tilde_e3t_b', ldstop = .FALSE. )
-            id4 = iom_varid( numror, 'tilde_e3t_n', ldstop = .FALSE. )
-            id5 = iom_varid( numror, 'hdiv_lf', ldstop = .FALSE. )
-            !
+         IF( ln_rstart ) THEN
+                   !* Read the restart file
+#if defined key_agrif
+            IF ( (.NOT.Agrif_root()).AND.(ln_init_chfrpar) ) THEN
+               ! skip reading restart if initialized from parent:
+               id1 = -1 ; id2 = -1 ; id3 = -1 ; id4 = -1 ; id5 = -1
+            ELSE
+#endif
+               CALL rst_read_open                  !  open the restart file if necessary
+               CALL iom_get( numror, jpdom_auto, 'sshn'   , ssh(:,:,Kmm), ldxios = lrxios    )
+               !
+               id1 = iom_varid( numror, 'e3t_b', ldstop = .FALSE. )
+               id2 = iom_varid( numror, 'e3t_n', ldstop = .FALSE. )
+               id3 = iom_varid( numror, 'tilde_e3t_b', ldstop = .FALSE. )
+               id4 = iom_varid( numror, 'tilde_e3t_n', ldstop = .FALSE. )
+               id5 = iom_varid( numror, 'hdiv_lf', ldstop = .FALSE. )
+#if defined key_agrif
+            ENDIF
+#endif
             !                             ! --------- !
             !                             ! all cases !
@@ -926 +938 @@
                ! is set up:
 !               CALL usr_def_istate( gdept_0, tmask, ts(:,:,:,:,Kbb), uu(:,:,:,Kbb), vv(:,:,:,Kbb), ssh(:,:,Kbb)  )
-!               !
+               !
 !               DO jk=1,jpk
 !                  e3t(:,:,jk,Kbb) =  e3t_0(:,:,jk) * ( ht_0(:,:) + ssh(:,:,Kbb) ) &
-!                     &            / ( ht_0(:,:) + 1._wp -ssmask(:,:) ) * tmask(:,:,jk)
+!                     &                             / ( ht_0(:,:) + 1._wp -ssmask(:,:) ) * tmask(:,:,jk) &
+!                     &             + e3t_0(:,:,jk) * ( 1._wp - tmask(:,:,jk) )
 !               END DO
 !               e3t(:,:,:,Kmm) = e3t(:,:,:,Kbb)
-                ssh(:,:,Kmm)=0._wp
-                e3t(:,:,:,Kmm)=e3t_0(:,:,:)
-                e3t(:,:,:,Kbb)=e3t_0(:,:,:)
+               ssh(:,:,Kmm)=0._wp
+               e3t(:,:,:,Kmm)=e3t_0(:,:,:)
+               e3t(:,:,:,Kbb)=e3t_0(:,:,:)
                !
             END IF           ! end of ll_wd edits
@@ -944 +957 @@
             END IF
          ENDIF
+
+#if defined key_agrif
+         IF ( .NOT.Agrif_root().AND.(ln_init_chfrpar) ) THEN
+            ! Interpolate initial ssh from parent:
+            CALL Agrif_istate_ssh( Kbb, Kmm )
+            !
+            DO jk = 1, jpk
+               e3t(:,:,jk,Kmm) =  e3t_0(:,:,jk) * ( ht_0(:,:) + ssh(:,:,Kmm)  ) &
+                 &                              / ( ht_0(:,:) + 1._wp - ssmask(:,:) ) * tmask(:,:,jk)   &
+                 &              + e3t_0(:,:,jk) * ( 1._wp - tmask(:,:,jk) )
+            END DO
+            e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm)
+         ENDIF
+#endif
          !
       ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file