Changeset 3616

Timestamp:

Feb 10, 2025, 2:28:11 PM (5 months ago)

Author:

jbclement

Message:

PEM:

• Improvement of the loops to update surface temperature in case of CO2 ice is disappearing.
• Improvement of messages printed on the terminal.

JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (12 diffs)
• read_data_PCM_mod.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -584 +584 @@
     - missing initialization for slope number (r3609);
     - taking into account that there may be no "startfi.nc" at the beginning of the PEM simulation.
+
+== 10/02/2025 == JBC
+- Improvement of the loops to update surface temperature in case of CO2 ice is disappearing.
+- Improvement of messages printed on the terminal.

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

@@ -231 +231 @@
 
 ! Some variables for the PEM run
-real, parameter   :: year_step = 1   ! Timestep for the pem
-real              :: i_myear_leg     ! Number of iteration
-real              :: n_myear_leg     ! Maximum number of iterations before stopping
-real              :: i_myear         ! Global number of Martian years of the chained simulations
-real              :: n_myear         ! Maximum number of Martian years of the chained simulations
-real              :: timestep        ! Timestep [s]
-character(100)    :: arg             ! To read command-line arguments program was invoked
-logical           :: timewall        ! Flag to use the time limit stopping criterion in case of a PEM job
-integer(kind = 8) :: cr              ! Number of clock ticks per second (count rate)
-integer(kind = 8) :: c1, c2          ! Counts of processor clock
-character(100)    :: chtimelimit     ! Time limit for the PEM job outputted by the SLURM command
-real              :: timelimit       ! Time limit for the PEM job in seconds
-real, parameter   :: antetime = 1200 ! Anticipation time to prevent reaching the time limit: 1200 s = 20 min by default
-integer           :: cstat, days, hours, minutes, seconds
-character(1)      :: sep
+real, parameter       :: year_step = 1   ! Timestep for the pem
+real                  :: i_myear_leg     ! Number of iteration
+real                  :: n_myear_leg     ! Maximum number of iterations before stopping
+real                  :: i_myear         ! Global number of Martian years of the chained simulations
+real                  :: n_myear         ! Maximum number of Martian years of the chained simulations
+real                  :: timestep        ! Timestep [s]
+character(100)        :: arg             ! To read command-line arguments program was invoked
+logical               :: timewall        ! Flag to use the time limit stopping criterion in case of a PEM job
+integer(kind = 8)     :: cr              ! Number of clock ticks per second (count rate)
+integer(kind = 8)     :: c1, c2          ! Counts of processor clock
+character(100)        :: chtimelimit     ! Time limit for the PEM job outputted by the SLURM command
+real                  :: timelimit       ! Time limit for the PEM job in seconds
+real, parameter       :: antetime = 1200 ! Anticipation time to prevent reaching the time limit: 1200 s = 20 min by default
+integer               :: cstat, days, hours, minutes, seconds
+character(1)          :: sep
+character(8)          :: date
+character(10)         :: time
+character(5)          :: zone
+integer, dimension(8) :: values
+character(128)        :: dir = ' '
+character(32)         :: logname = ' '
+character(32)         :: hostname = ' '
 
 #ifdef CPP_STD
@@ -358 +365 @@
 endif
 
+! Some user info
+call date_and_time(date,time,zone,values)
+call getcwd(dir)      ! Current directory
+call getlog(logname)  ! User name
+call hostnm(hostname) ! Machine/station name
+write(*,*)
+write(*,*) '######### PEM information #########'
+write(*,*) '+ User     : '//trim(logname)
+write(*,*) '+ Machine  : '//trim(hostname)
+write(*,*) '+ Directory: '//trim(dir)
+write(*,'(a,i2,a,i2,a,i4)') ' + Date     : ',values(3),'/',values(2),'/',values(1)
+write(*,'(a,i2,a,i2,a,i2,a)') ' + Time     : ',values(5),':',values(6),':',values(7)
+
 ! Parallel variables
 #ifndef CPP_STD
