Changeset 4195 for palm/trunk/SOURCE


Ignore:
Timestamp:
Aug 28, 2019 1:44:27 PM (5 years ago)
Author:
schwenkel
Message:

Bugfix for simple_corrector interpolation method in case of ocean runs and introduce particle advection interpolation into header

File:
1 edited

Legend:

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

    r4182 r4195  
    2525! -----------------
    2626! $Id$
     27! Bugfix for simple_corrector interpolation method in case of ocean runs and
     28! output particle advection interpolation method into header
     29!
     30! 4182 2019-08-22 15:20:23Z scharf
    2731! Corrected "Former revisions" section
    2832!
     
    622626
    623627    CHARACTER (LEN=40) ::  output_format       !< netcdf format
    624  
     628
    625629    INTEGER(iwp) ::  i               !<
    626630    INTEGER(iwp), INTENT(IN) ::  io  !< Unit of the output file
    627631
    628  
     632
    629633     IF ( humidity  .AND.  cloud_droplets )  THEN
    630634       WRITE ( io, 433 )
     
    643647!
    644648!--    Particle attributes
    645        WRITE ( io, 480 )  particle_advection_start, dt_prel, bc_par_lr, &
     649       WRITE ( io, 480 )  particle_advection_start, TRIM(particle_advection_interpolation), &
     650                          dt_prel, bc_par_lr, &
    646651                          bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
    647652                          end_time_prel
     
    708713            '       Particle advection is active (switched on at t = ', F7.1, &
    709714                    ' s)'/ &
     715            '       Interpolation of particle velocities is done by using ', A, &
     716                    ' method'/ &
    710717            '       Start of new particle generations every  ',F6.1,' s'/ &
    711718            '       Boundary conditions: left/right: ', A, ' north/south: ', A/&
     
    12221229
    12231230    TYPE(particle_type),TARGET ::  tmp_particle !< temporary particle used for initialization
     1231
    12241232
    12251233!
     
    32973305!
    32983306!--          z_direction
    3299              k_next = MAX( MIN( FLOOR( zp / (zw(kkw+1)-zw(kkw)) ), nzt ), 0)
     3307             k_next = MAX( MIN( FLOOR( zp / (zw(kkw+1)-zw(kkw)) + offset_ocean_nzt ), nzt ), 0)
    33003308             gamma = MAX( MIN( ( zp - zw(k_next) ) /                      &
    33013309                               ( zw(k_next+1) - zw(k_next) ), 1.0_wp ), 0.0_wp )
Note: See TracChangeset for help on using the changeset viewer.