Changeset 3972 for trunk/LMDZ.PLUTO/libf

Timestamp:

Nov 21, 2025, 2:57:40 PM (3 weeks ago)

Author:

tbertrand

Message:

PLUTO PCM:
Fixing bug in ch4clouds
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/ch4cloud.F90 (modified) (4 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/ch4cloud.F90

@@ -5 +5 @@
 
       use comgeomfi_h
-      use comcstfi_mod, only: pi, g, cpp
+      use comcstfi_mod, only: pi, g, cpp, r
       use tracer_h, only: igcm_ch4_gas, igcm_ch4_ice, & 
              rho_ch4_ice,lw_ch4,micro_indx
@@ -75 +75 @@
       REAL zqsat(ngrid,nlay)    ! saturation
       REAL zt(ngrid,nlay)       ! local value of temperature
+      !REAL rho(ngrid,nlay)       ! local value of atm density
       REAL nmix(ngrid,nlay)     ! Profil of Nmix
 
@@ -140 +141 @@
       zq(:,:,i_ice)=pq(:,:,i_ice)+pdq(:,:,i_ice)*ptimestep
 
+      !rho(:,:)=pplay(:,:)/(r*zt(:,:))
+      ! microphysics tracers are in particle/kg
       if (callmufi) then
         zq(:,:,micro_indx(1))=pq(:,:,micro_indx(1))+ &
@@ -162 +165 @@
       endif
 
-      pdqscloud(1:ngrid,1:nq)=0
-      pdqcloud(1:ngrid,1:nlay,1:nq)=0
-      pdtcloud(1:ngrid,1:nlay)=0
+      pdqscloud(1:ngrid,1:nq)=0.
+      pdqcloud(1:ngrid,1:nlay,1:nq)=0.
+      pdtcloud(1:ngrid,1:nlay)=0.
 
       ! Nmix profile
       if (callmufi) then
-         nmix(:,:)=zq(ig,l,micro_indx(1))+zq(ig,l,micro_indx(3))
+        nmix(:,:)=zq(:,:,micro_indx(1))+zq(:,:,micro_indx(3))
       else
-         nmix(:,:)=Nmix_ch4
+        nmix(:,:)=Nmix_ch4
       endif
-
 !    ----------------------------------------------
 !