Changeset 1025 for trunk/LMDZ.GENERIC/libf

Timestamp:

Aug 30, 2013, 3:47:25 PM (11 years ago)

Author:

sglmd

Message:

Helium case added

File:

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

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

@@ -77 +77 @@
                ! was out by a factor 2!
             elseif(igas.eq.igas_He)then
-               print*,'Helium not ready yet!'
-               tauconsti(igas) = 0.0
-               !call abort
+               tauconsti(igas) = (10.0/g)*0.00086*scalep/101325.0
             else
                print*,'Error in calc_rayleigh: Gas species not recognised!'
@@ -128 +126 @@
                      tauvari(igas) = 1.0/wl**4 
                   elseif(igas.eq.igas_He)then
-                     print*,'Helium not ready yet!'
-                     tauvari(igas) = 0.0
+                     tauvari(igas) = 1.0/wl**4 
                   else
                      print*,'Error in calc_rayleigh: Gas species not recognised!'