Changeset 3738 for trunk/LMDZ.COMMON/libf

Timestamp:

Apr 25, 2025, 6:30:25 PM (8 weeks ago)

Author:

jbclement

Message:

PEM:
Resulting change due to revision r3737 in the 1D Mars PCM.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

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

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

@@ -625 +625 @@
 Updating the average pressure at each iteration.
 
-== 1/4/2025 == EV
+== 01/04/2025 == EV
 updating paramaters to work with the layering
 
-== 2/4/2025 == EV
+== 02/04/2025 == EV
 updating the h2o and co2 ice values that will be used by the PCM according to the top layer in the stratum
 updating h2o ice recomp tendency to reduce the ice table loss rate due increase of depth and a decrease in the humidity around the ice.
 
-== 3/4/2025 == EV
+== 03/04/2025 == EV
 adding the growth of the lag layer
 changing the stopping critiria to take into account the top startum layer
 removing the case where if we have 0 ice the PEM stops.
 
-
-
+== 25/04/2025 == JBC
+Resulting change due to revision r3737 in the 1D Mars PCM.

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

@@ -111 +111 @@
     use comvert_mod,              only: ap, bp
     use writerestart1D_mod,       only: writerestart1D
-   ! use layering_mod,             only: layering, array2stratif, nb_str_max, layering_algo
 #endif
 
@@ -287 +286 @@
     real, dimension(:),     allocatable :: zqsat
     real, dimension(:,:,:), allocatable :: dq, dqdyn
-    real, dimension(nlayer)             :: play, w
+    real, dimension(nlayer)             :: play, w, q_def_h2o_vap
     real, dimension(nlayer + 1)         :: plev
 #else
@@ -323 +322 @@
         endif
     enddo
-    
+
     if (num_str) then ! This is a numeric sting so we considerer this is the job id
            ! Execute the system command
@@ -436 +435 @@
     call init_testphys1d('start1D.txt','startfi.nc',therestart1D,therestartfi,ngrid,nlayer,610.,nq,q,         &
                          time_0,pdyn(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w, &
-                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
+                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,q_def_h2o_vap)
     nsplit_phys = 1
     day_step = steps_per_sol
@@ -1302 +1301 @@
 
 write(*,*)
-write(*,*) '********* PEM finalization *********'
+write(*,*) '****** PEM final information *******'
 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."