Changeset 11741 for NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_user.F90

Timestamp:

2019-10-21T12:26:39+02:00 (4 years ago)

Author:

jchanut

Message:

#2222: correct definition of parent vertical grid on the child domain to perform vertical interpolation at boundaries. Use additionnal parent depths and number of levels arrays interpolated on the child grid domain to do so.
Correction of vertical interpolation of viscosity remains to be done as well as duplication of changes for passive tracers.

File:

• NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_user.F90 (modified) (7 diffs)

NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/agrif_user.F90

@@ -97 +97 @@
       USE oce
       USE lib_mpp
-      !
-      IMPLICIT NONE
-      !
+      USE lbclnk
+      !
+      IMPLICIT NONE
+      !
+      INTEGER :: ji, jj
       LOGICAL :: check_namelist
       CHARACTER(len=15) :: cl_check1, cl_check2, cl_check3, cl_check4 
+#if defined key_vertical
+      REAL(wp), DIMENSION(jpi,jpj) ::   zk   ! workspace
+#endif
       !!----------------------------------------------------------------------
 
@@ -110 +115 @@
       ! 2. First interpolations of potentially non zero fields
       !-------------------------------------------------------
+
+#if defined key_vertical
+      ! Build consistent parent bathymetry and number of levels
+      ! on the child grid 
+      Agrif_UseSpecialValue = .FALSE.
+      ht0_parent(:,:) = 0._wp
+      mbkt_parent(:,:) = 0
+      !
+      CALL Agrif_Bc_variable(ht0_id ,calledweight=1.,procname=interpht0 )
+      CALL Agrif_Bc_variable(mbkt_id,calledweight=1.,procname=interpmbkt)
+      !
+      DO jj = 1, jpjm1
+         DO ji = 1, jpim1
+            hu0_parent(ji,jj)  = MIN( ht0_parent(ji,jj), ht0_parent(ji+1,jj) ) 
+            hv0_parent(ji,jj)  = MIN( ht0_parent(ji,jj), ht0_parent(ji,jj+1) ) 
+            mbku_parent(ji,jj) = MIN(  mbkt_parent(ji+1,jj  ) , mbkt_parent(ji,jj)  )
+            mbkv_parent(ji,jj) = MIN(  mbkt_parent(ji  ,jj+1) , mbkt_parent(ji,jj)  )
+         END DO
+      END DO
+      !
+      CALL lbc_lnk( 'Agrif_InitValues_cont', hu0_parent, 'U', 1. )
+      CALL lbc_lnk( 'Agrif_InitValues_cont', hv0_parent, 'V', 1. )
+      zk(:,:) = REAL( mbku_parent(:,:), wp )   ;   CALL lbc_lnk( 'Agrif_InitValues_cont', zk, 'U', 1. )
+      mbku_parent(:,:) = MAX( NINT( zk(:,:) ), 1 )
+      zk(:,:) = REAL( mbkv_parent(:,:), wp )   ;   CALL lbc_lnk( 'Agrif_InitValues_cont', zk, 'V', 1. )
+      mbkv_parent(:,:) = MAX( NINT( zk(:,:) ), 1 )      
+#endif
+
       Agrif_SpecialValue    = 0._wp
       Agrif_UseSpecialValue = .TRUE.
@@ -225 +258 @@
             END IF
          ENDIF
+
+#if defined key_vertical
+    IF ( Agrif_Parent(jpk).GT.jpk ) THEN
+       CALL ctl_stop( ' With key_vertical, child grids must have jpk greater or equal to the parent value' )
+         ENDIF
+#endif
          !
       ENDIF
@@ -275 +314 @@
       CALL agrif_declare_variable((/1,2,0/),(/ind2,ind3,0/),(/'x','y','N'/),(/1,1,1/),(/nlci,nlcj,jpk/),umsk_id)
       CALL agrif_declare_variable((/2,1,0/),(/ind3,ind2,0/),(/'x','y','N'/),(/1,1,1/),(/nlci,nlcj,jpk/),vmsk_id)
+# if defined key_vertical
+      CALL agrif_declare_variable((/2,2/),(/ind3,ind3/),(/'x','y'/),(/1,1/),(/nlci,nlcj/),mbkt_id)
+      CALL agrif_declare_variable((/2,2/),(/ind3,ind3/),(/'x','y'/),(/1,1/),(/nlci,nlcj/),ht0_id)
+# endif
 
       CALL agrif_declare_variable((/2,2,0,0/),(/ind3,ind3,0,0/),(/'x','y','N','N'/),(/1,1,1,1/),(/nlci,nlcj,jpk,3/),scales_t_id)
@@ -311 +354 @@
       CALL Agrif_Set_bcinterp(ub2b_interp_id,interp1=Agrif_linear,interp2=AGRIF_ppm)
       CALL Agrif_Set_bcinterp(vb2b_interp_id,interp1=AGRIF_ppm,interp2=Agrif_linear)
-
+!
+! > Divergence conserving alternative:
+!      CALL Agrif_Set_bcinterp(sshn_id,interp=AGRIF_constant)
+!      CALL Agrif_Set_bcinterp(unb_id,interp1=Agrif_linear,interp2=AGRIF_constant)
+!      CALL Agrif_Set_bcinterp(vnb_id,interp1=AGRIF_constant,interp2=Agrif_linear)
+!      CALL Agrif_Set_bcinterp(ub2b_interp_id,interp1=Agrif_linear,interp2=AGRIF_constant)
+!      CALL Agrif_Set_bcinterp(vb2b_interp_id,interp1=AGRIF_constant,interp2=Agrif_linear)
+!<
 
       CALL Agrif_Set_bcinterp(un_sponge_id,interp1=Agrif_linear,interp2=AGRIF_ppm)
@@ -319 +369 @@
       CALL Agrif_Set_bcinterp(umsk_id,interp=AGRIF_constant)
       CALL Agrif_Set_bcinterp(vmsk_id,interp=AGRIF_constant)
+
+# if defined key_vertical
+      CALL Agrif_Set_bcinterp(mbkt_id,interp=AGRIF_constant)
+      CALL Agrif_Set_bcinterp(ht0_id ,interp=AGRIF_constant)
+# endif
 
       IF( ln_zdftke.OR.ln_zdfgls )   CALL Agrif_Set_bcinterp( avm_id, interp=AGRIF_linear )
@@ -341 +396 @@
       CALL Agrif_Set_bc( umsk_id, (/0,0/) )
       CALL Agrif_Set_bc( vmsk_id, (/0,0/) )
-
+# if defined key_vertical
+      ! extend the interpolation zone by 1 more point than necessary:
+      CALL Agrif_Set_bc(  mbkt_id, (/-nn_sponge_len*Agrif_irhox()-2,ind1/) )
+      CALL Agrif_Set_bc(  ht0_id,  (/-nn_sponge_len*Agrif_irhox()-2,ind1/) )
+# endif
 
       IF( ln_zdftke.OR.ln_zdfgls )   CALL Agrif_Set_bc( avm_id, (/0,ind1/) )