Changeset 716 for trunk/LMDZ.GENERIC/libf/phystd/calc_rayleigh.F90

Timestamp:

Jul 3, 2012, 8:09:54 PM (12 years ago)

Author:

rwordsworth

Message:

Mainly updates to radiative transfer and gas management scheme.
Most CIA data now read from standard HITRAN datafiles. For the H2O
continuum, two options have been added: the standard CKD continuum,
and the empirical formula in PPC (Pierrehumbert 2010). Use the toggle
'H2Ocont_simple' in callphys.def to choose.

Note to Martians: I've changed the default values of 'sedimentation' and
'co2cond' in inifis to false. Both these are defined in the standard deftank
callphys.def file, so there shouldn't be any compatibility problems.

File:

• trunk/LMDZ.GENERIC/libf/phystd/calc_rayleigh.F90 (modified) (2 diffs)

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

@@ -74 +74 @@
             elseif(gnom(igas).eq.'He_')then
                print*,'Helium not ready yet!'
+               tauconsti(igas) = 0.0
+               call abort
             else
                print*,'Error in calc_rayleigh: Gas species not recognised!'
@@ -117 +119 @@
                   elseif(gnom(igas).eq.'H2_')then
                      tauvari(igas) = 1.0/wl**4 
+                  elseif(gnom(igas).eq.'He_')then
+                     print*,'Helium not ready yet!'
+                     tauvari(igas) = 0.0
                   else
                      print*,'Error in calc_rayleigh: Gas species not recognised!'