- Timestamp:
- 2019-09-20T17:28:02+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/ORCA2_ICE_PISCES/EXPREF/context_nemo.xml
r11413 r11586 42 42 <axis id="iax_20C" long_name="20 degC isotherm" unit="degC" /> 43 43 <axis id="iax_28C" long_name="28 degC isotherm" unit="degC" /> 44 <!-- ABL vertical axis definition -->45 <axis id="ght_abl" long_name="ABL Vertical T levels" unit="m" positive="up" />46 <axis id="ghw_abl" long_name="ABL Vertical W levels" unit="m" positive="up" />47 44 </axis_definition> 48 45 -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/field_def_nemo-ice.xml
r11413 r11586 46 46 47 47 <!-- melt ponds --> 48 <field id="iceapnd" long_name="melt pond fraction" standard_name="sea_ice_meltpond_fraction" unit="%" /> 48 <field id="iceapnd" long_name="melt pond concentration" standard_name="sea_ice_meltpond_concentration" unit="" /> 49 <field id="icehpnd" long_name="melt pond depth" standard_name="sea_ice_meltpond_depth" unit="m" /> 49 50 <field id="icevpnd" long_name="melt pond volume" standard_name="sea_ice_meltpond_volume" unit="m" /> 50 51 … … 328 329 <field field_ref="icesalt" name="sisali" /> 329 330 <field field_ref="iceapnd" name="siapnd" /> 331 <field field_ref="icehpnd" name="sihpnd" /> 330 332 <field field_ref="icevpnd" name="sivpnd" /> 331 333 <field field_ref="iceage" name="siage" /> -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/field_def_nemo-oce.xml
r11413 r11586 237 237 238 238 <!-- SBC --> 239 239 240 <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> 240 241 … … 300 301 <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" /> 301 302 302 303 303 <!-- available if key_oasis3 + conservative method --> 304 304 <field id="rain" long_name="Liquid precipitation" standard_name="rainfall_flux" unit="kg/m2/s" /> … … 348 348 <!-- ABL --> 349 349 <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> 350 350 351 <!-- variables available with ABL on atmospheric T grid--> 351 352 <field_group id="grid_ABL3D" grid_ref="grid_TA_3D" > … … 354 355 <field id="t_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 355 356 <field id="q_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 357 <!-- debug (to be removed) --> 358 <field id="u_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 359 <field id="v_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 360 <field id="t_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 361 <field id="q_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 356 362 <field id="coeft" long_name="ABL nudging coefficient" standard_name="coeft" unit="" /> 357 363 </field_group> 364 358 365 <field_group id="grid_ABL2D" grid_ref="grid_TA_2D" > 359 366 <field id="pblh" long_name="ABL height" standard_name="abl_height" unit="m" /> 367 <field id="uz1_abl" long_name="ABL i-horizontal velocity" standard_name="abl_x_velocity" unit="m/s" /> 368 <field id="vz1_abl" long_name="ABL j-horizontal velocity" standard_name="abl_y_velocity" unit="m/s" /> 369 <field id="uvz1_abl" long_name="ABL wind speed module" standard_name="abl_wind_speed" unit="m/s" > sqrt( uz1_abl^2 + vz1_abl^2 ) </field> 370 <field id="tz1_abl" long_name="ABL potential temperature" standard_name="abl_theta" unit="K" /> 371 <field id="qz1_abl" long_name="ABL specific humidity" standard_name="abl_qspe" unit="kg/kg" /> 372 <field id="uz1_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 373 <field id="vz1_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 374 <field id="uvz1_dta" long_name="DTA wind speed module" standard_name="dta_wind_speed" unit="m/s" > sqrt( uz1_dta^2 + vz1_dta^2 ) </field> 375 <field id="tz1_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 376 <field id="qz1_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> 377 <!-- debug (to be removed) --> 378 <field id="uz1_geo" long_name="GEO i-horizontal velocity" standard_name="geo_x_velocity" unit="m/s" /> 379 <field id="vz1_geo" long_name="GEO j-horizontal velocity" standard_name="geo_y_velocity" unit="m/s" /> 380 <field id="uvz1_geo" long_name="GEO wind speed module" standard_name="geo_wind_speed" unit="m/s" > sqrt( uz1_geo^2 + vz1_geo^2 ) </field> 360 381 </field_group> 382 361 383 </field_group> <!-- ABL --> 384 362 385 363 386 <!-- U grid --> … … 375 398 <field id="uocet" long_name="ocean transport along i-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_U_3D" /> 376 399 <field id="uoces" long_name="ocean transport along i-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_U_3D" /> 400 <field id="ssuww" long_name="ocean surface wind work along i-axis" standard_name="surface_x_wind_work" unit="N/m*s" > utau * ssu </field> 377 401 378 402 <!-- u-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> … … 429 453 <field id="vocet" long_name="ocean transport along j-axis times temperature (CRS)" unit="degC*m/s" grid_ref="grid_V_3D" /> 430 454 <field id="voces" long_name="ocean transport along j-axis times salinity (CRS)" unit="1e-3*m/s" grid_ref="grid_V_3D" /> 455 <field id="ssvww" long_name="ocean surface wind work along j-axis" standard_name="surface_y_wind_work" unit="N/m*s" > vtau * ssv </field> 456 431 457 432 458 <!-- v-eddy diffusivity coefficients (available if ln_traldf_OFF=F) --> -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/namelist_ice_ref
r11413 r11586 70 70 ln_str_H79 = .true. ! ice strength param.: Hibler_79 => P = pstar*<h>*exp(-c_rhg*A) 71 71 rn_pstar = 2.0e+04 ! ice strength thickness parameter [N/m2] 72 rn_crhg = 20.0! ice strength conc. parameter (-)72 rn_crhg = 20.0 ! ice strength conc. parameter (-) 73 73 ! -- ice_rdgrft -- ! 74 74 rn_csrdg = 0.5 ! fraction of shearing energy contributing to ridging … … 177 177 &namthd_pnd ! Melt ponds 178 178 !------------------------------------------------------------------------------ 179 ln_pnd_H12 = .false. ! activate evolutive melt ponds (from Holland et al 2012) 180 ln_pnd_CST = .false. ! activate constant melt ponds 181 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 182 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 183 ln_pnd_alb = .false. ! melt ponds affect albedo or not 179 ln_pnd = .false. ! activate melt ponds or not 180 ln_pnd_H12 = .false. ! activate evolutive melt ponds (from Holland et al 2012) 181 ln_pnd_CST = .false. ! activate constant melt ponds 182 rn_apnd = 0.2 ! prescribed pond fraction, at Tsu=0 degC 183 rn_hpnd = 0.05 ! prescribed pond depth, at Tsu=0 degC 184 ln_pnd_alb = .false. ! melt ponds affect albedo or not 184 185 / 185 186 !------------------------------------------------------------------------------ … … 207 208 rn_hpd_ini_n = 0.05 ! initial pond depth (m), North 208 209 rn_hpd_ini_s = 0.05 ! " " South 209 ! ! if ln_iceini_file=T210 ! -- for ln_iceini_file = T 210 211 sn_hti = 'Ice_initialization' , -12 ,'hti' , .false. , .true., 'yearly' , '' , '', '' 211 212 sn_hts = 'Ice_initialization' , -12 ,'hts' , .false. , .true., 'yearly' , '' , '', '' … … 215 216 sn_tsu = 'Ice_initialization' , -12 ,'tsu' , .false. , .true., 'yearly' , '' , '', '' 216 217 sn_tms = 'NOT USED' , -12 ,'tms' , .false. , .true., 'yearly' , '' , '', '' 218 ! melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 217 219 sn_apd = 'NOT USED' , -12 ,'apd' , .false. , .true., 'yearly' , '' , '', '' 218 220 sn_hpd = 'NOT USED' , -12 ,'hpd' , .false. , .true., 'yearly' , '' , '', '' … … 232 234 &namdia ! Diagnostics 233 235 !------------------------------------------------------------------------------ 234 ln_icediachk = .false. ! check online the heat, mass & salt budgets (T) or not (F) 236 ln_icediachk = .false. ! check online the heat, mass & salt budgets at each time step 237 ! ! rate of ice spuriously gained/lost. For ex., rn_icechk=1. <=> 1mm/year, rn_icechk=0.1 <=> 1mm/10years 238 rn_icechk_cel = 1. ! check at any gridcell => stops the code if violated (and writes a file) 239 rn_icechk_glo = 0.1 ! check over the entire ice cover => only prints warnings 235 240 ln_icediahsb = .false. ! output the heat, mass & salt budgets (T) or not (F) 236 241 ln_icectl = .false. ! ice points output for debug (T or F) -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SHARED/namelist_ref
r11413 r11586 199 199 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 200 200 ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) 201 ln_abl = .false. ! ABL formulation (T => fill namsbc_abl ) 201 202 ! ! Type of coupling (Ocean/Ice/Atmosphere) : 202 203 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) … … 280 281 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 281 282 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 283 / 284 !----------------------------------------------------------------------- 285 &namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T) 286 !----------------------------------------------------------------------- 287 cn_dir = './' ! root directory for the location of the ABL grid file 288 cn_dom = 'dom_cfg_abl.nc' 289 ln_hpgls_frc = .false. 290 ln_geos_winds = .false. 291 nn_dyn_restore = 2 ! restoring option for dynamical ABL variables: = 0 no restoring 292 ! = 1 equatorial restoring 293 ! = 2 global restoring 294 rn_ldyn_min = 4.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 295 rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour] 296 rn_ltra_min = 4.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour] 297 rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour] 298 nn_amxl = 0 ! mixing length: = 0 Deardorff 80 length-scale 299 ! = 1 length-scale based on the distance to the PBL height 300 ! = 2 Bougeault & Lacarrere 89 length-scale 301 rn_Cm = 0.0667 ! 0.126 in MesoNH 302 rn_Ct = 0.1667 ! 0.143 in MesoNH 303 rn_Ce = 0.4 ! 0.4 in MesoNH 304 rn_Ceps = 0.7 ! 0.85 in MesoNH 305 rn_Rod = 0.15 ! c0 in RMCA17 mixing length formulation (not yet implemented) 306 rn_Ric = 0.139 ! Critical Richardson number (to compute PBL height and diffusivities) 282 307 / 283 308 !----------------------------------------------------------------------- … … 600 625 nn_ice_dta = 0 ! = 0, bdy data are equal to the initial state 601 626 ! ! = 1, bdy data are read in 'bdydata .nc' files 602 rn_ice_tem = 270. ! si3 only: arbitrary temperature of incoming sea ice603 rn_ice_sal = 10. ! si3 only: -- salinity --604 rn_ice_age = 30. ! si3 only: -- age --605 627 ! 606 628 ln_tra_dmp =.false. ! open boudaries conditions for tracers … … 632 654 bn_sal = 'amm12_bdyT_tra' , 24. , 'vosaline', .true. , .false., 'daily' , '' , '' , '' 633 655 !* for si3 634 ! bn_a_i = 'amm12_bdyT_ice' , 24. , 'ileadfra', .true. , .false., 'daily' , '' , '' , '' 635 ! bn_h_i = 'amm12_bdyT_ice' , 24. , 'iicethic', .true. , .false., 'daily' , '' , '' , '' 636 ! bn_h_s = 'amm12_bdyT_ice' , 24. , 'isnowthi', .true. , .false., 'daily' , '' , '' , '' 656 bn_a_i = 'amm12_bdyT_ice' , 24. , 'siconc' , .true. , .false., 'daily' , '' , '' , '' 657 bn_h_i = 'amm12_bdyT_ice' , 24. , 'sithic' , .true. , .false., 'daily' , '' , '' , '' 658 bn_h_s = 'amm12_bdyT_ice' , 24. , 'snthic' , .true. , .false., 'daily' , '' , '' , '' 659 bn_t_i = 'NOT USED' , 24. , 'sitemp' , .true. , .false., 'daily' , '' , '' , '' 660 bn_t_s = 'NOT USED' , 24. , 'sntemp' , .true. , .false., 'daily' , '' , '' , '' 661 bn_tsu = 'NOT USED' , 24. , 'sittop' , .true. , .false., 'daily' , '' , '' , '' 662 bn_s_i = 'NOT USED' , 24. , 'sisalt' , .true. , .false., 'daily' , '' , '' , '' 663 ! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd) 664 bn_aip = 'NOT USED' , 24. , 'siapnd' , .true. , .false., 'daily' , '' , '' , '' 665 bn_hip = 'NOT USED' , 24. , 'sihpnd' , .true. , .false., 'daily' , '' , '' , '' 666 ! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds 667 rn_ice_tem = 270. ! arbitrary temperature of incoming sea ice 668 rn_ice_sal = 10. ! -- salinity -- 669 rn_ice_age = 30. ! -- age -- 670 rn_ice_apnd = 0.2 ! -- pond fraction = a_ip/a_i -- 671 rn_ice_hpnd = 0.05 ! -- pond depth -- 637 672 / 638 673 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SPITZ12/EXPREF/namelist_cfg
r11413 r11586 170 170 nn_ice_dta = 1 ! = 0, bdy data are equal to the initial state 171 171 ! ! = 1, bdy data are read in 'bdydata .nc' files 172 rn_ice_tem = 267. ! si3 only: arbitrary temperature of incoming sea ice173 rn_ice_sal = 6. ! si3 only: -- salinity --174 rn_ice_age = 365. ! si3 only: -- age --175 !176 172 nn_rimwidth = 1 ! width of the relaxation zone 177 173 ln_vol = .false. ! total volume correction (see nn_volctl parameter) -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/SPITZ12/EXPREF/namelist_ice_cfg
r11263 r11586 80 80 &namthd_pnd ! Melt ponds 81 81 !------------------------------------------------------------------------------ 82 ln_pnd_H12 = .true. ! activate evolutive melt ponds (from Holland et al 2012) 83 ln_pnd_alb = .true. ! melt ponds affect albedo or not 82 ln_pnd = .true. ! activate melt ponds or not 83 ln_pnd_H12 = .true. ! activate evolutive melt ponds (from Holland et al 2012) 84 ln_pnd_alb = .true. ! melt ponds affect albedo or not 84 85 / 85 86 -
NEMO/branches/2019/dev_r11265_ASINTER-01_Guillaume_ABL1D/cfgs/ref_cfgs.txt
r11413 r11586 14 14 eORCA025_ICE OCE ICE 15 15 eORCA025_ICE_ABL OCE ICE ABL 16 eORCA025_SAS_ICE_ABL OCE SAS ICE ABL
Note: See TracChangeset
for help on using the changeset viewer.