source: trunk/LMDZ.COMMON/libf/evolution/evol_co2_ice_s_mod_slope.F90 @ 2841

Last change on this file since 2841 was 2835, checked in by romain.vande, 3 years ago

Mars PEM:
Introduction of the possibility to follow an orbital forcing.
Introduction of new control parameters.
Cleaning of the PEM (removing unused files, add comments and new files)

A file named run_PEM.def can be added to the run.def. It contains the following variables:

_ evol_orbit_pem: Boolean. Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
_ year_bp_ini: Integer. Number of year before present to start the pem run if evol_orbit_pem=.true. , default=0
_ Max_iter_pem: Integer. Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
_ dt_pem: Integer. Time step of the PEM in year, default=1
_ alpha_criterion: Real. Acceptance rate of sublimating ice surface change, default=0.2
_ soil_pem: Boolean. Do you want to run with subsurface physical processes in the PEM? default=.true.

RV
File size: 1.5 KB
Line 
1!
2! $Id $
3!
4SUBROUTINE evol_co2_ice_s_slope(qsurf,tendencies_co2_ice_phys,&
5                             iim_input,jjm_input,ngrid,cell_area,STOPPING,nslope)
6
7  USE temps_mod_evol, ONLY: dt_pem
8
9      IMPLICIT NONE
10
11!=======================================================================
12!
13!  Routine that compute the evolution of the water ice
14!
15!=======================================================================
16
17!   arguments:
18!   ----------
19
20!   INPUT
21
22  INTEGER, intent(in) :: iim_input,jjm_input, ngrid,nslope   ! # of grid points along longitude/latitude/ total
23  REAL, intent(in) ::  cell_area(ngrid)
24
25!   OUTPUT
26  REAL, INTENT(INOUT) ::  qsurf(ngrid,nslope)                ! physical point field : Previous and actual density of water ice
27  LOGICAL :: STOPPING
28  REAL, intent(inout) ::  tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
29
30!   local:
31!   ----
32
33  INTEGER :: i,j,ig0,islope                                  ! loop variable
34!  REAL :: not_budget, budget
35  REAL :: pos_tend, neg_tend, real_coefficient,negative_part
36  REAL ::  new_tendencies(ngrid)
37
38  STOPPING=.false.
39
40
41! Evolution of the water ice for each physical point
42  do i=1,ngrid
43    do islope=1,nslope
44      qsurf(i,islope)=qsurf(i,islope)+tendencies_co2_ice_phys(i,islope)*dt_pem
45      if (qsurf(i,islope).lt.0) then
46        qsurf(i,islope)=0.
47        tendencies_co2_ice_phys(i,islope)=0.
48      endif
49    enddo
50  enddo
51
52
53END SUBROUTINE evol_co2_ice_s_slope
Note: See TracBrowser for help on using the repository browser.