Changeset 12149 for NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW
- Timestamp:
- 2019-12-10T15:03:24+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/context_nemo.xml
r9930 r12149 6 6 <context id="nemo"> 7 7 <!-- $id$ --> 8 <variable_definition> 9 <!-- Year of time origin for NetCDF files; defaults to 1800 --> 10 <variable id="ref_year" type="int" > 1800 </variable> 11 <variable id="rau0" type="float" > 1026.0 </variable> 12 <variable id="cpocean" type="float" > 3991.86795711963 </variable> 13 <variable id="convSpsu" type="float" > 0.99530670233846 </variable> 14 <variable id="rhoic" type="float" > 917.0 </variable> 15 <variable id="rhosn" type="float" > 330.0 </variable> 16 <variable id="missval" type="float" > 1.e20 </variable> 17 </variable_definition> 8 18 <!-- Fields definition --> 9 19 <field_definition src="./field_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics --> … … 18 28 19 29 <axis_definition> 20 <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" />21 <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" />22 <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" />23 <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" />24 <axis id="nfloat" long_name="Float number" unit="-" />25 <axis id="icbcla" long_name="Iceberg class" unit="1"/>26 <axis id="ncatice" long_name="Ice category" unit="1"/>27 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC"/>28 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC"/>30 <axis id="deptht" long_name="Vertical T levels" unit="m" positive="down" /> 31 <axis id="depthu" long_name="Vertical U levels" unit="m" positive="down" /> 32 <axis id="depthv" long_name="Vertical V levels" unit="m" positive="down" /> 33 <axis id="depthw" long_name="Vertical W levels" unit="m" positive="down" /> 34 <axis id="nfloat" long_name="Float number" unit="-" /> 35 <axis id="icbcla" long_name="Iceberg class" unit="1" /> 36 <axis id="ncatice" long_name="Ice category" unit="1" /> 37 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC" /> 38 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC" /> 29 39 </axis_definition> 30 40 -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg
r10075 r12149 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_ubs_cfg
r10075 r12149 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_ubs_cfg
r10075 r12149 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT2_flux_ubs_cfg
r10075 r12149 211 211 ! ! coefficients 212 212 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 213 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)213 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 214 214 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 215 215 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg
r10075 r12149 34 34 cn_exp = "OVF_zps_FCT4_flux_ubs" ! experience name 35 35 nn_it000 = 1 ! first time step 36 !nn_itend = 6120 ! here 17h of simulation (=6120 time-step) 37 nn_itend = 5760 ! here 16h of simulation (=5760 time-step) abort after 5802 for zps: pb of physics conditions 36 nn_itend = 6120 ! here 17h of simulation (=6120 time-step) 38 37 nn_istate = 0 ! output the initial state (1) or not (0) 39 nn_stock = 1080 ! frequency of creation of a restart file (modulo referenced to 1)40 nn_write = 1080 ! frequency of write in the output file (modulo referenced to nn_it000)38 nn_stock = 6120 ! frequency of creation of a restart file (modulo referenced to 1) 39 nn_write = 6120 ! frequency of write in the output file (modulo referenced to nn_it000) 41 40 / 42 41 !----------------------------------------------------------------------- … … 255 254 &namzdf ! vertical physics manager (default: NO selection) 256 255 !----------------------------------------------------------------------- 256 ! ! adaptive-implicit vertical advection 257 ln_zad_Aimp = .true. ! Courant number dependent scheme (Shchepetkin 2015) 258 ! 257 259 ! ! type of vertical closure (required) 258 260 ln_zdfcst = .true. ! constant mixing … … 279 281 ! ! coefficients 280 282 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 281 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)283 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 282 284 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 283 285 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) … … 291 293 !! namdiu Cool skin and warm layer models (default: OFF) 292 294 !! namdiu Cool skin and warm layer models (default: OFF) 293 !! namflo float parameters ( "key_float")294 !! nam_diaharm Harmonic analysis of tidal constituents ( "key_diaharm")295 !! nam dct transports through some sections ("key_diadct")295 !! namflo float parameters (default: OFF) 296 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF) 297 !! nam_diadct transports through some sections (default: OFF) 296 298 !! nam_diatmb Top Middle Bottom Output (default: OFF) 297 299 !! nam_dia25h 25h Mean Output (default: OFF) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_vect_een_cfg
r10075 r12149 210 210 ! ! coefficients 211 211 rn_avm0 = 1.e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 212 rn_avt0 = 0.0 e! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)212 rn_avt0 = 0.0 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 213 213 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 214 214 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) -
NEMO/branches/2019/ENHANCE-03_closea/tests/OVERFLOW/MY_SRC/usrdef_nam.F90
r10074 r12149 39 39 CONTAINS 40 40 41 SUBROUTINE usr_def_nam( ldtxt, ldnam,cd_cfg, kk_cfg, kpi, kpj, kpk, kperio )41 SUBROUTINE usr_def_nam( cd_cfg, kk_cfg, kpi, kpj, kpk, kperio ) 42 42 !!---------------------------------------------------------------------- 43 43 !! *** ROUTINE dom_nam *** … … 51 51 !! ** input : - namusr_def namelist found in namelist_cfg 52 52 !!---------------------------------------------------------------------- 53 CHARACTER(len=*), DIMENSION(:), INTENT(out) :: ldtxt, ldnam ! stored print information54 53 CHARACTER(len=*) , INTENT(out) :: cd_cfg ! configuration name 55 54 INTEGER , INTENT(out) :: kk_cfg ! configuration resolution … … 57 56 INTEGER , INTENT(out) :: kperio ! lateral global domain b.c. 58 57 ! 59 INTEGER :: ios , ii! Local integer58 INTEGER :: ios ! Local integer 60 59 !! 61 60 NAMELIST/namusr_def/ ln_zco, ln_zps, ln_sco, rn_dx, rn_dz 62 61 !!---------------------------------------------------------------------- 63 62 ! 64 ii = 165 !66 63 REWIND( numnam_cfg ) ! Namelist namusr_def (exist in namelist_cfg only) 67 64 READ ( numnam_cfg, namusr_def, IOSTAT = ios, ERR = 902 ) 68 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' , .TRUE.)65 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namusr_def in configuration namelist' ) 69 66 ! 70 WRITE( ldnam(:), namusr_def )67 IF(lwm) WRITE( numond, namusr_def ) 71 68 ! 72 69 cd_cfg = 'OVERFLOW' ! name & resolution (not used) … … 79 76 ! 80 77 ! ! control print 81 WRITE( ldtxt(ii),*) ' ' ; ii = ii + 182 WRITE( ldtxt(ii),*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' ; ii = ii + 183 WRITE( ldtxt(ii),*) '~~~~~~~~~~~ ' ; ii = ii + 184 WRITE( ldtxt(ii),*) ' Namelist namusr_def : OVERFLOW test case' ; ii = ii + 185 WRITE( ldtxt(ii),*) ' type of vertical coordinate : ' ; ii = ii + 186 WRITE( ldtxt(ii),*) ' z-coordinate flag ln_zco = ', ln_zco ; ii = ii + 187 WRITE( ldtxt(ii),*) ' z-partial-step coordinate flag ln_zps = ', ln_zps ; ii = ii + 188 WRITE( ldtxt(ii),*) ' s-coordinate flag ln_sco = ', ln_sco ; ii = ii + 189 WRITE( ldtxt(ii),*) ' horizontal resolution rn_dx = ', rn_dx, ' meters' ; ii = ii + 190 WRITE( ldtxt(ii),*) ' vertical resolution rn_dz = ', rn_dz, ' meters' ; ii = ii + 191 WRITE( ldtxt(ii),*) ' OVERFLOW domain = 200 km x 3 grid-points x 2000 m' ; ii = ii + 192 WRITE( ldtxt(ii),*) ' resulting global domain size : jpiglo = ', kpi ; ii = ii + 193 WRITE( ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 194 WRITE( ldtxt(ii),*) ' jpkglo = ', kpk ; ii = ii + 178 WRITE(numout,*) ' ' 79 WRITE(numout,*) 'usr_def_nam : read the user defined namelist (namusr_def) in namelist_cfg' 80 WRITE(numout,*) '~~~~~~~~~~~ ' 81 WRITE(numout,*) ' Namelist namusr_def : OVERFLOW test case' 82 WRITE(numout,*) ' type of vertical coordinate : ' 83 WRITE(numout,*) ' z-coordinate flag ln_zco = ', ln_zco 84 WRITE(numout,*) ' z-partial-step coordinate flag ln_zps = ', ln_zps 85 WRITE(numout,*) ' s-coordinate flag ln_sco = ', ln_sco 86 WRITE(numout,*) ' horizontal resolution rn_dx = ', rn_dx, ' meters' 87 WRITE(numout,*) ' vertical resolution rn_dz = ', rn_dz, ' meters' 88 WRITE(numout,*) ' OVERFLOW domain = 200 km x 3 grid-points x 2000 m' 89 WRITE(numout,*) ' resulting global domain size : jpiglo = ', kpi 90 WRITE(numout,*) ' jpjglo = ', kpj 91 WRITE(numout,*) ' jpkglo = ', kpk 95 92 ! 96 93 ! ! Set the lateral boundary condition of the global domain 97 94 kperio = 0 ! OVERFLOW configuration : close basin 98 95 ! 99 WRITE( ldtxt(ii),*) ' ' ; ii = ii + 1100 WRITE( ldtxt(ii),*) ' Lateral boundary condition of the global domain' ; ii = ii + 1101 WRITE( ldtxt(ii),*) ' OVERFLOW : closed basin jperio = ', kperio ; ii = ii + 196 WRITE(numout,*) ' ' 97 WRITE(numout,*) ' Lateral boundary condition of the global domain' 98 WRITE(numout,*) ' OVERFLOW : closed basin jperio = ', kperio 102 99 ! 103 100 END SUBROUTINE usr_def_nam
Note: See TracChangeset
for help on using the changeset viewer.