Changeset 13286 for NEMO/trunk/tests
- Timestamp:
- 2020-07-09T17:48:29+02:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 38 edited
-
. (modified) (1 prop)
-
tests/BENCH/EXPREF/namelist_cfg_orca025_like (modified) (2 diffs)
-
tests/BENCH/EXPREF/namelist_cfg_orca12_like (modified) (2 diffs)
-
tests/BENCH/EXPREF/namelist_cfg_orca1_like (modified) (2 diffs)
-
tests/BENCH/MY_SRC/usrdef_hgr.F90 (modified) (3 diffs)
-
tests/BENCH/MY_SRC/usrdef_istate.F90 (modified) (4 diffs)
-
tests/BENCH/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/BENCH/MY_SRC/usrdef_sbc.F90 (modified) (2 diffs)
-
tests/BENCH/MY_SRC/usrdef_zgr.F90 (modified) (1 diff)
-
tests/CANAL/MY_SRC/domvvl.F90 (modified) (7 diffs)
-
tests/CANAL/MY_SRC/usrdef_hgr.F90 (modified) (3 diffs)
-
tests/CANAL/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/ICE_ADV1D/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/ICE_ADV1D/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/ICE_ADV2D/MY_SRC/usrdef_hgr.F90 (modified) (4 diffs)
-
tests/ICE_ADV2D/MY_SRC/usrdef_nam.F90 (modified) (3 diffs)
-
tests/ICE_AGRIF/EXPREF/AGRIF_FixedGrids.in (modified) (1 diff)
-
tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90 (modified) (3 diffs)
-
tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90 (modified) (3 diffs)
-
tests/ISOMIP+/MY_SRC/sbcfwb.F90 (modified) (1 diff)
-
tests/ISOMIP+/MY_SRC/tradmp.F90 (modified) (1 diff)
-
tests/ISOMIP/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/ISOMIP/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/ISOMIP/MY_SRC/usrdef_zgr.F90 (modified) (3 diffs)
-
tests/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/OVERFLOW/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/OVERFLOW/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/OVERFLOW/MY_SRC/usrdef_zgr.F90 (modified) (1 diff)
-
tests/STATION_ASF/MY_SRC/nemogcm.F90 (modified) (5 diffs)
-
tests/STATION_ASF/MY_SRC/usrdef_hgr.F90 (modified) (1 diff)
-
tests/STATION_ASF/MY_SRC/usrdef_nam.F90 (modified) (1 diff)
-
tests/VORTEX/MY_SRC/domvvl.F90 (modified) (7 diffs)
-
tests/VORTEX/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/VORTEX/MY_SRC/usrdef_nam.F90 (modified) (3 diffs)
-
tests/WAD/MY_SRC/usrdef_hgr.F90 (modified) (2 diffs)
-
tests/WAD/MY_SRC/usrdef_nam.F90 (modified) (2 diffs)
-
tests/WAD/MY_SRC/usrdef_zgr.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 2 2 ^/utils/build/makenemo@HEAD makenemo 3 3 ^/utils/build/mk@HEAD mk 4 ^/utils/tools /@HEADtools4 ^/utils/tools@HEAD tools 5 5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM … … 8 8 9 9 # SETTE 10 ^/utils/CI/ sette@12931sette10 ^/utils/CI/r12931_sette_ticket2366@HEAD sette
-
- Property svn:externals
-
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca025_like
r12530 r13286 15 15 &namusr_def ! User defined : BENCH configuration: Flat bottom, beta-plane 16 16 !----------------------------------------------------------------------- 17 nn_isize = 144 2! number of point in i-direction of global(local) domain if >0 (<0)18 nn_jsize = 120 7 !! 1050! number of point in j-direction of global(local) domain if >0 (<0)17 nn_isize = 1440 ! number of point in i-direction of global(local) domain if >0 (<0) 18 nn_jsize = 1206 !! 1049 ! number of point in j-direction of global(local) domain if >0 (<0) 19 19 nn_ksize = 75 ! total number of point in k-direction 20 20 nn_perio = 4 ! periodicity … … 30 30 &namctl ! Control prints (default: OFF) 31 31 !----------------------------------------------------------------------- 32 nn_print = 0 ! level of print (0 no extra print)33 32 ln_timing = .false. ! timing by routine write out in timing.output file 34 33 / -
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca12_like
r12530 r13286 15 15 &namusr_def ! User defined : BENCH configuration: Flat bottom, beta-plane 16 16 !----------------------------------------------------------------------- 17 nn_isize = 432 2! number of point in i-direction of global(local) domain if >0 (<0)18 nn_jsize = 314 7! number of point in j-direction of global(local) domain if >0 (<0)17 nn_isize = 4320 ! number of point in i-direction of global(local) domain if >0 (<0) 18 nn_jsize = 3146 ! number of point in j-direction of global(local) domain if >0 (<0) 19 19 nn_ksize = 75 ! total number of point in k-direction 20 20 nn_perio = 4 ! periodicity … … 30 30 &namctl ! Control prints (default: OFF) 31 31 !----------------------------------------------------------------------- 32 nn_print = 0 ! level of print (0 no extra print)33 32 ln_timing = .false. ! timing by routine write out in timing.output file 34 33 / -
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca1_like
r12530 r13286 15 15 &namusr_def ! User defined : BENCH configuration: Flat bottom, beta-plane 16 16 !----------------------------------------------------------------------- 17 nn_isize = 36 2! number of point in i-direction of global(local) domain if >0 (<0)18 nn_jsize = 33 2! number of point in j-direction of global(local) domain if >0 (<0)17 nn_isize = 360 ! number of point in i-direction of global(local) domain if >0 (<0) 18 nn_jsize = 331 ! number of point in j-direction of global(local) domain if >0 (<0) 19 19 nn_ksize = 75 ! total number of point in k-direction 20 20 nn_perio = 6 ! periodicity … … 30 30 &namctl ! Control prints (default: OFF) 31 31 !----------------------------------------------------------------------- 32 nn_print = 0 ! level of print (0 no extra print)33 32 ln_timing = .false. ! timing by routine write out in timing.output file 34 33 / -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90
r12740 r13286 61 61 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 62 62 ! 63 INTEGER :: ji, jj ! dummy loop indices63 INTEGER :: ji, jj ! dummy loop indices 64 64 REAL(wp) :: zres, zf0 65 REAL(wp) :: zti, z ui, ztj, zvj! local scalars65 REAL(wp) :: zti, ztj ! local scalars 66 66 !!------------------------------------------------------------------------------- 67 67 ! … … 72 72 IF(lwp) WRITE(numout,*) ' given by rn_dx and rn_dy' 73 73 ! 74 !75 74 ! Position coordinates (in grid points) 76 75 ! ========== 77 76 DO_2D_11_11 78 77 79 zti = REAL( ji - 1 + nimpp - 1, wp ) ; ztj = REAL( jj - 1 + njmpp - 1, wp )80 z ui = REAL( ji - 1 + nimpp - 1, wp ) + 0.5_wp ; zvj = REAL( jj - 1 + njmpp - 1, wp ) + 0.5_wp78 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 79 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 81 80 82 81 plamt(ji,jj) = zti 83 plamu(ji,jj) = z ui82 plamu(ji,jj) = zti + 0.5_wp 84 83 plamv(ji,jj) = zti 85 plamf(ji,jj) = z ui84 plamf(ji,jj) = zti + 0.5_wp 86 85 87 86 pphit(ji,jj) = ztj 88 pphiv(ji,jj) = zvj89 87 pphiu(ji,jj) = ztj 90 pphif(ji,jj) = zvj 88 pphiv(ji,jj) = ztj + 0.5_wp 89 pphif(ji,jj) = ztj + 0.5_wp 91 90 92 91 END_2D … … 109 108 kff = 1 ! indicate not to compute Coriolis parameter afterward 110 109 ! 111 zf0 110 zf0 = 2._wp * omega * SIN( rad * 45 ) ! constant coriolis factor corresponding to 45°N 112 111 pff_f(:,:) = zf0 113 112 pff_t(:,:) = zf0 -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90
r12794 r13286 57 57 REAL(wp) :: zfact 58 58 INTEGER :: ji, jj, jk 59 INTEGER :: igloi, igloj ! to be removed in the future, see comment bellow 59 60 !!---------------------------------------------------------------------- 60 61 ! … … 63 64 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ ' 64 65 ! 65 ! define unique value on each point. z2d ranging from 0.05 to -0.05 66 DO_2D_11_11 67 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig(ji) + (mjg(jj)-1) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) 66 ! define unique value on each point of the inner global domain. z2d ranging from 0.05 to -0.05 67 ! 68 ! WARNING: to keep compatibility with the trunk that was including periodocity into the input data, 69 ! we must define z2d as bellow. 70 ! Once we decide to forget trunk compatibility, we must simply define z2d as: 71 !!$ DO_2D_00_00 72 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 73 !!$ END_2D 74 igloi = Ni0glo + 2 * COUNT( (/ jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 /) ) 75 igloj = Nj0glo + 2 * COUNT( (/ jperio == 2 .OR. jperio == 7 /) ) + 1 * COUNT( (/ jperio >= 4 .AND. jperio <= 6 /) ) 76 DO_2D_00_00 77 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 68 78 END_2D 69 79 ! … … 71 81 pssh(:,:) = z2d(:,:) ! +/- 0.05 m 72 82 ! 73 DO jk = 1, jpk83 DO_3D_00_00( 1, jpkm1 ) 74 84 zfact = REAL(jk-1,wp) / REAL(jpk-1,wp) ! 0 to 1 to add a basic stratification 75 85 ! temperature choosen to lead to ~50% ice at the beginning if rn_thres_sst = 0.5 … … 80 90 pu(:,:,jk) = z2d(:,:) * 0.1_wp * umask(:,:,jk) ! +/- 0.005 m/s 81 91 pv(:,:,jk) = z2d(:,:) * 0.01_wp * vmask(:,:,jk) ! +/- 0.0005 m/s 82 ENDDO 92 END_3D 93 pts(:,:,jpk,:) = 0._wp 94 pu( :,:,jpk ) = 0._wp 95 pv( :,:,jpk ) = 0._wp 83 96 ! 84 97 CALL lbc_lnk('usrdef_istate', pssh, 'T', 1. ) ! apply boundary conditions -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_nam.F90
r12563 r13286 58 58 !! 59 59 NAMELIST/namusr_def/ nn_isize, nn_jsize, nn_ksize, nn_perio 60 NAMELIST/nammpp/ jpni, jpnj, ln_nnogather, ln_listonly60 NAMELIST/nammpp/ jpni, jpnj, nn_hls, ln_nnogather, ln_listonly 61 61 !!---------------------------------------------------------------------- 62 62 ! … … 77 77 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nammpp in configuration namelist' ) 78 78 79 kpi = ( -nn_isize - 2*nn_hls ) * jpni + 2*nn_hls80 kpj = ( -nn_jsize - 2*nn_hls ) * jpnj + 2*nn_hls79 kpi = -nn_isize * jpni 80 kpj = -nn_jsize * jpnj 81 81 ELSE 82 82 kpi = nn_isize -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90
r12740 r13286 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 bellow 101 102 !!--------------------------------------------------------------------- 102 103 #if defined key_si3 … … 104 105 ! 105 106 ! define unique value on each point. z2d ranging from 0.05 to -0.05 106 DO_2D_11_11 107 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( nimpp + ji - 1 + ( njmpp + jj - 2 ) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) 107 ! 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_00_00 112 !!$ z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 113 !!$ END_2D 114 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 DO_2D_00_00 117 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0_oldcmp(ji) + (mjg0_oldcmp(jj)-1) * igloi, wp ) / REAL( igloi * igloj, wp ) ) 108 118 END_2D 109 utau_ice(:,:) = 0.1_wp + 110 vtau_ice(:,:) = 0.1_wp + 119 utau_ice(:,:) = 0.1_wp + z2d(:,:) 120 vtau_ice(:,:) = 0.1_wp + z2d(:,:) 111 121 112 122 CALL lbc_lnk_multi( 'usrdef_sbc', utau_ice, 'U', -1., vtau_ice, 'V', -1. ) -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_zgr.F90
r12377 r13286 192 192 z2d(:,:) = REAL( jpkm1 , wp ) ! flat bottom 193 193 ! 194 IF( jperio == 3 .OR. jperio ==4 ) THEN ! add a small island in the upper corners to avoid model instabilities... 195 z2d(mi0( 1):mi1( 3),mj0(jpjglo-2):mj1(jpjglo)) = 0. 196 z2d(mi0(jpiglo-2):mi1(jpiglo),mj0(jpjglo-2):mj1(jpjglo)) = 0. 197 ENDIF 194 ! 195 ! BENCH should work without these 2 small islands on the 2 poles of the folding... 196 ! -> Comment out these lines if instabilities are too large... 197 ! 198 199 !!$ IF( jperio == 3 .OR. jperio == 4 ) THEN ! add a small island in the upper corners to avoid model instabilities... 200 !!$ z2d(mi0( nn_hls):mi1( nn_hls+2 ),mj0(jpjglo-nn_hls-1):mj1(jpjglo-nn_hls+1)) = 0. 201 !!$ z2d(mi0(jpiglo-nn_hls):mi1(MIN(jpiglo,jpiglo-nn_hls+2)),mj0(jpjglo-nn_hls-1):mj1(jpjglo-nn_hls+1)) = 0. 202 !!$ z2d(mi0(jpiglo/2 ):mi1( jpiglo/2 +2 ),mj0(jpjglo-nn_hls-1):mj1(jpjglo-nn_hls+1)) = 0. 203 !!$ ENDIF 204 !!$ ! 205 !!$ IF( jperio == 5 .OR. jperio == 6 ) THEN ! add a small island in the upper corners to avoid model instabilities... 206 !!$ z2d(mi0( nn_hls):mi1( nn_hls+1),mj0(jpjglo-nn_hls):mj1(jpjglo-nn_hls+1)) = 0. 207 !!$ z2d(mi0(jpiglo-nn_hls):mi1(jpiglo-nn_hls+1),mj0(jpjglo-nn_hls):mj1(jpjglo-nn_hls+1)) = 0. 208 !!$ z2d(mi0(jpiglo/2 ):mi1(jpiglo/2 +1),mj0(jpjglo-nn_hls):mj1(jpjglo-nn_hls+1)) = 0. 209 !!$ ENDIF 210 198 211 ! 199 212 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) -
NEMO/trunk/tests/CANAL/MY_SRC/domvvl.F90
r12740 r13286 450 450 ELSE 451 451 ijk_max = MAXLOC( ze3t(:,:,:) ) 452 ijk_max(1) = ijk_max(1) + nimpp - 1453 ijk_max(2) = ijk_max(2) + njmpp - 1452 ijk_max(1) = mig0_oldcmp(ijk_max(1)) 453 ijk_max(2) = mjg0_oldcmp(ijk_max(2)) 454 454 ijk_min = MINLOC( ze3t(:,:,:) ) 455 ijk_min(1) = ijk_min(1) + nimpp - 1456 ijk_min(2) = ijk_min(2) + njmpp - 1455 ijk_min(1) = mig0_oldcmp(ijk_min(1)) 456 ijk_min(2) = mjg0_oldcmp(ijk_min(2)) 457 457 ENDIF 458 458 IF (lwp) THEN … … 793 793 IF( ln_rstart ) THEN !* Read the restart file 794 794 CALL rst_read_open ! open the restart file if necessary 795 CALL iom_get( numror, jpdom_auto glo, 'sshn' , ssh(:,:,Kmm), ldxios = lrxios )795 CALL iom_get( numror, jpdom_auto, 'sshn' , ssh(:,:,Kmm), ldxios = lrxios ) 796 796 ! 797 797 id1 = iom_varid( numror, 'e3t_b', ldstop = .FALSE. ) … … 806 806 ! 807 807 IF( MIN( id1, id2 ) > 0 ) THEN ! all required arrays exist 808 CALL iom_get( numror, jpdom_auto glo, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios )809 CALL iom_get( numror, jpdom_auto glo, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios )808 CALL iom_get( numror, jpdom_auto, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios ) 809 CALL iom_get( numror, jpdom_auto, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios ) 810 810 ! needed to restart if land processor not computed 811 811 IF(lwp) write(numout,*) 'dom_vvl_rst : e3t(:,:,:,Kbb) and e3t(:,:,:,Kmm) found in restart files' … … 821 821 IF(lwp) write(numout,*) 'e3t_n set equal to e3t_b.' 822 822 IF(lwp) write(numout,*) 'l_1st_euler is forced to true' 823 CALL iom_get( numror, jpdom_auto glo, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios )823 CALL iom_get( numror, jpdom_auto, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios ) 824 824 e3t(:,:,:,Kmm) = e3t(:,:,:,Kbb) 825 825 l_1st_euler = .true. … … 828 828 IF(lwp) write(numout,*) 'e3t_b set equal to e3t_n.' 829 829 IF(lwp) write(numout,*) 'l_1st_euler is forced to true' 830 CALL iom_get( numror, jpdom_auto glo, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios )830 CALL iom_get( numror, jpdom_auto, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios ) 831 831 e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 832 832 l_1st_euler = .true. … … 853 853 ! ! ----------------------- ! 854 854 IF( MIN( id3, id4 ) > 0 ) THEN ! all required arrays exist 855 CALL iom_get( numror, jpdom_auto glo, 'tilde_e3t_b', tilde_e3t_b(:,:,:), ldxios = lrxios )856 CALL iom_get( numror, jpdom_auto glo, 'tilde_e3t_n', tilde_e3t_n(:,:,:), ldxios = lrxios )855 CALL iom_get( numror, jpdom_auto, 'tilde_e3t_b', tilde_e3t_b(:,:,:), ldxios = lrxios ) 856 CALL iom_get( numror, jpdom_auto, 'tilde_e3t_n', tilde_e3t_n(:,:,:), ldxios = lrxios ) 857 857 ELSE ! one at least array is missing 858 858 tilde_e3t_b(:,:,:) = 0.0_wp … … 863 863 ! ! ------------ ! 864 864 IF( id5 > 0 ) THEN ! required array exists 865 CALL iom_get( numror, jpdom_auto glo, 'hdiv_lf', hdiv_lf(:,:,:), ldxios = lrxios )865 CALL iom_get( numror, jpdom_auto, 'hdiv_lf', hdiv_lf(:,:,:), ldxios = lrxios ) 866 866 ELSE ! array is missing 867 867 hdiv_lf(:,:,:) = 0.0_wp -
NEMO/trunk/tests/CANAL/MY_SRC/usrdef_hgr.F90
r12740 r13286 63 63 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 64 64 ! 65 INTEGER :: ji, jj ! dummy loop indices65 INTEGER :: ji, jj ! dummy loop indices 66 66 REAL(wp) :: zphi0, zlam0, zbeta, zf0 67 REAL(wp) :: zti, z ui, ztj, zvj ! local scalars67 REAL(wp) :: zti, ztj ! local scalars 68 68 !!------------------------------------------------------------------------------- 69 69 ! … … 77 77 ! Position coordinates (in kilometers) 78 78 ! ========== 79 zlam0 = -REAL(NINT( jpiglo*rn_0xratio)-1, wp) * rn_dx80 zphi0 = -REAL(NINT( jpjglo*rn_0yratio)-1, wp) * rn_dy79 zlam0 = -REAL(NINT(Ni0glo*rn_0xratio)-1, wp) * rn_dx 80 zphi0 = -REAL(NINT(Nj0glo*rn_0yratio)-1, wp) * rn_dy 81 81 82 82 #if defined key_agrif … … 90 90 #endif 91 91 92 DO_2D_11_11 93 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 )94 z ui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp92 DO_2D_11_11 93 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 94 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 95 95 96 plamt(ji,jj) = zlam0 + rn_dx * zti97 plamu(ji,jj) = zlam0 + rn_dx * zui96 plamt(ji,jj) = zlam0 + rn_dx * zti 97 plamu(ji,jj) = zlam0 + rn_dx * ( zti + 0.5_wp ) 98 98 plamv(ji,jj) = plamt(ji,jj) 99 99 plamf(ji,jj) = plamu(ji,jj) 100 100 101 pphit(ji,jj) = zphi0 + rn_dy * ztj102 pphiv(ji,jj) = zphi0 + rn_dy * zvj101 pphit(ji,jj) = zphi0 + rn_dy * ztj 102 pphiv(ji,jj) = zphi0 + rn_dy * ( ztj + 0.5_wp ) 103 103 pphiu(ji,jj) = pphit(ji,jj) 104 104 pphif(ji,jj) = pphiv(ji,jj) -
NEMO/trunk/tests/CANAL/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain16 USE dom_oce 17 17 USE par_oce ! ocean space and time domain 18 18 USE phycst ! physical constants … … 106 106 kk_cfg = INT( rn_dx ) 107 107 ! 108 ! Global Domain size: EW_CANAL global domain is 1800 km x 1800 Km x 5000 m 109 kpi = NINT( rn_domszx / rn_dx ) + 1 110 kpj = NINT( rn_domszy / rn_dy ) + 3 111 kpk = NINT( rn_domszz / rn_dz ) + 1 112 #if defined key_agrif 113 IF( .NOT. Agrif_Root() ) THEN 114 kpi = nbcellsx + 2 + 2*nbghostcells 115 kpj = nbcellsy + 2 + 2*nbghostcells 108 IF( Agrif_Root() ) THEN ! Global Domain size: EW_CANAL global domain is 1800 km x 1800 Km x 5000 m 109 kpi = NINT( rn_domszx / rn_dx ) + 1 110 kpj = NINT( rn_domszy / rn_dy ) + 3 111 ELSE ! Global Domain size: add nbghostcells + 1 "land" point on each side 112 kpi = nbcellsx + nbghostcells_x + nbghostcells_x + 2 113 kpj = nbcellsy + nbghostcells_y_s + nbghostcells_y_n + 2 116 114 ENDIF 117 #endif 115 kpk = MAX( 2, NINT( rn_domszz / rn_dz ) + 1 ) 118 116 ! 119 117 zh = (kpk-1)*rn_dz -
NEMO/trunk/tests/ICE_ADV1D/MY_SRC/usrdef_hgr.F90
r12740 r13286 64 64 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 65 65 ! 66 INTEGER :: ji, jj ! dummy loop indices66 INTEGER :: ji, jj ! dummy loop indices 67 67 REAL(wp) :: zphi0, zlam0, zbeta, zf0 68 REAL(wp) :: zti, z ui, ztj, zvj ! local scalars68 REAL(wp) :: zti, ztj ! local scalars 69 69 !!------------------------------------------------------------------------------- 70 70 ! … … 75 75 76 76 ! ========== 77 zlam0 = - (jpiglo-1)/2* 1.e-3 * rn_dx78 zphi0 = - (jpjglo-1)/2* 1.e-3 * rn_dy77 zlam0 = -REAL( (Ni0glo-2)/2, wp) * 1.e-3 * rn_dx 78 zphi0 = -REAL( (Nj0glo-2)/2, wp) * 1.e-3 * rn_dy 79 79 80 80 DO_2D_11_11 81 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 )82 z ui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp81 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 82 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 83 83 84 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti85 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui84 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti 85 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * ( zti + 0.5_wp ) 86 86 plamv(ji,jj) = plamt(ji,jj) 87 87 plamf(ji,jj) = plamu(ji,jj) 88 88 89 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj90 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj89 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj 90 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * ( ztj + 0.5_wp ) 91 91 pphiu(ji,jj) = pphit(ji,jj) 92 92 pphif(ji,jj) = pphiv(ji,jj) -
NEMO/trunk/tests/ICE_ADV1D/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain17 16 USE par_oce ! ocean space and time domain 18 17 USE phycst ! physical constants … … 91 90 WRITE(numout,*) ' LX [km]: ', zlx 92 91 WRITE(numout,*) ' LY [km]: ', zly 93 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi94 WRITE(numout,*) ' jpjglo = ', kpj92 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 93 WRITE(numout,*) ' Nj0glo = ', kpj 95 94 WRITE(numout,*) ' jpkglo = ', kpk 96 95 WRITE(numout,*) ' Coriolis:', ln_corio -
NEMO/trunk/tests/ICE_ADV2D/MY_SRC/usrdef_hgr.F90
r12740 r13286 64 64 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 65 65 ! 66 INTEGER :: ji, jj ! dummy loop indices66 INTEGER :: ji, jj ! dummy loop indices 67 67 REAL(wp) :: zphi0, zlam0, zbeta, zf0 68 REAL(wp) :: zti, z ui, ztj, zvj ! local scalars68 REAL(wp) :: zti, ztj ! local scalars 69 69 !!------------------------------------------------------------------------------- 70 70 ! … … 76 76 77 77 ! ========== 78 zlam0 = - (jpiglo-1)/2* 1.e-3 * rn_dx79 zphi0 = - (jpjglo-1)/2 * 1.e-3 * rn_dy78 zlam0 = -REAL( (Ni0glo-2)/2, wp) * 1.e-3 * rn_dx 79 zphi0 = -REAL( (Nj0glo-2)/2, wp) * 1.e-3 * rn_dy 80 80 81 81 #if defined key_agrif … … 83 83 !clem zlam0 = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 84 84 !clem zphi0 = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 85 zlam0 = ( 0.5_wp - ( Agrif_parent(jpiglo) - 1 ) / 2) * 1.e-3 * Agrif_irhox() * rn_dx &85 zlam0 = ( 0.5_wp - REAL( (Agrif_parent(Ni0glo) - 2 ) / 2, wp ) ) * 1.e-3 * Agrif_irhox() * rn_dx & 86 86 & + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 87 zphi0 = ( 0.5_wp - ( Agrif_parent(jpjglo) - 1 ) / 2) * 1.e-3 * Agrif_irhoy() * rn_dy &87 zphi0 = ( 0.5_wp - REAL( (Agrif_parent(Nj0glo) - 2 ) / 2, wp ) ) * 1.e-3 * Agrif_irhoy() * rn_dy & 88 88 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 89 89 ENDIF … … 91 91 92 92 DO_2D_11_11 93 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 )94 z ui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp93 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 94 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 95 95 96 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti97 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui96 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti 97 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * ( zti + 0.5_wp ) 98 98 plamv(ji,jj) = plamt(ji,jj) 99 99 plamf(ji,jj) = plamu(ji,jj) 100 100 101 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj102 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj101 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj 102 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * ( ztj + 0.5_wp ) 103 103 pphiu(ji,jj) = pphit(ji,jj) 104 104 pphif(ji,jj) = pphiv(ji,jj) -
NEMO/trunk/tests/ICE_ADV2D/MY_SRC/usrdef_nam.F90
r12740 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp, Agrif_Root ! i- & j-indices of the local domain16 USE dom_oce 17 17 USE par_oce ! ocean space and time domain 18 18 USE phycst ! physical constants … … 82 82 kk_cfg = NINT( rn_dx ) 83 83 ! 84 IF( Agrif_Root() ) THEN ! Global Domain size: 84 IF( Agrif_Root() ) THEN ! Global Domain size: ICE_AGRIF domain is 300 km x 300 Km x 10 m 85 85 kpi = NINT( 300.e3 / rn_dx ) - 1 86 86 kpj = NINT( 300.e3 / rn_dy ) - 1 87 ELSE 88 kpi = nbcellsx + 2 + 2*nbghostcells89 kpj = nbcellsy + 2 + 2*nbghostcells87 ELSE ! Global Domain size: add nbghostcells + 1 "land" point on each side 88 kpi = nbcellsx + nbghostcells_x + nbghostcells_x + 2 89 kpj = nbcellsy + nbghostcells_y_s + nbghostcells_y_n + 2 90 90 ENDIF 91 kpk = 191 kpk = 2 92 92 ! 93 93 !! zlx = (kpi-2)*rn_dx*1.e-3 … … 110 110 WRITE(numout,*) ' LX [km]: ', zlx 111 111 WRITE(numout,*) ' LY [km]: ', zly 112 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi113 WRITE(numout,*) ' jpjglo = ', kpj112 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 113 WRITE(numout,*) ' Nj0glo = ', kpj 114 114 WRITE(numout,*) ' jpkglo = ', kpk 115 115 WRITE(numout,*) ' Coriolis:', ln_corio -
NEMO/trunk/tests/ICE_AGRIF/EXPREF/AGRIF_FixedGrids.in
r9159 r13286 1 1 1 2 3 4 63 34 633 3 32 33 62 33 62 3 3 3 3 3 0 -
NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90
r12740 r13286 64 64 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 65 65 ! 66 INTEGER :: ji, jj ! dummy loop indices66 INTEGER :: ji, jj ! dummy loop indices 67 67 REAL(wp) :: zphi0, zlam0, zbeta, zf0 68 REAL(wp) :: zti, z ui, ztj, zvj ! local scalars68 REAL(wp) :: zti, ztj ! local scalars 69 69 !!------------------------------------------------------------------------------- 70 70 ! … … 76 76 77 77 ! ========== 78 zlam0 = -(jpiglo-1)/2 * 1.e-3 * rn_dx79 zphi0 = -(jpjglo-1)/2 * 1.e-3 * rn_dy80 81 78 #if defined key_agrif 82 IF( .NOT. Agrif_Root() ) THEN 79 IF( Agrif_Root() ) THEN 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 86 #if defined key_agrif 87 ELSE 88 ! ! let lower left longitude and latitude from parent 83 89 !clem zlam0 = Agrif_Parent(zlam0) + (Agrif_ix())*Agrif_Parent(rn_dx) * 1.e-5 84 90 !clem zphi0 = Agrif_Parent(zphi0) + (Agrif_iy())*Agrif_Parent(rn_dy) * 1.e-5 85 zlam0 = ( 0.5_wp - ( Agrif_parent(jpiglo) - 1 ) / 2 ) * 1.e-3 * Agrif_irhox() * rn_dx & 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 & 86 95 & + ( Agrif_Ix() + nbghostcells - 1 ) * Agrif_irhox() * rn_dx * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dx * 1.e-3 87 zphi0 = ( 0.5_wp - ( Agrif_parent(jpjglo) - 1 ) / 2) * 1.e-3 * Agrif_irhoy() * rn_dy &96 zphi0 = ( 0.5_wp - REAL( ( Agrif_parent(Nj0glo)-1 ) / 2, wp) ) * 1.e-3 * Agrif_irhoy() * rn_dy & 88 97 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 89 98 ENDIF … … 91 100 92 101 DO_2D_11_11 93 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 )94 z ui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp102 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 103 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 95 104 96 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti97 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui105 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti 106 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * ( zti + 0.5_wp ) 98 107 plamv(ji,jj) = plamt(ji,jj) 99 108 plamf(ji,jj) = plamu(ji,jj) 100 109 101 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj102 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj110 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj 111 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * ( ztj + 0.5_wp ) 103 112 pphiu(ji,jj) = pphit(ji,jj) 104 113 pphif(ji,jj) = pphiv(ji,jj) -
NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90
r12597 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain16 USE dom_oce 17 17 USE par_oce ! ocean space and time domain 18 18 USE phycst ! physical constants … … 85 85 kpi = NINT( 300.e3 / rn_dx ) - 1 86 86 kpj = NINT( 300.e3 / rn_dy ) - 1 87 ELSE 88 kpi = nbcellsx + 2 + 2*nbghostcells 89 kpj = nbcellsy + 2 + 2*nbghostcells 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 89 ELSE ! Global Domain size: add nbghostcells + 1 "land" point on each side 90 kpi = nbcellsx + 2 * ( nbghostcells + 1 ) 91 kpj = nbcellsy + 2 * ( nbghostcells + 1 ) 92 !!$ kpi = nbcellsx + nbghostcells_x + nbghostcells_x + 2 93 !!$ kpj = nbcellsy + nbghostcells_y_s + nbghostcells_y_n + 2 90 94 ENDIF 91 95 kpk = 2 … … 110 114 WRITE(numout,*) ' LX [km]: ', zlx 111 115 WRITE(numout,*) ' LY [km]: ', zly 112 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi113 WRITE(numout,*) ' jpjglo = ', kpj116 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 117 WRITE(numout,*) ' Nj0glo = ', kpj 114 118 WRITE(numout,*) ' jpkglo = ', kpk 115 119 WRITE(numout,*) ' Coriolis:', ln_corio -
NEMO/trunk/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r12905 r13286 211 211 erp(:,:) = erp(:,:) + zerp_cor(:,:) 212 212 ! 213 IF( nprint == 1 .AND.lwp ) THEN ! control print213 IF( lwp ) THEN ! control print 214 214 IF( z_fwf < 0._wp ) THEN 215 215 WRITE(numout,*)' z_fwf < 0' -
NEMO/trunk/tests/ISOMIP+/MY_SRC/tradmp.F90
r12905 r13286 208 208 ! ! Read in mask from file 209 209 CALL iom_open ( cn_resto, imask) 210 CALL iom_get ( imask, jpdom_auto glo, 'resto', resto )210 CALL iom_get ( imask, jpdom_auto, 'resto', resto ) 211 211 CALL iom_close( imask ) 212 212 ENDIF -
NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_hgr.F90
r12740 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh for ISOMIP configuration 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain16 USE dom_oce 17 17 USE par_oce ! ocean space and time domain 18 18 USE phycst ! physical constants … … 79 79 DO_2D_11_11 80 80 ! ! longitude (west coast at lon=0°) 81 plamt(ji,jj) = rn_e1deg * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) )82 plamu(ji,jj) = rn_e1deg * ( REAL( ji-1 + nimpp-1 , wp ) )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 ) ) 83 83 plamv(ji,jj) = plamt(ji,jj) 84 84 plamf(ji,jj) = plamu(ji,jj) 85 85 ! ! latitude (south coast at lat= 81°) 86 pphit(ji,jj) = rn_e2deg * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) ) - 80._wp86 pphit(ji,jj) = rn_e2deg * ( - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp ) ) - 80._wp 87 87 pphiu(ji,jj) = pphit(ji,jj) 88 pphiv(ji,jj) = rn_e2deg * ( REAL( jj-1 + njmpp-1 , wp ) ) - 80_wp88 pphiv(ji,jj) = rn_e2deg * ( REAL( mjg0_oldcmp(jj)-1 , wp ) ) - 80_wp 89 89 pphif(ji,jj) = pphiv(ji,jj) 90 90 END_2D -
NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_nam.F90
r12377 r13286 15 15 !! usr_def_hgr : initialize the horizontal mesh 16 16 !!---------------------------------------------------------------------- 17 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain18 17 USE dom_oce , ONLY: ln_zco, ln_zps, ln_sco ! flag of type of coordinate 19 18 USE par_oce ! ocean space and time domain … … 95 94 WRITE(numout,*) ' vertical resolution rn_e3 = ', rn_e3 , ' meters' 96 95 WRITE(numout,*) ' ISOMIP domain = 15° x 10° x 900 m' 97 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi98 WRITE(numout,*) ' jpjglo = ', kpj96 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 97 WRITE(numout,*) ' Nj0glo = ', kpj 99 98 WRITE(numout,*) ' jpkglo = ', kpk 100 99 WRITE(numout,*) ' ' -
NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_zgr.F90
r12740 r13286 16 16 !!--------------------------------------------------------------------- 17 17 USE oce ! ocean variables 18 USE dom_oce , ONLY: mj0 , mj1 , nimpp , njmpp! ocean space and time domain19 USE dom_oce , ONLY: glamt , gphit 18 USE dom_oce , ONLY: mj0 , mj1 ! ocean space and time domain 19 USE dom_oce , ONLY: glamt , gphit ! ocean space and time domain 20 20 USE usrdef_nam ! User defined : namelist variables 21 21 ! … … 67 67 REAL(wp), DIMENSION(jpi,jpj) :: zht , zhu ! bottom depth 68 68 REAL(wp), DIMENSION(jpi,jpj) :: zhisf, zhisfu ! top depth 69 REAL(wp), DIMENSION(jpi,jpj) :: zmsk70 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2d workspace71 69 !!---------------------------------------------------------------------- 72 70 ! … … 87 85 ! !== isfdraft ==! 88 86 ! 89 ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=090 z2d(:,:) = 1._wp ! surface ocean is the 1st level91 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)92 zmsk(:,:) = NINT( z2d(:,:) )93 !94 !95 87 zht (:,:) = rbathy 96 88 zhisf(:,:) = 200._wp 97 ij0 = 1 ; ij1 = 4089 ij0 = 1 ; ij1 = 40+nn_hls 98 90 DO jj = mj0(ij0), mj1(ij1) 99 91 zhisf(:,jj)=700.0_wp-(gphit(:,jj)+80.0_wp)*125.0_wp 100 92 END DO 101 zhisf(:,:) = zhisf(:,:) * zmsk(:,:)102 93 ! 103 94 CALL zgr_z1d( pdept_1d, pdepw_1d, pe3t_1d , pe3w_1d ) ! Reference z-coordinate system -
NEMO/trunk/tests/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90
r12740 r13286 13 13 !! usr_def_hgr : initialize the horizontal mesh for LOCK_EXCHANGE configuration 14 14 !!---------------------------------------------------------------------- 15 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain15 USE dom_oce 16 16 USE par_oce ! ocean space and time domain 17 17 USE phycst ! physical constants … … 76 76 DO_2D_11_11 77 77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) )78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0_oldcmp(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) )83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0_oldcmp(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D -
NEMO/trunk/tests/LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain17 16 USE par_oce ! ocean space and time domain 18 17 USE phycst ! physical constants … … 85 84 WRITE(numout,*) ' vertical resolution rn_dz = ', rn_dz, ' meters' 86 85 WRITE(numout,*) ' LOCK_EXCHANGE domain = 64 km x 3 grid-points x 20 m' 87 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi88 WRITE(numout,*) ' jpjglo = ', kpj86 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 87 WRITE(numout,*) ' Nj0glo = ', kpj 89 88 WRITE(numout,*) ' jpkglo = ', kpk 90 89 WRITE(numout,*) ' ' -
NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_hgr.F90
r12740 r13286 13 13 !! usr_def_hgr : initialize the horizontal mesh for OVERFLOW configuration 14 14 !!---------------------------------------------------------------------- 15 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain15 USE dom_oce 16 16 USE par_oce ! ocean space and time domain 17 17 USE phycst ! physical constants … … 76 76 DO_2D_11_11 77 77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) )78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0_oldcmp(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) )83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0_oldcmp(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D -
NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain17 16 USE dom_oce , ONLY: ln_zco, ln_zps, ln_sco ! flag of type of coordinate 18 17 USE par_oce ! ocean space and time domain … … 86 85 WRITE(numout,*) ' vertical resolution rn_dz = ', rn_dz, ' meters' 87 86 WRITE(numout,*) ' OVERFLOW domain = 200 km x 3 grid-points x 2000 m' 88 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi89 WRITE(numout,*) ' jpjglo = ', kpj87 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 88 WRITE(numout,*) ' Nj0glo = ', kpj 90 89 WRITE(numout,*) ' jpkglo = ', kpk 91 90 ! -
NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_zgr.F90
r12740 r13286 15 15 !!--------------------------------------------------------------------- 16 16 USE oce ! ocean variables 17 USE dom_oce , ONLY: mi0, mi1 , nimpp, njmpp! ocean space and time domain18 USE dom_oce , ONLY: glamt 17 USE dom_oce , ONLY: mi0, mi1 ! ocean space and time domain 18 USE dom_oce , ONLY: glamt ! ocean space and time domain 19 19 USE usrdef_nam ! User defined : namelist variables 20 20 ! -
NEMO/trunk/tests/STATION_ASF/MY_SRC/nemogcm.F90
r13011 r13286 30 30 USE step_c1d ! Time stepping loop for the 1D configuration 31 31 ! 32 USE prtctl ! Print control 32 33 USE in_out_manager ! I/O manager 33 34 USE lib_mpp ! distributed memory computing … … 131 132 INTEGER :: ios, ilocal_comm ! local integers 132 133 !! 133 NAMELIST/namctl/ sn_cfctl, nn_print, nn_ictls, nn_ictle, & 134 & nn_isplt , nn_jsplt, nn_jctls, nn_jctle, & 135 & ln_timing, ln_diacfl 134 NAMELIST/namctl/ sn_cfctl, ln_timing, ln_diacfl, & 135 & nn_isplt, nn_jsplt, nn_ictls, nn_ictle, nn_jctls, nn_jctle 136 136 NAMELIST/namcfg/ ln_read_cfg, cn_domcfg, ln_closea, ln_write_cfg, cn_domcfg_out, ln_use_jattr 137 137 !!---------------------------------------------------------------------- … … 232 232 ! 233 233 IF( ln_read_cfg ) THEN ! Read sizes in domain configuration file 234 CALL domain_cfg ( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio )234 CALL domain_cfg ( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 235 235 ELSE ! user-defined namelist 236 CALL usr_def_nam( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio )236 CALL usr_def_nam( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 237 237 ENDIF 238 238 ! … … 306 306 WRITE(numout,*) ' sn_cfctl%procincr = ', sn_cfctl%procincr 307 307 WRITE(numout,*) ' sn_cfctl%ptimincr = ', sn_cfctl%ptimincr 308 WRITE(numout,*) ' level of print nn_print = ', nn_print309 WRITE(numout,*) ' Start i indice for SUM control nn_ictls = ', nn_ictls310 WRITE(numout,*) ' End i indice for SUM control nn_ictle = ', nn_ictle311 WRITE(numout,*) ' Start j indice for SUM control nn_jctls = ', nn_jctls312 WRITE(numout,*) ' End j indice for SUM control nn_jctle = ', nn_jctle313 WRITE(numout,*) ' number of proc. following i nn_isplt = ', nn_isplt314 WRITE(numout,*) ' number of proc. following j nn_jsplt = ', nn_jsplt315 308 WRITE(numout,*) ' timing by routine ln_timing = ', ln_timing 316 309 WRITE(numout,*) ' CFL diagnostics ln_diacfl = ', ln_diacfl 317 310 ENDIF 318 311 ! 319 nprint = nn_print ! convert DOCTOR namelist names into OLD names 320 nictls = nn_ictls 321 nictle = nn_ictle 322 njctls = nn_jctls 323 njctle = nn_jctle 324 isplt = nn_isplt 325 jsplt = nn_jsplt 326 312 IF( .NOT.ln_read_cfg ) ln_closea = .false. ! dealing possible only with a domcfg file 327 313 IF(lwp) THEN ! control print 328 314 WRITE(numout,*) … … 335 321 WRITE(numout,*) ' use file attribute if exists as i/p j-start ln_use_jattr = ', ln_use_jattr 336 322 ENDIF 337 IF( .NOT.ln_read_cfg ) ln_closea = .false. ! dealing possible only with a domcfg file338 !339 ! ! Parameter control340 !341 IF( sn_cfctl%l_prtctl .OR. sn_cfctl%l_prttrc ) THEN ! sub-domain area indices for the control prints342 IF( lk_mpp .AND. jpnij > 1 ) THEN343 isplt = jpni ; jsplt = jpnj ; ijsplt = jpni*jpnj ! the domain is forced to the real split domain344 ELSE345 IF( isplt == 1 .AND. jsplt == 1 ) THEN346 CALL ctl_warn( ' - isplt & jsplt are equal to 1', &347 & ' - the print control will be done over the whole domain' )348 ENDIF349 ijsplt = isplt * jsplt ! total number of processors ijsplt350 ENDIF351 IF(lwp) WRITE(numout,*)' - The total number of processors over which the'352 IF(lwp) WRITE(numout,*)' print control will be done is ijsplt : ', ijsplt353 !354 ! ! indices used for the SUM control355 IF( nictls+nictle+njctls+njctle == 0 ) THEN ! print control done over the default area356 lsp_area = .FALSE.357 ELSE ! print control done over a specific area358 lsp_area = .TRUE.359 IF( nictls < 1 .OR. nictls > jpiglo ) THEN360 CALL ctl_warn( ' - nictls must be 1<=nictls>=jpiglo, it is forced to 1' )361 nictls = 1362 ENDIF363 IF( nictle < 1 .OR. nictle > jpiglo ) THEN364 CALL ctl_warn( ' - nictle must be 1<=nictle>=jpiglo, it is forced to jpiglo' )365 nictle = jpiglo366 ENDIF367 IF( njctls < 1 .OR. njctls > jpjglo ) THEN368 CALL ctl_warn( ' - njctls must be 1<=njctls>=jpjglo, it is forced to 1' )369 njctls = 1370 ENDIF371 IF( njctle < 1 .OR. njctle > jpjglo ) THEN372 CALL ctl_warn( ' - njctle must be 1<=njctle>=jpjglo, it is forced to jpjglo' )373 njctle = jpjglo374 ENDIF375 ENDIF376 ENDIF377 323 ! 378 324 IF( 1._wp /= SIGN(1._wp,-0._wp) ) CALL ctl_stop( 'nemo_ctl: The intrinsec SIGN function follows f2003 standard.', & -
NEMO/trunk/tests/STATION_ASF/MY_SRC/usrdef_hgr.F90
r12629 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain17 16 USE c1d , ONLY: rn_lon1d, rn_lat1d ! ocean lon/lat define by namelist 18 17 USE par_oce ! ocean space and time domain -
NEMO/trunk/tests/STATION_ASF/MY_SRC/usrdef_nam.F90
r12629 r13286 15 15 !! usr_def_hgr : initialize the horizontal mesh 16 16 !!---------------------------------------------------------------------- 17 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain18 !!! USE dom_oce , ONLY: ln_zco, ln_zps, ln_sco ! flag of type of coordinate19 17 USE par_oce ! ocean space and time domain 20 18 USE phycst ! physical constants -
NEMO/trunk/tests/VORTEX/MY_SRC/domvvl.F90
r12740 r13286 450 450 ELSE 451 451 ijk_max = MAXLOC( ze3t(:,:,:) ) 452 ijk_max(1) = ijk_max(1) + nimpp - 1453 ijk_max(2) = ijk_max(2) + njmpp - 1452 ijk_max(1) = mig0_oldcmp(ijk_max(1)) 453 ijk_max(2) = mjg0_oldcmp(ijk_max(2)) 454 454 ijk_min = MINLOC( ze3t(:,:,:) ) 455 ijk_min(1) = ijk_min(1) + nimpp - 1456 ijk_min(2) = ijk_min(2) + njmpp - 1455 ijk_min(1) = mig0_oldcmp(ijk_min(1)) 456 ijk_min(2) = mjg0_oldcmp(ijk_min(2)) 457 457 ENDIF 458 458 IF (lwp) THEN … … 793 793 IF( ln_rstart ) THEN !* Read the restart file 794 794 CALL rst_read_open ! open the restart file if necessary 795 CALL iom_get( numror, jpdom_auto glo, 'sshn' , ssh(:,:,Kmm), ldxios = lrxios )795 CALL iom_get( numror, jpdom_auto, 'sshn' , ssh(:,:,Kmm), ldxios = lrxios ) 796 796 ! 797 797 id1 = iom_varid( numror, 'e3t_b', ldstop = .FALSE. ) … … 806 806 ! 807 807 IF( MIN( id1, id2 ) > 0 ) THEN ! all required arrays exist 808 CALL iom_get( numror, jpdom_auto glo, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios )809 CALL iom_get( numror, jpdom_auto glo, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios )808 CALL iom_get( numror, jpdom_auto, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios ) 809 CALL iom_get( numror, jpdom_auto, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios ) 810 810 ! needed to restart if land processor not computed 811 811 IF(lwp) write(numout,*) 'dom_vvl_rst : e3t(:,:,:,Kbb) and e3t(:,:,:,Kmm) found in restart files' … … 821 821 IF(lwp) write(numout,*) 'e3t_n set equal to e3t_b.' 822 822 IF(lwp) write(numout,*) 'l_1st_euler is forced to true' 823 CALL iom_get( numror, jpdom_auto glo, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios )823 CALL iom_get( numror, jpdom_auto, 'e3t_b', e3t(:,:,:,Kbb), ldxios = lrxios ) 824 824 e3t(:,:,:,Kmm) = e3t(:,:,:,Kbb) 825 825 l_1st_euler = .true. … … 828 828 IF(lwp) write(numout,*) 'e3t_b set equal to e3t_n.' 829 829 IF(lwp) write(numout,*) 'l_1st_euler is forced to true' 830 CALL iom_get( numror, jpdom_auto glo, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios )830 CALL iom_get( numror, jpdom_auto, 'e3t_n', e3t(:,:,:,Kmm), ldxios = lrxios ) 831 831 e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 832 832 l_1st_euler = .true. … … 853 853 ! ! ----------------------- ! 854 854 IF( MIN( id3, id4 ) > 0 ) THEN ! all required arrays exist 855 CALL iom_get( numror, jpdom_auto glo, 'tilde_e3t_b', tilde_e3t_b(:,:,:), ldxios = lrxios )856 CALL iom_get( numror, jpdom_auto glo, 'tilde_e3t_n', tilde_e3t_n(:,:,:), ldxios = lrxios )855 CALL iom_get( numror, jpdom_auto, 'tilde_e3t_b', tilde_e3t_b(:,:,:), ldxios = lrxios ) 856 CALL iom_get( numror, jpdom_auto, 'tilde_e3t_n', tilde_e3t_n(:,:,:), ldxios = lrxios ) 857 857 ELSE ! one at least array is missing 858 858 tilde_e3t_b(:,:,:) = 0.0_wp … … 863 863 ! ! ------------ ! 864 864 IF( id5 > 0 ) THEN ! required array exists 865 CALL iom_get( numror, jpdom_auto glo, 'hdiv_lf', hdiv_lf(:,:,:), ldxios = lrxios )865 CALL iom_get( numror, jpdom_auto, 'hdiv_lf', hdiv_lf(:,:,:), ldxios = lrxios ) 866 866 ELSE ! array is missing 867 867 hdiv_lf(:,:,:) = 0.0_wp -
NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_hgr.F90
r12740 r13286 63 63 REAL(wp), DIMENSION(:,:), INTENT(out) :: pe1e2u, pe1e2v ! u- & v-surfaces (if reduction in strait) [m2] 64 64 ! 65 INTEGER :: ji, jj ! dummy loop indices65 INTEGER :: ji, jj ! dummy loop indices 66 66 REAL(wp) :: zphi0, zlam0, zbeta, zf0 67 REAL(wp) :: zti, z ui, ztj, zvj ! local scalars67 REAL(wp) :: zti, ztj ! local scalars 68 68 !!------------------------------------------------------------------------------- 69 69 ! … … 77 77 ! Position coordinates (in kilometers) 78 78 ! ========== 79 zlam0 = -(jpiglo-1)/2 * 1.e-3 * rn_dx 80 zphi0 = -(jpjglo-1)/2 * 1.e-3 * rn_dy 81 79 #if defined key_agrif 80 IF( Agrif_Root() ) THEN 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 87 #if defined key_agrif 83 ! ! let lower left longitude and latitude from parent 84 IF (.NOT.Agrif_root()) THEN 85 zlam0 = (0.5_wp-(Agrif_parent(jpiglo)-1)/2)*1.e-3*Agrif_irhox()*rn_dx & 86 &+(Agrif_Ix()+nbghostcells-1)*Agrif_irhox()*rn_dx*1.e-3-(0.5_wp+nbghostcells)*rn_dx*1.e-3 87 zphi0 = (0.5_wp-(Agrif_parent(jpjglo)-1)/2)*1.e-3*Agrif_irhoy()*rn_dy & 88 &+(Agrif_Iy()+nbghostcells-1)*Agrif_irhoy()*rn_dy*1.e-3-(0.5_wp+nbghostcells)*rn_dy*1.e-3 88 ELSE 89 ! ! 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 & 94 & + ( 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 & 96 & + ( Agrif_Iy() + nbghostcells - 1 ) * Agrif_irhoy() * rn_dy * 1.e-3 - ( 0.5_wp + nbghostcells ) * rn_dy * 1.e-3 89 97 ENDIF 90 98 #endif 91 99 92 100 DO_2D_11_11 93 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 )94 z ui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp101 zti = REAL( mig0_oldcmp(ji) - 1, wp ) ! start at i=0 in the global grid without halos 102 ztj = REAL( mjg0_oldcmp(jj) - 1, wp ) ! start at j=0 in the global grid without halos 95 103 96 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti97 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui104 plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti 105 plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * ( zti + 0.5_wp ) 98 106 plamv(ji,jj) = plamt(ji,jj) 99 107 plamf(ji,jj) = plamu(ji,jj) 100 108 101 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj102 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj109 pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj 110 pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * ( ztj + 0.5_wp ) 103 111 pphiu(ji,jj) = pphit(ji,jj) 104 112 pphif(ji,jj) = pphiv(ji,jj) -
NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain16 USE dom_oce 17 17 USE par_oce ! ocean space and time domain 18 18 USE phycst ! physical constants … … 84 84 kpi = NINT( 1800.e3 / rn_dx ) + 3 85 85 kpj = NINT( 1800.e3 / rn_dy ) + 3 86 ELSE 87 kpi = nbcellsx + 2 + 2*nbghostcells 88 kpj = nbcellsy + 2 + 2*nbghostcells 86 ELSE ! Global Domain size: add nbghostcells + 1 "land" point on each side 87 kpi = nbcellsx + 2 * ( nbghostcells + 1 ) 88 kpj = nbcellsy + 2 * ( nbghostcells + 1 ) 89 !!$ kpi = nbcellsx + nbghostcells_x + nbghostcells_x + 2 90 !!$ kpj = nbcellsy + nbghostcells_y_s + nbghostcells_y_n + 2 89 91 ENDIF 90 92 kpk = NINT( 5000._wp / rn_dz ) + 1 … … 104 106 WRITE(numout,*) ' horizontal resolution rn_dy = ', rn_dy, ' m' 105 107 WRITE(numout,*) ' vertical resolution rn_dz = ', rn_dz, ' m' 108 WRITE(numout,*) ' resulting global domain size : Ni0glo = ', kpi 109 WRITE(numout,*) ' Nj0glo = ', kpj 110 WRITE(numout,*) ' jpkglo = ', kpk 106 111 WRITE(numout,*) ' VORTEX domain: ' 107 112 WRITE(numout,*) ' LX [km]: ', zlx -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_hgr.F90
r12740 r13286 13 13 !! usr_def_hgr : initialize the horizontal mesh for WAD_TEST_CASES configuration 14 14 !!---------------------------------------------------------------------- 15 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain15 USE dom_oce 16 16 USE par_oce ! ocean space and time domain 17 17 USE phycst ! physical constants … … 76 76 DO_2D_11_11 77 77 ! ! longitude 78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) )79 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) )78 plamt(ji,jj) = zfact * ( - 0.5 + REAL( mig0_oldcmp(ji)-1 , wp ) ) 79 plamu(ji,jj) = zfact * ( REAL( mig0_oldcmp(ji)-1 , wp ) ) 80 80 plamv(ji,jj) = plamt(ji,jj) 81 81 plamf(ji,jj) = plamu(ji,jj) 82 82 ! ! latitude 83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) )83 pphit(ji,jj) = zfact * ( - 0.5 + REAL( mjg0_oldcmp(jj)-1 , wp ) ) 84 84 pphiu(ji,jj) = pphit(ji,jj) 85 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) )85 pphiv(ji,jj) = zfact * ( REAL( mjg0_oldcmp(jj)-1 , wp ) ) 86 86 pphif(ji,jj) = pphiv(ji,jj) 87 87 END_2D -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_nam.F90
r12377 r13286 14 14 !! usr_def_hgr : initialize the horizontal mesh 15 15 !!---------------------------------------------------------------------- 16 USE dom_oce , ONLY: nimpp , njmpp ! i- & j-indices of the local domain17 16 USE par_oce ! ocean space and time domain 18 17 USE phycst ! physical constants … … 77 76 ! ! Set the lateral boundary condition of the global domain 78 77 kperio = 0 ! WAD_TEST_CASES configuration : closed domain 79 IF( nn_wad_test == 8 ) kperio = 7 ! North-South cyclic test 78 IF( nn_wad_test == 8 ) THEN 79 kperio = 7 ! North-South cyclic test 80 kpi = kpi - 2 ! no closed boundary 81 kpj = kpj - 2 ! no closed boundary 82 ENDIF 80 83 ! 81 84 ! ! control print -
NEMO/trunk/tests/WAD/MY_SRC/usrdef_zgr.F90
r12740 r13286 15 15 !!--------------------------------------------------------------------- 16 16 USE oce ! ocean variables 17 USE dom_oce , ONLY: ht_0, mi0, mi1, nimpp, njmpp, & 18 & mj0, mj1, glamt, gphit ! ocean space and time domain 17 USE dom_oce , ONLY: ht_0, mi0, mi1, mj0, mj1, glamt, gphit ! ocean space and time domain 19 18 USE usrdef_nam ! User defined : namelist variables 20 19 USE wet_dry , ONLY: rn_wdmin1, rn_wdmin2, rn_wdld ! Wetting and drying
Note: See TracChangeset
for help on using the changeset viewer.
