Changeset 6906
- Timestamp:
- 2016-09-02T15:06:17+02:00 (8 years ago)
- Location:
- branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG
- Files:
-
- 1 added
- 9 edited
-
LOCK_EXCHANGE/EXP00/namelist_cfg (modified) (12 diffs)
-
LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90 (modified) (4 diffs)
-
LOCK_EXCHANGE/MY_SRC/usrdef_istate.F90 (modified) (3 diffs)
-
LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90 (modified) (4 diffs)
-
LOCK_EXCHANGE/MY_SRC/usrdef_zgr.F90 (added)
-
LOCK_EXCHANGE/cpp_LOCK_EXCHANGE.fcm (modified) (1 diff)
-
OVERFLOW/EXP00/namelist_cfg (modified) (4 diffs)
-
OVERFLOW/MY_SRC/usrdef_istate.F90 (modified) (2 diffs)
-
OVERFLOW/MY_SRC/usrdef_nam.F90 (modified) (4 diffs)
-
OVERFLOW/MY_SRC/usrdef_zgr.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/EXP00/namelist_cfg
r6901 r6906 6 6 !----------------------------------------------------------------------- 7 7 ln_zco = .true. ! z-coordinate flag 8 ln_zps = .false. ! z-partial-step coordinate flag9 ln_sco = .false. ! s-coordinate flag10 8 rn_dx = 500. ! horizontal resolution in meters 11 9 rn_dz = 1. ! vertical resolution in meters … … 25 23 &namcfg ! parameters of the configuration 26 24 !----------------------------------------------------------------------- 27 ln_read_cfg = . true. ! (=T) read the domain configuration in 'domain_cfg.nc" file25 ln_read_cfg = .false. ! (=T) read the domain configuration in 'domain_cfg.nc" file 28 26 ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules 29 ln_write_cfg= .false. 27 ln_write_cfg= .false. ! (=T) create the domain configuration file 30 28 ! 31 29 cp_cfg = "lock" ! name of the configuration … … 49 47 / 50 48 !----------------------------------------------------------------------- 51 &namcrs ! Grid coarsening for dynamics output and/or52 ! passive tracer coarsened online simulations53 !-----------------------------------------------------------------------54 /55 !-----------------------------------------------------------------------56 49 &namtsd ! data : Temperature & Salinity 57 50 !----------------------------------------------------------------------- 58 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 59 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 60 sn_tem = 'output.init_0000' , 1 ,'votemper' , .true. , .true. , 'yearly' , '' , '' , '' 61 sn_sal = 'output.init_0000' , 1 ,'vosaline' , .true. , .true. , 'yearly' , '' , '' , '' 62 ! 63 cn_dir = './' ! root directory for the location of the runoff files 64 ln_tsd_init = .true. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 51 ln_tsd_init = .false. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 65 52 ln_tsd_tradmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) 66 53 / … … 78 65 / 79 66 !----------------------------------------------------------------------- 80 &namtra_qsr ! penetrative solar radiation81 !-----------------------------------------------------------------------82 /83 !-----------------------------------------------------------------------84 &namsbc_rnf ! runoffs namelist surface boundary condition85 !-----------------------------------------------------------------------86 /87 !-----------------------------------------------------------------------88 &namsbc_ssr ! surface boundary condition : sea surface restoring89 !-----------------------------------------------------------------------90 /91 !-----------------------------------------------------------------------92 &namsbc_alb ! albedo parameters93 !-----------------------------------------------------------------------94 /95 !-----------------------------------------------------------------------96 &namberg ! iceberg parameters97 !-----------------------------------------------------------------------98 /99 !-----------------------------------------------------------------------100 67 &namlbc ! lateral momentum boundary condition 101 68 !----------------------------------------------------------------------- 102 ! rn_shlat = 0. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat103 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat69 ! ! free slip ! partial slip ! no slip ! strong slip 70 rn_shlat = 0. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 104 71 / 105 72 !----------------------------------------------------------------------- … … 112 79 &nambbc ! bottom temperature boundary condition (default: NO) 113 80 !----------------------------------------------------------------------- 114 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom115 81 / 116 82 !----------------------------------------------------------------------- … … 122 88 !----------------------------------------------------------------------- 123 89 nn_eos = 1 ! type of equation of state and Brunt-Vaisala frequency 124 125 126 127 ! ! S-EOS coefficients (ln_seos=T):128 ! 129 rn_a0 = 2.e-1! thermal expension coefficient (nn_eos= 1)130 rn_b0 = 8.e-1! saline expension coefficient (nn_eos= 1)131 rn_lambda1 = 0. ! cabbeling coeff in T^2 (=0 for linear eos)132 rn_lambda2 = 0. ! cabbeling coeff in S^2 (=0 for linear eos)133 rn_mu1 = 0. ! thermobaric coeff. in T (=0 for linear eos)134 rn_mu2 = 0. ! thermobaric coeff. in S (=0 for linear eos)135 rn_nu = 0. ! cabbeling coeff in T*S (=0 for linear eos)90 ! =-1, TEOS-10 91 ! = 0, EOS-80 92 ! = 1, S-EOS (simplified eos) 93 ! ! S-EOS coefficients (nn_eos=1): 94 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 95 rn_a0 = 0.2 ! thermal expension coefficient (nn_eos= 1) 96 rn_b0 = 0. ! saline expension coefficient (nn_eos= 1) 97 rn_lambda1 = 0. ! cabbeling coeff in T^2 (=0 for linear eos) 98 rn_lambda2 = 0. ! cabbeling coeff in S^2 (=0 for linear eos) 99 rn_mu1 = 0. ! thermobaric coeff. in T (=0 for linear eos) 100 rn_mu2 = 0. ! thermobaric coeff. in S (=0 for linear eos) 101 rn_nu = 0. ! cabbeling coeff in T*S (=0 for linear eos) 136 102 / 137 103 !----------------------------------------------------------------------- … … 153 119 / 154 120 !----------------------------------------------------------------------- 155 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) 156 !----------------------------------------------------------------------- 157 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation 158 / 159 !---------------------------------------------------------------------------------- 121 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO) 122 !----------------------------------------------------------------------- 123 / 124 !----------------------------------------------------------------------- 160 125 &namtra_ldf ! lateral diffusion scheme for tracers 161 126 !---------------------------------------------------------------------------------- … … 203 168 / 204 169 !----------------------------------------------------------------------- 205 &namtra_dmp ! tracer: T & S newtonian damping 170 &namtra_ldfeiv ! eddy induced velocity param. (default: NO) 171 !----------------------------------------------------------------------- 172 / 173 !----------------------------------------------------------------------- 174 &namtra_dmp ! tracer: T & S newtonian damping (default: YES) 206 175 !----------------------------------------------------------------------- 207 176 ln_tradmp = .false. ! add a damping termn (T) or not (F) … … 220 189 !----------------------------------------------------------------------- 221 190 ln_vvl_zstar = .true. ! zstar vertical coordinate 222 !----------------------------------------------------------------------- 223 &namdyn_vor ! option of physics/algorithm (not control by CPP keys) 191 / 192 !----------------------------------------------------------------------- 193 &namdyn_vor ! option of physics/algorithm 224 194 !----------------------------------------------------------------------- 225 195 ln_dynvor_ene = .false. ! enstrophy conserving scheme … … 252 222 ! ! Type of the operator : 253 223 ! ! no diffusion: set ln_dynldf_lap=..._blp=F 254 ln_dynldf_lap = . true.! laplacian operator224 ln_dynldf_lap = .false. ! laplacian operator 255 225 ln_dynldf_blp = .false. ! bilaplacian operator 256 226 ! ! Direction of action : … … 276 246 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 277 247 rn_avt0 = 0. ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 278 ln_zdfevd = .false. ! enhanced vertical diffusion (evd) (T) or not (F) 279 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) 280 / 281 !----------------------------------------------------------------------- 282 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 283 !----------------------------------------------------------------------- 284 / 285 !----------------------------------------------------------------------- 286 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") 287 !----------------------------------------------------------------------- 288 / 289 !----------------------------------------------------------------------- 290 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") 291 !----------------------------------------------------------------------- 248 ln_zdfevd = .false. ! enhanced vertical diffusion (evd) 249 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm 292 250 / 293 251 !----------------------------------------------------------------------- … … 299 257 !----------------------------------------------------------------------- 300 258 / 301 !-----------------------------------------------------------------------302 &namptr ! Poleward Transport Diagnostic303 !-----------------------------------------------------------------------304 /305 !-----------------------------------------------------------------------306 &namhsb ! Heat and salt budgets307 !-----------------------------------------------------------------------308 /309 !-----------------------------------------------------------------------310 &namobs ! observation usage311 !-----------------------------------------------------------------------312 /313 !-----------------------------------------------------------------------314 &nam_asminc ! assimilation increments ('key_asminc')315 !-----------------------------------------------------------------------316 / -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90
r6901 r6906 10 10 !!---------------------------------------------------------------------- 11 11 12 !!---------------------------------------------------------------------- 12 !!---------------------------------------------------------------------------------- 13 13 !! usr_def_hgr : initialize the horizontal mesh for LOCK_EXCHANGE configuration 14 !!---------------------------------------------------------------------- 14 !!---------------------------------------------------------------------------------- 15 15 USE dom_oce , ONLY: nimpp, njmpp ! ocean space and time domain 16 16 USE par_oce ! ocean space and time domain … … 46 46 !! 47 47 !! ** Method : set all intent(out) argument to a proper value 48 !! LOCK_EXCHANGE configuration : uniform grid spacing ( 1km)48 !! LOCK_EXCHANGE configuration : uniform grid spacing (500m) 49 49 !! without Coriolis force (f=0) 50 50 !! … … 71 71 IF(lwp) WRITE(numout,*) 72 72 IF(lwp) WRITE(numout,*) 'usr_def_hgr : LOCK_EXCHANGE configuration bassin' 73 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ uniform 1kmgrid spacing WITHOUT Coriolis force (f=0)'73 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ uniform grid spacing WITHOUT Coriolis force (f=0)' 74 74 ! 75 75 ! !== grid point position ==! (in kilometers) … … 77 77 DO jj = 1, jpj 78 78 DO ji = 1, jpi ! longitude 79 p glamt(ji,jj) = zfact * ( - 0.5 + ( REAL( ji-1 + nimpp-1 ) , wp ) )80 p glamu(ji,jj) = zfact * ( - 0.5 + ( REAL( ji-1 + nimpp-1 ) + 0.5, wp ) )81 p glamv(ji,jj) = pglamt(ji,jj)82 p glamf(ji,jj) = pglamu(ji,jj)79 plamt(ji,jj) = zfact * ( - 0.5 + REAL( ji-1 + nimpp-1 , wp ) ) 80 plamu(ji,jj) = zfact * ( REAL( ji-1 + nimpp-1 , wp ) ) 81 plamv(ji,jj) = plamt(ji,jj) 82 plamf(ji,jj) = plamu(ji,jj) 83 83 ! ! latitude 84 p gphit(ji,jj) = zfact * ( - 0.5 + ( REAL( jj-1 + njmpp-1 ), wp ) )85 p gphiu(ji,jj) = pgphit(ji,jj)86 p gphiv(ji,jj) = zfact * ( - 0.5 + ( REAL( jj-1 + njmpp-1 ) + 0.5, wp ) )87 p gphif(ji,jj) = pgphiv(ji,jj)84 pphit(ji,jj) = zfact * ( - 0.5 + REAL( jj-1 + njmpp-1 , wp ) ) 85 pphiu(ji,jj) = pphit(ji,jj) 86 pphiv(ji,jj) = zfact * ( REAL( jj-1 + njmpp-1 , wp ) ) 87 pphif(ji,jj) = pphiv(ji,jj) 88 88 END DO 89 89 END DO -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_istate.F90
r6901 r6906 14 14 !!---------------------------------------------------------------------- 15 15 USE par_oce ! ocean space and time domain 16 USE dom_oce , ONLY : glamt 16 17 USE phycst ! physical constants 17 18 ! … … 49 50 ! 50 51 INTEGER :: jk ! dummy loop indices 52 REAL(wp) :: zdam ! location of dam [Km] 51 53 !!---------------------------------------------------------------------- 52 54 ! 53 55 IF(lwp) WRITE(numout,*) 54 IF(lwp) WRITE(numout,*) 'usr_def_istate : LOCK_EXCHANGE configuration, analytical definition of initial state ' 55 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ Ocean at rest, with a dam on T profile, and uniform S profile' 56 IF(lwp) WRITE(numout,*) 'usr_def_istate : LOCK_EXCHANGE configuration, analytical definition of initial state' 57 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ Ocean at rest, with a constant salinity (not used as rho=F(T) ' 58 IF(lwp) WRITE(numout,*) ' and a vertical density front with a 5 kg/m3 difference located at glam=32km' 59 IF(lwp) WRITE(numout,*) ' (i.e. a temperature difference of 25 degrees with rn_a0 = 0.2' 60 ! 61 ! rn_a0 = 2.e-1 ! thermal expension coefficient (nn_eos= 1) 62 ! rho = rau0 - rn_a0 * (T-10) 63 ! delta_T = 25 degrees ==>> delta_rho = 25 * rn_a0 = 5 kg/m3 56 64 ! 57 65 pu (:,:,:) = 0._wp ! ocean at rest … … 60 68 ! 61 69 ! ! T & S profiles 62 pts(:,:,:,jp_tem) = 30._wp * tmask(:,:,:) 70 zdam = 32. ! in kilometers 71 pts(:,:,:,jp_tem) = 30._wp * ptmask(:,:,:) 63 72 DO jk = 1, jpkm1 64 !SF ADD rn_dx2 65 WHERE( gphit(:,:) <= zdam ) pts(:,:,jk,jp_tem) = 5._wp * tmask(:,:,jk) 73 WHERE( glamt(:,:) <= zdam ) pts(:,:,jk,jp_tem) = 5._wp * ptmask(:,:,jk) 66 74 END DO 67 75 ! 68 pts(:,:,:,jp_sal) = 35. *tmask(:,:,:)76 pts(:,:,:,jp_sal) = 35._wp * ptmask(:,:,:) 69 77 ! 70 78 END SUBROUTINE usr_def_istate -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_nam.F90
r6901 r6906 55 55 INTEGER , INTENT(out) :: kperio ! lateral global domain b.c. 56 56 ! 57 INTEGER :: ios ! Local integer57 INTEGER :: ios, ii ! Local integer 58 58 !! 59 59 NAMELIST/namusr_def/ ln_zco, ln_zps, ln_sco, rn_dx, rn_dz 60 60 !!---------------------------------------------------------------------- 61 ! 62 ii = 1 61 63 ! 62 64 REWIND( numnam_cfg ) ! Namelist namusr_def (exist in namelist_cfg only) … … 64 66 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist', .TRUE. ) 65 67 ! 66 !!gm This does not work... I don't know how to write namusr_def in "output.namelist.dyn" 67 WRITE( ldnam(ii), namusr_def ) 68 !!gm 68 WRITE( ldnam(:), namusr_def ) 69 69 ! 70 70 ! Global Domain size: LOCK_EXCHANGE domain is 64 km x 3 grid-points x 20 m 71 71 kpi = INT( 64.e3 / rn_dx ) + 2 72 72 kpj = 3 73 kpk = INT( (2.e-1 ) * 1.e3/ rn_dz ) + 173 kpk = INT( 20. / rn_dz ) + 1 74 74 ! 75 75 ! ! control print … … 77 77 WRITE(ldtxt(ii),*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' ; ii = ii + 1 78 78 WRITE(ldtxt(ii),*) '~~~~~~~~~~~ ' ; ii = ii + 1 79 WRITE(ldtxt(ii),*) ' Namelist namusr_def : LOCK_EXCHANGE test case' 79 WRITE(ldtxt(ii),*) ' Namelist namusr_def : LOCK_EXCHANGE test case' ; ii = ii + 1 80 80 WRITE(ldtxt(ii),*) ' type of vertical coordinate : ' ; ii = ii + 1 81 81 WRITE(ldtxt(ii),*) ' z-coordinate flag ln_zco = ', ln_zco ; ii = ii + 1 … … 84 84 WRITE(ldtxt(ii),*) ' horizontal resolution rn_dx = ', rn_dx, ' meters' ; ii = ii + 1 85 85 WRITE(ldtxt(ii),*) ' vertical resolution rn_dz = ', rn_dz, ' meters' ; ii = ii + 1 86 WRITE(ldtxt(ii),*) ' LOCK_EXCHANGE domain = 64 km x 3 grid-points x 20 m' 86 WRITE(ldtxt(ii),*) ' LOCK_EXCHANGE domain = 64 km x 3 grid-points x 20 m' ; ii = ii + 1 87 87 WRITE(ldtxt(ii),*) ' resulting global domain size : jpiglo = ', kpi ; ii = ii + 1 88 88 WRITE(ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 1 -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/cpp_LOCK_EXCHANGE.fcm
r6901 r6906 1 bld::tool::fppkeys key_zdfcst key_iomput key_mpp_mpi key_mpp_rep 1 bld::tool::fppkeys key_zdfcst key_iomput key_mpp_mpi key_mpp_rep key_nosignedzero -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/OVERFLOW/EXP00/namelist_cfg
r6905 r6906 6 6 !----------------------------------------------------------------------- 7 7 ! ! type of vertical coordinate 8 ln_zco = .false. 8 ln_zco = .false. ! z-coordinate 9 9 ln_zps = .false. ! z-partial-step coordinate 10 ln_sco = .true. ! s-coordinate10 ln_sco = .true. ! s-coordinate 11 11 rn_dx = 1000. ! horizontal resolution [meters] 12 12 rn_dz = 20. ! vertical resolution [meters] … … 19 19 cn_exp = "overfl-FCT2-flux-ubs-ens" ! experience name 20 20 nn_it000 = 1 ! first time step 21 nn_itend = 1 80 ! for 34h of simulation (=12240)21 nn_itend = 12240 ! for 34h of simulation (=12240) 22 22 nn_istate = 1 ! output the initial state (1) or not (0) 23 23 nn_stock = 12240 ! frequency of creation of a restart file (modulo referenced to 1) … … 29 29 ln_read_cfg = .false. ! (=T) read the domain configuration in 'domain_cfg.nc" file 30 30 ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules 31 ln_write_cfg= . true. !.false. ! (=T) create the domain configuration file31 ln_write_cfg= .false. ! (=T) create the domain configuration file 32 32 ! 33 33 cp_cfg = "overflow" ! name of the configuration … … 51 51 / 52 52 !----------------------------------------------------------------------- 53 &namcrs ! Grid coarsening for dynamics output and/or54 ! passive tracer coarsened online simulations55 !-----------------------------------------------------------------------56 /57 !-----------------------------------------------------------------------58 53 &namtsd ! data : Temperature & Salinity 59 54 !----------------------------------------------------------------------- 60 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !61 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !62 sn_tem = 'output.init_0000' , 1 ,'votemper' , .true. , .true. , 'yearly' , '' , '' , ''63 sn_sal = 'output.init_0000' , 1 ,'vosaline' , .true. , .true. , 'yearly' , '' , '' , ''64 !65 cn_dir = './' ! root directory for the location of the runoff files66 55 ln_tsd_init = .false. ! Initialisation of ocean T & S with T &S input data (T) or not (F) 67 56 ln_tsd_tradmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/OVERFLOW/MY_SRC/usrdef_istate.F90
r6905 r6906 50 50 ! 51 51 INTEGER :: jk ! dummy loop indices 52 REAL(wp) :: zdam ! location of dam [ m]52 REAL(wp) :: zdam ! location of dam [Km] 53 53 !!---------------------------------------------------------------------- 54 54 ! … … 56 56 IF(lwp) WRITE(numout,*) 'usr_def_istate : OVERFLOW configuration, analytical definition of initial state' 57 57 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ Ocean at rest, with a constant salinity (not used as rho=F(T) ' 58 IF(lwp) WRITE(numout,*) ' and a vertical density front with a 2 kg/m3 difference located'59 58 IF(lwp) WRITE(numout,*) ' and a vertical density front with a 2 kg/m3 difference located at glam=20km' 60 59 IF(lwp) WRITE(numout,*) ' (i.e. a temperature difference of 10 degrees with rn_a0 = 0.2' -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/OVERFLOW/MY_SRC/usrdef_nam.F90
r6905 r6906 47 47 !! ** Method : read in namusr_def containing all the user specific namelist parameter 48 48 !! 49 !! Here GYREconfiguration49 !! Here OVERFLOW configuration 50 50 !! 51 51 !! ** input : - namusr_def namelist found in namelist_cfg … … 68 68 WRITE( ldnam(:), namusr_def ) 69 69 ! 70 ! Global Domain size: OVERFLOW domain is 70 ! Global Domain size: OVERFLOW domain is 200 km x 3 grid-points x 2000 m 71 71 kpi = INT( 200.e3 / rn_dx ) + 2 72 72 kpj = 3 … … 78 78 WRITE(ldtxt(ii),*) '~~~~~~~~~~~ ' ; ii = ii + 1 79 79 WRITE(ldtxt(ii),*) ' Namelist namusr_def : OVERFLOW test case' ; ii = ii + 1 80 WRITE(ldtxt(ii),*) ' type of vertical coordinate : ' 80 WRITE(ldtxt(ii),*) ' type of vertical coordinate : ' ; ii = ii + 1 81 81 WRITE(ldtxt(ii),*) ' z-coordinate flag ln_zco = ', ln_zco ; ii = ii + 1 82 82 WRITE(ldtxt(ii),*) ' z-partial-step coordinate flag ln_zps = ', ln_zps ; ii = ii + 1 … … 84 84 WRITE(ldtxt(ii),*) ' horizontal resolution rn_dx = ', rn_dx, ' meters' ; ii = ii + 1 85 85 WRITE(ldtxt(ii),*) ' vertical resolution rn_dz = ', rn_dz, ' meters' ; ii = ii + 1 86 WRITE(ldtxt(ii),*) ' OVERFLOW domain = 3 grid-points x 200 km x 2000 m'; ii = ii + 186 WRITE(ldtxt(ii),*) ' OVERFLOW domain = 200 km x 3 grid-points x 2000 m' ; ii = ii + 1 87 87 WRITE(ldtxt(ii),*) ' resulting global domain size : jpiglo = ', kpi ; ii = ii + 1 88 88 WRITE(ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 1 -
branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/OVERFLOW/MY_SRC/usrdef_zgr.F90
r6905 r6906 7 7 !! 8 8 !!============================================================================== 9 !! History : 4.0 ! 2016-08 ( S. Flavoni) Original code9 !! History : 4.0 ! 2016-08 (G. Madec) Original code 10 10 !!---------------------------------------------------------------------- 11 11 … … 154 154 ! 155 155 IF ( ln_zps ) THEN !== zps-coordinate ==! (partial bottom-steps) 156 !157 ! CALL ctl_stop( 'STOP', ' zps coordinate not yet coded' )158 156 ! 159 157 ze3min = 0.1_wp * rn_dz
Note: See TracChangeset
for help on using the changeset viewer.