@@ -382 +402 @@
 !    I_a Read the "run.def"
 !------------------------
+write(*,*)
+write(*,*) '######### PEM initialization #########'
 write(*,*) '> Reading "run.def" (PEM)'
 #ifndef CPP_1D
@@ -706 +728 @@
 !    II_a Update pressure, ice and tracers
 !------------------------
+write(*,*)
+write(*,*) '######### PEM cycle #########'
 i_myear_leg = 0
 stopPEM = 0
@@ -723 +747 @@
 do while (i_myear_leg < n_myear_leg .and. i_myear < n_myear)
 ! II.a.1. Compute updated global pressure
+    write(*,'(a,f10.2)') ' #### Iteration of the PEM leg (Martian years): ', i_myear_leg + 1
     write(*,*) "> Updating the surface pressure"
     ps_avg_global_old = ps_avg_global_new
@@ -826 +851 @@
             if (is_co2ice_ini(ig,islope) .and. co2_ice(ig,islope) < 1.e-10 .and. .not. co2ice_disappeared(ig,islope)) then
                 co2ice_disappeared(ig,islope) = .true.
-                if (latitude_deg(ig) > 0.) then
-                    outer1: do ig_loop = ig,ngrid
-                        do islope_loop = islope,iflat,-1
+                if (latitude_deg(ig) > 0.) then ! North hemisphere
+                    outer1: do ig_loop = ig,ngrid ! Go towards equator
+                        do islope_loop = islope - 1,1,-1 ! Go over the slopes (backward numbering - should be equator-ward)
+                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
+                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
+                                exit outer1
+                            endif
+                        enddo
+                        do islope_loop = islope + 1,nslope ! Go over the slopes (forward numbering - should be pole-ward)
                             if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
                                 tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
@@ -835 +866 @@
                         enddo
                     enddo outer1
-                else
-                    outer2: do ig_loop = ig,1,-1
-                        do islope_loop = islope,iflat
+                else ! South hemisphere
+                    outer2: do ig_loop = ig,1,-1 ! Go towards equator
+                        do islope_loop = islope + 1,nslope ! Go over the slopes (forward numbering - should be equator-ward)
+                            if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
+                                tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
+                                exit outer2
+                            endif
+                        enddo
+                        do islope_loop = islope - 1,1,-1 ! Go over the slopes (backward numbering - should be pole-ward)
                             if (.not. is_co2ice_ini(ig_loop,islope_loop) .and. co2_ice(ig_loop,islope_loop) < 1.e-10) then
                                 tsurf_avg(ig,islope) = tsurf_avg(ig_loop,islope_loop)
@@ -845 +882 @@
                     enddo outer2
                 endif
-            else if (co2_ice(ig,islope) > 1.e-10 .and. d_co2ice(ig,islope) > 1.e-10) then ! Put tsurf as tcond CO2
+            else if (co2_ice(ig,islope) > 1.e-10 .and. d_co2ice(ig,islope) > 1.e-10) then ! Put tsurf as Tcond CO2
                 call computeTcondCO2(timelen,ngrid,nslope,vmr_co2_PEM_phys,ps_timeseries,ps_avg_global_ini,ps_avg_global_new,tsurf_avg)
             endif
@@ -976 +1013 @@
 !------------------------
     call recomp_tend_co2(ngrid,nslope,timelen,d_co2ice,d_co2ice_ini,co2_ice,emis,vmr_co2_PCM,vmr_co2_PEM_phys,ps_timeseries,ps_avg_global_old,ps_avg_global_new)
+    
     do ig = 1,ngrid
         do islope = 1,nslope
@@ -999 +1037 @@
         select case (stopPEM)
             case(1)
-                write(*,*) "STOPPING because surface of h2o ice sublimating is too low:", stopPEM, "See message above."
+                write(*,'(a,i0,a)') " #### STOPPING because surface of h2o ice sublimating is too low:", stopPEM, "See message above."
             case(2)
