Changeset 3589 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Jan 20, 2025, 7:05:11 PM (3 days ago)

Author:

jbclement

Message:

PEM:
The initialization of constants between dynamics and physics was messed up which led to uninitialized variables.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (4 diffs)

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

@@ -84 +84 @@
     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, mugaz, r
+    use comcstfi_h,         only: mugaz
     use surfini_mod,        only: surfini
-    use comconst_mod,       only: kappa, cpp
+    use comconst_mod,       only: pi, rad, g, r, cpp, kappa
 #else
     use tracer_h,           only: noms, igcm_h2o_ice, igcm_co2 ! Tracer names
@@ -405 +405 @@
     call gr_dyn_fi(1,iip1,jjp1,ngridmx,ps_start_dyn,ps_start)
 
-    call iniconst !new
-    call inigeom
+    call iniconst ! Initialization of dynamical constants (comconst_mod)
+    call inigeom ! Initialization of geometry
 
     allocate(ap(nlayer + 1))
@@ -417 +417 @@
     status = nf90_close(ncid)
 
+    ! Initialization of physics constants (comcstfi_h)
     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
@@ -512 +513 @@
     if (noms(nnq) == "co2") igcm_co2 = nnq
 enddo
-r = 8.314511*1000./mugaz
 
 !------------------------