Changeset 1819 for trunk/LMDZ.TITAN/libf/muphytitan/mm_haze.f90

Timestamp:

Nov 8, 2017, 4:45:50 PM (7 years ago)

Author:

jvatant

Message:

Commit various bugs leading to negative tracer tendencies
Also deleted some useless cfg files
--JVO

File:

• trunk/LMDZ.TITAN/libf/muphytitan/mm_haze.f90 (modified) (1 diff)

trunk/LMDZ.TITAN/libf/muphytitan/mm_haze.f90

@@ -743 +743 @@
         IF ((ratio <= 0.9_mm_wp .OR. ratio >= 1.1_mm_wp) .AND. wth(i) /= 0._mm_wp) THEN
           wdt = wth(i)*dt
-          zcorf(i) = dzb/wdt * (exp(-wdt*log(ratio)/dzb)-1._mm_wp) / (1._mm_wp-ratio)
+          ! bugfix: max exponential arg to 30)
+          zcorf(i) = dzb/wdt * (exp(MIN(30._mm_wp,-wdt*log(ratio)/dzb))-1._mm_wp) / (1._mm_wp-ratio)
+          !zcorf(i) = dzb/wdt * (exp(-wdt*log(ratio)/dzb)-1._mm_wp) / (1._mm_wp-ratio)
         ENDIF
       ENDIF