Changeset 14223 for NEMO/trunk/tests
- Timestamp:
- 2020-12-19T11:22:45+01:00 (3 years ago)
- Location:
- NEMO/trunk/tests
- Files:
-
- 13 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90
r13295 r14223 74 74 ! Position coordinates (in grid points) 75 75 ! ========== 76 DO_2D( 1, 1, 1, 1)76 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 77 77 78 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos79 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos78 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 79 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 80 80 81 81 plamt(ji,jj) = zti -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90
r14053 r14223 53 53 REAL(wp), DIMENSION(jpi,jpj,jpk) , INTENT( out) :: pu ! i-component of the velocity [m/s] 54 54 REAL(wp), DIMENSION(jpi,jpj,jpk) , INTENT( out) :: pv ! j-component of the velocity [m/s] 55 !!st REAL(wp), DIMENSION(jpi,jpj) , INTENT( out) :: pssh ! sea-surface height56 55 ! 57 56 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace 58 57 REAL(wp) :: zfact 59 58 INTEGER :: ji, jj, jk 60 INTEGER :: igloi, igloj ! to be removed in the future, see comment bellow61 59 !!---------------------------------------------------------------------- 62 60 ! … … 67 65 ! define unique value on each point of the inner global domain. z2d ranging from 0.05 to -0.05 68 66 ! 69 ! WARNING: to keep compatibility with the trunk that was including periodocity into the input data, 70 ! we must define z2d as bellow. 71 ! Once we decide to forget trunk compatibility, we must simply define z2d as: 72 !!$ DO_2D( 0, 0, 0, 0 ) 73 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 74 !!$ END_2D 75 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 76 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 77 DO_2D( 0, 0, 0, 0 ) 78 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 67 DO_2D( 0, 0, 0, 0 ) ! +/- 0.05 68 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 79 69 END_2D 80 !81 ! sea level:82 !!st pssh(:,:) = z2d(:,:) ! +/- 0.05 m83 70 ! 84 71 DO_3D( 0, 0, 0, 0, 1, jpkm1 ) 85 72 zfact = REAL(jk-1,wp) / REAL(jpk-1,wp) ! 0 to 1 to add a basic stratification 86 73 ! temperature choosen to lead to ~50% ice at the beginning if rn_thres_sst = 0.5 87 pts(:,:,jk,jp_tem) = 20._wp*z2d(:,:) - 1._wp - 0.5_wp * zfact ! -1 to -1.5 +/- 1.0 degG74 pts(:,:,jk,jp_tem) = 20._wp*z2d(:,:) - 1._wp - 0.5_wp * zfact ! -1 to -1.5 +/- 1.0 degG 88 75 ! salinity: 89 pts(:,:,jk,jp_sal) = 30._wp + 1._wp * zfact + z2d(:,:) ! 30 to 31 +/- 0.05 psu76 pts(:,:,jk,jp_sal) = 30._wp + 1._wp * zfact + z2d(:,:) ! 30 to 31 +/- 0.05 psu 90 77 ! velocities: 91 78 pu(:,:,jk) = z2d(:,:) * 0.1_wp * umask(:,:,jk) ! +/- 0.005 m/s … … 96 83 pv( :,:,jpk ) = 0._wp 97 84 ! 98 !!st CALL lbc_lnk('usrdef_istate', pssh, 'T', 1. ) ! apply boundary conditions99 85 CALL lbc_lnk('usrdef_istate', pts, 'T', 1. ) ! apply boundary conditions 100 86 CALL lbc_lnk('usrdef_istate', pu, 'U', -1. ) ! apply boundary conditions … … 117 103 ! 118 104 INTEGER :: ji, jj 119 INTEGER :: igloi, igloj ! to be removed in the future, see usr_def_istate comment120 105 !!---------------------------------------------------------------------- 121 106 ! … … 123 108 IF(lwp) WRITE(numout,*) 'usr_def_istate_ssh : BENCH configuration, analytical definition of initial ssh' 124 109 ! 125 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 126 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 127 ! sea level: +/- 0.05 m 128 DO_2D( 0, 0, 0, 0 ) 129 pssh(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 110 DO_2D( 0, 0, 0, 0 ) ! sea level: +/- 0.05 m 111 pssh(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 130 112 END_2D 131 113 ! 132 CALL lbc_lnk('usrdef_istate', pssh, 'T', 1. ) 114 CALL lbc_lnk('usrdef_istate', pssh, 'T', 1. ) ! apply boundary conditions 133 115 ! 134 116 END SUBROUTINE usr_def_istate_ssh -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90
r13472 r14223 99 99 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace 100 100 INTEGER :: ji, jj 101 INTEGER :: igloi, igloj ! to be removed in the future, see comment bellow102 101 !!--------------------------------------------------------------------- 103 102 #if defined key_si3 … … 106 105 ! define unique value on each point. z2d ranging from 0.05 to -0.05 107 106 ! 108 ! WARNING: to keep compatibility with the trunk that was including periodocity into the input data,109 ! we must define z2d as bellow.110 ! Once we decide to forget trunk compatibility, we must simply define z2d as:111 !!$ DO_2D( 0, 0, 0, 0 )112 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) )113 !!$ END_2D114 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) )115 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) )116 107 DO_2D( 0, 0, 0, 0 ) 117 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0 _oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) )108 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 118 109 END_2D 119 110 utau_ice(:,:) = 0.1_wp + z2d(:,:) -
NEMO/trunk/tests/CANAL/MY_SRC/usrdef_hgr.F90
r13295 r14223 77 77 ! Position coordinates (in kilometers) 78 78 ! ========== 79 zlam0 = -REAL(N INT(Ni0glo*rn_0xratio)-1, wp)* rn_dx80 zphi0 = -REAL(N INT(Nj0glo*rn_0yratio)-1, wp) * rn_dy79 zlam0 = -REAL(Ni0glo, wp) * rn_0xratio * rn_dx 80 zphi0 = -REAL(Nj0glo, wp) * rn_0yratio * rn_dy 81 81 82 82 #if defined key_agrif … … 90 90 #endif 91 91 92 DO_2D( 1, 1, 1, 1)93 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos94 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos92 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 93 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 94 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 95 95 96 96 plamt(ji,jj) = zlam0 + rn_dx * zti -
NEMO/trunk/tests/ICE_ADV1D/MY_SRC/usrdef_hgr.F90
r13295 r14223 75 75 76 76 ! ========== 77 zlam0 = -REAL( (Ni0glo-2)/2, wp)* 1.e-3 * rn_dx78 zphi0 = -REAL( (Nj0glo-2)/2, wp)* 1.e-3 * rn_dy77 zlam0 = -REAL(Ni0glo, wp) * 0.5 * 1.e-3 * rn_dx 78 zphi0 = -REAL(Nj0glo, wp) * 0.5 * 1.e-3 * rn_dy 79 79 80 DO_2D( 1, 1, 1, 1)81 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos82 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos80 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 81 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 82 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 83 83 84 84 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti -
NEMO/trunk/tests/ICE_ADV2D/MY_SRC/usrdef_hgr.F90
r13295 r14223 76 76 77 77 ! ========== 78 zlam0 = -REAL( (Ni0glo-2)/2, wp)* 1.e-3 * rn_dx79 zphi0 = -REAL( (Nj0glo-2)/2, wp) * 1.e-3 * rn_dy78 zlam0 = -REAL(Ni0glo, wp) * 0.5 * 1.e-3 * rn_dx 79 zphi0 = -REAL(Nj0glo, wp) * 0.5 * 1.e-3 * rn_dy 80 80 81 81 #if defined key_agrif … … 90 90 #endif 91 91 92 DO_2D( 1, 1, 1, 1)93 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos94 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos92 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 93 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 94 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 95 95 96 96 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti -
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 ) -
NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90
r13286 r14223 83 83 ! 84 84 IF( Agrif_Root() ) THEN ! Global Domain size: ICE_AGRIF domain is 300 km x 300 Km x 10 m 85 kpi = NINT( 300.e3 / rn_dx ) - 1 86 kpj = NINT( 300.e3 / rn_dy ) - 1 87 kpi = kpi - 2 ! for compatibility with old version (because kerio=7) --> to be removed 88 kpj = kpj - 2 ! for compatibility with old version (because kerio=7) --> to be removed 85 kpi = NINT( 300.e3 / rn_dx ) - 3 86 kpj = NINT( 300.e3 / rn_dy ) - 3 89 87 ELSE ! Global Domain size: add nbghostcells + 1 "land" point on each side 90 88 kpi = nbcellsx + 2 * ( nbghostcells + 1 ) -
NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_hgr.F90
r13295 r14223 77 77 ! 78 78 ! !== grid point position ==! (in degrees) 79 DO_2D( 1, 1, 1, 1)79 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 80 80 ! ! longitude (west coast at lon=0°) 81 plamt(ji,jj) = rn_e1deg * ( - 0.5 + REAL( mig0 _oldcmp(ji)-1 , wp ) )82 plamu(ji,jj) = rn_e1deg * ( REAL( mig0 _oldcmp(ji)-1 , wp ) )81 plamt(ji,jj) = rn_e1deg * ( - 0.5 + REAL( mig0(ji)-1 , wp ) ) 82 plamu(ji,jj) = rn_e1deg * ( REAL( mig0(ji)-1 , wp ) ) 83 83 plamv(ji,jj) = plamt(ji,jj) 84 84 plamf(ji,jj) = plamu(ji,jj) 85 ! ! latitude (south coast at lat= 81°)86 pphit(ji,jj) = rn_e2deg * ( - 0.5 + REAL( mjg0 _oldcmp(jj)-1 , wp ) ) - 80._wp85 ! ! latitude (south coast at lat=-80°) 86 pphit(ji,jj) = rn_e2deg * ( - 0.5 + REAL( mjg0(jj)-1 , wp ) ) - 80._wp 87 87 pphiu(ji,jj) = pphit(ji,jj) 88 pphiv(ji,jj) = rn_e2deg * ( REAL( mjg0 _oldcmp(jj)-1 , wp ) ) - 80_wp88 pphiv(ji,jj) = rn_e2deg * ( REAL( mjg0(jj)-1 , wp ) ) - 80._wp 89 89 pphif(ji,jj) = pphiv(ji,jj) 90 90 END_2D 91 91 ! 92 92 ! !== Horizontal scale factors ==! (in meters) 93 DO_2D( 1, 1, 1, 1)93 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 94 94 ! ! e1 (zonal) 95 95 pe1t(ji,jj) = ra * rad * COS( rad * pphit(ji,jj) ) * rn_e1deg -
NEMO/trunk/tests/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90
r13295 r14223 74 74 ! !== grid point position ==! (in kilometers) 75 75 zfact = rn_dx * 1.e-3 ! conversion in km 76 DO_2D( 1, 1, 1, 1)77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0 _oldcmp(ji)-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( mig0 _oldcmp(ji)-1 , wp ) )76 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 77 ! ! longitude (west coast at lon=0°) 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0 _oldcmp(jj)-1 , wp ) )82 ! ! latitude (south coast at lat= 0°) 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( mjg0 _oldcmp(jj)-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D -
NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_hgr.F90
r13295 r14223 74 74 ! !== grid point position ==! (in kilometers) 75 75 zfact = rn_dx * 1.e-3 ! conversion in km 76 DO_2D( 1, 1, 1, 1)77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0 _oldcmp(ji)-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( mig0 _oldcmp(ji)-1 , wp ) )76 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 77 ! ! longitude (west coast at lon=0°) 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0 _oldcmp(jj)-1 , wp ) )82 ! ! latitude (south coast at lat= 0°) 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( mjg0 _oldcmp(jj)-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D -
NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_hgr.F90
r13295 r14223 80 80 IF( Agrif_Root() ) THEN 81 81 #endif 82 ! Compatibility WITH old version: 83 ! jperio = 0 => Ni0glo = jpigo_old_version 84 ! => jpiglo-1 replaced by Ni0glo-1 85 zlam0 = -REAL( (Ni0glo-1)/2, wp) * 1.e-3 * rn_dx 86 zphi0 = -REAL( (Nj0glo-1)/2, wp) * 1.e-3 * rn_dy 82 zlam0 = -REAL(Ni0glo, wp) * 0.5 * 1.e-3 * rn_dx 83 zphi0 = -REAL(Nj0glo, wp) * 0.5 * 1.e-3 * rn_dy 87 84 #if defined key_agrif 88 85 ELSE 89 86 ! ! let lower left longitude and latitude from parent 90 ! Compatibility WITH old version: 91 ! jperio = 0 => Ni0glo = jpigo_old_version 92 ! => Agrif_parent(jpiglo)-1 replaced by Agrif_parent(Ni0glo)-1 93 zlam0 = ( 0.5_wp - REAL( ( Agrif_parent(Ni0glo)-1 ) / 2, wp) ) * 1.e-3 * Agrif_irhox() * rn_dx & 87 zlam0 = ( 0.5_wp - REAL(Ni0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhox() * rn_dx & 94 88 & + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 95 zphi0 = ( 0.5_wp - REAL( ( Agrif_parent(Nj0glo)-1 ) / 2, wp)) * 1.e-3 * Agrif_irhoy() * rn_dy &89 zphi0 = ( 0.5_wp - REAL(Nj0glo, wp) * 0.5 ) * 1.e-3 * Agrif_irhoy() * rn_dy & 96 90 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 97 91 ENDIF 98 92 #endif 99 93 100 DO_2D( 1, 1, 1, 1)101 zti = REAL( mig0 _oldcmp(ji) - 1, wp ) ! start at i=0in the global grid without halos102 ztj = REAL( mjg0 _oldcmp(jj) - 1, wp ) ! start at j=0in the global grid without halos94 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 95 zti = REAL( mig0(ji), wp ) - 0.5_wp ! start at i=0.5 in the global grid without halos 96 ztj = REAL( mjg0(jj), wp ) - 0.5_wp ! start at j=0.5 in the global grid without halos 103 97 104 98 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_hgr.F90
r13295 r14223 74 74 ! !== grid point position ==! (in kilometers) 75 75 zfact = rn_dx * 1.e-3 ! conversion in km 76 DO_2D( 1, 1, 1, 1)77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0 _oldcmp(ji)-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( mig0 _oldcmp(ji)-1 , wp ) )76 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 77 ! ! longitude (west coast at lon=0°) 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0 _oldcmp(jj)-1 , wp ) )82 ! ! latitude (south coast at lat= 0°) 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( mjg0 _oldcmp(jj)-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D
Note: See TracChangeset
for help on using the changeset viewer.