- Timestamp:
- 2020-06-26T10:26:32+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement
- Files:
-
- 18 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@12931 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg
r12489 r13159 81 81 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 82 82 ! Sea-ice : 83 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 84 ! except in AGRIF zoom where it has to be specified 83 nn_ice = 2 ! =0 no ice boundary condition 84 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 85 ! ! =2 or 3 for SI3 and CICE, respectively 85 86 ! Misc. options of sbc : 86 87 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg
r12489 r13159 78 78 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 79 79 ! Sea-ice : 80 nn_ice = 2 ! =0 Use SI3 model 80 nn_ice = 2 ! =0 no ice boundary condition 81 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 82 ! ! =2 or 3 for SI3 and CICE, respectively 81 83 ! Misc. options of sbc : 82 84 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg
r12489 r13159 78 78 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 79 79 ! Sea-ice : 80 nn_ice = 2 ! =0 Use SI3 model 80 nn_ice = 2 ! =0 no ice boundary condition 81 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 82 ! ! =2 or 3 for SI3 and CICE, respectively 81 83 ! Misc. options of sbc : 82 84 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg
r12489 r13159 81 81 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 82 82 ! Sea-ice : 83 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 84 ! except in AGRIF zoom where it has to be specified 83 nn_ice = 2 ! =0 no ice boundary condition 84 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 85 ! ! =2 or 3 for SI3 and CICE, respectively 85 86 ! Misc. options of sbc : 86 87 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/C1D_PAPA/EXPREF/file_def_nemo-oce.xml
r9799 r13159 53 53 <file id="file4" name_suffix="_grid_W" description="ocean W grid variables" > 54 54 <field field_ref="e3w" /> 55 <field field_ref="woce" name="wo" />56 55 <field field_ref="avt" name="difvho" /> 57 56 </file> -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/C1D_PAPA/EXPREF/namelist_cfg
r12489 r13159 49 49 &namdom ! time and space domain 50 50 !----------------------------------------------------------------------- 51 ln_linssh = .true. ! =T linear free surface ==>> model level are fixed in time 52 ! 51 53 rn_Dt = 360. ! time step for the dynamics and tracer 52 54 / … … 358 360 &namdyn_spg ! surface pressure gradient (default: NO selection) 359 361 !----------------------------------------------------------------------- 360 ln_dynspg_ts = .true. ! split-explicit free surface361 ln_bt_fw = .false. ! Forward integration of barotropic Eqs.362 ln_bt_av = .true. ! Time filtering of barotropic variables363 362 / 364 363 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg
r12749 r13159 22 22 cn_exp = "ORCA2" ! experience name 23 23 nn_it000 = 1 ! first time step 24 nn_itend = 5840 ! last time step (std 5840)24 nn_itend = 16 ! last time step (std 5475) 25 25 nn_date0 = 20130101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 26 26 nn_write = 4 ! frequency of write in the output file (modulo referenced to nn_it000) … … 31 31 &namdom ! time and space domain 32 32 !----------------------------------------------------------------------- 33 rn_Dt 33 rn_Dt = 5400. ! time step for the dynamics and tracer 34 34 / 35 35 !----------------------------------------------------------------------- … … 56 56 sn_sal = 'data_1m_salinity_nomask' , -1. ,'vosaline', .true. , .true. , 'yearly' , '' , '' , '' 57 57 / 58 59 58 !!====================================================================== 60 59 !! *** Surface Boundary Condition namelists *** !! … … 70 69 !! namsbc_rnf river runoffs (ln_rnf =T) 71 70 !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T) 71 !! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr ) 72 !! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T) 72 73 !! namsbc_wave external fields from wave model (ln_wave =T) 73 74 !! namberg iceberg floats (ln_icebergs=T) … … 78 79 !----------------------------------------------------------------------- 79 80 nn_fsbc = 1 ! frequency of SBC module call 80 ! ! (controlsea-ice & iceberg model call)81 ! Type of air-sea fluxes 81 ! (also = the frequency of sea-ice & iceberg model call) 82 ! Type of air-sea fluxes 82 83 ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) 83 84 ln_abl = .true. ! ABL formulation (T => fill namsbc_abl ) … … 86 87 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 87 88 ! ! =2 or 3 for SI3 and CICE, respectively 88 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 89 ! ! =F levitating ice (no pressure, mass and salt exchanges) 90 ! Misc. options of sbc : 89 ! Misc. options of sbc : 91 90 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 91 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 92 92 ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave 93 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)94 nn_fwb = 2 ! FreshWater Budget: =0 unchecked95 ! ! =1 global mean of e-p-r set to zero at each time step96 ! ! =2 annual global mean of e-p-r set to zero97 93 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 94 nn_fwb = 2 ! FreshWater Budget: 95 ! ! =2 annual global mean of e-p-r set to zero 98 96 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) 99 97 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 100 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 98 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 101 99 nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 102 100 ! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] … … 114 112 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 115 113 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 116 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 114 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31) 115 rn_zqt = 10. ! Air temperature & humidity reference height (m) 116 rn_zu = 10. ! Wind vector reference height (m) 117 117 ! 118 ln_humi_sph = .true. ! humidity "sn_humi" is specific humidity [kg/kg] 119 ln_tpot = .false. !!GS: compute potential temperature or not 118 ! Skin is ONLY available in ECMWF and COARE algorithms: 119 ln_skin_cs = .false. ! use the cool-skin parameterization => set nn_fsbc=1 and ln_dm2dc=.true.! 120 ln_skin_wl = .false. ! use the warm-layer " => set nn_fsbc=1 and ln_dm2dc=.true.! 121 ! 122 ln_humi_sph = .true. ! humidity specified below in "sn_humi" is specific humidity [kg/kg] if .true. 123 ln_humi_dpt = .false. ! humidity specified below in "sn_humi" is dew-point temperature [K] if .true. 124 ln_humi_rlh = .false. ! humidity specified below in "sn_humi" is relative humidity [%] if .true. 120 125 ! 121 cn_dir = './'! root directory for the bulk data location126 cn_dir = './' ! root directory for the bulk data location 122 127 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! 123 128 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! … … 142 147 cn_dir = './' ! root directory for the location of the ABL grid file 143 148 cn_dom = 'dom_cfg_abl_L25Z10.nc' 144 145 ln_rstart_abl = .false.146 149 ln_hpgls_frc = .true. 147 150 ln_geos_winds = .false. 148 151 nn_dyn_restore = 1 149 150 rn_ldyn_min = 4.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 152 rn_ldyn_min = 7.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 151 153 rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour] 152 rn_ltra_min = 4.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour]154 rn_ltra_min = 7.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour] 153 155 rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour] 154 156 / … … 200 202 &namberg ! iceberg parameters (default: OFF) 201 203 !----------------------------------------------------------------------- 202 ln_icebergs = .true. ! activate iceberg floats (force =F with "key_agrif")203 204 cn_dir = './' ! root directory for the location of drag coefficient files205 !______!___________!___________________!______________!______________!_________!___________!__________!__________!_______________!206 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !207 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename !208 sn_icb = 'calving', -1. , 'calving' , .true. , .true. , 'yearly' , '' , '' , ''209 204 / 210 205 !!====================================================================== -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r12501 r13159 80 80 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 81 81 ! Sea-ice : 82 nn_ice = 2 ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 83 ! except in AGRIF zoom where it has to be specified 82 nn_ice = 2 ! =0 no ice boundary condition 83 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 84 ! ! =2 or 3 for SI3 and CICE, respectively 84 85 ! Misc. options of sbc : 85 86 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) … … 120 121 / 121 122 !----------------------------------------------------------------------- 123 &namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T) 124 !----------------------------------------------------------------------- 125 / 126 !----------------------------------------------------------------------- 122 127 &namtra_qsr ! penetrative solar radiation (ln_traqsr =T) 123 128 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_top_cfg
r12377 r13159 20 20 ! 21 21 ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) 22 ln_trcbc = . true.! Enables Boundary conditions22 ln_trcbc = .false. ! Enables Boundary conditions 23 23 ! ! ! ! ! ! 24 24 ! ! name ! title of the field ! units ! init ! sbc ! cbc ! obc ! -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_top_cfg
r12377 r13159 20 20 ! 21 21 ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) 22 ln_trcbc = . true.! Enables Boundary conditions22 ln_trcbc = .false. ! Enables Boundary conditions 23 23 ! ! ! ! ! ! 24 24 ! ! name ! title of the field ! units ! init ! sbc ! cbc ! obc ! -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg
r12377 r13159 59 59 nn_fsbc = 1 ! frequency of SBC module call 60 60 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 61 nn_ice = 2 ! =2 sea-ice model ("key_SI3" or "key_cice") 61 nn_ice = 2 ! =0 no ice boundary condition 62 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 63 ! ! =2 or 3 for SI3 and CICE, respectively 62 64 / 63 65 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/SHARED/field_def_nemo-oce.xml
r12749 r13159 1 <?xml version="1.0"?> 1 <?xml version="1.0"?> 2 2 <!-- $id$ --> 3 3 … … 16 16 Configuration of multiple-linear-regression analysis (diamlr) 17 17 ===================================================================================================== 18 18 19 19 This field group configures diamlr for tidal harmonic analysis of field 20 20 ssh: in addition to a regressor for fitting the mean value (diamlr_r101), … … 73 73 74 74 </field_group> 75 76 <!-- 75 76 <!-- 77 77 ============================================================================================================ 78 78 = definition of all existing variables = … … 101 101 </field_group> 102 102 103 <!-- 103 <!-- 104 104 ============================================================================================================ 105 105 Physical ocean model variables … … 108 108 109 109 <!-- T grid --> 110 110 111 111 <field_group id="grid_T" grid_ref="grid_T_2D" > 112 112 <field id="e3t" long_name="T-cell thickness" standard_name="cell_thickness" unit="m" grid_ref="grid_T_3D" /> 113 113 <field id="e3ts" long_name="T-cell thickness" field_ref="e3t" standard_name="cell_thickness" unit="m" grid_ref="grid_T_SFC"/> 114 114 <field id="e3t_0" long_name="Initial T-cell thickness" standard_name="ref_cell_thickness" unit="m" grid_ref="grid_T_3D" /> 115 <field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/> 115 <field id="e3tb" long_name="bottom T-cell thickness" standard_name="bottom_cell_thickness" unit="m" grid_ref="grid_T_2D"/> 116 116 <field id="e3t_300" field_ref="e3t" grid_ref="grid_T_zoom_300" detect_missing_value="true" /> 117 117 <field id="e3t_vsum300" field_ref="e3t_300" grid_ref="grid_T_vsum" detect_missing_value="true" /> 118 118 <field id="masscello" long_name="Sea Water Mass per unit area" standard_name="sea_water_mass_per_unit_area" unit="kg/m2" grid_ref="grid_T_3D"/> 119 <field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/> 119 <field id="volcello" long_name="Ocean Volume" standard_name="ocean_volume" unit="m3" grid_ref="grid_T_3D"/> 120 120 <field id="toce" long_name="temperature" standard_name="sea_water_potential_temperature" unit="degC" grid_ref="grid_T_3D"/> 121 121 <field id="toce_e3t" long_name="temperature (thickness weighted)" unit="degC" grid_ref="grid_T_3D" > toce * e3t </field > … … 146 146 <field id="sst_cs" long_name="Delta SST of cool skin" unit="degC" /> 147 147 <field id="temp_3m" long_name="temperature at 3m" unit="degC" /> 148 148 149 149 <field id="sss" long_name="sea surface salinity" standard_name="sea_surface_salinity" unit="1e-3" /> 150 150 <field id="sss2" long_name="square of sea surface salinity" unit="1e-6" > sss * sss </field > … … 152 152 <field id="sssmin" long_name="min of sea surface salinity" field_ref="sss" operation="minimum" /> 153 153 <field id="sbs" long_name="sea bottom salinity" unit="0.001" /> 154 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 155 156 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 154 <field id="somint" long_name="vertical integral of salinity times density" standard_name="integral_wrt_depth_of_product_of_density_and_salinity" unit="(kg m2) x (1e-3)" /> 155 156 <field id="taubot" long_name="bottom stress module" unit="N/m2" /> 157 157 158 158 <!-- Case EOS = TEOS-10 : output potential temperature --> … … 295 295 <field id="us_y" long_name="j component of Stokes drift" unit="m/s" /> 296 296 </field_group> 297 297 298 298 <!-- SBC --> 299 299 <field_group id="SBC" > <!-- time step automaticaly defined based on nn_fsbc --> … … 311 311 <field id="precip" long_name="Total precipitation" standard_name="precipitation_flux" unit="kg/m2/s" /> 312 312 <field id="wclosea" long_name="closed sea empmr correction" standard_name="closea_empmr" unit="kg/m2/s" /> 313 313 314 314 <field id="qt" long_name="Net Downward Heat Flux" standard_name="surface_downward_heat_flux_in_sea_water" unit="W/m2" /> 315 315 <field id="qns" long_name="non solar Downward Heat Flux" unit="W/m2" /> … … 321 321 <field id="taum" long_name="wind stress module" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 322 322 <field id="wspd" long_name="wind speed module" standard_name="wind_speed" unit="m/s" /> 323 323 324 324 <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 325 325 <field id="ssh_ib" long_name="Inverse barometer sea surface height" standard_name="sea_surface_height_correction_due_to_air_pressure_at_low_frequency" unit="m" /> … … 369 369 <field id="taum_oce" long_name="wind stress module over open ocean" standard_name="magnitude_of_surface_downward_stress" unit="N/m2" /> 370 370 371 <!-- variables computed by the bulk parameterization algorithms (ln_blk) --> 372 <field id="Cd_oce" long_name="Drag coefficient over open ocean" standard_name="drag_coefficient_water" unit="" /> 373 <field id="Ce_oce" long_name="Evaporaion coefficient over open ocean" standard_name="evap_coefficient_water" unit="" /> 374 <field id="Ch_oce" long_name="Sensible heat coefficient over open ocean" standard_name="sensible_heat_coefficient_water" unit="" /> 375 <field id="theta_zt" long_name="Potential air temperature at z=zt" standard_name="potential_air_temperature_at_zt" unit="degC" /> 376 <field id="q_zt" long_name="Specific air humidity at z=zt" standard_name="specific_air_humidity_at_zt" unit="kg/kg" /> 377 <field id="theta_zu" long_name="Potential air temperature at z=zu" standard_name="potential_air_temperature_at_zu" unit="degC" /> 378 <field id="q_zu" long_name="Specific air humidity at z=zu" standard_name="specific_air_humidity_at_zu" unit="kg/kg" /> 379 <field id="ssq" long_name="Saturation specific humidity of air at z=0" standard_name="surface_air_saturation_spec_humidity" unit="kg/kg" /> 380 <field id="wspd_blk" long_name="Bulk wind speed at z=zu" standard_name="bulk_wind_speed_at_zu" unit="m/s" /> 381 <!-- ln_blk + key_si3 --> 382 <field id="Cd_ice" long_name="Drag coefficient over ice" standard_name="drag_coefficient_ice" unit="" /> 383 <field id="Ce_ice" long_name="Evaporaion coefficient over ice" standard_name="evap_coefficient_ice" unit="" /> 384 <field id="Ch_ice" long_name="Sensible heat coefficient over ice" standard_name="sensible_heat_coefficient_ice" unit="" /> 385 371 386 <!-- available key_oasis3 --> 372 387 <field id="snow_ao_cea" long_name="Snow over ice-free ocean (cell average)" standard_name="snowfall_flux" unit="kg/m2/s" /> … … 405 420 <!-- ice field (nn_ice=1) --> 406 421 <field id="ice_cover" long_name="Ice fraction" standard_name="sea_ice_area_fraction" unit="1" /> 407 422 408 423 <!-- dilution --> 409 424 <field id="emp_x_sst" long_name="Concentration/Dilution term on SST" unit="kg*degC/m2/s" /> 410 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 425 <field id="emp_x_sss" long_name="Concentration/Dilution term on SSS" unit="kg*1e-3/m2/s" /> 411 426 <field id="rnf_x_sst" long_name="Runoff term on SST" unit="kg*degC/m2/s" /> 412 427 <field id="rnf_x_sss" long_name="Runoff term on SSS" unit="kg*1e-3/m2/s" /> 413 428 414 429 <!-- sbcssm variables --> 415 430 <field id="sst_m" unit="degC" /> … … 422 437 423 438 </field_group> 424 439 425 440 426 441 </field_group> <!-- SBC --> 427 442 428 443 <!-- ABL --> 429 444 <field_group id="ABL" > <!-- time step automaticaly defined based on nn_fsbc --> … … 458 473 <field id="uz1_dta" long_name="DTA i-horizontal velocity" standard_name="dta_x_velocity" unit="m/s" /> 459 474 <field id="vz1_dta" long_name="DTA j-horizontal velocity" standard_name="dta_y_velocity" unit="m/s" /> 460 <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> 475 <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> 461 476 <field id="tz1_dta" long_name="DTA potential temperature" standard_name="dta_theta" unit="K" /> 462 477 <field id="qz1_dta" long_name="DTA specific humidity" standard_name="dta_qspe" unit="kg/kg" /> … … 464 479 <field id="uz1_geo" long_name="GEO i-horizontal velocity" standard_name="geo_x_velocity" unit="m/s" /> 465 480 <field id="vz1_geo" long_name="GEO j-horizontal velocity" standard_name="geo_y_velocity" unit="m/s" /> 466 <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> 481 <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> 467 482 </field_group> 468 483 469 484 </field_group> <!-- ABL --> 470 485 471 486 472 487 <!-- U grid --> 473 488 474 489 <field_group id="grid_U" grid_ref="grid_U_2D"> 475 490 <field id="e2u" long_name="U-cell width in meridional direction" standard_name="cell_width" unit="m" /> … … 480 495 <field id="uoce_e3u" long_name="ocean current along i-axis (thickness weighted)" unit="m/s" grid_ref="grid_U_3D" > uoce * e3u </field> 481 496 <field id="uoce_e3u_vsum" long_name="ocean current along i-axis * e3u summed on the vertical" field_ref="uoce_e3u" unit="m3/s" grid_ref="grid_U_vsum"/> 482 <field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field> 497 <field id="uocetr_vsum" long_name="ocean transport along i-axis summed on the vertical" field_ref="e2u" unit="m3/s"> this * uoce_e3u_vsum </field> 483 498 484 499 <field id="uocetr_vsum_op" long_name="ocean current along i-axis * e3u * e2u summed on the vertical" read_access="true" freq_op="1mo" field_ref="e2u" unit="m3/s"> @uocetr_vsum </field> 485 <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" /> 500 <field id="uocetr_vsum_cumul" long_name="ocean current along i-axis * e3u * e2u cumulated from southwest point" freq_offset="_reset_" operation="instant" freq_op="1mo" unit="m3/s" /> 486 501 <field id="msftbarot" long_name="ocean_barotropic_mass_streamfunction" unit="kg s-1" > uocetr_vsum_cumul * $rau0 </field> 487 502 … … 536 551 <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis" standard_name="ocean_salt_x_transport_due_to_diffusion" unit="1e-3*kg/s" /> 537 552 </field_group> 538 553 539 554 <!-- V grid --> 540 555 541 556 <field_group id="grid_V" grid_ref="grid_V_2D"> 542 557 <field id="e1v" long_name="V-cell width in longitudinal direction" standard_name="cell_width" unit="m" /> … … 595 610 <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis" standard_name="ocean_salt_y_transport_due_to_diffusion" unit="1e-3*kg/s" /> 596 611 </field_group> 597 612 598 613 <!-- W grid --> 599 614 600 615 <field_group id="grid_W" grid_ref="grid_W_3D"> 601 616 <field id="e3w" long_name="W-cell thickness" standard_name="cell_thickness" unit="m" /> 602 617 <field id="woce" long_name="ocean vertical velocity" standard_name="upward_sea_water_velocity" unit="m/s" /> 603 <field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field> 618 <field id="woce_e3w" long_name="ocean vertical velocity * e3w" unit="m2/s" > woce * e3w </field> 604 619 <field id="wocetr_eff" long_name="effective ocean vertical transport" unit="m3/s" /> 605 620 … … 611 626 612 627 <field id="avt" long_name="vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 613 <field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field> 628 <field id="avt_e3w" long_name="vertical heat diffusivity * e3w" unit="m3/s" > avt * e3w </field> 614 629 <field id="logavt" long_name="logarithm of vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 615 630 <field id="avm" long_name="vertical eddy viscosity" standard_name="ocean_vertical_momentum_diffusivity" unit="m2/s" /> 616 <field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field> 631 <field id="avm_e3w" long_name="vertical eddy viscosity * e3w" unit="m3/s" > avm * e3w </field> 617 632 618 633 <!-- avs: /= avt with ln_zdfddm=T --> 619 634 <field id="avs" long_name="salt vertical eddy diffusivity" standard_name="ocean_vertical_salt_diffusivity" unit="m2/s" /> 620 <field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field> 635 <field id="avs_e3w" long_name="vertical salt diffusivity * e3w" unit="m3/s" > avs * e3w </field> 621 636 <field id="logavs" long_name="logarithm of salt vertical eddy diffusivity" standard_name="ocean_vertical_heat_diffusivity" unit="m2/s" /> 622 637 623 638 <!-- avt_evd and avm_evd: available with ln_zdfevd --> 624 639 <field id="avt_evd" long_name="convective enhancement of vertical diffusivity" standard_name="ocean_vertical_tracer_diffusivity_due_to_convection" unit="m2/s" /> 625 <field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field> 640 <field id="avt_evd_e3w" long_name="convective enhancement to vertical diffusivity * e3w " unit="m3/s" > avt_evd * e3w </field> 626 641 <field id="avm_evd" long_name="convective enhancement of vertical viscosity" standard_name="ocean_vertical_momentum_diffusivity_due_to_convection" unit="m2/s" /> 627 642 … … 636 651 <field id="wstokes" long_name="Stokes Drift vertical velocity" standard_name="upward_StokesDrift_velocity" unit="m/s" /> 637 652 638 <!-- variables available with diaar5 --> 653 <!-- variables available with diaar5 --> 639 654 <field id="w_masstr" long_name="vertical mass transport" standard_name="upward_ocean_mass_transport" unit="kg/s" /> 640 655 <field id="w_masstr2" long_name="square of vertical mass transport" standard_name="square_of_upward_ocean_mass_transport" unit="kg2/s2" /> 641 656 642 657 </field_group> 643 658 644 659 <!-- F grid --> 645 660 <!-- AGRIF sponge --> … … 696 711 </field_group> 697 712 698 713 699 714 <!-- variables available with ln_floats --> 700 715 … … 711 726 <!-- variables available with iceberg trajectories --> 712 727 713 <field_group id="icbvar" domain_ref="grid_T" > 728 <field_group id="icbvar" domain_ref="grid_T" > 714 729 <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" /> 715 730 <field id="berg_melt_hcflx" long_name="icb heat flux to ocean due to melting heat content" unit="J/m2/s" /> … … 729 744 </field_group> 730 745 731 <!-- Poleward transport : ptr --> 732 <field_group id="diaptr" > 746 <!-- Poleward transport : ptr --> 747 <field_group id="diaptr" > 733 748 <field id="zomsf" long_name="Overturning Stream-Function : All basins" unit="Sv" grid_ref="grid_znl_W_3D" /> 734 749 <field id="zotem" long_name="Zonal Mean Temperature : All basins" unit="degree_C" grid_ref="grid_znl_T_3D" /> … … 738 753 <field id="sopstove" long_name="Overturning Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 739 754 <field id="sophtbtr" long_name="Barotropic Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 740 <field id="sopstbtr" long_name="Barotropic Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 755 <field id="sopstbtr" long_name="Barotropic Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 741 756 <field id="sophtadv" long_name="Advective Heat Transport: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 742 757 <field id="sopstadv" long_name="Advective Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> … … 744 759 <field id="sopstldf" long_name="Diffusive Salt Transport: All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 745 760 <field id="sophtvtr" long_name="Heat Transport : All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 746 <field id="sopstvtr" long_name="Salt Transport : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 761 <field id="sopstvtr" long_name="Salt Transport : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> 747 762 <field id="sophteiv" long_name="Heat Transport from mesoscale eddy advection: All basins" unit="PW" grid_ref="grid_znl_T_2D" /> 748 763 <field id="sopsteiv" long_name="Salt Transport from mesoscale eddy advection : All basins" unit="Giga g/s" grid_ref="grid_znl_T_2D" /> … … 760 775 761 776 762 <!-- 777 <!-- 763 778 ============================================================================================================ 764 779 Physical ocean model trend diagnostics : temperature, KE, PE, momentum … … 901 916 <field id="ketrd_zdf" long_name="ke-trend: vertical diffusion" unit="W/s^3" /> 902 917 <field id="ketrd_tau" long_name="ke-trend: wind stress " unit="W/s^3" grid_ref="grid_T_2D" /> 903 <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" /> 904 <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" /> 905 <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" /> 918 <field id="ketrd_bfr" long_name="ke-trend: bottom friction (explicit)" unit="W/s^3" /> 919 <field id="ketrd_bfri" long_name="ke-trend: bottom friction (implicit)" unit="W/s^3" /> 920 <field id="ketrd_atf" long_name="ke-trend: asselin time filter trend" unit="W/s^3" /> 906 921 <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" /> 907 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 922 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 908 923 909 924 <!-- variables available when explicit lateral mixing is used (ln_dynldf_OFF=F) --> 910 <field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 911 <field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 925 <field id="dispkexyfo" long_name="KE-trend: lateral mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_xy_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 926 <field id="dispkevfo" long_name="KE-trend: vertical mixing induced dissipation" standard_name="ocean_kinetic_energy_dissipation_per_unit_area_due_to_vertical_friction" unit="W/m^2" grid_ref="grid_T_2D" /> 912 927 <!-- variables available with ln_traadv_eiv=T and ln_diaeiv=T --> 913 <field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" /> 928 <field id="eketrd_eiv" long_name="EKE-trend due to parameterized eddy advection" standard_name="tendency_of_ocean_eddy_kinetic_energy_content_due_to_parameterized_eddy_advection" unit="W/m^2" grid_ref="grid_T_2D" /> 914 929 915 930 <!-- variables available with ln_PE_trd --> … … 928 943 <field id="petrd_bbc" long_name="pe-trend: geothermal heating" unit="W/m^3" /> 929 944 <field id="petrd_atf" long_name="pe-trend: asselin time filter" unit="W/m^3" /> 930 <field id="PEanom" long_name="potential energy anomaly" unit="1" /> 931 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" /> 932 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 945 <field id="PEanom" long_name="potential energy anomaly" unit="1" /> 946 <field id="alphaPE" long_name="partial deriv. of PEanom wrt T" unit="degC-1" /> 947 <field id="betaPE" long_name="partial deriv. of PEanom wrt S" unit="1e3" /> 933 948 </field_group> 934 949 … … 947 962 <field id="utrd_zdf" long_name="i-trend: vertical diffusion" unit="m/s^2" /> 948 963 <field id="utrd_tau" long_name="i-trend: wind stress " unit="m/s^2" grid_ref="grid_U_2D" /> 949 <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" /> 950 <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" /> 951 <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" /> 952 <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" /> 964 <field id="utrd_bfr" long_name="i-trend: bottom friction (explicit)" unit="m/s^2" /> 965 <field id="utrd_bfri" long_name="i-trend: bottom friction (implicit)" unit="m/s^2" /> 966 <field id="utrd_tot" long_name="i-trend: total momentum trend before atf" unit="m/s^2" /> 967 <field id="utrd_atf" long_name="i-trend: asselin time filter trend" unit="m/s^2" /> 953 968 </field_group> 954 969 … … 967 982 <field id="vtrd_zdf" long_name="j-trend: vertical diffusion" unit="m/s^2" /> 968 983 <field id="vtrd_tau" long_name="j-trend: wind stress " unit="m/s^2" grid_ref="grid_V_2D" /> 969 <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" /> 970 <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" /> 971 <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" /> 972 <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" /> 984 <field id="vtrd_bfr" long_name="j-trend: bottom friction (explicit)" unit="m/s^2" /> 985 <field id="vtrd_bfri" long_name="j-trend: bottom friction (implicit)" unit="m/s^2" /> 986 <field id="vtrd_tot" long_name="j-trend: total momentum trend before atf" unit="m/s^2" /> 987 <field id="vtrd_atf" long_name="j-trend: asselin time filter trend" unit="m/s^2" /> 973 988 </field_group> 974 989 975 990 976 <!-- 991 <!-- 977 992 ============================================================================================================ 978 993 Definitions for iodef_demo.xml … … 992 1007 <field field_ref="strd_zdfp_li" name="osaltdiff" /> 993 1008 </field_group> 994 1009 995 1010 <field_group id="mooring" > 996 1011 <field field_ref="toce" name="thetao" long_name="sea_water_potential_temperature" /> … … 1001 1016 <field field_ref="avt" name="difvho" long_name="ocean_vertical_heat_diffusivity" /> 1002 1017 <field field_ref="avm" name="difvmo" long_name="ocean_vertical_momentum_diffusivity" /> 1003 1018 1004 1019 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1005 1020 <field field_ref="sst2" name="tossq" long_name="square_of_sea_surface_temperature" /> … … 1049 1064 <field field_ref="BLT" name="blt" long_name="Barrier Layer Thickness" /> 1050 1065 </field_group> 1051 1066 1052 1067 <field_group id="groupU" > 1053 1068 <field field_ref="uoce" name="uo" long_name="sea_water_x_velocity" /> 1054 1069 <field field_ref="utau" name="tauuo" long_name="surface_downward_x_stress" /> 1055 1070 </field_group> 1056 1071 1057 1072 <field_group id="groupV" > 1058 1073 <field field_ref="voce" name="vo" long_name="sea_water_y_velocity" /> 1059 1074 <field field_ref="vtau" name="tauvo" long_name="surface_downward_y_stress" /> 1060 1075 </field_group> 1061 1076 1062 1077 <field_group id="groupW" > 1063 1078 <field field_ref="woce" name="wo" long_name="ocean vertical velocity" /> … … 1102 1117 </field_group> 1103 1118 1104 <!-- 1119 <!-- 1105 1120 ============================================================================================================ 1106 1121 --> 1107 <!-- output variables for my configuration (example) --> 1108 1122 <!-- output variables for my configuration (example) --> 1123 1109 1124 <field_group id="myvarOCE" > 1110 <!-- grid T --> 1125 <!-- grid T --> 1111 1126 <field field_ref="e3t" name="e3t" long_name="vertical scale factor" /> 1112 1127 <field field_ref="sst" name="tos" long_name="sea_surface_temperature" /> 1113 1128 <field field_ref="sss" name="sos" long_name="sea_surface_salinity" /> 1114 1129 <field field_ref="ssh" name="zos" long_name="sea_surface_height_above_geoid" /> 1115 1116 <!-- grid U --> 1130 1131 <!-- grid U --> 1117 1132 <field field_ref="e3u" name="e3u" long_name="vertical scale factor" /> 1118 1133 <field field_ref="ssu" name="uos" long_name="sea_surface_x_velocity" /> 1119 1120 <!-- grid V --> 1134 1135 <!-- grid V --> 1121 1136 <field field_ref="e3v" name="e3v" long_name="vertical scale factor" /> 1122 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 1123 </field_group> 1137 <field field_ref="ssv" name="vos" long_name="sea_surface_y_velocity" /> 1138 </field_group> 1124 1139 1125 1140 </field_definition> -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/SHARED/namelist_pisces_ref
r12377 r13159 352 352 ! 353 353 cn_dir = './' ! root directory for the location of the dynamical files 354 ln_ironsed = . true. ! boolean for Fe input from sediments355 ln_ironice = . true. ! boolean for Fe input from sea ice356 ln_hydrofe = . true. ! boolean for from hydrothermal vents354 ln_ironsed = .false. ! boolean for Fe input from sediments 355 ln_ironice = .false. ! boolean for Fe input from sea ice 356 ln_hydrofe = .false. ! boolean for from hydrothermal vents 357 357 sedfeinput = 2.e-9 ! Coastal release of Iron 358 358 distcoast = 5.e3 ! Distance off the coast for Iron from sediments -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/SHARED/namelist_ref
r12749 r13159 259 259 !----------------------------------------------------------------------- 260 260 ! ! bulk algorithm : 261 ln_NCAR = . false.! "NCAR" algorithm (Large and Yeager 2008)261 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 262 262 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 263 263 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) … … 271 271 rn_pfac = 1. ! multipl. factor for precipitation (total & snow) 272 272 rn_efac = 1. ! multipl. factor for evaporation (0. or 1.) 273 rn_vfac = 1. ! multipl. factor for ocean & ice velocity273 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 274 274 ! ! used to calculate the wind stress 275 275 ! ! (0. => absolute or 1. => relative winds) … … 277 277 ln_skin_wl = .false. ! use the warm-layer parameterization 278 278 ! ! ==> only available in ECMWF and COARE algorithms 279 ln_humi_sph = . false.! humidity "sn_humi" is specific humidity [kg/kg]279 ln_humi_sph = .true. ! humidity "sn_humi" is specific humidity [kg/kg] 280 280 ln_humi_dpt = .false. ! humidity "sn_humi" is dew-point temperature [K] 281 281 ln_humi_rlh = .false. ! humidity "sn_humi" is relative humidity [%] … … 1138 1138 ! ! = 3 as =2 with distinct dissipative an mixing length scale 1139 1139 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) 1140 nn_mxlice = 0 ! type of scaling under sea-ice 1141 ! = 0 no scaling under sea-ice 1142 ! = 1 scaling with constant sea-ice thickness 1143 ! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model ) 1144 ! = 3 scaling with maximum sea-ice thickness 1145 rn_mxlice = 10. ! max constant ice thickness value when scaling under sea-ice ( nn_mxlice=1) 1140 1146 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value 1141 1147 ln_drg = .false. ! top/bottom friction added as boundary condition of TKE … … 1398 1404 &namctl ! Control prints (default: OFF) 1399 1405 !----------------------------------------------------------------------- 1400 sn_cfctl%l_glochk = .FALSE. ! Range sanity checks are local (F) or global (T). Set T for debugging only 1401 sn_cfctl%l_allon = .FALSE. ! IF T activate all options. If F deactivate all unless l_config is T 1402 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1403 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1404 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1405 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1406 sn_cfctl%l_layout = .FALSE. ! 1407 sn_cfctl%l_prtctl = .FALSE. ! 1408 sn_cfctl%l_prttrc = .FALSE. ! 1409 sn_cfctl%l_oasout = .FALSE. ! 1410 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1411 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 1412 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 1413 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 1414 nn_print = 0 ! level of print (0 no extra print) 1415 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1416 nn_ictle = 0 ! end i indice of control sum multi processor runs 1417 nn_jctls = 0 ! start j indice of control over a subdomain) 1418 nn_jctle = 0 ! end j indice of control 1419 nn_isplt = 1 ! number of processors in i-direction 1420 nn_jsplt = 1 ! number of processors in j-direction 1421 ln_timing = .false. ! timing by routine write out in timing.output file 1422 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 1406 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1407 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1408 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1409 sn_cfctl%l_layout = .FALSE. ! 1410 sn_cfctl%l_prtctl = .FALSE. ! 1411 sn_cfctl%l_prttrc = .FALSE. ! 1412 sn_cfctl%l_oasout = .FALSE. ! 1413 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1414 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 1415 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 1416 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 1417 nn_print = 0 ! level of print (0 no extra print) 1418 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 1419 nn_ictle = 0 ! end i indice of control sum multi processor runs 1420 nn_jctls = 0 ! start j indice of control over a subdomain) 1421 nn_jctle = 0 ! end j indice of control 1422 nn_isplt = 1 ! number of processors in i-direction 1423 nn_jsplt = 1 ! number of processors in j-direction 1424 ln_timing = .false. ! timing by routine write out in timing.output file 1425 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 1423 1426 / 1424 1427 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/WED025/EXPREF/file_def_nemo-ice.xml
r11844 r13159 78 78 </file> 79 79 80 <file id="file22" name_suffix="_SBC_scalar" description="scalar variables" enabled=".true." >81 <!-- global contents -->82 <field field_ref="ibgvol_tot" grid_ref="grid_1point" name="ibgvol_tot" />83 <field field_ref="sbgvol_tot" grid_ref="grid_1point" name="sbgvol_tot" />84 <field field_ref="ibgarea_tot" grid_ref="grid_1point" name="ibgarea_tot" />85 <field field_ref="ibgsalt_tot" grid_ref="grid_1point" name="ibgsalt_tot" />86 <field field_ref="ibgheat_tot" grid_ref="grid_1point" name="ibgheat_tot" />87 <field field_ref="sbgheat_tot" grid_ref="grid_1point" name="sbgheat_tot" />88 89 <!-- global drifts (conservation checks) -->90 <field field_ref="ibgvolume" grid_ref="grid_1point" name="ibgvolume" />91 <field field_ref="ibgsaltco" grid_ref="grid_1point" name="ibgsaltco" />92 <field field_ref="ibgheatco" grid_ref="grid_1point" name="ibgheatco" />93 <field field_ref="ibgheatfx" grid_ref="grid_1point" name="ibgheatfx" />94 95 <!-- global forcings -->96 <field field_ref="ibgfrcvoltop" grid_ref="grid_1point" name="ibgfrcvoltop" />97 <field field_ref="ibgfrcvolbot" grid_ref="grid_1point" name="ibgfrcvolbot" />98 <field field_ref="ibgfrctemtop" grid_ref="grid_1point" name="ibgfrctemtop" />99 <field field_ref="ibgfrctembot" grid_ref="grid_1point" name="ibgfrctembot" />100 <field field_ref="ibgfrcsal" grid_ref="grid_1point" name="ibgfrcsal" />101 <field field_ref="ibgfrchfxtop" grid_ref="grid_1point" name="ibgfrchfxtop" />102 <field field_ref="ibgfrchfxbot" grid_ref="grid_1point" name="ibgfrchfxbot" />103 </file>104 105 80 </file_group> 106 81 -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/WED025/EXPREF/namelist_cfg
r12489 r13159 5 5 !! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 6 6 !! namsbc_sas, namtra_qsr, namsbc_rnf, 7 !! nam sbc_isf, namsbc_iscpl, namsbc_apr,7 !! namisf, namsbc_apr, 8 8 !! namsbc_ssr, namsbc_wave, namberg) 9 9 !! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) … … 38 38 nn_it000 = 1 ! first time step 39 39 nn_itend = 26280 ! last time step (std 5475) 40 nn_date0 = 19760301 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)40 nn_date0 = 20000101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 41 41 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 42 42 nn_rstctl = 2 ! restart control ==> activated only if ln_rstart=T … … 61 61 ln_tsd_init = .true. ! ocean initialisation 62 62 ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) 63 63 64 64 cn_dir = './' ! root directory for the T-S data location 65 !___________!_____________________ ____!___________________!___________!_____________!________!___________!__________________!__________!_______________!66 ! ! file name 67 ! ! 68 sn_tem = ' dta_temp_WED025' , -12 , 'votemper', .true., .true. , 'yearly' , '' , '' , ''69 sn_sal = ' dta_sal_WED025' , -12 , 'vosaline', .true., .true. , 'yearly' , '' , '' , ''65 !___________!_____________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 66 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 67 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 68 sn_tem = 'WED025_init_JRA_200001.nc', -12 , 'votemper', .false. , .true. , 'yearly' , '' , '' , '' 69 sn_sal = 'WED025_init_JRA_200001.nc', -12 , 'vosaline', .false. , .true. , 'yearly' , '' , '' , '' 70 70 / 71 71 !----------------------------------------------------------------------- … … 116 116 ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk ) 117 117 ! Sea-ice : 118 nn_ice = 2 ! =0 no ice boundary condition 118 nn_ice = 2 ! =0 no ice boundary condition 119 119 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 120 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 121 ! ! except in AGRIF zoom where it has to be specified 120 ! ! =2 or 3 for SI3 and CICE, respectively 122 121 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 123 122 ! ! =F levitating ice (no pressure, mass and salt exchanges) 124 123 ! Misc. options of sbc : 125 124 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr) 126 ln_dm2dc = . true.! daily mean to diurnal cycle on short wave125 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 127 126 ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 128 127 nn_fwb = 0 ! FreshWater Budget: =0 unchecked … … 141 140 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 142 141 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 143 ln_COARE_3p 5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)144 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)142 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 143 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 145 144 146 145 cn_dir = './' ! root directory for the bulk data location … … 148 147 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 149 148 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 150 sn_wndi = 'u10_ core' , 6 , 'U_10_MOD', .true. , .false. , 'yearly' , 'weights_bicubic_core.nc' , 'Uwnd' , ''151 sn_wndj = 'v10_ core' , 6 , 'V_10_MOD', .true. , .false. , 'yearly' , 'weights_bicubic_core.nc' , 'Vwnd' , ''152 sn_qsr = ' qsw_core' , 24 , 'SWDN_MOD', .false. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''153 sn_qlw = ' qlw_core' , 24 , 'LWDN_MOD', .false. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''154 sn_tair = 't10_ core' , 6 , 'T_10_MOD', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''155 sn_humi = 'q10_ core' , 6 , 'Q_10_MOD', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''156 sn_prec = 'precip_ core' , -1 , 'TPRECIP', .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''157 sn_snow = 'snow_ core' , -1 , 'SNOW' , .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''158 sn_slp = 'slp_ core' , 6 , 'SLP' , .true. , .false. , 'yearly' , 'weights_bilin_core.nc' , '' , ''149 sn_wndi = 'u10_JRA' , 3 , 'uas_10m' , .true. , .false. , 'yearly' , 'weights_bicubic_JRA.nc' , 'Uwnd' , '' 150 sn_wndj = 'v10_JRA' , 3 , 'vas_10m' , .true. , .false. , 'yearly' , 'weights_bicubic_JRA.nc' , 'Vwnd' , '' 151 sn_qsr = 'rsds_JRA' , 3 , 'rsds' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 152 sn_qlw = 'rlds_JRA' , 3 , 'rlds' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 153 sn_tair = 't10_JRA' , 3 , 'tas_10m' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 154 sn_humi = 'q10_JRA' , 3 , 'huss_10m', .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 155 sn_prec = 'precip_JRA' , 3 , 'prto' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 156 sn_snow = 'snow_JRA' , 3 , 'prsn' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 157 sn_slp = 'slp_JRA' , 3 , 'psl' , .true. , .false. , 'yearly' , 'weights_bilin_JRA.nc' , '' , '' 159 158 / 160 159 !----------------------------------------------------------------------- … … 201 200 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 202 201 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 203 sn_rnf = ' runoff_WED025' , -1 , 'runoff' , .true. , .false., 'yearly' , '' , '' , ''202 sn_rnf = 'WED025_icb' , -1 , 'runoff' , .true. , .false., 'yearly' , '' , '' , '' 204 203 / 205 204 !----------------------------------------------------------------------- … … 221 220 cn_isfcav_mlt = '3eq' ! ice shelf melting formulation (spe/2eq/3eq/oasis) 222 221 ! ! spe = fwfisf is read from a forcing field 223 ! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 )224 ! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 201 5)222 ! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 for a short description) 223 ! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 225 224 ! ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 226 225 ! ! cn_isfcav_mlt = 2eq or 3eq cases: 227 226 cn_gammablk = 'vel' ! scheme to compute gammat/s (spe,ad15,hj99) 228 ! ! ad15 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 229 ! ! hj99 = velocity and stability dependent Gamma (Holland et al. 1999) 230 rn_gammat0 = 1.4e-2 ! gammat coefficient used in blk formula 231 rn_gammas0 = 4.e-4 ! gammas coefficient used in blk formula 227 ! ! spe = constant transfert velocity (rn_gammat0, rn_gammas0) 228 ! ! vel = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 229 ! ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 230 rn_gammat0 = 1.4e-2 ! gammat coefficient used in spe, vel and vel_stab gamma computation method 231 rn_gammas0 = 4.0e-4 ! gammas coefficient used in spe, vel and vel_stab gamma computation method 232 232 ! 233 233 rn_htbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) … … 255 255 sn_isfpar_zmin = 'isfmlt_par', -12. , 'sozisfmin' , .false. , .true. , 'yearly' , '' , '' , '' 256 256 !* 'spe' and 'oasis' case 257 sn_isfpar_fwf = 'isfmlt_par' , -12. , 257 sn_isfpar_fwf = 'isfmlt_par' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 258 258 !* 'bg03' case 259 sn_isfpar_Leff = 'isfmlt_par', 0. , 259 sn_isfpar_Leff = 'isfmlt_par', 0. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 260 260 ! 261 261 ! ---------------- ice sheet coupling ------------------------------- … … 300 300 ln_tide = .true. ! Activate tides 301 301 ln_tide_pot = .false. ! use tidal potential forcing 302 clname(1) = 'M2' ! name of constituent - all tidal components must be set in namelist_cfg303 clname(2) = 'S2'304 clname(3) = 'K1'305 clname(4) = 'O1'302 sn_tide_cnames(1) = 'M2' ! name of constituent - all tidal components must be set in namelist_cfg 303 sn_tide_cnames(2) = 'S2' 304 sn_tide_cnames(3) = 'K1' 305 sn_tide_cnames(4) = 'O1' 306 306 / 307 307 !----------------------------------------------------------------------- … … 340 340 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 341 341 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 342 bn_ssh = ' bdyT_ssh_WED025' , -1 , 'sossheig' , .true. , .false., 'yearly' , '' , '' , ''343 bn_u2d = ' bdyU_u2d_WED025' , -1 , 'vobtcrtx' , .true. , .false., 'yearly' , '' , '' , ''344 bn_v2d = ' bdyV_u2d_WED025' , -1 , 'vobtcrty' , .true. , .false., 'yearly' , '' , '' , ''345 bn_u3d = ' bdyU_u3d_WED025' , -1 , 'vozocrtx' , .true. , .false., 'yearly' , '' , '' , ''346 bn_v3d = ' bdyV_u3d_WED025' , -1 , 'vomecrty' , .true. , .false., 'yearly' , '' , '' , ''347 bn_tem = ' bdyT_tra_WED025' , -1 , 'votemper' , .true. , .false., 'yearly' , '' , '' , ''348 bn_sal = ' bdyT_tra_WED025' , -1 , 'vosaline' , .true. , .false., 'yearly' , '' , '' , ''342 bn_ssh = 'WED025_bdyT_ssh' , -1 , 'sossheig' , .true. , .false., 'yearly' , '' , '' , '' 343 bn_u2d = 'WED025_bdyU_u2d' , -1 , 'vobtcrtx' , .true. , .false., 'yearly' , '' , '' , '' 344 bn_v2d = 'WED025_bdyV_u2d' , -1 , 'vobtcrty' , .true. , .false., 'yearly' , '' , '' , '' 345 bn_u3d = 'WED025_bdyU_u3d' , -1 , 'vozocrtx' , .true. , .false., 'yearly' , '' , '' , '' 346 bn_v3d = 'WED025_bdyV_u3d' , -1 , 'vomecrty' , .true. , .false., 'yearly' , '' , '' , '' 347 bn_tem = 'WED025_bdyT_tra' , -1 , 'votemper' , .true. , .false., 'yearly' , '' , '' , '' 348 bn_sal = 'WED025_bdyT_tra' , -1 , 'vosaline' , .true. , .false., 'yearly' , '' , '' , '' 349 349 !* for si3 350 bn_a_i = ' bdyT_ice_WED025' , -1 , 'ileadfra' , .true. , .false., 'yearly' , '' , '' , ''351 bn_h_i = ' bdyT_ice_WED025' , -1 , 'iicethic' , .true. , .false., 'yearly' , '' , '' , ''352 bn_h_s = ' bdyT_ice_WED025' , -1 , 'isnowthi' , .true. , .false., 'yearly' , '' , '' , ''350 bn_a_i = 'WED025_bdyT_ice' , -1 , 'ileadfra' , .true. , .false., 'yearly' , '' , '' , '' 351 bn_h_i = 'WED025_bdyT_ice' , -1 , 'iicethic' , .true. , .false., 'yearly' , '' , '' , '' 352 bn_h_s = 'WED025_bdyT_ice' , -1 , 'isnowthi' , .true. , .false., 'yearly' , '' , '' , '' 353 353 / 354 354 !----------------------------------------------------------------------- 355 355 &nambdy_tide ! tidal forcing at open boundaries (default: OFF) 356 356 !----------------------------------------------------------------------- 357 filtide = ' bdytide_WED025_' ! file name root of tidal forcing files357 filtide = 'WED025_bdytide_' ! file name root of tidal forcing files 358 358 / 359 359 … … 658 658 &namctl ! Control prints (default: OFF) 659 659 !----------------------------------------------------------------------- 660 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 661 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 662 sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 663 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 664 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 665 sn_cfctl%l_layout = .FALSE. ! 666 sn_cfctl%l_mppout = .FALSE. ! 667 sn_cfctl%l_mpptop = .FALSE. ! 668 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 669 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 670 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 671 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 672 nn_print = 0 ! level of print (0 no extra print) 673 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus 674 nn_ictle = 0 ! end i indice of control sum multi processor runs 675 nn_jctls = 0 ! start j indice of control over a subdomain) 676 nn_jctle = 0 ! end j indice of control 677 nn_isplt = 1 ! number of processors in i-direction 678 nn_jsplt = 1 ! number of processors in j-direction 679 ln_timing = .true. ! timing by routine write out in timing.output file 680 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 660 sn_cfctl%l_runstat = .true. ! switches and which areas produce reports with the proc integer settings. 661 ln_timing = .true. ! timing by routine write out in timing.output file 681 662 / 682 663 !----------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/cfgs/WED025/EXPREF/namelist_ice_cfg
r11487 r13159 26 26 &namitd ! Ice discretization 27 27 !------------------------------------------------------------------------------ 28 ln_cat_hfn = .true. ! ice categories are defined by a function following rn_himean**(-0.05) 29 rn_himean = 2.0 ! expected domain-average ice thickness (m) 30 rn_himin = 0.01 ! minimum ice thickness (m) used in remapping 28 31 / 29 32 !------------------------------------------------------------------------------ 30 33 &namdyn ! Ice dynamics 31 34 !------------------------------------------------------------------------------ 35 ln_landfast_L16 = .true. ! landfast: parameterization from Lemieux 2016 32 36 / 33 37 !------------------------------------------------------------------------------ … … 42 46 &namdyn_adv ! Ice advection 43 47 !------------------------------------------------------------------------------ 48 ln_adv_Pra = .false. ! Advection scheme (Prather) 49 ln_adv_UMx = .true. ! Advection scheme (Ultimate-Macho) 50 nn_UMx = 5 ! order of the scheme for UMx (1-5 ; 20=centered 2nd order) 44 51 / 45 52 !------------------------------------------------------------------------------ … … 62 69 &namthd_do ! Ice growth in open water 63 70 !------------------------------------------------------------------------------ 71 rn_hinew = 0.02 ! thickness for new ice formation in open water (m), must be larger than rn_himin 72 ln_frazil = .true. ! Frazil ice parameterization (ice collection as a function of wind) 64 73 / 65 74 !------------------------------------------------------------------------------ … … 70 79 &namthd_pnd ! Melt ponds 71 80 !------------------------------------------------------------------------------ 81 ln_pnd = .true. ! activate melt ponds or not 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 72 84 / 85 73 86 !------------------------------------------------------------------------------ 74 87 &namini ! Ice initialization 75 88 !------------------------------------------------------------------------------ 89 ln_iceini = .true. ! activate ice initialization (T) or not (F) 90 ln_iceini_file = .true. ! netcdf file provided for initialization (T) or not (F) 91 ! -- for ln_iceini_file = T 92 sn_hti = 'WED025_init_JRA_200001.nc', -12 ,'icethic_cea', .false. , .true., 'yearly' , '' , '', '' 93 sn_hts = 'WED025_init_JRA_200001.nc', -12 ,'icesnow_cea', .false. , .true., 'yearly' , '' , '', '' 94 sn_ati = 'WED025_init_JRA_200001.nc', -12 ,'ice_cover' , .false. , .true., 'yearly' , '' , '', '' 95 sn_smi = 'NOT USED' , -12 ,'smi' , .false. , .true., 'yearly' , '' , '', '' 96 sn_tmi = 'NOT USED' , -12 ,'tmi' , .false. , .true., 'yearly' , '' , '', '' 97 sn_tsu = 'NOT USED' , -12 ,'tsu' , .false. , .true., 'yearly' , '' , '', '' 98 sn_tms = 'NOT USED' , -12 ,'tms' , .false. , .true., 'yearly' , '' , '', '' 99 ! melt ponds (be careful, sn_apd is the pond concentration (not fraction), so it differs from rn_apd) 100 sn_apd = 'NOT USED' , -12 ,'apd' , .false. , .true., 'yearly' , '' , '', '' 101 sn_hpd = 'NOT USED' , -12 ,'hpd' , .false. , .true., 'yearly' , '' , '', '' 102 cn_dir='./' 76 103 / 77 104 !------------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.