Model help:GIPL: Difference between revisions
No edit summary |
No edit summary |
||
Line 127: | Line 127: | ||
|width="20%"|<span>The upper boundary file</span> | |width="20%"|<span>The upper boundary file</span> | ||
|width="60%"|<span>Time series of the temperature, usually air temperature. The default file include daily temperature distribution. First number correspond to number of variables. The right column correspond to number of day. The left column correspond to daily temperature.</span> | |width="60%"|<span>Time series of the temperature, usually air temperature. The default file include daily temperature distribution. First number correspond to number of variables. The right column correspond to number of day. The left column correspond to daily temperature.</span> | ||
|width="20%"|<span >-</span> | |width="20%"|<span >days/^oC</span> | ||
|-valign="top" | |||
|width="20%"|<span>The snow depth file</span> | |||
|width="60%"|<span>Time series of the snow depth. The default file include daily snow depth. First number correspond to number of variables. The right column correspond to number of day. The left column correspond to daily snow depth.</span> | |||
|width="20%"|<span >days/meters</span> | |||
|-valign="top" | |||
|width="20%"|<span>The snow conductivity file</span> | |||
|width="60%"|<span>Time series of the snow conductivity. The default file include daily snow conductivity. First number correspond to number of variables. The right column correspond to number of day. The left column correspond to daily snow conductivity.</span> | |||
|width="20%"|<span >days/W/(mK)</span> | |||
|-valign="top" | |||
|width="20%"|<span>The grid file</span> | |||
|width="60%"|<span>The default file include number of grid points from 1.5m above the surface to 33 m depth below and id number of the output points.</span> | |||
|width="20%"|<span >meters</span> | |||
|-valign="top" | |||
|width="20%"|<span>Initial temperature distribution file</span> | |||
|width="60%"|<span>The first number in the file correspond to point id, the second to the number of temperature points</span> | |||
|width="20%"|<span >meters</span> | |||
|} | |} | ||
Line 141: | Line 162: | ||
==Main equations== | ==Main equations== | ||
< | 1) Heat Flow Equation | ||
::::{| | |||
|width=1400px|<math> \frac{\partial H(x,t)}{\partial \tau}=\texttt{div}(k(x,t)\nabla t(x,\tau)) </math> | |||
|width=50p=x align="right"|(1) | |||
|} | |||
::::{| | |||
|width=1400px|<math> H(x,t)=\int\limits_0^tC(x,s)ds+L\Theta(x,t) </math> | |||
|width=50p=x align="right"|(2) | |||
|} | |||
::::{| | |||
|width=1400px|<math> \frac{\partial t(l_2,\tau)}{\partial x}=g</math> | |||
|width=50p=x align="right"|(3) | |||
|} | |||
::::{| | |||
|width=1400px|<math> t(x,0)=t_0(x)</math> | |||
|width=50p=x align="right"|(4) | |||
|} | |||
::::{| | |||
|width=1400px|<math> \Theta(x,t)=\eta(x)\cdot\begin{cases} 1 , & t\ge t_* \\ a|t|^{-b}, & t< t_* \end{cases} | |||
</math> | |||
|width=50p=x align="right"|(5) | |||
|} | |||
==Notes== | ==Notes== |
Revision as of 16:04, 17 August 2011
GIPL
GIPL is a heat flow with phase change model use to model permafrost thermal state and active layer depth, can be used for site specific and regional ground temperature distribution modeling.
Model introduction
GIPL(Geophysical Institute Permafrost Laboratory) is an implicit finite difference transient one-dimensional heat flow model. The model simulates ground temperature dynamics and the depth of the active layer by solving non-linear heat equation with phase change numerically. The model employs the Enthalpy method which does not require explicit treatment of the freeze/thaw moving boundary. In this model the process of freezing or thawing is occurring in accordance with unfrozen water content and soil thermal properties, and depends on the degree of soil saturation.
Model parameters
Uses ports
This will be something that the CSDMS facility will add
Provides ports
This will be something that the CSDMS facility will add
Main equations
1) Heat Flow Equation
<math> \frac{\partial H(x,t)}{\partial \tau}=\texttt{div}(k(x,t)\nabla t(x,\tau)) </math> (1)
<math> H(x,t)=\int\limits_0^tC(x,s)ds+L\Theta(x,t) </math> (2)
<math> \frac{\partial t(l_2,\tau)}{\partial x}=g</math> (3)
<math> t(x,0)=t_0(x)</math> (4)
<math> \Theta(x,t)=\eta(x)\cdot\begin{cases} 1 , & t\ge t_* \\ a|t|^{-b}, & t< t_* \end{cases} </math>
(5)
Notes
Any notes, comments, you want to share with the user
Numerical scheme
Examples
An example run with input parameters, BLD files, as well as a figure / movie of the output
Follow the next steps to include images / movies of simulations:
- Upload file: http://csdms.colorado.edu/wiki/Special:Upload
- Create link to the file on your page: [[Image:<file name>]].
See also: Help:Images or Help:Movies
Developer(s)
Name of the module developer(s)
References
Key papers
Links
Any link, eg. to the model questionnaire, etc.