import os
import sys
import numpy as np
import pencil as pc
# parameters unradable
ssat_const=1.63e-2
npswarm = 819200000.0
rho0 = 1.06
#
param = pc.read_param(quiet=True)
param2 = pc.read_param(quiet=True, param2=True)
pdim = pc.read_pdim()
dim = pc.read_dim()
#est=2.53e11*np.exp(-5420./293.25)
est=param2.c1*np.exp(-param2.c2/param.tt_const)
qvs=est/(param2.rv*rho0*param.tt_const)
supersaturation = ssat_const/qvs-1
print('supersaturation=', supersaturation)
#
condensationRate = 4*np.pi*param2.rhopmat/rho0*param2.g_condensation/(param.xyz1[0]*2/dim.nxgrid)**3*(supersaturation*param.a0_initdist*(pdim.npar/dim.nxgrid**3))
print('condensationRate =', condensationRate)