GLM output files

11 views
Skip to first unread message

Elizabeth Kenneally

unread,
Mar 17, 2026, 10:30:58 PMMar 17
to Aquatic Ecosystem MOdelling Network (AEMON) discussion forum
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.
!/

Reply all
Reply to author
Forward
0 new messages