Changeset 1795 for trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

Timestamp:

Oct 12, 2017, 12:26:18 PM (7 years ago)

Author:

jvatant

Message:

Making Titan's hazy again, part II
+ Added calmufi and inimufi routines that interface YAMMS model
+ Major changes of the tracer gestion in tracer_h (new query by name)
+ Update the deftank
JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

@@ -10 +10 @@
 
   use radinc_h, only: ini_radinc_h
-  use datafile_mod, only: datadir
+  use datafile_mod, only: datadir,config_mufi
   use comdiurn_h, only: sinlat, coslat, sinlon, coslon
   use comgeomfi_h, only: totarea, totarea_planet
@@ -319 +319 @@
      call getin_p("ichim",ichim)
      write(*,*)" ichim = ",ichim
+
+! Microphysics
+
+     write(*,*) "Run with or without microphysics?"
+     callmufi=.false. ! default value
+     call getin_p("callmufi",callmufi)
+     write(*,*)" callmufi = ",callmufi
+
+     ! sanity check
+     if (callmufi.and.(.not.tracer)) then
+       print*,"You are running microphysics without tracer"
+       print*,"Please start again with tracer =.true."
+       stop
+     endif
+
+     write(*,*) "Compute clouds?"
+     callclouds=.false. ! default value
+     call getin_p("callclouds",callclouds)
+     write(*,*)" callclouds = ",callclouds
+
+     ! sanity check
+     if (callclouds.and.(.not.callmufi)) then
+       print*,"You are trying to make clouds without microphysics !"
+       print*,"Please start again with callmufi =.true."
+       stop
+     endif
+
+     write(*,*) "Fractal dimension ?"
+     df_mufi=2.0 ! default value
+     call getin_p("df_mufi",df_mufi)
+
+     write(*,*) "Monomer radius (m) ?"
+     rm_mufi=6.66e-08 ! default value
+     call getin_p("rm_mufi",rm_mufi)
+
+     write(*,*) "Aerosol density (kg.m-3)?"
+     rho_aer_mufi=1.e3 ! default value
+     call getin_p("rho_aer_mufi",rho_aer_mufi)
+
+     write(*,*) "Pressure level of aer. production (Pa) ?"
+     p_prod=1.0 ! default value
+     call getin_p("p_prod",p_prod)
+     
+     write(*,*) "Aerosol production rate (kg.m-2.s-1) ?"
+     tx_prod=3.5e-13 ! default value
+     call getin_p("tx_prod",tx_prod)
+
+     write(*,*) "Equivalent radius production (m) ?"
+     rc_prod=2.0e-8 ! default value
+     call getin_p("rc_prod",rc_prod)
+
+     write(*,*) "Radius of air (nitrogen) molecule (m) ?"
+     air_rad=1.75e-10 ! default value
+     call getin_p("air_rad",air_rad)
+
+     write(*,*) "Path to microphys. config file ?"
+     config_mufi='datagcm/microphysics/config.cfg' ! default value
+     call getin_p("config_mufi",config_mufi)
 
 ! Soil model