source: trunk/LMDZ.COMMON/libf/evolution/evol_co2_ice_s_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.5 KB
Line 
1MODULE evol_co2_ice_s_mod
2
3implicit none
4
5!=======================================================================
6contains
7!=======================================================================
8
9SUBROUTINE evol_co2_ice_s(qsurf,tendencies_co2_ice_phys,iim_input,jjm_input,ngrid,cell_area,nslope)
10
11use time_evol_mod, only: dt_pem
12
13implicit none
14
15!=======================================================================
16!
17!  Routine that compute the evolution of the CO2 ice
18!
19!=======================================================================
20!   arguments:
21!   ----------
22!   INPUT
23integer,                intent(in) :: iim_input, jjm_input, ngrid, nslope ! # of grid points along longitude/latitude/ total
24REAL, dimension(ngrid), intent(in) :: cell_area
25
26!   OUTPUT
27real, dimension(ngrid,nslope), intent(inout) :: qsurf                   ! physical point field: Previous and actual density of water ice
28real, dimension(ngrid,nslope), intent(inout) :: tendencies_co2_ice_phys ! physical point field: Evolution of perennial ice over one year
29
30!   local:
31!   ----
32integer :: i, j, ig0, islope ! loop variable
33!=======================================================================
34! Evolution of the CO2 ice for each physical point
35do i=1,ngrid
36    do islope=1,nslope
37        qsurf(i,islope) = qsurf(i,islope) + tendencies_co2_ice_phys(i,islope)*dt_pem
38        if (qsurf(i,islope) < 0) then
39            qsurf(i,islope) = 0.
40            tendencies_co2_ice_phys(i,islope) = 0.
41        endif
42    enddo
43enddo
44
45END SUBROUTINE evol_co2_ice_s
46
47END MODULE evol_co2_ice_s_mod
Note: See TracBrowser for help on using the repository browser.