Ignore:
Timestamp:
Jun 8, 2017 2:59:01 PM (4 years ago)
Author:
schwenkel
Message:

Implemented splitting and merging algorithm

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/lpm_init.f90

    r2233 r2263  
    2525! -----------------
    2626! $Id$
     27! Implemented splitting and merging algorithm
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
     
    175178                block_offset, block_offset_def, collision_kernel,              &
    176179                curvature_solution_effects,                                    &
    177                 density_ratio, grid_particles,                                 &
     180                density_ratio, grid_particles, isf,i_splitting_mode,           &
    178181                initial_weighting_factor, ibc_par_b, ibc_par_lr, ibc_par_ns,   &
    179182                ibc_par_t, iran_part, log_z_z0,                                &
     
    190193                radius, random_start_position, read_particles_from_restartfile,&
    191194                seed_follows_topography, sgs_wf_part, sort_count,              &
    192                 total_number_of_particles,                                     &
     195                splitting_function, splitting_mode, total_number_of_particles,  &
    193196                use_sgs_for_particles,                                         &
    194197                write_particle_statistics, uniform_particles, zero_particle,   &
     
    473476         
    474477    END SELECT
     478    SELECT CASE ( splitting_mode )
     479   
     480       CASE ( 'const' )
     481          i_splitting_mode = 1
     482
     483       CASE ( 'cl_av' )
     484          i_splitting_mode = 2
     485
     486       CASE ( 'gb_av' )
     487          i_splitting_mode = 3
     488         
     489       CASE DEFAULT
     490          WRITE( message_string, * )  'unknown splitting condition ',   &
     491                                       'splitting_mode = "', TRIM( splitting_mode ), '"'
     492          CALL message( 'lpm_init', 'PA0146', 1, 2, 0, 6, 0 )
     493         
     494    END SELECT
     495    SELECT CASE ( splitting_function )
     496   
     497       CASE ( 'gamma' )
     498          isf = 1
     499
     500       CASE ( 'log' )
     501          isf = 2
     502
     503       CASE ( 'exp' )
     504          isf = 3
     505         
     506       CASE DEFAULT
     507          WRITE( message_string, * )  'unknown splitting function ',   &
     508                                       'splitting_function = "', TRIM( splitting_function ), '"'
     509          CALL message( 'lpm_init', 'PA0147', 1, 2, 0, 6, 0 )
     510         
     511    END SELECT
     512   
    475513
    476514!
Note: See TracChangeset for help on using the changeset viewer.