Ignore:
Timestamp:
Apr 7, 2016 7:49:42 AM (5 years ago)
Author:
hoffmann
Message:

changes in LPM and bulk cloud microphysics

File:
1 edited

Legend:

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

    r1818 r1822  
    1919! Current revisions:
    2020! -----------------
    21 !
    22 ! 
     21! Tails removed. icloud_scheme replaced by microphysics_*
     22!
    2323! Former revisions:
    2424! -----------------
     
    275275               density_ratio, dissipation_classes, dt_min_part, dt_prel,       &
    276276               dt_write_particle_data, end_time_prel,                          &
    277                maximum_number_of_tailpoints, maximum_tailpoint_age,            &
    278                minimum_tailpoint_distance, number_of_particle_groups,          &
    279                particle_advection, particle_advection_start,                   &
     277               number_of_particle_groups, particle_advection,                  &
     278               particle_advection_start,                                       &
    280279               particles_per_point, pdx, pdy, pdz,  psb, psl, psn, psr, pss,   &
    281280               pst, radius, radius_classes, random_start_position,             &
    282281               seed_follows_topography,                                        &
    283                total_number_of_particles, use_particle_tails,                  &
    284                use_sgs_for_particles, total_number_of_tails,                   &
     282               total_number_of_particles, use_sgs_for_particles,               &
    285283               vertical_particle_advection, write_particle_statistics
    286284       
     
    17001698                                   slicer_range_limits_dvrp(:,m)
    17011699             ENDIF
    1702           ELSEIF ( mode_dvrp(i)(1:9) == 'particles' )  THEN
    1703              WRITE ( io, 363 )  dvrp_psize
    1704              IF ( particle_dvrpsize /= 'none' )  THEN
    1705                 WRITE ( io, 364 )  'size', TRIM( particle_dvrpsize ), &
    1706                                    dvrpsize_interval
    1707              ENDIF
    1708              IF ( particle_color /= 'none' )  THEN
    1709                 WRITE ( io, 364 )  'color', TRIM( particle_color ), &
    1710                                    color_interval
    1711              ENDIF
    17121700          ENDIF
    17131701          i = i + 1
     
    17821770       WRITE ( io, 415 )
    17831771       WRITE ( io, 416 ) surface_pressure, r_d, rho_surface, cp, l_v
    1784        IF ( icloud_scheme == 0 )  THEN
     1772       IF ( microphysics_seifert )  THEN
    17851773          WRITE ( io, 510 ) 1.0E-6_wp * nc_const
    1786           IF ( precipitation )  WRITE ( io, 511 ) c_sedimentation
     1774          WRITE ( io, 511 ) c_sedimentation
    17871775       ENDIF
    17881776    ENDIF
     
    17961784       WRITE ( io, 432 )
    17971785       IF ( cloud_top_radiation )  WRITE ( io, 132 )
    1798        IF ( icloud_scheme == 1 )  THEN
    1799           IF ( precipitation )  WRITE ( io, 133 )
    1800        ELSEIF ( icloud_scheme == 0 )  THEN
     1786       IF ( microphysics_kessler )  THEN
     1787          WRITE ( io, 133 )
     1788       ELSEIF ( microphysics_seifert )  THEN
    18011789          IF ( drizzle )  WRITE ( io, 506 )
    1802           IF ( precipitation )  THEN
    1803              WRITE ( io, 505 )
    1804              IF ( turbulence )  WRITE ( io, 507 )
    1805              IF ( ventilation_effect )  WRITE ( io, 508 )
    1806              IF ( limiter_sedimentation )  WRITE ( io, 509 )
    1807           ENDIF
     1790          WRITE ( io, 505 )
     1791          IF ( turbulence )  WRITE ( io, 507 )
     1792          IF ( ventilation_effect )  WRITE ( io, 508 )
     1793          IF ( limiter_sedimentation )  WRITE ( io, 509 )
    18081794       ENDIF
    18091795    ELSEIF ( humidity  .AND.  cloud_droplets )  THEN
     
    18721858       IF ( particles_per_point > 1 )  WRITE ( io, 489 )  particles_per_point
    18731859       WRITE ( io, 495 )  total_number_of_particles
    1874        IF ( use_particle_tails  .AND.  maximum_number_of_tailpoints /= 0 )  THEN
    1875           WRITE ( io, 483 )  maximum_number_of_tailpoints
    1876           IF ( minimum_tailpoint_distance /= 0 )  THEN
    1877              WRITE ( io, 484 )  total_number_of_tails,      &
    1878                                 minimum_tailpoint_distance, &
    1879                                 maximum_tailpoint_age
    1880           ENDIF
    1881        ENDIF
    18821860       IF ( dt_write_particle_data /= 9999999.9_wp )  THEN
    18831861          WRITE ( io, 485 )  dt_write_particle_data
     
    22152193362 FORMAT (/'       Slicer plane ',A/ &
    22162194            '       Slicer limits: [',F6.2,',',F6.2,']')
    2217 363 FORMAT (/'       Particles'/ &
    2218             '          particle size:  ',F7.2,' m')
    2219 364 FORMAT ('          particle ',A,' controlled by "',A,'" with interval [', &
    2220                        F6.2,',',F6.2,']')
    22212195365 FORMAT (/'       Groundplate color: (',F4.2,',',F4.2,',',F4.2,') (R,G,B)'/ &
    22222196            '       Superelevation along (x,y,z): (',F4.1,',',F4.1,',',F4.1, &
     
    23772351481 FORMAT ('       Particles have random start positions'/)
    23782352482 FORMAT ('          Particles are advected only horizontally'/)
    2379 483 FORMAT ('       Particles have tails with a maximum of ',I3,' points')
    2380 484 FORMAT ('            Number of tails of the total domain: ',I10/ &
    2381             '            Minimum distance between tailpoints: ',F8.2,' m'/ &
    2382             '            Maximum age of the end of the tail:  ',F8.2,' s')
    23832353485 FORMAT ('       Particle data are written on file every ', F9.1, ' s')
    23842354486 FORMAT ('       Particle statistics are written on file'/)
Note: See TracChangeset for help on using the changeset viewer.