source: trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_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.9 KB
Line 
1!
2! $Id $
3!
4SUBROUTINE recomp_tend_co2_slope(tendencies_co2_ice_phys,tendencies_co2_ice_phys_ini,vmr_co2_gcm,vmr_co2_pem,ps_GCM_2,global_ave_press_GCM,global_ave_press_new,timelen,ngrid,nslope)
5
6      IMPLICIT NONE
7
8!=======================================================================
9!
10!  Routine that compute the evolution of the tendencie for co2 ice
11!
12!=======================================================================
13
14!   arguments:
15!   ----------
16
17!   INPUT
18  INTEGER, intent(in) :: timelen,ngrid,nslope
19  REAL, INTENT(in) ::  vmr_co2_gcm(ngrid,timelen)                ! physical point field : Volume mixing ratio of co2 in the first layer
20  REAL, INTENT(in) ::  vmr_co2_pem(ngrid,timelen)                ! physical point field : Volume mixing ratio of co2 in the first layer
21  REAL, intent(in) :: ps_GCM_2(ngrid,timelen)                 ! physical point field : Surface pressure in the GCM
22  REAL, intent(in) :: global_ave_press_GCM
23  REAL, intent(in) :: global_ave_press_new
24  REAL, intent(in) ::  tendencies_co2_ice_phys_ini(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
25
26!   OUTPUT
27  REAL, intent(inout) ::  tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
28
29!   local:
30!   ----
31
32  INTEGER :: i,t,islope
33  REAL :: eps, sigma, L, beta, alpha, coef, ave
34
35  eps=0.95
36  sigma=5.678E-8
37  L=5.71*10**5
38  beta=3182.48
39  alpha=23.3494
40
41  coef=669*24*3600*eps*sigma/L
42
43! Evolution of the water ice for each physical point
44  do i=1,ngrid
45    do islope=1,nslope
46      ave=0.
47      do t=1,timelen
48        ave=ave+(beta/(alpha-log(vmr_co2_gcm(i,t)*ps_GCM_2(i,t)/100)))**4  &
49              -(beta/(alpha-log(vmr_co2_pem(i,t)*ps_GCM_2(i,t)*global_ave_press_GCM/global_ave_press_new/100)))**4
50      enddo
51      tendencies_co2_ice_phys(i,islope)=tendencies_co2_ice_phys_ini(i,islope)-coef*ave/timelen
52    enddo
53  enddo
54
55
56END SUBROUTINE recomp_tend_co2_slope
Note: See TracBrowser for help on using the repository browser.