Changeset 1053 for palm/trunk/SOURCE/read_var_list.f90

Timestamp:

Nov 13, 2012 5:11:03 PM (11 years ago)

Author:

hoffmann

Message:

two-moment cloud physics implemented

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/read_var_list.f90 (modified) (12 diffs)

palm/trunk/SOURCE/read_var_list.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! necessary expansions according to the two new prognostic equations (nr, qr) 
+! of the two-moment cloud physics scheme:
+! +bc_*_b, +bc_*_t, +bc_*_t_val, *_init, *_surface, *_surface_initial_change,
+! +*_vertical_gradient, +*_vertical_gradient_level, *_vertical_gradient_level_ind,
+! +surface_waterflux_*
+!
+! in addition, steering parameters parameters of the two-moment cloud physics 
+! scheme:   
+! +cloud_scheme, +drizzle, +mu_constant, +mu_constant_value, +ventilation_effect
 !
 ! Former revisions:
@@ -169 +177 @@
 !-- Make version number check first
     READ ( 13 )  version_on_file
-    binary_version = '3.5'
+    binary_version = '3.6'
     IF ( TRIM( version_on_file ) /= TRIM( binary_version ) )  THEN
        WRITE( message_string, * ) 'version mismatch concerning control ', &
@@ -235 +243 @@
        ALLOCATE( lad(0:nz+1), ug(0:nz+1), u_init(0:nz+1), vg(0:nz+1),    &
                  v_init(0:nz+1), pt_init(0:nz+1), q_init(0:nz+1),        &
-                 sa_init(0:nz+1),                                        &
+                 qr_init(0:nz+1), nr_init(0:nz+1), sa_init(0:nz+1),      &
                  hom(0:nz+1,2,pr_palm+max_pr_user,0:statistic_regions),  &
                  hom_sum(0:nz+1,pr_palm+max_pr_user,0:statistic_regions) )
@@ -268 +276 @@
           CASE ( 'bc_ns' )
              READ ( 13 )  bc_ns
+          CASE ( 'bc_nr_b' )
+             READ ( 13 )  bc_nr_b
+          CASE ( 'bc_nr_t' )
+             READ ( 13 )  bc_nr_t
+          CASE ( 'bc_nr_t_val' )
+             READ ( 13 )  bc_nr_t_val
           CASE ( 'bc_p_b' )
              READ ( 13 )  bc_p_b
@@ -284 +298 @@
           CASE ( 'bc_q_t_val' )
              READ ( 13 )  bc_q_t_val
+          CASE ( 'bc_qr_b' )
+             READ ( 13 )  bc_qr_b
+          CASE ( 'bc_qr_t' )
+             READ ( 13 )  bc_qr_t
+          CASE ( 'bc_qr_t_val' )
+             READ ( 13 )  bc_qr_t_val
           CASE ( 'bc_s_b' )
              READ ( 13 )  bc_s_b
@@ -326 +346 @@
           CASE ( 'cloud_physics' )
              READ ( 13 )  cloud_physics
+          CASE ( 'cloud_scheme' )
+             READ ( 13 )  cloud_scheme
           CASE ( 'collective_wait' )
              READ ( 13 )  collective_wait
@@ -356 +378 @@
           CASE ( 'drag_coefficient' )
              READ ( 13 )  drag_coefficient
+          CASE ( 'drizzle' )
+             READ ( 13 )  drizzle
           CASE ( 'dt_pr_1d' )
              READ ( 13 )  dt_pr_1d
@@ -438 +462 @@
           CASE ( 'momentum_advec' )
              READ ( 13 )  momentum_advec
+          CASE ( 'mu_constant' )
+             READ ( 13 )  mu_constant
+          CASE ( 'mu_constant_value' )
+             READ ( 13 )  mu_constant_value
           CASE ( 'netcdf_precision' )
              READ ( 13 )  netcdf_precision
@@ -444 +472 @@
           CASE ( 'ngsrb' )
              READ ( 13 )  ngsrb
+          CASE ( 'nr_init' )
+             READ ( 13 )  nr_init
+          CASE ( 'nr_surface' )
+             READ ( 13 )  nr_surface
+          CASE ( 'nr_surface_initial_change' )
+             READ ( 13 )  nr_surface_initial_change
+          CASE ( 'nr_vertical_gradient' )
+             READ ( 13 )  nr_vertical_gradient
+          CASE ( 'nr_vertical_gradient_level' )
+             READ ( 13 )  nr_vertical_gradient_level
+          CASE ( 'nr_vertical_gradient_level_ind' )
+             READ ( 13 )  nr_vertical_gradient_level_ind
           CASE ( 'nsor' )
              READ ( 13 )  nsor
@@ -510 +550 @@
           CASE ( 'q_vertical_gradient_level_ind' )
              READ ( 13 )  q_vertical_gradient_level_ind
+          CASE ( 'qr_init' )
+             READ ( 13 )  qr_init
+          CASE ( 'qr_surface' )
+             READ ( 13 )  qr_surface
+          CASE ( 'qr_surface_initial_change' )
+             READ ( 13 )  qr_surface_initial_change
+          CASE ( 'qr_vertical_gradient' )
+             READ ( 13 )  qr_vertical_gradient
+          CASE ( 'qr_vertical_gradient_level' )
+             READ ( 13 )  qr_vertical_gradient_level
+          CASE ( 'qr_vertical_gradient_level_ind' )
+             READ ( 13 )  qr_vertical_gradient_level_ind
           CASE ( 'radiation' )
              READ ( 13 )  radiation
@@ -555 +607 @@
              READ ( 13 )  surface_scalarflux             
           CASE ( 'surface_waterflux' )
-             READ ( 13 )  surface_waterflux             
+             READ ( 13 )  surface_waterflux      
+          CASE ( 'surface_waterflux_nr' )
+             READ ( 13 )  surface_waterflux_nr    
+          CASE ( 'surface_waterflux_qr' )
+             READ ( 13 )  surface_waterflux_qr    
           CASE ( 's_surface' )
              READ ( 13 )  s_surface
@@ -656 +712 @@
           CASE ( 'v_max_ijk' )
              READ ( 13 )  v_max_ijk
+          CASE ( 'ventilation_effect' )
+             READ ( 13 )  ventilation_effect
           CASE ( 'vg' )
              READ ( 13 )  vg