- Timestamp:
- Jan 31, 2023, 9:52:30 PM (22 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/evolution/criterion_ice_stop_mod_slope.F90
r2835 r2885 2 2 ! $Id $ 3 3 ! 4 SUBROUTINE criterion_ice_stop_slope(cell_area,ini_surf,qsurf,STOPPING, ngrid,initial_h2o_ice,global_ave_press_GCM,global_ave_press_new,nslope)4 SUBROUTINE criterion_ice_stop_slope(cell_area,ini_surf,qsurf,STOPPING,STOPPING_ps,ngrid,initial_h2o_ice,global_ave_press_GCM,global_ave_press_new,nslope) 5 5 6 6 USE temps_mod_evol, ONLY: alpha_criterion … … 28 28 29 29 ! OUTPUT 30 LOGICAL, intent(out) :: STOPPING ! Logical : is the criterion reached?30 LOGICAL, intent(out) :: STOPPING,STOPPING_ps ! Logical : is the criterion reached? 31 31 32 32 ! local: … … 39 39 ! initialisation to false 40 40 STOPPING=.FALSE. 41 STOPPING_ps=.FALSE. 41 42 42 43 ! computation of the actual surface … … 67 68 if(global_ave_press_new.LT.global_ave_press_GCM*(0.9) .OR. & 68 69 global_ave_press_new.GT.global_ave_press_GCM*(1.1)) then 69 STOPPING =.TRUE.70 STOPPING_ps=.TRUE. 70 71 print *, "Reason of stopping : The global pressure reach the threshold:" 71 72 print *, "Current global pressure=", global_ave_press_new
Note: See TracChangeset
for help on using the changeset viewer.