Changeset 3442 for trunk/LMDZ.COMMON/libf

Timestamp:

Sep 27, 2024, 5:29:39 PM (3 months ago)

Author:

jbclement

Message:

PEM:
Correction: the potential temperature teta is now updated due to the change of surface pressure at the end of the PEM to write the "restart.nc" file.
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

@@ -430 +430 @@
 == 20/09/2024 == JBC
 Small corrections for the launching script (run numbers must be integer) and the PEM stopping criteria algorithm related to r3430.
+
+== 27/09/2024 == JBC
+Correction: the potential temperature teta is now updated due to the change of surface pressure at the end of the PEM to write the "restart.nc" file.

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

@@ -83 +83 @@
     use mod_phys_lmdz_para, only: is_parallel, is_sequential, is_mpi_root, is_omp_root, is_master
     use planete_h,          only: aphelie, periheli, year_day, peri_day, obliquit, iniorbit
-    use comcstfi_h,         only: pi, rad, g, cpp, mugaz, r
+    use comcstfi_h,         only: pi, rad, g, mugaz, r
     use surfini_mod,        only: surfini
+    use comconst_mod,       only: kappa, cpp
 #else
     use tracer_h,           only: noms, igcm_h2o_ice, igcm_co2 ! Tracer names
@@ -92 +93 @@
     use aerosol_mod,        only: iniaerosol
     use planete_mod,        only: apoastr, periastr, year_day, peri_day, obliquit
-    use comcstfi_mod,       only: pi, rad, g, cpp, mugaz, r
+    use comcstfi_mod,       only: pi, rad, g, mugaz, r
 #endif
 
@@ -138 +139 @@
 real, dimension(:,:,:), allocatable :: q             ! champs advectes
 real, dimension(ip1jmp1)            :: ps            ! pression au sol
+real, dimension(ip1jmp1)            :: ps0           ! pression au sol initiale
 real, dimension(:),     allocatable :: ps_start_PCM  ! (ngrid) surface pressure
 real, dimension(:,:),   allocatable :: ps_timeseries ! (ngrid x timelen) instantaneous surface pressure
@@ -409 +411 @@
     ps_start_PCM(1) = ps(1)
 #endif
+
+! Save initial surface pressure
+ps0 = ps
 
 ! In the PCM, these values are given to the physic by the dynamic.
@@ -1140 +1145 @@
 allocate(p(ip1jmp1,nlayer + 1))
 #ifndef CPP_1D
-    call pression (ip1jmp1,ap,bp,ps,p)
+    ! Correction on teta due to surface pressure changes
+    do l = 1,nlayer
+        do i = 1,ip1jmp1
+            teta(i,l)= teta(i,l)*(1 + (ps(i) - ps0(i))*year_iter/ps(i))**kappa
+        enddo
+    enddo
+    ! Correction on atmospheric pressure
+    call pression(ip1jmp1,ap,bp,ps,p)
+    ! Correction on the mass of atmosphere
     call massdair(p,masse)
     call dynredem0("restart.nc",day_ini,phis)