source: trunk/LMDZ.COMMON/libf/evolution/compute_tendencies_mod_slope.F90 @ 2835

Last change on this file since 2835 was 2835, checked in by romain.vande, 2 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: 2.1 KB
RevLine 
[2779]1!
2! $Id $
3!
4SUBROUTINE compute_tendencies_slope(tendencies_h2o_ice,min_h2o_ice_Y1,&
5     min_h2o_ice_Y2,iim_input,jjm_input,ngrid,tendencies_h2o_ice_phys,nslope)
6
7      IMPLICIT NONE
8
9!=======================================================================
10!
11!  Compute the tendencies of the evolution of water ice over the years
12!
13!=======================================================================
14
15!   arguments:
16!   ----------
17
18!   INPUT
19
20     INTEGER, intent(in) :: iim_input,jjm_input,ngrid  ,nslope                           ! # of grid points along longitude/latitude/ total
21     REAL, intent(in) , dimension(iim_input+1,jjm_input+1,nslope):: min_h2o_ice_Y1       ! LON x LAT field : minimum of water ice at each point for the first year
22     REAL, intent(in) , dimension(iim_input+1,jjm_input+1,nslope):: min_h2o_ice_Y2       ! LON x LAT field : minimum of water ice at each point for the second year
23
24!   OUTPUT
25     REAL, intent(out) , dimension(iim_input+1,jjm_input+1,nslope) :: tendencies_h2o_ice ! LON x LAT field : difference between the minima = evolution of perenial ice
26     REAL, intent(out) , dimension(ngrid,nslope)   :: tendencies_h2o_ice_phys            ! physical point field : difference between the minima = evolution of perenial ice
27
28!   local:
29!   ------
30
31     INTEGER :: i,j,ig0,islope                                                           ! loop variable
32
33!=======================================================================
34
35!  We compute the difference
36
37  DO j=1,jjm_input+1
38    DO i = 1, iim_input
39       DO islope = 1, nslope
40         tendencies_h2o_ice(i,j,islope)=min_h2o_ice_Y2(i,j,islope)-min_h2o_ice_Y1(i,j,islope)
41       enddo
42    ENDDO
43  ENDDO
44
45!  If the difference is too small; there is no evolution
46  DO j=1,jjm_input+1
47    DO i = 1, iim_input
48       DO islope = 1, nslope
49         if(abs(tendencies_h2o_ice(i,j,islope)).LT.1.0E-10) then
50            tendencies_h2o_ice(i,j,islope)=0.
51         endif
52       enddo
53    ENDDO
54  ENDDO
55
[2835]56  DO islope = 1,nslope
57    CALL gr_dyn_fi(1,iim_input+1,jjm_input+1,ngrid,tendencies_h2o_ice(:,:,islope),tendencies_h2o_ice_phys(:,islope))
[2779]58  ENDDO
59
60END SUBROUTINE compute_tendencies_slope
61
Note: See TracBrowser for help on using the repository browser.