- Timestamp:
- 08/01/14 14:20:46 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/gfluxi.F
r222 r265 104 104 ! -- same results for most thin atmospheres 105 105 ! -- and stabilizes integrations 106 NT = int(TLEV(2*L+1)*NTfac) - NTstar+1106 ! NT = int(TLEV(2*L+1)*NTfac) - NTstar+1 107 107 !! For deep, opaque, thick first layers (e.g. Saturn) 108 108 !! what is below works much better, not unstable, ... 109 109 !! ... and actually fully accurate because 1st layer temp (JL) 110 !NT = int(TLEV(2*L)*NTfac) - NTstar+1110 NT = int(TLEV(2*L)*NTfac) - NTstar+1 111 111 !! (or this one yields same results 112 112 !NT = int( (TLEV(2*L)+TLEV(2*L+1))*0.5*NTfac ) - NTstar+1
Note: See TracChangeset
for help on using the changeset viewer.