Changeset 995 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jun 21, 2013, 1:07:34 AM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. Updated saturn 1D files. Added comments around the gluxi problem (will be further investigated). removed abort because Helium rayleigh scattering is missing

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• calc_rayleigh.F90 (modified) (1 diff)
• gfluxi.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/calc_rayleigh.F90

@@ -77 +77 @@
                print*,'Helium not ready yet!'
                tauconsti(igas) = 0.0
-               call abort
+               !call abort
             else
                print*,'Error in calc_rayleigh: Gas species not recognised!'

trunk/LMDZ.GENERIC/libf/phystd/gfluxi.F

@@ -101 +101 @@
       LAMDA(L) = ALPHA(L)*(1.0D0-W0(L)*COSBAR(L))/UBARI
 
+      ! Tsurf is used for 1st layer source function
+      ! -- same results for most thin atmospheres
+      ! -- and stabilizes integrations
       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
+      !! (or this one yields same results
+      !NT    = int( (TLEV(2*L)+TLEV(2*L+1))*0.5*NTfac ) - NTstar+1
+
       NT2   = int(TLEV(2*L)*NTfac)   - NTstar+1
       B1(L) = (PLANCKIR(NW,NT)-PLANCKIR(NW,NT2))/DTAU(L)