Changeset 3585 for trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

Timestamp:

Jan 20, 2025, 10:17:54 AM (9 hours ago)

Author:

debatzbr

Message:

Connecting microphysics to radiative transfer + miscellaneous cleans

File:

• trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90

@@ -12 +12 @@
   use radinc_h, only: ini_radinc_h, naerkind
   use radcommon_h, only: ini_radcommon_h
-  use radii_mod, only: radfixed, Nmix_n2
-  use datafile_mod, only: datadir,config_mufi,hazeprop_file,hazerad_file,hazemmr_file,hazedens_file
+  use radii_mod, only: radfixed
+  use datafile_mod, only: datadir,hazeprop_file,hazerad_file,hazemmr_file,hazedens_file, &
+                          config_mufi, aersprop_file, aerfprop_file
   use comdiurn_h, only: sinlat, coslat, sinlon, coslon
   use comgeomfi_h, only: totarea, totarea_planet
@@ -687 +688 @@
      call getin_p("config_mufi",config_mufi)
      if (is_master) write(*,*)" config_mufi = ",config_mufi
+
+     if (is_master) write(*,*) "Spherical aerosol optical properties datafile"
+     aersprop_file="optprop_rannou_r2-200nm_nu003.dat"  ! default file
+     call getin_p("aersprop_file",aersprop_file)
+     if (is_master) write(*,*) trim(rname)//" aersprop_file = ",trim(aersprop_file)
+
+     if (is_master) write(*,*) "Fractal aerosol optical properties datafile"
+     aerfprop_file="optprop_rannou_fractal_r010nm_N1_1e4_d2.dat"  ! default file
+     call getin_p("aerfprop_file",aerfprop_file)
+     if (is_master) write(*,*) trim(rname)//" aerfprop_file = ",trim(aerfprop_file)
  
      if (is_master) write(*,*) "Use haze production from CH4 photolysis or production rate?"
@@ -1285 +1296 @@
        call abort_physic(rname, 'if haze are on, tracers must be on!', 1)
      endif
+     if ((callmufi).and.(haze)) then
+       call abort_physic(rname, 'if haze are on, microphysics should be desactivated!', 1)
+     endif
+     if ((haze).and.(naerkind.gt.1)) then
+      call abort_physic(rname, 'if haze are on, naerkind must be equal to 1!', 1)
+     endif
+     if ((callmufi).and..not.(naerkind.gt.1)) then
+      call abort_physic(rname, 'if microphysics is on, naerkind must be > 1!', 1)
+     endif
+     if (.not.(callmufi.or.haze).and.(optichaze)) then
+      call abort_physic(rname, 'if microphysics and haze are off, optichaze must be desactivated!', 1)
+     endif
+     if ((callmufi.and.call_haze_prod_pCH4).and..not.(methane)) then
+      call abort_physic(rname, 'if haze production from CH4 photolysis is on, methane must be activated!', 1)
+     endif
      if (haze_proffix.and.sedimentation) then
          call abort_physic(rname, 'if haze profile is set, sedimentation must be deactivated', 1)