source: trunk/LMDZ.COMMON/libf/evolution/info_run_PEM.F90 @ 2962

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
RevLine 
[2886]1SUBROUTINE info_run_PEM(year_iter, criterion_stop)
2
[2888]3   IMPLICIT NONE
[2886]4
5  INTEGER, intent(in) :: year_iter, criterion_stop                  ! # of year and reason to stop
6
[2888]7  logical :: exist
[2886]8
9
[2888]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)
[2886]18
19END SUBROUTINE info_run_PEM
Note: See TracBrowser for help on using the repository browser.