Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/dyn1d/rcm1d.F

Timestamp:

Mar 19, 2026, 2:35:46 PM (10 days ago)

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/dyn1d/rcm1d.F (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/dyn1d/rcm1d.F

@@ -19 +19 @@
      &         obliquit,nres,z0,coefvis,coefir,
      &         timeperi,e_elips,p_elips
-      use comcstfi_mod, only: pi, cpp, rad, g, r,
-     &                        mugaz, rcp, omeg
+      use comcstfi_mod, only: pi, cppd_ref, rad, g, rd_ref,
+     &                        mugaz_ref, rcp_ref, omeg
       use time_phylmdz_mod, only: daysec, dtphys, nday,
      &                            diagfi_output_rate
@@ -651 +651 @@
 !!! - physical constants: nevermind, things are done allright below
 !!! - physical frequency: nevermind, in inifis this is a simple print
-      cpp=-9999. ! dummy init for inifis, will be rewrite later on
-      r=-9999.   ! dummy init for inifis, will be rewrite later on
+      cppd_ref=-9999. ! dummy init for inifis, will be rewrite later on
+      rd_ref=-9999.   ! dummy init for inifis, will be rewrite later on
       CALL inifis(1,llm,nq,day0,daysec,nday,dtphys,
-     .            latitude,longitude,cell_area,rad,g,r,cpp)
+     .            latitude,longitude,cell_area,rad,g,rd_ref,cppd_ref)
 
       nsoil=nsoilmx
@@ -664 +664 @@
 !!! We check everything is OK.
       PRINT *,"CHECK"
-      PRINT *,"--> mugaz = ",mugaz
-      PRINT *,"--> cpp = ",cpp
-      r = 8.314463E+0 * 1000.E+0 / mugaz
-      rcp = r / cpp
+      PRINT *,"--> mugaz_ref = ",mugaz_ref
+      PRINT *,"--> cppd_ref = ",cppd_ref
+      rd_ref = 8.314463E+0 * 1000.E+0 / mugaz_ref
+      rcp_ref = rd_ref / cppd_ref
 
 c output spectrum?
@@ -1006 +1006 @@
 c  profil de temperature au premier appel
 c  --------------------------------------
-      pks=psurf**rcp
+      pks=psurf**rcp_ref
 
       if (restart) then
@@ -1182 +1182 @@
 !                 s(ilayer)=(play(ilayer)/psurf)**rcp
 !              else
-          s(ilayer)=(aps(ilayer)/psurf+bps(ilayer))**rcp
+          s(ilayer)=(aps(ilayer)/psurf+bps(ilayer))**rcp_ref
 !              endif
               !s(ilayer)=(aps(ilayer)/psurf+bps(ilayer))**rcp
-          h(ilayer)=cpp*temp(ilayer)/(pks*s(ilayer))
+          h(ilayer)=cppd_ref*temp(ilayer)/(pks*s(ilayer))
        ENDDO