Changeset 4164 for trunk/LMDZ.GENERIC/libf/phygeneric/pcm_thermodynamics_mod.F90

Timestamp:

Apr 2, 2026, 9:15:00 AM (12 days ago)

Author:

emillour

Message:

Genric PCM:
Some changes following the "thermodynamics update":

• bug fix on the size of pplev
• renaming the module to something more explicit (thermo_mod.F90 => pcm_thermodynamics_mod.F90) and the Exner&potential temperature update routine "thermodynamics" => "thermodynamics_update"
• some optimizations in "thermodynamics_update" : in the 'thermo_uni_ideal' case seting r(),cpp() and rcp() can be done at first call only.

EM

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/pcm_thermodynamics_mod.F90 (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/thermo_mod.F90) (6 diffs)

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

@@ -1 @@
-      module thermo_mod
+      module pcm_thermodynamics_mod
 
       use comcstfi_mod, only: cppd_ref,cppv_ref,rd_ref, rcp_ref, mugaz_ref
@@ -12 +12 @@
       contains
       
-      subroutine thermodynamics(thermo_phy, pplay, pplev, t, ngrid, nlayer, nq, q, iq, zh, zpopsk)
-        CHARACTER(64), INTENT(IN) :: thermo_phy               ! flag
-        REAL, INTENT(IN)    :: pplay(ngrid,nlayer)      ! pressure (Pa)
-        REAL, INTENT(IN)    :: pplev(ngrid,nlayer)      ! pressure (Pa)
+      subroutine thermodynamics_update(thermo_phy, pplay, pplev, t, ngrid, nlayer, nq, q, iq, zh, zpopsk)
+        CHARACTER(LEN=*), INTENT(IN) :: thermo_phy               ! flag
+        REAL, INTENT(IN)    :: pplay(ngrid,nlayer)    ! pressure (Pa) at mid-layer
+        REAL, INTENT(IN)    :: pplev(ngrid,nlayer+1)  ! pressure (Pa) at layer interfaces
         REAL, INTENT(IN)    :: t(ngrid,nlayer)      ! temperature (K)
         INTEGER, INTENT(IN) :: ngrid, nlayer,nq     ! Number of cells, vertical layers and tracers
@@ -27 +27 @@
         logical, save :: firstcall=.true.
 !$OMP THREADPRIVATE(firstcall)
+        character(len=80),parameter :: myname = "thermodynamics_update" 
         
         if (firstcall) then
@@ -32 +33 @@
           ALLOCATE(cpp(ngrid,nlayer))
           ALLOCATE(rcp(ngrid,nlayer))
+          
+          SELECT CASE (TRIM(thermo_phy))
+            CASE('thermo_uni_ideal')
+              ! Ideal gas, homogeneous
+              r(:,:) = rd_ref
+              cpp(:,:) = cppd_ref
+              rcp(:,:) = rcp_ref
+            CASE DEFAULT
+            write(*,*) 'Bad selector for thermodynamics mod: <', TRIM(thermo_phy), '>'
+            call abort_physic(trim(myname),'Bad selector for thermodynamics mod!',1)
+          END SELECT
+          
           firstcall=.false.
-        endif
+        endif ! of if (firstcall)
         
         SELECT CASE (TRIM(thermo_phy))
@@ -39 +52 @@
          CASE('thermo_uni_ideal')
             ! Ideal gas
-            r(:,:) = rd_ref
-            cpp(:,:) = cppd_ref
-            rcp(:,:) = rcp_ref
             do l=1,nlayer
               do ig=1,ngrid
@@ -51 +61 @@
             write(*,*) 'Bad selector for thermodynamics mod: <', TRIM(thermo_phy), '>'
             write(*,*) 'Option is <thermo_uni_ideal>'
-            call abort
+            call abort_physic(trim(myname),'Bad selector for thermodynamics mod!',1)
         END SELECT
         
-        end subroutine thermodynamics
+        end subroutine thermodynamics_update
         
-        end module thermo_mod
+        end module pcm_thermodynamics_mod