Changeset 4185 for trunk/LMDZ.VENUS

Timestamp:

Apr 13, 2026, 4:04:43 PM (4 weeks ago)

Author:

flefevre

Message:

Bug fix in mean Cp calculation used in the physics

File:

• trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F (modified) (4 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/concentrations2.F

@@ -210 +210 @@
       end do
 
-
 !     update mass mixing ratio tracers
 
@@ -216 +215 @@
          do ig = 1,klon
             do i = 1,nqmx
-!               iq = niq(i) 
                zq(ig,l,i) = max(1.e-30, tr_seri(ig,l,i))
             end do
@@ -232 +230 @@
          do ig = 1,klon
             do i = 1,nqmx-nmicro
-c               iq = niq(i) 
                mmean(ig,l) = mmean(ig,l) + zq(ig,l,i)/M_tr(i)
             end do
@@ -248 +245 @@
 
       do l = 1,klev
-          do ig = 1,klon
-
-            ntot = pplay(ig,l)/(RKBOL*zt(ig,l))*1.e-6  ! in #/cm3 
+         do ig = 1,klon
+
+            ntot = pplay(ig,l)/(RKBOL*zt(ig,l))*1.e-6       ! in #/cm3 
             rho(ig,l) = (ntot * mmean(ig,l))/RNAVO*1.e3     ! in kg/m3
-
-c            write(*,*),'Air density: ',ig, l, rho(0,l)            
-
-!!  WARNING -> Cp here below doesn't depend on T (cpdet)
 
             do i = 1,nbq
                iq = niq(i)
+               if (iq == i_co2) then
+                 cpnew(ig,l) = cpnew(ig,l) + zq(ig,l,iq)*cpdet(zt(ig,l))
+               else
+                 cpnew(ig,l) = cpnew(ig,l) + zq(ig,l,iq)*cpi(i)
+               end if
                ni(i) = ntot*zq(ig,l,iq)*mmean(ig,l)/M_tr(iq)
-!! On a une super formule pour calculer cp_co2 sur Venus
-               if (iq == i_co2) then
-                  cpnew(ig,l) = cpnew(ig,l) + ni(i)*cpdet(zt(ig,l))
-               else
-                  cpnew(ig,l) = cpnew(ig,l) + ni(i)*cpi(i)
-               end if
                akknew(ig,l) = akknew(ig,l) + ni(i)*aki(i)
             end do 
  
-            cpnew(ig,l) = cpnew(ig,l)/ntot
             akknew(ig,l)= akknew(ig,l)/ntot
 
-
-          end do 
-       end do
+         end do 
+      end do
 
       return