Changeset 995

Timestamp:

Jun 21, 2013, 1:07:34 AM (11 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

Files:

• deftank/saturn1d/callphys.def (modified) (2 diffs)
• deftank/saturn1d/rcm1d.def (modified) (2 diffs)
• libf/phystd/calc_rayleigh.F90 (modified) (1 diff)
• libf/phystd/gfluxi.F (modified) (1 diff)

trunk/LMDZ.GENERIC/deftank/saturn1d/callphys.def

@@ -6 +6 @@
 diurnal   = .false. #.true.
 # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
-season    = .false.
+season    = .true. 
 # Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart
 tlocked   = .false.
@@ -32 +32 @@
 callrad    = .true.
 # the rad. transfer is computed every "iradia" physical timestep
-iradia     = 120 180 60 20 1 160
+iradia     = 160
 # call multilayer correlated-k radiative transfer ?
 corrk      = .true.

trunk/LMDZ.GENERIC/deftank/saturn1d/rcm1d.def

@@ -18 +18 @@
 year_day       = 24430. # Read System IIIw
 ## MIN DIST STAR-PLANET in AU [periastron]
-periastr       = 9.04
+periastr       = 9.02152
 ## MAX DIST STAR-PLANET in AU [apoastron]
-apoastr        = 10.12
+apoastr        = 10.05448
 ## DATE OF PERIASTRON in days
-peri_day       = 18700.
+peri_day       = 19280.
 ## OBLIQUITY in deg
 obliquit       = 26.73
 ## SURFACE PRESSURE in Pa
 psurf          = 1.e6
+
 
 ## Time integration parameters
@@ -37 +38 @@
 day_step       = 1
 # Number of sols to run 
-ndt            = 293160
+ndt            = 342020
 # Number of steps between each writing in diagfi
 ecritphy=5000

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)