Changeset 2831 for trunk/LMDZ.GENERIC/libf/phystd/suaer_corrk.F90

Timestamp:

Nov 23, 2022, 4:41:34 PM (2 years ago)

Author:

emillour

Message:

Generic PCM:
Add the possibility to include Venus-like aerosols (triggered by option
aerovenus=.true. in callphys.def); baseline is to use 5 distinct scatterers
but each may be turned on/off (via aerovenus1, aerovenus2, aerovenus2p,
aerovenus3, aerovenusUV flags which may be specified in callphys.def).
GG

File:

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

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

@@ -232 +232 @@
 ! added by SG
        endif
+
+! VENUS CLOUDS
+
+       if (iaer.eq.iaero_venus1) then
+         print*, 'naerkind= venus1', iaer
+
+!     visible
+         file_id(iaer,1) = 'optprop_h2so4vis_n50.dat'
+         lamrefvis(iaer)=1.5E-6   ! no idea, must find
+!     infrared
+         file_id(iaer,2) = 'optprop_h2so4ir_n50.dat'
+         lamrefir(iaer)=9.3E-6 ! no idea, must find
+! added by SL
+       endif
+
+       if (iaer.eq.iaero_venus2) then
+         print*, 'naerkind= venus2', iaer
+
+!     visible
+         file_id(iaer,1) = 'optprop_h2so4vis_n50.dat'
+         lamrefvis(iaer)=1.5E-6   ! no idea, must find
+!     infrared
+         file_id(iaer,2) = 'optprop_h2so4ir_n50.dat'
+         lamrefir(iaer)=9.3E-6 ! no idea, must find
+! added by SL
+       endif
+
+       if (iaer.eq.iaero_venus2p) then
+         print*, 'naerkind= venus2p', iaer
+
+!     visible
+         file_id(iaer,1) = 'optprop_h2so4vis_n50.dat'
+         lamrefvis(iaer)=1.5E-6   ! no idea, must find
+!     infrared
+         file_id(iaer,2) = 'optprop_h2so4ir_n50.dat'
+         lamrefir(iaer)=9.3E-6 ! no idea, must find
+! added by SL
+       endif
+
+       if (iaer.eq.iaero_venus3) then
+         print*, 'naerkind= venus3', iaer
+
+!     visible
+         file_id(iaer,1) = 'optprop_h2so4vis_n50.dat'
+         lamrefvis(iaer)=1.5E-6   ! no idea, must find
+!     infrared
+         file_id(iaer,2) = 'optprop_h2so4ir_n50.dat'
+         lamrefir(iaer)=9.3E-6 ! no idea, must find
+! added by SL
+       endif
+
+       if (iaer.eq.iaero_venusUV) then
+         print*, 'naerkind= venusUV', iaer
+
+!     visible
+         file_id(iaer,1) = 'optprop_venusUVvis.dat'
+         lamrefvis(iaer)=3.5E-7   ! Haus et al. 2015
+!     infrared
+         file_id(iaer,2) = 'optprop_venusUVir.dat'
+         lamrefir(iaer)=9.3E-6 ! not used anyway
+! added by SL
+       endif
+
+! END VENUS CLOUDS
+       
 ! the following was added by LT
        do ia=1,aerogeneric ! Read Radiative Generic Condensable Species data
@@ -265 +330 @@
          endif
        enddo ! ia=1,aerogeneric
-      enddo
+      enddo ! of do iaer=1,naerkind
       
 !------------------------------------------------------------------