Changeset 849 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Mar 15, 2012 10:35:09 AM (12 years ago)

Author:

raasch

Message:

Changed:

---

Original routine advec_particles split into several new subroutines and renamed
lpm.
init_particles renamed lpm_init
user_advec_particles renamed user_lpm_advec,
particle_boundary_conds renamed lpm_boundary_conds,
set_particle_attributes renamed lpm_set_attributes,
user_init_particles renamed user_lpm_init,
user_particle_attributes renamed user_lpm_set_attributes
(Makefile, lpm_droplet_collision, lpm_droplet_condensation, init_3d_model, modules, palm, read_var_list, time_integration, write_var_list, deleted: advec_particles, init_particles, particle_boundary_conds, set_particle_attributes, user_advec_particles, user_init_particles, user_particle_attributes, new: lpm, lpm_advec, lpm_boundary_conds, lpm_calc_liquid_water_content, lpm_data_output_particles, lpm_droplet_collision, lpm_drollet_condensation, lpm_exchange_horiz, lpm_extend_particle_array, lpm_extend_tails, lpm_extend_tail_array, lpm_init, lpm_init_sgs_tke, lpm_pack_arrays, lpm_read_restart_file, lpm_release_set, lpm_set_attributes, lpm_sort_arrays, lpm_write_exchange_statistics, lpm_write_restart_file, user_lpm_advec, user_lpm_init, user_lpm_set_attributes

File:

• palm/trunk/SOURCE/lpm_init.f90 (moved) (moved from palm/trunk/SOURCE/init_particles.f90) (15 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -1 @@
- SUBROUTINE init_particles
+ SUBROUTINE lpm_init
 
 !------------------------------------------------------------------------------!
 ! Current revisions:
 ! -----------------
-! 
+! routine renamed: init_particles -> lpm_init
+! de_dx, de_dy, de_dz are allocated here (instead of automatic arrays in
+! advec_particles),
+! sort_particles renamed lpm_sort_arrays, user_init_particles renamed lpm_init
 !
 ! Former revisions:
@@ -65 +68 @@
 ! ------------
 ! This routine initializes a set of particles and their attributes (position,
-! radius, ..). Advection of these particles is carried out by advec_particles,
-! plotting is done in data_output_dvrp.
+! radius, ..) which are used by the Lagrangian particle model (see lpm).
 !------------------------------------------------------------------------------!
 
@@ -83 +85 @@
     IMPLICIT NONE
 
-    CHARACTER (LEN=10) ::  particle_binary_version, version_on_file
-
     INTEGER ::  i, j, n, nn
 #if defined( __parallel )
@@ -124 +124 @@
                                   '&number_of_particle_groups reset to ', &
                                   max_number_of_particle_groups
-       CALL message( 'init_particles', 'PA0213', 0, 1, 0, 6, 0 )
+       CALL message( 'lpm_init', 'PA0213', 0, 1, 0, 6, 0 )
        number_of_particle_groups = max_number_of_particle_groups
     ENDIF
@@ -154 +154 @@
 
 !
+!-- Allocate arrays required for calculating particle SGS velocities
+    IF ( use_sgs_for_particles )  THEN
+       ALLOCATE( de_dx(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
+                 de_dy(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
+                 de_dz(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+    ENDIF
+
+!
 !-- Initialize collision kernels
     IF ( collision_kernel /= 'none' )  CALL init_kernels
@@ -159 +167 @@
 !
 !-- For the first model run of a possible job chain initialize the
-!-- particles, otherwise read the particle data from file.
+!-- particles, otherwise read the particle data from restart file.
     IF ( TRIM( initializing_actions ) == 'read_restart_data'  &
          .AND.  read_particles_from_restartfile )  THEN
 
-!
-!--    Read particle data from previous model run.
-!--    First open the input unit.
-       IF ( myid_char == '' )  THEN
-          OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN'//myid_char, &
-                     FORM='UNFORMATTED' )
-       ELSE
-          OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN/'//myid_char, &
-                     FORM='UNFORMATTED' )
-       ENDIF
-
-!
-!--    First compare the version numbers
-       READ ( 90 )  version_on_file
-       particle_binary_version = '3.0'
-       IF ( TRIM( version_on_file ) /= TRIM( particle_binary_version ) )  THEN
-          message_string = 'version mismatch concerning data from prior ' // &
-                           'run &version on file    = "' //                  &
-                                         TRIM( version_on_file ) //          &
-                           '&version in program = "' //                      &
-                                         TRIM( particle_binary_version ) // '"'
-          CALL message( 'init_particles', 'PA0214', 1, 2, 0, 6, 0 )
-       ENDIF
-
-!
-!--    Read some particle parameters and the size of the particle arrays,
-!--    allocate them and read their contents.
-       READ ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,                  &
-                    maximum_number_of_particles, maximum_number_of_tailpoints, &
-                    maximum_number_of_tails, number_of_initial_particles,      &
-                    number_of_particles, number_of_particle_groups,            &
-                    number_of_tails, particle_groups, time_prel,               &
-                    time_write_particle_data, uniform_particles
-
-       IF ( number_of_initial_particles /= 0 )  THEN
-          ALLOCATE( initial_particles(1:number_of_initial_particles) )
-          READ ( 90 )  initial_particles
-       ENDIF
-
-       ALLOCATE( prt_count(nzb:nzt+1,nysg:nyng,nxlg:nxrg),       &
-                 prt_start_index(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
-                 particle_mask(maximum_number_of_particles),     &
-                 part_1(maximum_number_of_particles),            &
-                 part_2(maximum_number_of_particles) )
-
-       part_1 = particle_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
-                               0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
-                               0.0, 0, 0, 0, 0 )
-
-       part_2 = particle_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
-                               0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
-                               0.0, 0, 0, 0, 0 )
-
-       sort_count = 0
-
-       particles => part_1
-
-       READ ( 90 )  prt_count, prt_start_index
-       READ ( 90 )  particles
-
-       IF ( use_particle_tails )  THEN
-          ALLOCATE( particle_tail_coordinates(maximum_number_of_tailpoints,5, &
-                    maximum_number_of_tails),                                 &
-                    new_tail_id(maximum_number_of_tails),                     &
-                    tail_mask(maximum_number_of_tails) )
-          READ ( 90 )  particle_tail_coordinates
-       ENDIF
-
-       CLOSE ( 90 )
+       CALL lpm_read_restart_file
 
     ELSE
@@ -278 +218 @@
              WRITE( message_string, * ) 'particle group #', i, 'has a', &
                                         'density ratio /= 0 but radius = 0'
-             CALL message( 'init_particles', 'PA0215', 1, 2, 0, 6, 0 )
+             CALL message( 'lpm_init', 'PA0215', 1, 2, 0, 6, 0 )
           ENDIF
           particle_groups(i)%density_ratio = density_ratio(i)
@@ -318 +258 @@
                                       maximum_number_of_particles, ') on PE ', &
                                              myid
-                               CALL message( 'init_particles', 'PA0216', &
-                                                                 2, 2, -1, 6, 1 )
+                               CALL message( 'lpm_init', 'PA0216', 2, 2, -1, 6,&
+                                             1 )
                             ENDIF
                             particles(n)%x             = pos_x
@@ -410 +350 @@
 !
 !--    User modification of initial particles
-       CALL user_init_particles
+       CALL user_lpm_init
 
 !
@@ -456 +396 @@
 !--    Sort particles in the sequence the gridboxes are stored in the memory.
 !--    Only required if cloud droplets are used.
-       IF ( cloud_droplets )  CALL sort_particles
+       IF ( cloud_droplets )  CALL lpm_sort_arrays
 
 !
@@ -575 +515 @@
           WRITE( message_string, * )  'unknown boundary condition ',   &
                                        'bc_par_b = "', TRIM( bc_par_b ), '"'
-          CALL message( 'init_particles', 'PA0217', 1, 2, 0, 6, 0 )
+          CALL message( 'lpm_init', 'PA0217', 1, 2, 0, 6, 0 )
           
     END SELECT
@@ -589 +529 @@
           WRITE( message_string, * ) 'unknown boundary condition ',   &
                                      'bc_par_t = "', TRIM( bc_par_t ), '"'
-          CALL message( 'init_particles', 'PA0218', 1, 2, 0, 6, 0 )
+          CALL message( 'lpm_init', 'PA0218', 1, 2, 0, 6, 0 )
           
     END SELECT
@@ -606 +546 @@
           WRITE( message_string, * ) 'unknown boundary condition ',   &
                                      'bc_par_lr = "', TRIM( bc_par_lr ), '"'
-          CALL message( 'init_particles', 'PA0219', 1, 2, 0, 6, 0 )
+          CALL message( 'lpm_init', 'PA0219', 1, 2, 0, 6, 0 )
           
     END SELECT
@@ -623 +563 @@
           WRITE( message_string, * ) 'unknown boundary condition ',   &
                                      'bc_par_ns = "', TRIM( bc_par_ns ), '"'
-          CALL message( 'init_particles', 'PA0220', 1, 2, 0, 6, 0 )
+          CALL message( 'lpm_init', 'PA0220', 1, 2, 0, 6, 0 )
           
     END SELECT
@@ -630 +570 @@
 8000 FORMAT (I6,1X,F7.2,4X,I6,71X,I6)
 
- END SUBROUTINE init_particles
+ END SUBROUTINE lpm_init