Changeset 1035 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Sep 11, 2013, 12:23:22 PM (11 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

File:

• trunk/LMDZ.MARS/libf/aeronomars/concentrations.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