Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/calc_cpp_mugaz.F90

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/calc_cpp_mugaz.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/calc_cpp_mugaz.F90

@@ -19 +19 @@
 
       use gases_h
-      use comcstfi_mod, only: cpp, mugaz
+      use comcstfi_mod, only: cppd_ref, mugaz_ref
       use callkeys_mod, only: cpp_mugaz_mode
       use mod_phys_lmdz_para, only : is_master
@@ -155 +155 @@
       endif
 
-      if(((abs(1.-cpp/cpp_c).gt.1.e-6) .or.  &
-           (abs(1.-mugaz/mugaz_c).gt.1.e-6)).and. is_master ) then
+      if(((abs(1.-cppd_ref/cpp_c).gt.1.e-6) .or.  &
+           (abs(1.-mugaz_ref/mugaz_c).gt.1.e-6)).and. is_master ) then
          ! Ehouarn: tolerate a small mismatch between computed/stored values
-         print*,'--> Values do not match with the predefined one ! (',cpp,',',mugaz,')'
+         print*,'--> Values do not match with the predefined one ! (',cppd_ref,',',mugaz_ref,')'
          print*,'    Because cp varies with temperature and that some gases may not appear in gases.def,'
          print*,'    a small discrepancy might be completely normal.'
@@ -168 +168 @@
       if (cpp_mugaz_mode==2) then
           if (is_master) print*,'*** cpp_mugaz_mode==2, so setting cpp & mugaz to computations in calc_cpp_mugaz.'
-          mugaz = mugaz_c
-          cpp = cpp_c
+          mugaz_ref = mugaz_c
+          cppd_ref = cpp_c
       else
           if (is_master) print*,'*** Leaving cpp & mugaz equal to predefined values'