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

Last change on this file since 3136 was 3130, checked in by jbclement, 2 years 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
RevLine 
[3076]1MODULE compute_tendencies_slope_mod
2
3implicit none
4
5!=======================================================================
6contains
7!=======================================================================
8
[3070]9SUBROUTINE compute_tendencies_slope(ngrid,nslope,min_ice_Y1,min_ice_Y2,tendencies_ice)
[2779]10
[3070]11implicit none
[2779]12
13!=======================================================================
14!
15!  Compute the tendencies of the evolution of water ice over the years
16!
17!=======================================================================
18
19!   arguments:
20!   ----------
21!   INPUT
[3070]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
[2779]25
26!   OUTPUT
[3130]27real, dimension(ngrid,nslope), intent(out) :: tendencies_ice ! physical point field : difference between the minima = evolution of perennial ice
[2779]28!=======================================================================
29
30!  We compute the difference
[3070]31tendencies_ice = min_ice_Y2 - min_ice_Y1
[2779]32
33!  If the difference is too small; there is no evolution
[3070]34where (abs(tendencies_ice) < 1.e-10) tendencies_ice = 0.
[2779]35
36END SUBROUTINE compute_tendencies_slope
37
[3076]38END MODULE compute_tendencies_slope_mod
39
Note: See TracBrowser for help on using the repository browser.