|
Last change
on this file since 2962 was
2888,
checked in by llange, 22 months ago
|
|
PEM
- Following r-2886-5, debugging some issues ( conservation of H2O, water_reservoir not initialized, info_pem which was not working)
- Cleaning of some redundant routines (e.g., stopping criterion of the PEM) and variables renamed (e.g., alpha_criterion now transofrmed in ice_criterion)
- Ice table subroutine is now in a module and recalled "compute_ice_table_equilibrium" (v.s. dynamic ice table, efforts ongoing)
- Abort_pem introduced to avoid just stopping the pem with raw "stop" in the codes
- conf_pem has been improved so that values are not hard coded (e.g., geothermal flux, mass of water_reservoir)
- Regolith thermal properties updated in a cleaner way
(Not atomic commit, sorry)
LL & RV
|
|
File size:
510 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE info_run_PEM(year_iter, criterion_stop) |
|---|
| 2 | |
|---|
| 3 | IMPLICIT NONE |
|---|
| 4 | |
|---|
| 5 | INTEGER, intent(in) :: year_iter, criterion_stop ! # of year and reason to stop |
|---|
| 6 | |
|---|
| 7 | logical :: exist |
|---|
| 8 | |
|---|
| 9 | |
|---|
| 10 | inquire(file='info_run_PEM.txt', exist=exist) |
|---|
| 11 | if (exist) then |
|---|
| 12 | open(12, file='info_run_PEM.txt', status="old", position="append", action="write") |
|---|
| 13 | else |
|---|
| 14 | open(12, file='info_run_PEM.txt', status="new", action="write") |
|---|
| 15 | end if |
|---|
| 16 | write(12, *) year_iter, criterion_stop |
|---|
| 17 | close(12) |
|---|
| 18 | |
|---|
| 19 | END SUBROUTINE info_run_PEM |
|---|
Note: See
TracBrowser
for help on using the repository browser.