-                write(*,*) "STOPPING because tendencies on h2o ice = 0:", stopPEM, "See message above."
+                write(*,'(a,i0,a)') " #### STOPPING because tendencies on h2o ice = 0:", stopPEM, "See message above."
             case(3)
-                write(*,*) "STOPPING because surface of co2 ice sublimating is too low:", stopPEM, "See message above."
+                write(*,'(a,i0,a)') " #### STOPPING because surface of co2 ice sublimating is too low:", stopPEM, "See message above."
             case(4)
-                write(*,*) "STOPPING because surface global pressure changed too much:", stopPEM, "See message above."
+                write(*,'(a,i0,a)') " #### STOPPING because surface global pressure changed too much:", stopPEM, "See message above."
             case(5)
-                write(*,*) "STOPPING because maximum number of iterations is reached (possibly due to orbital parameters):", stopPEM
+                write(*,'(a,i0)') " #### STOPPING because maximum number of iterations is reached (possibly due to orbital parameters):", stopPEM
             case(6)
-                write(*,*) "STOPPING because maximum number of Martian years to be simulated is reached:", stopPEM
+                write(*,'(a,i0)') " #### STOPPING because maximum number of Martian years to be simulated is reached:", stopPEM
             case(7)
-                write(*,*) "STOPPING because the time limit for the PEM job will be reached soon:", stopPEM
+                write(*,'(a,i0)') " #### STOPPING because the time limit for the PEM job will be reached soon:", stopPEM
             case(8)
-                write(*,*) "STOPPING because the layering algorithm met an hasty end:", stopPEM
+                write(*,'(a,i0)') " #### STOPPING because the layering algorithm met an hasty end:", stopPEM
             case default
-                write(*,*) "STOPPING with unknown stopping criterion code:", stopPEM
+                write(*,'(a,i0)') " #### STOPPING with unknown stopping criterion code:", stopPEM
         end select
         exit
     else
-        write(*,*) "The PEM can continue!"
-        write(*,*) "Number of iterations of the PEM  : i_myear_leg =", i_myear_leg
-        write(*,*) "Number of simulated Martian years: i_myear     =", i_myear
+        write(*,'(a,f10.2,a)') ' #### The chained simulation has run for ',i_myear,' Martian years.'
+        write(*,*) '#### The PEM can continue!'
+        write(*,*) '####'
     endif
 enddo ! big time iteration loop of the pem
@@ -1046 +1084 @@
 !------------------------
 ! III_a.1 Ice update for start file
+write(*,*)
+write(*,*) '######### PEM finalization #########'
 write(*,*) '> Reconstructing perennial ice and frost for the PCM'
 watercap = 0.
@@ -1251 +1291 @@
 call info_PEM(i_myear_leg,stopPEM,i_myear,n_myear)
 
-write(*,*) '************************************************************************'
-write(*,*) "> The PEM has run for", i_myear_leg, "Martian years."
-write(*,*) "> The chained simulation has run for", i_myear, "Martian years =", i_myear*convert_years, "Earth years."
-write(*,*) "> The reached date is now", (year_bp_ini + i_myear)*convert_years, "Earth years."
-write(*,*) "> PEM: so far, so good!"
-write(*,*) '************************************************************************'
+write(*,*)
+write(*,*) '######### PEM finalization #########'
+write(*,'(a,f10.2,a)') " + The PEM leg has run for ", i_myear_leg, " Martian years."
+write(*,'(a,f10.2,a,f10.2,a)') " + The chained simulation has run for ", i_myear, " Martian years =", i_myear*convert_years, "Earth years."
+write(*,'(a,f10.2,a)') " + The reached date is now ", (year_bp_ini + i_myear)*convert_years, " Earth years."
+write(*,*) "+ PEM: so far, so good!"
+write(*,*) '####################################'
 
 if (layering_algo) then

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

@@ -304 +304 @@
 ! Close the NetCDF file
 call error_msg(nf90_close(fID),"close",filename2)
-write(*,*) '"> '//filename2//'" processed!'
+write(*,*) '> "'//filename2//'" processed!'
 
 END SUBROUTINE read_data_PCM