Changeset 13141 for NEMO/branches/2020/dev_r12973_AGRIF_CMEMS/src/NST/agrif_user.F90

Timestamp:

2020-06-22T18:27:34+02:00 (4 years ago)

Author:

jchanut

Message:

#2129, corrections/add ons to initial state interpolation with AGRIF
1) add namelist flag for child grid initial state interpolation - ice not considered yet
2) provide depths and not thicknesses as inputs to vertical linear interpolation
3) extend initial state interpolation to a restart scenario for parent grid (warning should be added in that case in order to prevent users doing this at each model restart...)
The online interpolation seems to work fine in the VORTEX case (provided 0. is not considered as a special value in the initial velocity field, i.e. ln_spc_dyn=F)

File:

• NEMO/branches/2020/dev_r12973_AGRIF_CMEMS/src/NST/agrif_user.F90 (modified) (5 diffs)

NEMO/branches/2020/dev_r12973_AGRIF_CMEMS/src/NST/agrif_user.F90

@@ -52 +52 @@
       INTEGER :: jn
 
+      IF(lwp) WRITE(numout,*) ' '
+      IF(lwp) WRITE(numout,*) 'AGRIF: interp child initial state from parent'
+      IF(lwp) WRITE(numout,*) ' '
+
       l_ini_child = .TRUE.
       Agrif_SpecialValue    = 0._wp
       Agrif_UseSpecialValue = .TRUE.
+      CALL Agrif_childgrid_to_parentgrid()
+      l_ini_child = .TRUE.
+      CALL Agrif_parentgrid_to_childgrid()  
       uu(:,:,:,:) = 0.  ;  vv(:,:,:,:) = 0.   ;  ts(:,:,:,:,:) = 0.
        
@@ -61 +68 @@
       ! Brutal fix to pas 1x1 refinment. 
   !    IF(Agrif_Irhox() == 1) THEN
-         CALL Agrif_Init_Variable(tsini_id, procname=agrif_initts)
+  !       CALL Agrif_Init_Variable(tsini_id, procname=agrif_initts)
   !    ELSE
-   !      CALL Agrif_Init_Variable(tsini_id, procname=interptsn)
+      CALL Agrif_Init_Variable(tsini_id, procname=interptsn)
 
   !    ENDIF
+! just for VORTEX because Parent velocities can actually be exactly zero
+!      Agrif_UseSpecialValue = .FALSE.
       Agrif_UseSpecialValue = ln_spc_dyn
       use_sign_north = .TRUE.
       sign_north = -1.
-  !    CALL Agrif_Init_Variable(uini_id , procname=interpun )
-  !    CALL Agrif_Init_Variable(vini_id , procname=interpvn )
-       use_sign_north = .FALSE.
+      CALL Agrif_Init_Variable(uini_id , procname=interpun )
+      CALL Agrif_Init_Variable(vini_id , procname=interpvn )
+      use_sign_north = .FALSE.
 
       Agrif_UseSpecialValue = .FALSE.            !
       l_ini_child = .FALSE.
+      CALL Agrif_childgrid_to_parentgrid()
+      l_ini_child = .FALSE.
+      CALL Agrif_parentgrid_to_childgrid() 
+
       Krhs_a = Kaa ; Kmm_a = Kmm
 
@@ -281 +294 @@
       mbkv_parent(:,:) = MAX( NINT( zk(:,:) ), 1 )   
 
-
-      CALL Agrif_Init_Variable(sshini_id, procname=agrif_initssh)
-      CALL lbc_lnk( 'Agrif_Init_Domain', ssh(:,:,Kbb), 'T', 1. )
-      DO jk = 1, jpk
-            e3t(:,:,jk,Kbb) =  e3t_0(:,:,jk) * ( ht_0(:,:) + ssh(:,:,Kbb)  ) &
-      &                            / ( ht_0(:,:) + 1._wp - ssmask(:,:) ) * tmask(:,:,jk)   &
-                     &              + e3t_0(:,:,jk) * ( 1._wp - tmask(:,:,jk) )
-      END DO
+      IF ( ln_init_chfrpar ) THEN 
+         CALL Agrif_Init_Variable(sshini_id, procname=agrif_initssh)
+         CALL lbc_lnk( 'Agrif_Init_Domain', ssh(:,:,Kbb), 'T', 1. )
+         DO jk = 1, jpk
+               e3t(:,:,jk,Kbb) =  e3t_0(:,:,jk) * ( ht_0(:,:) + ssh(:,:,Kbb)  ) &
+                        &             / ( ht_0(:,:) + 1._wp - ssmask(:,:) ) * tmask(:,:,jk)   &
+                        &              + e3t_0(:,:,jk) * ( 1._wp - tmask(:,:,jk) )
+         END DO
+      ENDIF
 
       ! check if masks and bathymetries match
@@ -850 +864 @@
       !
       INTEGER  ::   ios                 ! Local integer output status for namelist read
-      NAMELIST/namagrif/ ln_agrif_2way, rn_sponge_tra, rn_sponge_dyn, rn_trelax_tra, rn_trelax_dyn, &
+      NAMELIST/namagrif/ ln_agrif_2way, ln_init_chfrpar, rn_sponge_tra, rn_sponge_dyn, rn_trelax_tra, rn_trelax_dyn, &
                        & ln_spc_dyn, ln_chk_bathy, ln_bry_south
       !!--------------------------------------------------------------------------------------
@@ -866 +880 @@
          WRITE(numout,*) '   Namelist namagrif : set AGRIF parameters'
          WRITE(numout,*) '      Two way nesting activated ln_agrif_2way         = ', ln_agrif_2way
+         WRITE(numout,*) '      child initial state from parent ln_init_chfrpar = ', ln_init_chfrpar
          WRITE(numout,*) '      sponge coefficient for tracers    rn_sponge_tra = ', rn_sponge_tra, ' m^2/s'
          WRITE(numout,*) '      sponge coefficient for dynamics   rn_sponge_tra = ', rn_sponge_dyn, ' m^2/s'