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

Last change on this file since 824 was 824, checked in by raasch, 10 years ago

preliminary checkin of new curvature/solution effects on droplet growth

  • Property svn:keywords set to Id
File size: 2.4 KB
Line 
1 SUBROUTINE init_cloud_physics
2
3!------------------------------------------------------------------------------!
4! Current revisions:
5! -----------------
6! calculation of b_cond replaced by calculation of bfactor
7!
8! Former revisions:
9! ------------------
10! $Id: init_cloud_physics.f90 824 2012-02-17 09:09:57Z raasch $
11!
12! 221 2009-01-12 15:32:23Z raasch
13! Bugfix: abort in case that absolute temperature is below zero
14!
15! 95 2007-06-02 16:48:38Z raasch
16! hydro_press renamed hyp
17!
18! February 2007
19! RCS Log replace by Id keyword, revision history cleaned up
20!
21! Revision 1.5  2005/06/26 19:55:58  raasch
22! Initialization of cloud droplet constants, gas_constant renamed r_d,
23! latent_heat renamed l_v
24!
25! Revision 1.1  2000/04/13 14:37:22  schroeter
26! Initial revision
27!
28!
29! Description:
30! ------------
31! Initialization of parameters for handling cloud-physics
32!------------------------------------------------------------------------------!
33
34    USE arrays_3d
35    USE cloud_parameters
36    USE constants
37    USE control_parameters
38    USE grid_variables
39    USE indices
40
41    IMPLICIT NONE
42
43    INTEGER ::  k
44    REAL    ::  t_surface
45
46    ALLOCATE( hyp(nzb:nzt+1), pt_d_t(nzb:nzt+1), t_d_pt(nzb:nzt+1) )
47
48!
49!-- Calculate frequently used parameters
50    l_d_cp = l_v / cp
51    l_d_r  = l_v / r_d
52    l_d_rv = l_v / r_v
53!
54!-- Calculate factor used in equation for droplet growth by condensation
55    bfactor = 3.0 * vanthoff * mass_of_solute * molecular_weight_of_water &
56              / ( 4.0 * pi * rho_l * molecular_weight_of_solute )
57!
58!-- Calculate:
59!-- pt / t : ratio of potential and actual temperature (pt_d_t)
60!-- t / pt : ratio of actual and potential temperature (t_d_pt)
61!-- p_0(z) : vertical profile of the hydrostatic pressure (hyp)
62    t_surface = pt_surface * ( surface_pressure / 1000.0 )**0.286
63    DO  k = nzb, nzt+1
64!
65!--    Check temperature in case of too large domain height
66       IF ( ( t_surface - g/cp * zu(k) ) < 0.0 )  THEN
67          WRITE( message_string, * )  'absolute temperature < 0.0 at zu(', k, &
68                                      ') = ', zu(k)
69          CALL message( 'init_cloud_physics', 'PA0142', 1, 2, 0, 6, 0 )
70       ENDIF
71       hyp(k)    = surface_pressure * 100.0 * &
72                   ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
73       pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
74       t_d_pt(k) = 1.0 / pt_d_t(k)       
75    ENDDO
76
77!
78!-- Compute reference density
79    rho_surface = surface_pressure * 100.0 / ( r_d * t_surface )
80
81
82 END SUBROUTINE init_cloud_physics
Note: See TracBrowser for help on using the repository browser.