Changeset 4195 for palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

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:

• palm/trunk/SOURCE/lagrangian_particle_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/lagrangian_particle_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix for simple_corrector interpolation method in case of ocean runs and
+! output particle advection interpolation method into header
+! 
+! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
 ! 
@@ -622 +626 @@
 
     CHARACTER (LEN=40) ::  output_format       !< netcdf format
- 
+
     INTEGER(iwp) ::  i               !<
     INTEGER(iwp), INTENT(IN) ::  io  !< Unit of the output file
 
- 
+
      IF ( humidity  .AND.  cloud_droplets )  THEN
        WRITE ( io, 433 )
@@ -643 +647 @@
 !
 !--    Particle attributes
-       WRITE ( io, 480 )  particle_advection_start, dt_prel, bc_par_lr, &
+       WRITE ( io, 480 )  particle_advection_start, TRIM(particle_advection_interpolation), &
+                          dt_prel, bc_par_lr, &
                           bc_par_ns, bc_par_b, bc_par_t, particle_maximum_age, &
                           end_time_prel
@@ -708 +713 @@
             '       Particle advection is active (switched on at t = ', F7.1, &
                     ' s)'/ &
+            '       Interpolation of particle velocities is done by using ', A, &
+                    ' method'/ &
             '       Start of new particle generations every  ',F6.1,' s'/ &
             '       Boundary conditions: left/right: ', A, ' north/south: ', A/&
@@ -1222 +1229 @@
 
     TYPE(particle_type),TARGET ::  tmp_particle !< temporary particle used for initialization
+
 
 !
@@ -3297 +3305 @@
 !
 !--          z_direction
-             k_next = MAX( MIN( FLOOR( zp / (zw(kkw+1)-zw(kkw)) ), nzt ), 0)
+             k_next = MAX( MIN( FLOOR( zp / (zw(kkw+1)-zw(kkw)) + offset_ocean_nzt ), nzt ), 0)
              gamma = MAX( MIN( ( zp - zw(k_next) ) /                      &
                                ( zw(k_next+1) - zw(k_next) ), 1.0_wp ), 0.0_wp )