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

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

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc.F (modified) (1 diff)

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