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

Timestamp:

Sep 26, 2017, 12:43:42 PM (7 years ago)

Author:

jvatant

Message:

Get rid of all the old-generic dummy aerosol scheme ( just left scatterers_h for compilation )
as the new microphysics for Titan will be plugged in
-> even removed sedimentation ( will be done in the microphysical model )
--JVO

File:

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

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

@@ -9 +9 @@
              prad,pg,pr,pcpp)
 
-  use radinc_h, only: ini_radinc_h, naerkind
+  use radinc_h, only: ini_radinc_h
   use datafile_mod, only: datadir
   use comdiurn_h, only: sinlat, coslat, sinlon, coslon
@@ -385 +385 @@
      endif
 
-     write(*,*)"Default planetary temperature?"
-     tplanet=215.0
-     call getin_p("tplanet",tplanet)
-     write(*,*)" tplanet = ",tplanet
-
      write(*,*)"Which star?"
      startype=1 ! default value = Sol
@@ -399 +394 @@
      call getin_p("Fat1AU",Fat1AU)
      write(*,*)" Fat1AU = ",Fat1AU
-
-
-! TRACERS:
-
-!         write(*,*)"Number of radiatively active aerosols:"
-!         naerkind=0. ! default value
-!         call getin_p("naerkind",naerkind)
-!         write(*,*)" naerkind = ",naerkind
-
- 
-!=================================
-
-     write(*,*)"Radiatively active two-layer aersols?"
-     aeroback2lay=.false.     ! default value
-     call getin_p("aeroback2lay",aeroback2lay)
-     write(*,*)" aeroback2lay = ",aeroback2lay
-
-     write(*,*)"TWOLAY AEROSOL: total optical depth ", &
-                    "in the tropospheric layer (visible)"
-     obs_tau_col_tropo=8.D0
-     call getin_p("obs_tau_col_tropo",obs_tau_col_tropo)
-     write(*,*)" obs_tau_col_tropo = ",obs_tau_col_tropo
-
-     write(*,*)"TWOLAY AEROSOL: total optical depth ", &
-                    "in the stratospheric layer (visible)"
-     obs_tau_col_strato=0.08D0
-     call getin_p("obs_tau_col_strato",obs_tau_col_strato)
-     write(*,*)" obs_tau_col_strato = ",obs_tau_col_strato
-
-     write(*,*)"TWOLAY AEROSOL: pres_bottom_tropo? in pa"
-     pres_bottom_tropo=66000.0
-     call getin_p("pres_bottom_tropo",pres_bottom_tropo)
-     write(*,*)" pres_bottom_tropo = ",pres_bottom_tropo
-
-     write(*,*)"TWOLAY AEROSOL: pres_top_tropo? in pa"
-     pres_top_tropo=18000.0
-     call getin_p("pres_top_tropo",pres_top_tropo)
-     write(*,*)" pres_top_tropo = ",pres_top_tropo
-
-     write(*,*)"TWOLAY AEROSOL: pres_bottom_strato? in pa"
-     pres_bottom_strato=2000.0
-     call getin_p("pres_bottom_strato",pres_bottom_strato)
-     write(*,*)" pres_bottom_strato = ",pres_bottom_strato
-
-     write(*,*)"TWOLAY AEROSOL: pres_top_strato? in pa"
-     pres_top_strato=100.0
-     call getin_p("pres_top_strato",pres_top_strato)
-     write(*,*)" pres_top_strato = ",pres_top_strato
-
-     write(*,*)"TWOLAY AEROSOL: particle size in the ", &
-                    "tropospheric layer, in meters"
-     size_tropo=2.e-6
-     call getin_p("size_tropo",size_tropo)
-     write(*,*)" size_tropo = ",size_tropo
-
-     write(*,*)"TWOLAY AEROSOL: particle size in the ", &
-                    "stratospheric layer, in meters"
-     size_strato=1.e-7
-     call getin_p("size_strato",size_strato)
-     write(*,*)" size_strato = ",size_strato
-
-!=================================
-
-     write(*,*) "Gravitationnal sedimentation ?"
-     sedimentation=.false. ! default value
-     call getin_p("sedimentation",sedimentation)
-     write(*,*) " sedimentation = ",sedimentation       
 
      write(*,*) "Does user want to force cpp and mugaz?"