Changeset 265 for codes/icosagcm/branches

Timestamp:

08/01/14 14:20:46 (10 years ago)

Author:

ymipsl

Message:

Fix for specific saturn physics

YM+AS

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/gfluxi.F (modified) (1 diff)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/gfluxi.F

@@ -104 +104 @@
       ! -- same results for most thin atmospheres
       ! -- and stabilizes integrations
-      NT    = int(TLEV(2*L+1)*NTfac) - NTstar+1
+      ! NT    = int(TLEV(2*L+1)*NTfac) - NTstar+1
       !! For deep, opaque, thick first layers (e.g. Saturn)
       !! what is below works much better, not unstable, ...
       !! ... and actually fully accurate because 1st layer temp (JL) 
-      !NT    = int(TLEV(2*L)*NTfac) - NTstar+1
+      NT    = int(TLEV(2*L)*NTfac) - NTstar+1
       !! (or this one yields same results
       !NT    = int( (TLEV(2*L)+TLEV(2*L+1))*0.5*NTfac ) - NTstar+1