Changeset 14223 for NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90
- Timestamp:
- 2020-12-19T11:22:45+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90
r13295 r14223 79 79 IF( Agrif_Root() ) THEN 80 80 #endif 81 ! Compatibility WITH old version: 82 ! jperio = 7 => Ni0glo = jpigo_old_version - 2 83 ! => jpiglo-1 replaced by Ni0glo+1 84 zlam0 = -REAL( (Ni0glo+1)/2, wp) * 1.e-3 * rn_dx 85 zphi0 = -REAL( (Nj0glo+1)/2, wp) * 1.e-3 * rn_dy ! +1 for compatibility with old version --> to be replaced by -1 as before 81 zlam0 = -REAL(Ni0glo, wp) * 0.5 * 1.e-3 * rn_dx 82 zphi0 = -REAL(Nj0glo, wp) * 0.5 * 1.e-3 * rn_dy 86 83 #if defined key_agrif 87 84 ELSE … … 89 86 !clem zlam0 = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 90 87 !clem zphi0 = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 91 ! Compatibility WITH old version: 92 ! jperio = 0 => Ni0glo = jpigo_old_version 93 ! => Agrif_parent(jpiglo)-1 replaced by Agrif_parent(Ni0glo)-1 94 zlam0 = ( 0.5_wp - REAL( ( Agrif_parent(Ni0glo)-1 ) / 2, wp) ) * 1.e-3 * Agrif_irhox() * rn_dx & 88 zlam0 = ( 0.5_wp - REAL(Ni0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhox() * rn_dx & 95 89 & + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 96 zphi0 = ( 0.5_wp - REAL( ( Agrif_parent(Nj0glo)-1 ) / 2, wp)) * 1.e-3 * Agrif_irhoy() * rn_dy &90 zphi0 = ( 0.5_wp - REAL(Nj0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhoy() * rn_dy & 97 91 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 98 92 ENDIF 99 93 #endif 100 94 101 DO_2D( 1, 1, 1, 1)102 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos103 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos104 95 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 96 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 97 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 98 105 99 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti 106 100 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * ( zti + 0.5_wp )
Note: See TracChangeset
for help on using the changeset viewer.