!------------------------------------------------------------------------------- ! Configuracion general del modelo !------------------------------------------------------------------------------- ! ! sim_name [string] titulo de la simulacion ! max_layers [integer] maximo numero de capas ! min_layer_vol [real] volumen minimo de la capa (m3 * 1000) ! min_layer_thick [real] minimo espesor de la capa (m) ! max_layer_thick [real] maximo espesor de la capa (m) ! Kw [real] atenuacion luminica de fondo (background light attenuation (m**-1) ! coef_inf_entrain [real] entrainment coefficient for inflows ! coef_mix_conv [real] eficiencia de mezclado- mezcla convectiva ! coef_wind_stir [real] eficiencia de mezclado- mezcla por viento ! coef_mix_turb [real] eficiencia de mezclado- efectos turbulentos de regimen variable ! coef_mix_shear [real] eficiencia de mezclado- generacion de esfuerzos de corte ! coef_mix_KH [real] eficiencia de mezclado- remolinos turbulentos hipolimneticos Kelvin-Helmholtz ! coef_mix_hyp [real] eficiencia de mezclado- turbulencia hipolimnetica !------------------------------------------------------------------------------- &glm_setup sim_name = 'GLM Simulation' max_layers = 500 min_layer_vol = 0.025 min_layer_thick = 0.50 max_layer_thick = 1.500 Kw = 0.6 coef_inf_entrain = 0. coef_mix_conv = 0.125 coef_wind_stir = 0.23 coef_mix_shear = 0.20 coef_mix_turb = 0.51 coef_mix_KH = 0.30 coef_mix_hyp = 0.5 / !------------------------------------------------------------------------------- ! Configuracion fabm !------------------------------------------------------------------------------- ! ode_method [integer] ! ODE esquema numerico para la dinamica de fuentes y depositions ! 1: primer-order explicito (no positivo) ! 2: segundo-order explicito Runge-Kutta (no positivo) ! 3: cuarto-order explicit Runge-Kutta (no positivo) ! 4: Patankar (primer-order, no conservativo) ! 5: Patankar-RK (segundo-order, no conservativo) ! 6: Patankar-RK (no funciona, no conservativo) ! 7: Patankar modificado (primer-order, conservat., posit.) ! 8: Patankar-RK modificado (segundo-order, conservat., posit.) ! 9: Patankar-RK modificado (no funciona, conservat., ! posit.) ! 10: Extended Modified Patankar (1st-order, conservat., ! posit.) ! 11: Extended Modified Patankar-RK (2nd-order, conservat., ! posit.) ! This variable is used only if bio_calc = True ! split_factor [integer, minimum = 1] ! number of biogeochemical time steps per physical time step ! bioshade_feedback [bool] ! retroalimentacion de bio-turbidez a la equacion de temperatura ######## ! repair_state [bool] ! Opcion de FABM para reparar las variables de estado que tienen that have -ve's####### ! fabm_nml_file [string] ! nombre del archivo fabm.nml para pasar a FABM ! multi_ben [bool] ! opcion especifica de GLM para FABM para hacer flujos benticos solo en la capa profunda, o en ! los costados de todas las capas (.true.) !------------------------------------------------------------------------------- &fabm_setup ode_method = 1 split_factor = 1 bioshade_feedback = .true. repair_state = .true. fabm_nml_file = 'fabm.nml' multi_ben = .true. / !------------------------------------------------------------------------------- !Detalles del lago !------------------------------------------------------------------------------- ! ! name [string] ! nombre del lago ! latitude [float, minimum = -90, maximum = 90, unit = deg North] ! latitud ! longitude [float, minimum = -360, maximum = 360, unit = deg East] ! longitud ! base_elev [float] ! elevacion de la base (m) ! crest_elev [float] ! elevacion de la cresta (m) ! bsn_len [float] ! largo del cuenco en la cresta (cuando el lago esta lleno) (m) ! bsn_wid [float] ! ancho del cuenco en la cresta (cuando el lago esta lleno) (m) ! bsn_vals [integer] ! numero de puntos de profundidad en la relacion altura-area ! H [float] ! altura (m) (comma separated list, len=bsn_vals) ! A [float] ! area (m2 * 1000) (comma separated list, len=bsn_vals) ! !------------------------------------------------------------------------------- &morphometry lake_name = 'AwesomeLake' latitude = 32 longitude = 35 base_elev = -252.9, crest_elev = -203.9, bsn_len = 21000, bsn_wid = 13000, bsn_vals = 45 ! H (m) A (m2 * 1000) H = -252.9, -251.9, -250.9, -249.9, -248.9, -247.9, -246.9, -245.9, -244.9, -243.9, -242.9, -241.9, -240.9, -239.9, -238.9, -237.9, -236.9, -235.9, -234.9, -233.9, -232.9, -231.9, -230.9, -229.9, -228.9, -227.9, -226.9, -225.9, -224.9, -223.9, -222.9, -221.9, -220.9, -219.9, -218.9, -217.9, -216.9, -215.9, -214.9, -213.9, -212.9, -211.9, -208.9, -207.9, -203.9 A = 0, 9250, 15200, 17875, 21975, 26625, 31700, 33950, 38250, 41100, 46800, 51675, 55725, 60200, 64675, 69600, 74475, 79850, 85400, 90975, 96400, 102000, 107000, 113000, 118000, 123000, 128000, 132000, 136000, 139000, 143000, 146000, 148000, 150000, 151000, 153000, 155000, 157000, 158000, 160000, 161000, 162000, 167000, 170000, 173000 V = 5425, 17650, 34187.5, 54112.5, 78412.5, 108000, 140000, 177000, 216000, 260000, 309000, 363000, 421000, 483000, 551000, 623000, 700000, 782000, 871000, 964000, 1060000, 1170000, 1280000, 1390000, 1510000, 1640000, 1770000, 1900000, 2040000, 2180000, 2330000, 2470000, 2620000, 2770000, 2930000, 3080000, 3230000, 3390000, 3550000, 3710000, 3870000, 4040000, 4300000 / !------------------------------------------------------------------------------- ! Duracion de la corrida !------------------------------------------------------------------------------- ! ! timefmt [integer] ! method to specify start and duration of model run ! 1: duration computed from number of time steps, MaxN (bogus start ! date used) [no longer implemented!!] ! 2: duration computed from given start and stop dates (number of time ! steps MaxN computed) ! 3: duration computed from number of time steps, MaxN (start date as ! specified, stop date computed) ! start [string, format = "yyyy-mm-dd hh:mm:ss"] ! nominal start date ! This variable is used only if timefmt != 1 ! stop [string, format = "yyyy-mm-dd hh:mm:ss"] ! nominal stop date ! This variable is used only if timefmt = 2 ! dt [float, minimum = 0.001, maximum = 86400, unit = s] ! Time step for integration ! numb_days [number of days to run the simulation ] ! This variable is used only if timefmt != 2 ! !------------------------------------------------------------------------------- &time timefmt = 2 start = '1997-01-02 00:00:00' !stop = '1999-05-02 00:00:00' stop = '2001-09-30 00:00:00' dt = 3600.0 num_days = 760 / !------------------------------------------------------------------------------- ! Formato para los archivos de salida (resultados) y nombres de archivos (format for output and filename(s) !------------------------------------------------------------------------------- ! ! out_dir [string] ! ruta para sacar/producir directorio (establecer permisos) ! out_fn [string] ! nombre del archivo de salida netcdf ! out_lkn [logical] ! producir el numero de lago en un archivo############# ! nsave [integer, minimum = 1, maximum = 86400] ! guardar resultados 'nsave' en la frequencia establecida############## ! csv_lake_fname [string] ! nombre del lago.csv. Informacion resumida de las simulacion diaria del lago. ##### ! csv_point_nlevs [integer] ! numero de profundidades a las cuales dejar un archivo de series de tiempo csv###### ! csv_point_at [real] ! altura desde el sedimento para los archivos csv de puntos (lista separada por comas)####### ! csv_point_fname [string] ! nombre de el/los archivos de salida csv (lista separada pro comas) ! csv_point_nvars [integer] ! numero de vaiables para sacar/producir en csv ! csv_point_vars [string] ! lista de nombres de variables de los archivos de resultados (output)-el orden es importante ! !------------------------------------------------------------------------------- &output out_dir = '.' out_fn = 'output' out_lkn = .false. nsave = 12 ! csv_lake_fname = 'lake' ! csv_point_nlevs = 1 csv_point_fname = 'WQ_' csv_point_at = 32. csv_point_nvars = 3 csv_point_vars = 'temp', 'salt', 'aed_oxygen_oxy', / !------------------------------------------------------------------------------- ! Perfiles de condiciones iniciales !------------------------------------------------------------------------------- ! ! lake_depth [float] profundidad incial del lago (m) ! num_depths [integer] numero de profundidad provisto por perfiles iniciales ! the_depths [float] las profundidades de los puntos de perfiles iniciales (m) ! the_temps [float] Temperatura de los puntos de pefiles iniciales. ! the_sals [float] la salinidad de los puntos de perfiles iniciales(psu) ! num_wq_vars [integer] numero de variable non GLM (e.g FABM) para ser inicializado ! wq_names [string] nombre de variables non GLM (e.g. FABM) para inicializarse####### ! wq_init_vals [float] arreglo de variables (fila = variables; columna = profundidad.) ! !------------------------------------------------------------------------------- &init_profiles lake_depth = 41.0 num_depths = 5 the_depths = 1,10,20,30,40,46 ! GLM the_temps = 18.00,18.00,18.00, 18., 18,18. the_sals = 0.5, 0.5, 0.5, 0.5, 0.5, 0.5 ! FABM num_wq_vars = 6 wq_names = 'aed_organic_matter_don', 'aed_organic_matter_pon', 'aed_organic_matter_dop', 'aed_organic_matter_pop', 'aed_organic_matter_doc', 'aed_organic_matter_poc' wq_init_vals = 1.1, 1.2, 1.3, 1.2, 1.3, 2.1, 2.2, 2.3, 1.2, 1.3, 3.1, 3.2, 3.3, 1.2, 1.3, 4.1, 4.2, 4.3, 1.2, 1.3, 5.1, 5.2, 5.3, 1.2, 1.3, 6.1, 6.2, 6.3, 1.2, 1.3 / !------------------------------------------------------------------------------- ! meteorologia !------------------------------------------------------------------------------- ! ! met_sw [bool] switch to include surface meteorological forcing ! lw_type [string] type of longwave data supplied (LW_IN/LW_CC/LW_NET) ! rain_sw [bool] include rainfall nutrient composition ! snow_sw [bool] include snowfall (m/d) ! atm_stab [bool] account for non-neutral atmospheric stability (not possible yet} ! meteo_fl [string] name of file with meteorology input data ! wind_factor [float] wind multiplication factor (-) ! ce [float] bulk aerodynamic coefficient for latent heat transfer ! ch [float] bulk aerodynamic coefficient for sensible heat transfer ! coef_wind_drag [float] bulk aerodynamic coefficient for transfer of momentum ! !------------------------------------------------------------------------------- &meteorology met_sw = .true. lw_type = 'LW_IN' rain_sw = .false. snow_sw = .false. atm_stab = .false. subdaily = .true. meteo_fl = 'met_hourly.csv' wind_factor = 1.0 ce = 0.0013 ch = 0.0013 coef_wind_drag = 0.0013 / !------------------------------------------------------------------------------- !Caudales de entrada !------------------------------------------------------------------------------- ! ! num_inflows [integer] numberos de rios tributarios(0+) ! names_of_strms [string] nombre de las entrada (comma separated list) ! strm_hf_angle [float] angulo de la mitad del rio (grados) ############################## ! strmbd_slope [float] streambed slope (degrees)################## ! strmbd_drag [float] streambed drag coefficient (-)###################### ! inflow_factor [float] multipicador de tasa de caudal de entrada(-) ! inflow_fl [string] nombres de los archivsod de caudales de entrada (lista con coma para separar) ! inflow_varnum [integer] numero de columnas a ser leidas (exlcluyendo las afechas) ! inflow_vars [string] nombre de las variables de las columnas de los archivos de tributarios ! Esto debe ser una lista perarada por comas y debe incluir ! incluir CAUDAL, SALINITY & TEMPERATURA (para GLM), y ! opcionalmente puede incluir los nombres de las variables de FABM ! !------------------------------------------------------------------------------- &inflow num_inflows = 2 names_of_strms = 'Riv1','Riv2' strm_hf_angle = 65.0, 65.0 strmbd_slope = 2.0, 2.0 strmbd_drag = 0.0160, 0.0160 inflow_factor = 1.0, 1.0 inflow_fl = 'inflow_1.csv', 'inflow_2.csv' inflow_varnum = 3 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' / !------------------------------------------------------------------------------- ! outflows !------------------------------------------------------------------------------- ! ! num_outlet [integer] numero de salidas ! flt_off_sw [bool] cambios en la toma de agua flotante. ! outl_elvs [float] alturas de las salida (comma separated list) ! bsn_len_outl [float] basin length at outlets (m)########################## ! bsn_wid_outl [float] basin width at outlets (m)########################### ! outflow_fl [string] archivos de datos de los caudales de salida ! outflow_factor [float] multiplicador de tasa de caudal de salida(-) ! !------------------------------------------------------------------------------- &outflow num_outlet = 1 flt_off_sw = .false. outl_elvs = -215.5 bsn_len_outl = 799 bsn_wid_outl = 399 outflow_fl = 'outflow.csv' outflow_factor = 0.8 /