Changeset 3654 for trunk/LMDZ.GENERIC/libf/phystd/setspv.F90

Timestamp:

Feb 27, 2025, 11:35:32 AM (4 months ago)

Author:

gmilcareck

Message:

Generic PCM:
Rewriting of the rayleigh scattering module. Before the commit, the Rayleigh scattering was separated into a constant part and a wavelength-dependent part.
However, in the constant part, because everything was added up/multiplied together before being added to the module,
we lose the information about where these numbers come from and if we want to add a new molecule or update a molecule... well... good luck.
Now, each molecule depends only on 2 physical parameters: the refractive index and the King Factor.
If you want add another molecule, you just need to add these 2 parameters and the references (please).
For water, the rayleigh scattering now depends on temperature and pressure.
GM

File:

• trunk/LMDZ.GENERIC/libf/phystd/setspv.F90 (modified) (3 diffs)

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

@@ -5 +5 @@
 !     Purpose
 !     -------
-!     Set up spectral intervals, stellar spectrum and Rayleigh 
-!     opacity in the shortwave. 
+!     Set up spectral intervals and stellar spectrum in the shortwave. 
 !     
 !     Authors
@@ -25 +24 @@
       use radinc_h,    only: L_NSPECTV, corrkdir, banddir
       use radcommon_h, only: BWNV,BLAMV,WNOV,DWNV,WAVEV, &
-                             STELLARF,TAURAY
+                             STELLARF
       use datafile_mod, only: datadir
-      use callkeys_mod, only: Fat1AU,rayleigh
+      use callkeys_mod, only: Fat1AU
 
       implicit none
@@ -133 +132 @@
       print*,' '
 
-
-!=======================================================================
-!     Set up the wavelength independent part of the Rayleigh scattering.
-!     The pressure dependent part will be computed elsewhere (OPTCV).
-!     WAVEV is in microns.  There is no Rayleigh scattering in the IR.
-
-      if(rayleigh) then
-         call calc_rayleigh
-      else
-         print*,'setspv: No Rayleigh scattering, check for NaN in output!'
-         do N=1,L_NSPECTV
-            TAURAY(N) = 1E-16
-         end do
-      endif
-
       RETURN
     END subroutine setspv