Warning: Undefined array key 1 in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 78

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 431

Warning: Undefined array key "DOKU_PREFS" in /home/np29546/public_html/elmerice/wiki/inc/common.php on line 2082

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/common.php on line 2090

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
solvers:surfaceboundaryenthalpy [Elmer/Ice Wiki]

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Last revision Both sides next revision
solvers:surfaceboundaryenthalpy [2020/05/15 13:44]
adriengilbert
solvers:surfaceboundaryenthalpy [2020/05/15 15:48]
adriengilbert
Line 4: Line 4:
  
   * **Solver Fortran File:** ''SurfaceBoundaryEnthalpy.f90''   * **Solver Fortran File:** ''SurfaceBoundaryEnthalpy.f90''
-  * **Solver Name:** ''SurfBoundarySolver''+  * **Solver Name:** ''SurfEnthBoundarySolver''
   * **Required Output Variable(s):** ''Surf Enth'', ''Mass Balance'', ''Densi'', ''Firn'', ''Melting'', ''Refreeze'', ''Accu'', ''Rad_Fact'', ''Rain'', ''PotRad''   * **Required Output Variable(s):** ''Surf Enth'', ''Mass Balance'', ''Densi'', ''Firn'', ''Melting'', ''Refreeze'', ''Accu'', ''Rad_Fact'', ''Rain'', ''PotRad''
   * **Required Input Variable(s):** ''Depth'' and ''SurfGrad'' from the FlowDepth Solver   * **Required Input Variable(s):** ''Depth'' and ''SurfGrad'' from the FlowDepth Solver
Line 15: Line 15:
 ==== General Description ==== ==== General Description ====
  
-SurfBoundarySolver is a pseudo-solver which compute surface mass balance and Dirichlet surface boundary condition for the Enthalpy solver. It takes into account firn heating processes by solving vertical melt-water percolation and refreezing.+SurfEnthBoundarySolver is a pseudo-solver which compute surface mass balance and Dirichlet surface boundary condition for the Enthalpy solver. It takes into account firn heating processes by solving vertical melt-water percolation and refreezing.
  
 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 following variables: 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 following variables:
Line 35: Line 35:
  
 ==== SIF contents ==== ==== SIF contents ====
 +
 +The parameters are set in the constant section of the sif file :
 +
 +<code>
 +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        ! Latent heat of fusion
 +
 +  AirTemperatureFile = File "YourTempFile.dat"  ! Contain daily temperature data
 +  PrecipFile = File "YourPrecipFile.dat"  ! Contain daily precipitation data (optional)
 +    
 +  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/2.0 ! Melting factor for snow from radiation
 +  Deg_jour = real 0.0114 ! Melting factor from air temperature
 +  
 +  seuil_precip = real 2.0 !Rain/Snow air temperature threshold (degree C)
 +  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,Y_output1D) (optional)
 +  X_output1D = real x_coordinate
 +  Y_output1D = real y_coordinate
 +  
 +
 +End
 +</code>
  
 The solver needs output from the FlowDepth Solver : The solver needs output from the FlowDepth Solver :
Line 61: Line 105:
  
   Equation = SurfBoundary   Equation = SurfBoundary
-  Variable = Mass Balance+  Variable = Surf Enth
   Variable DOFs = 1   Variable DOFs = 1
-  procedure =  "ElmerIceSolvers" "SurfBoundarySolver"+  procedure =  "ElmerIceSolvers" "SurfEnthBoundarySolver"
  
-  Exported Variable 1 = String "Surf Enth"+  Exported Variable 1 = String "Mass Balance"
   Exported Variable 1 DOFs = 1   Exported Variable 1 DOFs = 1
  
Line 109: Line 153:
  
 ==== Examples ==== ==== Examples ====
-An example solving for the enthalpy within the Rika Samba Glacier using the SurfBoundarySolver can be found in ''[ELMER_TRUNK]/elmerice/Tests/SurfaceBoundaryEnth''.+An example solving for the enthalpy within the Rika Samba Glacier using the SurfEnthBoundarySolver can be found in ''[ELMER_TRUNK]/elmerice/examples/Test_SurfaceBoundaryEnth''.
  
 ==== Reference ==== ==== Reference ====
solvers/surfaceboundaryenthalpy.txt ยท Last modified: 2020/05/18 11:06 by adriengilbert
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0