Changeset 1035 for trunk/LMDZ.MARS/libf

Timestamp:

Sep 11, 2013, 12:23:22 PM (12 years ago)

Author:

emillour

Message:

Mars GCM:

• Bug fix in vdifc.F: in some cases, some elements of pdqsdif() were not given any value. In all cases, it is safer to start with clean initialization of output tendencies to zero.
• Bug fix in concentration.F: error in cpi() and aki() indexes led to wrong computation of atmospheric conductivity and specific heat.

EM

Location:

trunk/LMDZ.MARS/libf

Files:

• aeronomars/concentrations.F (modified) (1 diff)
• phymars/vdifc.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/aeronomars/concentrations.F

@@ -271 +271 @@
                iq = niq(i) 
                ni(iq) = ntot*zq(ig,l,iq)*mmean(ig,l)/mmol(iq)
-               cpnew(ig,l) = cpnew(ig,l) + ni(iq)*cpi(iq)
-               akknew(ig,l) = akknew(ig,l) + ni(iq)*aki(iq)
+               cpnew(ig,l) = cpnew(ig,l) + ni(iq)*cpi(i)
+               akknew(ig,l) = akknew(ig,l) + ni(iq)*aki(i)
             end do 
             cpnew(ig,l) = cpnew(ig,l)/ntot

trunk/LMDZ.MARS/libf/phymars/vdifc.F

@@ -206 +206 @@
 
       nlev=nlay+1
+
+      ! initialize output tendencies to zero:
+      pdudif(1:ngrid,1:nlay)=0
+      pdvdif(1:ngrid,1:nlay)=0
+      pdhdif(1:ngrid,1:nlay)=0
+      pdtsrf(1:ngrid)=0
+      pdqdif(1:ngrid,1:nlay,1:nq)=0
+      pdqsdif(1:ngrid,1:nq)=0
 
 c    ** calcul de rho*dz et dt*rho/dz=dt*rho**2 g/dp