Changeset 3050 for trunk/LMDZ.COMMON/libf/evolution/criterion_pem_stop_mod.F90

Timestamp:

Sep 26, 2023, 3:57:38 PM (18 months ago)

Author:

jbclement

Message:

Mars PEM:
Minor changes concerning the form of the code in the PEM.
JBC

File:

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

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

@@ -63 +63 @@
      present_surf.GT.ini_surf*(1+water_ice_criterion)) then
     STOPPING=.TRUE. 
-    print *, "Reason of stopping : The surface of water ice sublimating reach the threshold:"
-    print *, "Current surface of water ice sublimating=", present_surf
-    print *, "Initial surface of water ice sublimating=", ini_surf
-    print *, "Percentage of change accepted=", water_ice_criterion*100
-    print *, "present_surf<ini_surf*(1-water_ice_criterion)", (present_surf.LT.ini_surf*(1-water_ice_criterion))
+    write(*,*) "Reason of stopping : The surface of water ice sublimating reach the threshold:"
+    write(*,*) "Current surface of water ice sublimating=", present_surf
+    write(*,*) "Initial surface of water ice sublimating=", ini_surf
+    write(*,*) "Percentage of change accepted=", water_ice_criterion*100
+    write(*,*) "present_surf<ini_surf*(1-water_ice_criterion)", (present_surf.LT.ini_surf*(1-water_ice_criterion))
   endif
 
@@ -129 +129 @@
      present_surf.GT.ini_surf*(1+co2_ice_criterion)) then
     STOPPING_ice=.TRUE. 
-    print *, "Reason of stopping : The surface of co2 ice sublimating reach the threshold:"
-    print *, "Current surface of co2 ice sublimating=", present_surf
-    print *, "Initial surface of co2 ice sublimating=", ini_surf
-    print *, "Percentage of change accepted=", co2_ice_criterion*100
-    print *, "present_surf<ini_surf*(1-co2_ice_criterion)", (present_surf.LT.ini_surf*(1-co2_ice_criterion))
+    write(*,*) "Reason of stopping : The surface of co2 ice sublimating reach the threshold:"
+    write(*,*) "Current surface of co2 ice sublimating=", present_surf
+    write(*,*) "Initial surface of co2 ice sublimating=", ini_surf
+    write(*,*) "Percentage of change accepted=", co2_ice_criterion*100
+    write(*,*) "present_surf<ini_surf*(1-co2_ice_criterion)", (present_surf.LT.ini_surf*(1-co2_ice_criterion))
   endif
 
@@ -143 +143 @@
      global_ave_press_new.GT.global_ave_press_GCM*(1+ps_criterion)) then
     STOPPING_ps=.TRUE. 
-    print *, "Reason of stopping : The global pressure reach the threshold:"
-    print *, "Current global pressure=", global_ave_press_new
-    print *, "GCM global pressure=", global_ave_press_GCM
-    print *, "Percentage of change accepted=", ps_criterion*100
-    print *, "global_ave_press_new<global_ave_press_GCM*(ps_criterion)", (global_ave_press_new.LT.global_ave_press_GCM*(1-ps_criterion))
+    write(*,*) "Reason of stopping : The global pressure reach the threshold:"
+    write(*,*) "Current global pressure=", global_ave_press_new
+    write(*,*) "GCM global pressure=", global_ave_press_GCM
+    write(*,*) "Percentage of change accepted=", ps_criterion*100
+    write(*,*) "global_ave_press_new<global_ave_press_GCM*(ps_criterion)", (global_ave_press_new.LT.global_ave_press_GCM*(1-ps_criterion))
   endif