Changeset 15300

Timestamp:

2021-09-29T09:25:34+02:00 (3 years ago)

Author:

smasson

Message:

trunk: improve BENCH for reproducibility tests

File:

• NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90 (modified) (4 diffs)

NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90

@@ -62 +62 @@
       !
       INTEGER  ::   ji, jj         ! dummy loop indices
+      REAL(wp), DIMENSION(jpi,jpj) ::   z2d   ! 2D workspace
       REAL(wp) ::   zres, zf0
       REAL(wp) ::   zti, ztj       ! local scalars
@@ -72 +73 @@
       IF(lwp) WRITE(numout,*) '          given by rn_dx and rn_dy' 
       !
+      ! define unique value on each point of the inner global domain. z2d ranging from 0.05 to -0.05
+      !
+      DO_2D( 0, 0, 0, 0 )   !  +/- 0.5
+         z2d(ji,jj) = 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp )
+      END_2D
+      !
       ! Position coordinates (in grid points)
       !                          ==========
-      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
+      DO_2D( 0, 0, 0, 0 )
          
-         zti = REAL( mig0(ji), wp ) - 0.5_wp  ! start at i=0.5 in the global grid without halos
-         ztj = REAL( mjg0(jj), wp ) - 0.5_wp  ! start at j=0.5 in the global grid without halos
+         zti = REAL( mig0(ji), wp ) - 0.5_wp   ! start at i=0.5 in the global grid without halos
+         ztj = REAL( mjg0(jj), wp ) - 0.5_wp   ! start at j=0.5 in the global grid without halos
          
-         plamt(ji,jj) = zti
-         plamu(ji,jj) = zti + 0.5_wp
-         plamv(ji,jj) = zti
-         plamf(ji,jj) = zti + 0.5_wp
+         plamt(ji,jj) =   zti            * (1. + 1.0e-5 * z2d(ji,jj) )
+         plamu(ji,jj) = ( zti + 0.5_wp ) * (1. + 2.0e-5 * z2d(ji,jj) )
+         plamv(ji,jj) =   zti            * (1. + 3.0e-5 * z2d(ji,jj) )
+         plamf(ji,jj) = ( zti + 0.5_wp ) * (1. + 4.0e-5 * z2d(ji,jj) )
          
-         pphit(ji,jj) = ztj
-         pphiu(ji,jj) = ztj
-         pphiv(ji,jj) = ztj + 0.5_wp
-         pphif(ji,jj) = ztj + 0.5_wp
+         pphit(ji,jj) =   ztj            * (1. + 1.5e-5 * z2d(ji,jj) )         
+         pphiu(ji,jj) =   ztj            * (1. + 2.5e-5 * z2d(ji,jj) )         
+         pphiv(ji,jj) = ( ztj + 0.5_wp ) * (1. + 3.5e-5 * z2d(ji,jj) )
+         pphif(ji,jj) = ( ztj + 0.5_wp ) * (1. + 4.5e-5 * z2d(ji,jj) )
 
       END_2D
@@ -93 +100 @@
       ! Horizontal scale factors (in meters)
       !                              ======
-      zres = 1.e+5   !  100km
-      pe1t(:,:) = zres  ;   pe2t(:,:) = zres 
-      pe1u(:,:) = zres  ;   pe2u(:,:) = zres
-      pe1v(:,:) = zres  ;   pe2v(:,:) = zres
-      pe1f(:,:) = zres  ;   pe2f(:,:) = zres 
-
+      DO_2D( 0, 0, 0, 0 )
+         zres = 1.e+5   !  100km
+         pe1t(ji,jj) = zres * (1. + 1.0e-5 * z2d(ji,jj) )  ;   pe2t(ji,jj) = zres * (1. + 1.5e-5 * z2d(ji,jj) )
+         pe1u(ji,jj) = zres * (1. + 2.0e-5 * z2d(ji,jj) )  ;   pe2u(ji,jj) = zres * (1. + 2.5e-5 * z2d(ji,jj) )
+         pe1v(ji,jj) = zres * (1. + 3.0e-5 * z2d(ji,jj) )  ;   pe2v(ji,jj) = zres * (1. + 3.5e-5 * z2d(ji,jj) )
+         pe1f(ji,jj) = zres * (1. + 4.0e-5 * z2d(ji,jj) )  ;   pe2f(ji,jj) = zres * (1. + 4.5e-5 * z2d(ji,jj) )
+      END_2D
       !                             ! NO reduction of grid size in some straits 
       ke1e2u_v = 0                  !    ==>> u_ & v_surfaces will be computed in dom_hgr routine
@@ -109 +117 @@
       !
       zf0 = 2._wp * omega * SIN( rad * 45 )   ! constant coriolis factor corresponding to 45°N
-      pff_f(:,:) = zf0
-      pff_t(:,:) = zf0
+      DO_2D( 0, 0, 0, 0 )
+         pff_f(ji,jj) = zf0 * (1. + 1.0e-5 * z2d(ji,jj) )
+         pff_t(ji,jj) = zf0 * (1. + 2.0e-5 * z2d(ji,jj) )
+      END_2D
+      !
+      ! calls lbc_lnk done in dom_hgr
       !
    END SUBROUTINE usr_def_hgr