Changeset 824 for palm/trunk/SOURCE/modules.f90

Timestamp:

Feb 17, 2012 9:09:57 AM (12 years ago)

Author:

raasch

Message:

preliminary checkin of new curvature/solution effects on droplet growth

File:

• palm/trunk/SOURCE/modules.f90 (modified) (6 diffs)

palm/trunk/SOURCE/modules.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! +bfactor, curvature_solution_effects, eps_ros, molecular_weight_of_water,
+! vanthoff, -b_cond in cloud_parameters,
+! dopts_num increased to 27, particle attributes speed_x|y|z_sgs renamed
+! rvar1|2|3
 !
 ! Former revisions:
@@ -332 +335 @@
     REAL, DIMENSION(:,:,:,:), ALLOCATABLE ::  rif_wall
 
-    REAL, DIMENSION(:,:,:), ALLOCATABLE :: var_x, var_y, var_z, gamma_x,        &
+    REAL, DIMENSION(:,:,:), ALLOCATABLE :: var_x, var_y, var_z, gamma_x,       &
                                            gamma_y, gamma_z
 
@@ -369 +372 @@
 !------------------------------------------------------------------------------!
 
-    REAL  ::  b_cond, cp = 1005.0, diff_coeff_l = 0.23E-4,                     &
-              effective_coll_efficiency, l_d_cp, l_d_r, l_d_rv, l_v = 2.5E+06, &
-              mass_of_solute, molecular_weight_of_solute,                      &
+    LOGICAL ::  curvature_solution_effects = .FALSE.
+
+    REAL  ::  bfactor, cp = 1005.0, diff_coeff_l = 0.23E-4,                    &
+              effective_coll_efficiency,                                       &
+              eps_ros = 1.0E-4,             &  ! accuracy of Rosenbrock method
+              l_d_cp, l_d_r, l_d_rv, l_v = 2.5E+06,                            &
+              mass_of_solute = 1.0E-17,              & ! soluted NaCl in kg
+              molecular_weight_of_solute = 0.05844,  & ! mol. mass NaCl (kg/mol)
+              molecular_weight_of_water = 0.01801528,& ! mol. mass H2O (kg/mol)
               prec_time_const = 0.001, ql_crit = 0.0005, rho_l = 1.0E3,        &
-              r_d = 287.0, r_v = 461.51, thermal_conductivity_l = 2.43E-2
+              r_d = 287.0, r_v = 461.51, thermal_conductivity_l = 2.43E-2,     &
+              vanthoff = 2.0                        ! van't Hoff factor (NaCl)
 
     REAL, DIMENSION(:), ALLOCATABLE   ::  pt_d_t, t_d_pt
@@ -1037 +1047 @@
 #endif
 
-    INTEGER, PARAMETER ::  dopr_norm_num = 7, dopts_num = 26, dots_max = 100
+    INTEGER, PARAMETER ::  dopr_norm_num = 7, dopts_num = 27, dots_max = 100
 
     INTEGER ::  dots_num = 23
@@ -1052 +1062 @@
           (/ 'tnpt   ', 'x_     ', 'y_     ', 'z_     ', 'z_abs  ', 'u      ', &
              'v      ', 'w      ', 'u"     ', 'v"     ', 'w"     ', 'npt_up ', &
-             'w_up   ', 'w_down ', 'npt_max', 'npt_min', 'x*2    ', 'y*2    ', &
-             'z*2    ', 'u*2    ', 'v*2    ', 'w*2    ', 'u"2    ', 'v"2    ', &
-             'w"2    ', 'npt*2  ' /)
+             'w_up   ', 'w_down ', 'radius ', 'npt_max', 'npt_min', 'x*2    ', &
+             'y*2    ', 'z*2    ', 'u*2    ', 'v*2    ', 'w*2    ', 'u"2    ', &
+             'v"2    ', 'w"2    ', 'npt*2  ' /)
 
     CHARACTER (LEN=7), DIMENSION(dopts_num) :: dopts_unit = &
           (/ 'number ', 'm      ', 'm      ', 'm      ', 'm      ', 'm/s    ', &
              'm/s    ', 'm/s    ', 'm/s    ', 'm/s    ', 'm/s    ', 'number ', &
-             'm/s    ', 'm/s    ', 'number ', 'number ', 'm2     ', 'm2     ', &
-             'm2     ', 'm2/s2  ', 'm2/s2  ', 'm2/s2  ', 'm2/s2  ', 'm2/s2  ', &
-             'm2/s2  ', 'number2' /)
+             'm/s    ', 'm/s    ', 'm      ', 'number ', 'number ', 'm2     ', &
+             'm2     ', 'm2     ', 'm2/s2  ', 'm2/s2  ', 'm2/s2  ', 'm2/s2  ', &
+             'm2/s2  ', 'm2/s2  ', 'number2' /)
 
     CHARACTER (LEN=7), DIMENSION(dots_max) :: dots_label = &
@@ -1230 +1240 @@
        SEQUENCE
        REAL    ::  age, age_m, dt_sum, dvrp_psize, e_m, origin_x, origin_y, &
-                   origin_z, radius, speed_x, speed_x_sgs, speed_y,         &
-                   speed_y_sgs, speed_z, speed_z_sgs, weight_factor, x, y, z
+                   origin_z, radius, rvar1, rvar2, rvar3, speed_x, speed_y, &
+                   speed_z, weight_factor, x, y, z
        INTEGER ::  color, group, tailpoints, tail_id
     END TYPE particle_type