This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:surfaceboundaryenthalpy [2020/05/15 13:33] adriengilbert |
solvers:surfaceboundaryenthalpy [2020/05/18 11:06] (current) adriengilbert |
||
---|---|---|---|
Line 4: | Line 4: | ||
* **Solver Fortran File:** '' | * **Solver Fortran File:** '' | ||
- | * **Solver Name:** '' | + | * **Solver Name:** '' |
- | * **Required Output Variable(s): | + | * **Required Output Variable(s): |
+ | * **Optional Output Variable(s): | ||
* **Required Input Variable(s): | * **Required Input Variable(s): | ||
* **Input Data:** Daily air temperature timeserie | * **Input Data:** Daily air temperature timeserie | ||
* **Optional Input Data: | * **Optional Input Data: | ||
+ | |||
+ | For vertically structured 3D mesh only. Works in serial and parallel. | ||
==== General Description ==== | ==== General Description ==== | ||
- | SurfBoundarySolver | + | SurfEnthBoundarySolver |
- | The solver uses the provided air temperature (and precipitation) daily record to compute the associated mean surface characteristic of the glacier over the time period covered by the provided data time-serie. It outputs | + | The solver uses the provided air temperature (and precipitation) daily record to compute the associated mean surface characteristic of the glacier over the time period covered by the provided data time-serie. It can output |
- '' | - '' | ||
Line 33: | Line 36: | ||
==== SIF contents ==== | ==== SIF contents ==== | ||
+ | |||
+ | The parameters are set in the constant section of the sif file : | ||
+ | |||
+ | < | ||
+ | Constants | ||
+ | |||
+ | rho_surf = real 350.0 ! Snow surface density | ||
+ | rho_ice = real 917.0 ! Ice density | ||
+ | rho_w = real 1000.0 ! Water density | ||
+ | Sr = real 0.005 ! Residual water saturation in Snow/Firn | ||
+ | T_ref_enthalpy = real 200.0 ! Use to compute Surf Enth (see Enthalpy solver) | ||
+ | L_heat = real 334000.0 | ||
+ | |||
+ | AirTemperatureFile = File " | ||
+ | PrecipFile = File " | ||
+ | | ||
+ | Precip = real 0.300 !Mean annual precipitation if PrecipFile not provided | ||
+ | TempCorrec= real -0.12 !Possibility of shifting temperature to get steady state mass balance for example (optional) | ||
+ | PrecipCorrec = real 1.0 !Possibility of adding a correcting factor on precipitation if PrecipFile provided (optional) | ||
+ | | ||
+ | GradTemp = real 0.0065 !Air temperature Lapse Rate (K m^-1) | ||
+ | GradPrecip= real 0.001 !Precipitation Lapse Rate (% m^-1) | ||
+ | z_temp = real 5310.0 !Elevation of temperature measurement from AirTemperatureFile (m) | ||
+ | z_precip = real 5310.0 !Elevation of Precipitation measurement (m) | ||
+ | | ||
+ | RadFact_ice = real 0.0000925 ! Melting factor for ice from radiation | ||
+ | RadFact_snow = real $0.0000925/ | ||
+ | Deg_jour = real 0.0114 ! Melting factor from air temperature | ||
+ | | ||
+ | seuil_precip = real 2.0 !Rain/ | ||
+ | seuil_fonte = real 0.0 !Melting air temperature threshold (degree C) | ||
+ | | ||
+ | firn_param = real 30.0 ! Firn densification factor (yr) | ||
+ | super_ice = real 0.15 ! Superimposed ice factor | ||
+ | | ||
+ | Latitude = real 28.82 !Latitude (degree) to compute Potential Solar Radiation | ||
+ | | ||
+ | !Possibility to export 1D profile simulation at one node of coordinate (X_output1D, | ||
+ | X_output1D = real x_coordinate | ||
+ | Y_output1D = real y_coordinate | ||
+ | | ||
+ | |||
+ | End | ||
+ | </ | ||
The solver needs output from the FlowDepth Solver : | The solver needs output from the FlowDepth Solver : | ||
Line 59: | Line 106: | ||
Equation = SurfBoundary | Equation = SurfBoundary | ||
- | Variable = Mass Balance | + | Variable = Surf Enth |
Variable DOFs = 1 | Variable DOFs = 1 | ||
- | procedure = " | + | procedure = " |
- | Exported Variable 1 = String "Surf Enth" | + | |
+ | |||
+ | | ||
Exported Variable 1 DOFs = 1 | Exported Variable 1 DOFs = 1 | ||
Line 92: | Line 141: | ||
End | End | ||
</ | </ | ||
+ | |||
+ | Boundary Condition to setup Dirichlet condition for the Enthalpy Solver: | ||
+ | |||
+ | < | ||
+ | ! Upper Surface | ||
+ | Boundary Condition 2 | ||
+ | Target Boundaries = 2 | ||
+ | |||
+ | Depth = real 0.0 | ||
+ | Enthalpy_h = Equals Surf Enth | ||
+ | |||
+ | End | ||
+ | </ | ||
+ | |||
+ | ==== Examples ==== | ||
+ | An example solving for the enthalpy within the Rika Samba Glacier using the SurfEnthBoundarySolver can be found in '' | ||
+ | |||
+ | ==== Reference ==== | ||
+ | |||
+ | Gilbert, A., Sinisalo, A., Gurung, T. R., Fujita, K., Maharjan, S. B., Sherpa, T. C., & Fukuda, T. (2020). The influence of water percolation through crevasses on the thermal regime of a Himalayan mountain glacier. The Cryosphere, 14(4), 1273–1288. https:// | ||