Changeset 1035

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

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/README

@@ -1892 +1892 @@
 - Cleaned and commented version of thermal plume model with automatic arrays
 - Checked: exact same results than before modifications
+
+== 09/11/2013 == EM
+- 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.

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

trunk/LMDZ.MARS/makegcm_ifort

@@ -384 +384 @@
         else if $LINUX then
            ## for ifort
-             set optim=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
-             set optim90=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
-             set optimtru90=" -g -no-ftz -traceback -ftrapuv -fp-stack-check -check"
+             set optim=" -g -O0 -fpe-all=0 -traceback -ftrapuv -fp-stack-check -check all -debug"
+             set optim90=" -g -O0 -fpe-all=0 -traceback -ftrapuv -fp-stack-check -check all -debug"
+             set optimtru90=" -g -O0 -fpe-all=0 -traceback -ftrapuv -fp-stack-check -check all -debug"
         else 
            echo "pas d option debug predefinie pour cette machine"