source: palm/trunk/SOURCE/init_cloud_physics.f90 @ 3

Last change on this file since 3 was 3, checked in by raasch, 15 years ago

RCS Log replace by Id keyword, revision history cleaned up

File size: 2.0 KB
Line 
1 SUBROUTINE init_cloud_physics
2
3!------------------------------------------------------------------------------!
4! Actual revisions:
5! -----------------
6!
7!
8! Former revisions:
9! ------------------
10! $Id$
11! RCS Log replace by Id keyword, revision history cleaned up
12!
13! Revision 1.5  2005/06/26 19:55:58  raasch
14! Initialization of cloud droplet constants, gas_constant renamed r_d,
15! latent_heat renamed l_v
16!
17! Revision 1.1  2000/04/13 14:37:22  schroeter
18! Initial revision
19!
20!
21! Description:
22! ------------
23! Initialization of parameters for handling cloud-physics
24!------------------------------------------------------------------------------!
25
26    USE arrays_3d
27    USE cloud_parameters
28    USE control_parameters
29    USE grid_variables
30    USE indices
31
32    IMPLICIT NONE
33
34    INTEGER ::  k
35    REAL    ::  t_surface
36
37    ALLOCATE( hydro_press(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
38
39!
40!-- Compute frequently used parameters
41    l_d_cp = l_v / cp
42    l_d_r  = l_v / r_d
43    l_d_rv = l_v / r_v
44
45!
46!-- Constant b in equation for droplet growth by condensation / evaporation.
47!-- Factor 1E-3 is needed because formula is in cgs units
48    mass_of_solute = 1.0E-17            ! in kg
49    molecular_weight_of_solute = 58.5   ! NaCl
50    b_cond = 4.3 * 2.0 * mass_of_solute / molecular_weight_of_solute * 1.0E-6
51
52!
53!-- Calculate:
54!-- pt / t : ratio of potential and actual temperature (pt_d_t)
55!-- t / pt : ratio of actual and potential temperature (t_d_pt)
56!-- p_0(z) : vertical profile of the hydrostatic pressure (hydro_press)
57    t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
58    DO  k = nzb, nzt+1
59       hydro_press(k) = surface_pressure * 100.0 * &
60                        ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
61       pt_d_t(k)      = ( 100000.0 / hydro_press(k) )**0.286
62       t_d_pt(k)      = 1.0 / pt_d_t(k)       
63    ENDDO
64
65!
66!-- Compute reference density
67    rho_surface = surface_pressure * 100.0 / ( r_d * t_surface )
68
69
70 END SUBROUTINE init_cloud_physics
Note: See TracBrowser for help on using the repository browser.