source: trunk/LMDZ.COMMON/libf/evolution/compute_tendencies_slope_mod.F90 @ 3130

Last change on this file since 3130 was 3130, checked in by jbclement, 19 months ago

PEM:
The perennial co2 ice is now taken into account with co2 frost (qsurf) to compute the tendency and to make the update + Rework of how co2 frost is converted to perennial co2 ice at the end of the PEM run + Correction of the value of 'threshold_co2_frost2perennial' to correspond to 10 m + Perennial co2 ice is now handled outside 'paleoclimate' in "phyetat0_mod.F90" of the Mars PCM + Some cleanings.

/!\ Commit for the PEM management of co2 ice before a rework of ice management in the PEM!
JBC
File size: 1.4 KB
Line 
1MODULE compute_tendencies_slope_mod
2
3implicit none
4
5!=======================================================================
6contains
7!=======================================================================
8
9SUBROUTINE compute_tendencies_slope(ngrid,nslope,min_ice_Y1,min_ice_Y2,tendencies_ice)
10
11implicit none
12
13!=======================================================================
14!
15!  Compute the tendencies of the evolution of water ice over the years
16!
17!=======================================================================
18
19!   arguments:
20!   ----------
21!   INPUT
22integer,                       intent(in) :: ngrid, nslope ! # of grid points along longitude/latitude/ total
23real, dimension(ngrid,nslope), intent(in) :: min_ice_Y1    ! LON x LAT field : minimum of water ice at each point for the first year
24real, dimension(ngrid,nslope), intent(in) :: min_ice_Y2    ! LON x LAT field : minimum of water ice at each point for the second year
25
26!   OUTPUT
27real, dimension(ngrid,nslope), intent(out) :: tendencies_ice ! physical point field : difference between the minima = evolution of perennial ice
28!=======================================================================
29
30!  We compute the difference
31tendencies_ice = min_ice_Y2 - min_ice_Y1
32
33!  If the difference is too small; there is no evolution
34where (abs(tendencies_ice) < 1.e-10) tendencies_ice = 0.
35
36END SUBROUTINE compute_tendencies_slope
37
38END MODULE compute_tendencies_slope_mod
39
Note: See TracBrowser for help on using the repository browser.