Changeset 2857 for trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90

Timestamp:

Dec 29, 2022, 10:58:39 AM (2 years ago)

Author:

llange

Message:

PEM
Fixing bug when updating tracer mass mixing ratio (ccn_number and dust_number were previously set to 1)
Tiny edits in evol_h2o_ice and criterion_ice_stop_mod_water
LL

File:

• trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod_slope.F90 (modified) (3 diffs)

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

@@ -26 +26 @@
 !   OUTPUT
   REAL, INTENT(INOUT) ::  qsurf(ngrid,nslope)                ! physical point field : Previous and actual density of water ice
-  LOGICAL :: STOPPING
+  LOGICAL, INTENT(INOUT) :: STOPPING
   REAL, intent(inout) ::  tendencies_h2o_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year
 
@@ -80 +80 @@
     print *, "Tendencies on ice increasing=", pos_tend
     print *, "This can be due to the absence of water ice in the PCM run!!"
-      call criterion_ice_stop_water_slope(cell_area,1,qsurf(:,:)*0.,STOPPING,ngrid,cell_area)
+      call criterion_ice_stop_water_slope(cell_area,1.,qsurf(:,:)*0.,STOPPING,ngrid,cell_area)
       do i=1,ngrid
          do islope=1,nslope
@@ -101 +101 @@
     enddo
   enddo
+  
 
   real_coefficient=negative_part/pos_tend
+  if(pos_tend.eq.0) real_coefficient = 0.
 
   do i=1,ngrid