Changeset 2831 for trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.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/inifis_mod.F90 (modified) (1 diff)

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

@@ -687 +687 @@
      endif
 
+     if (is_master) write(*,*)"Radiatively active Venus clouds ?"
+     aerovenus=.false. ! default value
+     call getin_p("aerovenus",aerovenus)
+     if (aerovenus) then
+       aerovenus1=.true.     ! default value
+       aerovenus2=.true.     ! default value
+       aerovenus2p=.true.     ! default value
+       aerovenus3=.true.     ! default value
+       aerovenusUV=.true.     ! default value
+     else
+       aerovenus1=.false.     ! default value
+       aerovenus2=.false.     ! default value
+       aerovenus2p=.false.     ! default value
+       aerovenus3=.false.     ! default value
+       aerovenusUV=.false.     ! default value
+     endif
+     ! in case the user wants to specifically set/unset sub-options
+     call getin_p("aerovenus1",aerovenus1)
+     if (is_master) write(*,*)" aerovenus1 = ",aerovenus1
+     call getin_p("aerovenus2",aerovenus2)
+     if (is_master) write(*,*)" aerovenus2 = ",aerovenus2
+     call getin_p("aerovenus2p",aerovenus2p)
+     if (is_master) write(*,*)" aerovenus2p= ",aerovenus2p
+     call getin_p("aerovenus3",aerovenus3)
+     if (is_master) write(*,*)" aerovenus3 = ",aerovenus3
+     call getin_p("aerovenusUV",aerovenusUV)
+     if (is_master) write(*,*)" aerovenusUV= ",aerovenusUV
+     ! Sanity check: if any of the aerovenus* is set to true
+     ! then aeronovenus should also be set to true
+     if ((.not.aerovenus).and.(aerovenus1.or.aerovenus2.or.aerovenus2p.or.&
+                               aerovenus3.or.aerovenusUV)) then
+      if(is_master) then
+       write(*,*)" Error, if you set some of the aerovenus* to true"
+       write(*,*)" then flag aerovenus should be set to true as well!"
+      endif
+      call abort_physic("inifis"," aerovenus* flags mismatch!",1)
+     endif
+     
      if (is_master) write(*,*)trim(rname)//&
        ": TWOLAY AEROSOL: total optical depth "//&