Changeset 2794 for trunk/LMDZ.COMMON/libf/evolution/evol_co2_ice_s_mod_slope.F90

Timestamp:

Sep 9, 2022, 4:02:51 PM (2 years ago)

Author:

llange

Message:

MARS PEM:

• Add a PEMETAT0 that read "startfi_pem.nc"
• Add the soil in the model: soil temperature, thermal properties, ice table
• Add a routine that compute CO2 + H2O adsorption
• Minor corrections in PEM.F90

LL

File:

• trunk/LMDZ.COMMON/libf/evolution/evol_co2_ice_s_mod_slope.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/evol_co2_ice_s_mod_slope.F90

@@ -20 +20 @@
 !   INPUT
 
-  INTEGER, intent(in) :: iim_input,jjm_input, ngrid,nslope      ! # of grid points along longitude/latitude/ total
-  REAL, intent(in) ::  cell_area(ngrid)                         ! Area of each cell
+  INTEGER, intent(in) :: iim_input,jjm_input, ngrid,nslope   ! # of grid points along longitude/latitude/ total
+!  REAL, intent(in) ::  tendencies_h2o_ice_phys(ngrid) ! physical point field : Evolution of perenial ice over one year
+  REAL, intent(in) ::  cell_area(ngrid)
 
 !   OUTPUT
-  REAL, INTENT(INOUT) ::  qsurf(ngrid,nslope)                   ! physical point field : Previous and actual density of water ice
-  LOGICAL :: STOPPING                                           ! Not used
+  REAL, INTENT(INOUT) ::  qsurf(ngrid,nslope)                ! physical point field : Previous and actual density of water ice
+  LOGICAL :: STOPPING
   REAL, intent(inout) ::  tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
 
@@ -33 +34 @@
 
   INTEGER :: i,j,ig0,islope                                  ! loop variable
+!  REAL :: not_budget, budget
+  REAL :: pos_tend, neg_tend, real_coefficient,negative_part
+  REAL ::  new_tendencies(ngrid)
 
   STOPPING=.false.
 
 
-! Evolution of the co2 ice for each physical point
+! Evolution of the water ice for each physical point
   do i=1,ngrid
     do islope=1,nslope