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

Timestamp:

Nov 20, 2017, 3:19:07 PM (7 years ago)

Author:

jvatant

Message:

Preliminary modifs for the optical coupling of haze
+ Moved inits of setspi/v before init of mufi
+ Added access to tarcers in optci/v
+ Some coherence in call to directories
JVO

File:

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

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

@@ -10 +10 @@
 
   use radinc_h, only: ini_radinc_h
-  use datafile_mod, only: datadir,config_mufi
+  use datafile_mod
   use comdiurn_h, only: sinlat, coslat, sinlon, coslon
   use comgeomfi_h, only: totarea, totarea_planet
@@ -220 +220 @@
      call getin_p("corrk",corrk)
      write(*,*) " corrk = ",corrk
+     
+     if (corrk) then
+       ! default path is set in datadir
+       write(*,*) "callcorrk: Correlated-k data base folder:",trim(datadir)
+       call getin_p("corrkdir",corrkdir)
+       write(*,*) " corrkdir = ",corrkdir
+     endif
+     
+     if (corrk .and. ngrid.eq.1) then
+       write(*,*) "simulate global averaged conditions ?"
+       global1d = .false. ! default value
+       call getin_p("global1d",global1d)
+       write(*,*) " global1d = ",global1d
+       
+       ! Test of incompatibility : if global1d is true, there should not be any diurnal cycle.
+       if (global1d.and.diurnal) then
+          write(*,*) "if global1d is true, diurnal must be set to false"
+          stop
+       endif
+
+       if (global1d) then
+          write(*,*) "Solar Zenith angle (deg.) ?"
+          write(*,*) "(assumed for averaged solar flux S/4)"
+          szangle=60.0  ! default value
+          call getin_p("szangle",szangle)
+          write(*,*) " szangle = ",szangle
+       endif
+     endif
 
      write(*,*) "prohibit calculations outside corrk T grid?"
@@ -349 +377 @@
      df_mufi=2.0 ! default value
      call getin_p("df_mufi",df_mufi)
+     write(*,*)" df_mufi = ",df_mufi
 
      write(*,*) "Monomer radius (m) ?"
      rm_mufi=6.66e-08 ! default value
      call getin_p("rm_mufi",rm_mufi)
+     write(*,*)" 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(*,*)" 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(*,*)" 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(*,*)" tx_prod = ",tx_prod
 
      write(*,*) "Equivalent radius production (m) ?"
      rc_prod=2.0e-8 ! default value
      call getin_p("rc_prod",rc_prod)
+     write(*,*)" rhc_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(*,*)" 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)
+     write(*,*)" config_mufi = ",config_mufi
 
 ! Soil model