Venus PCM - Microphysical scheme

From Planets
Jump to: navigation, search
 ! This is a work in progress to describe the microphysical scheme used in the Venus PCM ! 


How to activate the microphysical scheme

In the physiq.def file

In the physiq.def file, you need to modify the flag :

cl_scheme = 1 to cl_scheme = 2

A new set of flags are now useful to you :

debug_cloud = n -> get you a full set of information on what is happening in the microphysical scheme.  ! Be careful it is intended for 1D for now, you can expect a lot of data in 3D (to be improved) !


A set of flags to activate or not the individual processes of the microphysical scheme

mode_merging = y

hom_nuclea = y

het_nuclea = y

cond_evap =y

coag = y


You are able to modify the standard deviation of the size distribution of the mode 1 and mode 2 droplets

sigma_mode1 = 1.56

sigma_mode2 = 1.29


Finally you have a set of flags and free parameters to change of the aerosols production (which serves as possible condensation nuclei for the droplets)

First concerning the properties of the aerosols

prod_aer = y If you want to include the aerosols or not

rad_aer = 0.125e-6 (m) The mean radius of the aerosols produced

sigma_aer = 1.56 The standard deviation of the size distribution of the aerosols

rho_aer = 1800 (kg/m3) The density of the aerosols

And then concerning the spatial distribution of the aerosols production. We're making the assumption that the vertical production follows a gaussian shape.

alt_aer = 62000 (m) Altitude of the peak of production of the aerosols

sigmaz_aer = 2000 (m) Vertical standard deviation of the gaussian distribution

tx_prod_aer = 2.D-10 (/m2/s) Vertically integrated production rate of the aerosols


In the traceurs.def file

You need to add a set of new traceurs ! You need to add it at the end of the traceurs list (Until the rework of the traceurs in the Venus PCM) !

	
10 10 m0_aer
10 10 m3_aer
10 10 m0_mode1drop
10 10 m0_mode1ccn
10 10 m3_mode1sa
10 10 m3_mode1wv
10 10 m3_mode1ccn
10 10 m0_mode2drop
10 10 m0_mode2ccn
10 10 m3_mode2sa
10 10 m3_mode2wv
10 10 m3_mode2ccn