Changeset 3628 for trunk/LMDZ.COMMON

Timestamp:

Feb 18, 2025, 3:31:38 PM (5 days ago)

Author:

jbclement

Message:

PEM:
Bug correction about pressure rewriting in the "start.nc".
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

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

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

@@ -594 +594 @@
 == 12/02/2025 == JBC
 Few small adjustments in the code.
+
+== 18/02/2025 == JBC
+Bug correction about pressure rewriting in the "start.nc".

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

@@ -291 +291 @@
     real, dimension(:),   allocatable :: ap ! Coefficient ap read in start_name and written in restart
     real, dimension(:),   allocatable :: bp ! Coefficient bp read in start_name and written in restart
-    real, dimension(:,:), allocatable :: p  ! Grid points x Atmosphere: pressure to recompute and write in restart (ngrid,llmp1)
+    real, dimension(:,:), allocatable :: p  ! Grid points x Atmosphere: pressure to recompute and write in restart (ip1jmp1,llmp1)
 #endif
 
@@ -469 +469 @@
     call iniphysiq(iim,jjm,llm,(jjm - 1)*iim + 2,comm_lmdz,daysec,day_ini,dtphys/nsplit_phys,rlatu,rlatv,rlonu,rlonv,aire,cu,cv,rad,g,r,cpp,iflag_phys)
 #else
-    ps_start0 = pdyn
+    ps_start0(1) = pdyn(1)
 #endif
 deallocate(pdyn)
@@ -1039 +1039 @@
         select case (stopPEM)
             case(1)
-                write(*,'(a,i0,a)') " #### 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(*,'(a,i0,a)') " #### 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(*,'(a,i0,a)') " #### 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(*,'(a,i0,a)') " #### 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(*,'(a,i0)') " #### 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(*,'(a,i0)') " #### 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(*,'(a,i0)') " #### 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(*,'(a,i0)') " #### 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(*,'(a,i0)') " #### STOPPING with unknown stopping criterion code:", stopPEM
+                write(*,'(a,i0)') " #### STOPPING with unknown stopping criterion code: ", stopPEM
         end select
         exit
@@ -1251 +1251 @@
     call massdair(p,masse)
     call dynredem0("restart.nc",day_ini,phis)
-    call dynredem1("restart.nc",time_0,vcov,ucov,teta,q,masse,ps_start)
+    call dynredem1("restart.nc",time_0,vcov,ucov,teta,q,masse,pdyn)
     deallocate(ap,bp,p,pdyn)
 #else