Changeset 1263

Timestamp:

May 15, 2014, 5:48:33 PM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS possible floating exception in newcondens makes debug mode to stop. no impact on results in normal mode because a correction was proposed right afterwards. fixed the problem by extending the if loop to be really operational as intended.

File:

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

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

@@ -290 +290 @@
             zt(ig,l)=pt(ig,l)+ pdt(ig,l)*ptimestep
 !            ztcond(ig,l)=1./(bcond-acond*log(.0095*pplay(ig,l)))
-            ztcond(ig,l)=
+            if (pplay(ig,l).ge.1e-4) then
+              ztcond(ig,l)=
      &         1./(bcond-acond*log(.01*vmr_co2(ig,l)*pplay(ig,l)))
-            if (pplay(ig,l).lt.1e-4) ztcond(ig,l)=0.0 !mars Monica
+            else
+              ztcond(ig,l)=0.0 !mars Monica
+            endif
          ENDDO
       ENDDO