Model:GIPL

From CSDMS
Revision as of 07:38, 2 September 2010 by WikiSysop (talk | contribs) (repaired layout)

Contact

Name Elchin Jafarov
Type of contact
Institute / Organization Univ. of Alaska Fairbanks
Postal address 1
Postal address 2
Town / City Fairbanks
Postal code 99775
State Alaska
Country US"US" is not in the list (Afghanistan, Albania, Algeria, Andorra, Angola, Antigua and Barbuda, Argentina, Armenia, Australia, Austria, ...) of allowed values for the "Country" property.
Email address eejafarov@alaska.edu
Phone
Fax



GIPL


Metadata

Summary

Also known as
Model type Single
Model part of larger framework
Note on status model
Date note status model

Technical specs

Supported platforms
Unix, Linux, Windows
Other platform
Programming language

Fortran90, Matlab

Other program language
Code optimized Single Processor, Parallel

Computing"Parallel </br>Computing" is not in the list (Single Processor, Multiple Processors) of allowed values for the "Code optimized" property.

Multiple processors implemented
Nr of distributed processors
Nr of shared processors
Start year development 2000
Does model development still take place? Yes
If above answer is no, provide end year model development
Code development status
When did you indicate the 'code development status'?
Model availability As executable"As executable" is not in the list (As code, As teaching tool) of allowed values for the "Model availability" property.
Source code availability
(Or provide future intension)
Through owner"Through owner" is not in the list (Through web repository, Through CSDMS repository) of allowed values for the "Source code availability" property.
Source web address
Source csdms web address
Program license type Other
Program license type other
Memory requirements
Typical run time it takes less than a minite to run the serial model for one with daily time interval


In/Output

Describe input parameters Upper Boundary (Air temperature)

Lower Boundary (Temperature gradient) Initial conditions (Temperature distribution at initial time) Thermo-physical properties

Input format ASCII
Other input format
Describe output parameters Temperature distribution with depth

Active Layer Depth Freezing/Thawing day

Output format ASCII
Other output format netcdf, GIS
Pre-processing software needed? Yes
Describe pre-processing software For spatial case one can developed its own pre-processing in order to put the input dataset in the format readable for GIPL.
Post-processing software needed? Yes
Describe post-processing software To generate netcdf or GIS outputs one can write its own converter for that.
Visualization software needed? Yes
If above answer is yes ESRI, Matlab
Other visualization software Matlab, Microsoft Excel (for serial); Matlab, ARCGIS, ncview (for spatial model)


Process

Describe processes represented by the model Main purpose of the model is to calculate subsurface temperature profile, active layer depth and freeze-up day.
Describe key physical parameters and equations Thermal capacities and conductivities prescribed for each subsurface layer, volumetric water content and unfrozen water coefficients.
Describe length scale and resolution constraints
Describe time scale and resolution constraints
Describe any numerical limitations and issues


Testing

Describe available calibration data sets We have tested the model for different permafrost observation sites for Alaska(USA) and Siberia(Russia). Typically, the model results show good correlation with measured data (if observations are accurate).
Upload calibration data sets if available: Media:Sample.zip
Describe available test data sets
Upload test data sets if available:
Describe ideal data for testing


Other

Do you have current or future plans for collaborating with other researchers?
Is there a manual available? No
Upload manual if available:
Model website if any
Model forum / discussion board
Comments


Introduction

GIPL(Geophysical Institute Permafrost Laboratory) is an implicit finite difference one-dimensional heat flow numerical model. The model was developed by V.Romanovsky and G. Tipenko at University of Alaska Fairbanks. First time the model was introduced at the 2004 AGU Conference Tipenko G. (2004). The model uses coarse vertical resolution grid which preserves the latent-heat effects in the phase transition zone, even under conditions of rapid or abrupt changes in the temperature fields. The air temperature is a major driving force for the GIPL model. The model includes upper boundary condition, constant geothermal heat flux at the lower boundary (typically from 500 to 1000 m) and temperature distribution at initial time. The other inputs are precipitation, prescribed water content and thermal properties of the multilayered soil column. As an output the model produces temperature distributions at different depths, active layer thickness and calculates time of freeze up. S. Marchenko and others in 2008 extend the model by developing pre-processing procedures which convert GIS format input data into GIPL format. E. Jafarov parallelized the GIPL code in order to run the model on supercomputers with finer grid resolution. First run of parallel UAF-GIPL2.0 model was in November 2009. The model is still under constant testing and development. However the preliminary results are available in the form of netcdf files. The preliminary results includes temperatures at different depths and active layer thickness. The detailed paper about the parallel model is under development.


Issues

Does not include convective heat transfer.

Visualization

File:Example.jpg

Input Files

There are 6 input files.

input.txt contains input parameters and input points. thermo.txt contains thermal and unfrozen water properties. grid.txt contains used for calculations grid. bound.txt includes upper boundary temperatures, which are usually air temperature. snow.txt and rsnow.txt contains snow depth and conductivity.

Output Files

The are two main output files.

results.txt contains results which are depend on time step prescribed by user in input.txt file (usually timestep = 1 day).

mean.txt contains results which are depend on mean time step, also assigned in input.txt


References

Tipenko G, Marchenko S, Romanovsky S, Groshev V, Sazonova T. 2004. Spatially distributed model of permafrost dynamics in Alaska. Eos Transactions AGU, 85(47): Fall Meet. Suppl., Abstract C12A-02.

Marchenko SS, Romanovsky VE, Tipenko G. 2008. Numerical modeling of spatial permafrost dynamics in Alaska. In Proceedings of the Ninth International Conference on Permafrost, 29 June–3 July 2008, Fairbanks, Alaska. Institute of Northern Engineering, University of Alaska, Fairbanks.