Changeset 2794 for trunk/LMDZ.COMMON/libf/evolution/criterion_ice_stop_mod_water_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/criterion_ice_stop_mod_water_slope.F90 (modified) (3 diffs)

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

@@ -5 +5 @@
 
   USE temps_mod_evol, ONLY: alpha_criterion
-  use comslope_mod, ONLY: subslope_dist,nslope
+          use comslope_mod, ONLY: subslope_dist,nslope
 
       IMPLICIT NONE
@@ -19 +19 @@
 
 !   INPUT
-  INTEGER, intent(in) :: ngrid                         ! # of grid physical grid points 
-  REAL,    intent(in) :: cell_area(ngrid)              ! physical point field : Area of the cells
+  INTEGER, intent(in) :: ngrid                  ! # of grid physical grid points 
+  REAL,    intent(in) :: cell_area(ngrid)       ! physical point field : Area of the cells
   REAL,    intent(in) ::  qsurf(ngrid,nslope)          ! physical point field : Actual density of water ice
-  REAL,    intent(in) :: ini_surf                      ! Surface of sublimating ice in the GCM output
-  REAL,    intent(in) :: initial_h2o_ice(ngrid,nslope) ! Boolean to check if it is a sublimating area
+  REAL,    intent(in) :: ini_surf
+  REAL,    intent(in) :: initial_h2o_ice(ngrid,nslope)
 
 
 !   OUTPUT
-  LOGICAL, intent(out) :: STOPPING                     ! Logical : is the criterion reached?
+  LOGICAL, intent(out) :: STOPPING              ! Logical : is the criterion reached?
 
 !   local:
 !   -----
-  INTEGER :: i,islope                                  ! Loop
-  REAL :: present_surf                                 ! Initial/Actual surface of water ice
+  INTEGER :: i,islope                    ! Loop
+  REAL :: present_surf  ! Initial/Actual surface of water ice
 
 !=======================================================================
@@ -44 +44 @@
     do islope=1, nslope
       if (initial_h2o_ice(i,islope).GT.0.5 .and. qsurf(i,islope).GT.0.) then
+!         print *, "i", i
+!         print *, "initial_h2o_ice(i,islope)", initial_h2o_ice(i,islope)
+!         print *, "qsurf(i,islope)", qsurf(i,islope)
+!         print *, "cell_area(i)", cell_area(i)
+!         print *, "present_surf",present_surf
          present_surf=present_surf+cell_area(i)*subslope_dist(i,islope)
       endif
     enddo
   enddo
+
+!  print *, "initial_h2o_ice", initial_h2o_ice
+!  print *, "qsurf", qsurf
+
+!  print *, "present_surf", present_surf
+!  print *, "ini_surf", ini_surf
+!  print *, "ini_surf*0.8", ini_surf*(1-alpha_criterion)
   
 !   check of the criterion