Hello all,
I am new to GLM and completing my thesis on it looking at a wetland. I have some observed data which I am looking to compare my results with to calibrate the model. I am having some issues getting my output files I produce in smaller increments.
I understand the standard 'lake' file is daily, but I was wondering if there was any way to get the outflow in 6 minute increments and the lake level. Currently my outlet, outflow and lake files are all daily (dt = 360, nsave =1). Below is my code, can anyone assist? Thank you in advance!
!-------------------------------------------------------------------------------
! general model setup
!-------------------------------------------------------------------------------
&glm_setup
sim_name = 'GLM Simulation'
max_layers = 50
min_layer_vol = 0.01
min_layer_thick = 0.05
max_layer_thick = 0.25
density_model = 1
/
&mixing
surface_mixing = 1
coef_mix_conv = 0.02
coef_wind_stir = 0.35
coef_mix_shear = 0.2
coef_mix_turb = 0.51
coef_mix_KH = 0.3
deep_mixing = 2
coef_mix_hyp = 0.5
diff = 0.
/
&morphometry
lake_name = 'Wetland'
latitude = -37.997
longitude = 145.289
base_elev = 24.25,
crest_elev = 25.30,
bsn_vals = 23,
! H (m) A (m2)
A = 116.916,1750.97,1908.504,2068.312,2267.867,2738.112,3058.18,3389.616,3962.041,4620.513,5281.193,5977.344,6805.159,7688.316,8595.711,9527.993,11388.835,14283.868,17177.248,19904.432,25843.674,28619.239,30311.492
H = 24.25, 24.3, 24.35, 24.4, 24.45, 24.5, 24.55, 24.6, 24.65, 24.7, 24.75, 24.8, 24.85, 24.9, 24.95, 25, 25.05, 25.1, 25.15, 25.2, 25.25, 25.3, 25.35
/
&time
timefmt = 3
start = '2025-05-20 17:30:00'
!stop = '2025-05-31 17:30:00'
dt = 60.0
num_days = 11
/
&output
out_dir = 'output'
out_fn = 'output'
nsave = 6
! General summary file
csv_lake_fname = 'lake'
! Depth specific outputs
csv_point_nlevs = 1
csv_point_fname = 'WQ_'
csv_point_at = 2.
csv_point_nvars = 6
csv_point_vars = 'temp',
'salt',
'TRC_SS1'
'OXY_oxy',
'Uorb'
'Taub'
'SIL_rsi',
'NIT_amm',
'NIT_nit',
'PHS_frp',
'OGM_don',
'OGM_pon',
'OGM_dop',
'OGM_pop',
'OGM_doc',
'OGM_poc',
'PHY_green',
'PHY_crypto',
'PHY_diatom'
! Combined outlet file & overflow
!csv_outlet_allinone = .true.
csv_outlet_fname = 'outlet_'
csv_outlet_nvars = 17
csv_outlet_vars = 'flow',
'temp',
'salt',
'OXY_oxy',
'SIL_rsi',
'NIT_amm',
'NIT_nit',
'PHS_frp',
'OGM_don',
'OGM_pon',
'OGM_dop',
'OGM_pop',
'OGM_doc',
'OGM_poc',
'PHY_green',
'PHY_crypto',
'PHY_diatom'
csv_ovrflw_fname = "overflow"
/
&init_profiles
lake_depth = 1.1
num_depths = 2
the_depths = 0.5,0.7
! GLM
the_temps = 5.10,5.10
the_sals = 0, 0
! WQ
num_wq_vars = 2
wq_names = 'OXY_oxy',
'TRC_ss1'
wq_init_vals = 200.1, 200.2
6.2, 6.2
/
&meteorology
met_sw = .TRUE.
lw_type = 'LW_CC'
rain_sw = .false.
atm_stab = .false.
!catchrain = .true.
!rad_mode = 1
!albedo_mode = 1
cloud_mode = 1
!-- BC file details
subdaily = .true.
meteo_fl = 'bcs/met_200525_310525.csv'
wind_factor = 1.0
sw_factor = 1.0
lw_factor = 1.0
at_factor = 1.0
rh_factor = 1.0
rain_factor = 1
!-- Parameters
ce = 0.0013
ch = 0.0013
cd = 0.0013
!rain_threshold = 0.001
!runoff_coef = 0.0
/
&light
! light_mode = 0
! n_bands = 4
! light_extc = 1.0, 0.5, 2.0, 4.0
! energy_frac = 0.51, 0.45, 0.035, 0.005
! Benthic_Imin = 0.10
Kw = 0.8
/
&bird_model
! AP = 973
! Oz = 0.279
! WatVap = 1.1
! AOD500 = 0.033
! AOD380 = 0.038
! Albedo = 0.2
/
&inflow
num_inflows = 1
names_of_strms = 'Riv1'
!subm_flag = .false.,.false.
strm_hf_angle = 65.0
strmbd_slope = 2.0
strmbd_drag = 0.0160
inflow_factor = 0.5
inflow_fl = 'bcs/inflow_200525_310525.csv'
inflow_varnum = 2
inflow_vars = 'FLOW',
!'TEMP',
'SALT',
'aed_oxygen_oxy',
'aed_silica_rsi',
'aed_nitrogen_amm',
'aed_nitrogen_nit',
'aed_phosphorus_frp',
'aed_organic_matter_don',
'aed_organic_matter_pon',
'aed_organic_matter_dop',
'aed_organic_matter_pop',
'aed_organic_matter_doc',
'aed_organic_matter_poc',
'aed_phytoplankton_green',
'aed_phytoplankton_crypto',
'aed_phytoplankton_diatom'
coef_inf_entrain = 0.
! time_fmt = 'YYYY-MM-DD hh:mm:ss'
/
!-------------------------------------------------------------------------------
&outflow
num_outlet = 1
!flt_off_sw = .false.
outlet_type = 1
outl_elvs = 25.30
bsn_len_outl = 0.35
bsn_wid_outl = 0.5
outflow_fl = 'bcs/outflow_200525_310525.csv'
outflow_factor = 1.0
!seepage = .true.
!seepage_rate = 0.01
/
!&sediment
! sed_heat_Ksoil = 0.0
! sed_temp_depth = 0.2
! sed_temp_mean = 5,10,20
! sed_temp_amplitude = 6,8,10
! sed_temp_peak_doy = 80, 70, 60
! benthic_mode = 2
! n_zones = 3
! zone_heights = 10., 40., 90.
! sed_reflectivity = 0.1, 0.01, 0.01
! sed_roughness = 0.1, 0.01, 0.01
!/
!&snowice
! snow_albedo_factor = 1.0
! snow_rho_max = 300
! snow_rho_min = 50
!/
!&debugging
!disable_evap = .true.
!/