Ignore:
Timestamp:
Feb 22, 2012 11:21:19 AM (12 years ago)
Author:
raasch
Message:

Bugfix: array diss can be used only in case of Wang kernel

File:
1 edited

Legend:

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

    r832 r835  
    44! Current revisions:
    55! ------------------
    6 !
     6! Bugfix: array diss can be used only in case of Wang kernel
    77!
    88! Former revisions:
     
    180180                distance, drdt, drdt_ini, drdt_m, dt_ros, dt_ros_last, &
    181181                dt_ros_next, dt_ros_sum, dt_ros_sum_ini, d2rdt2, d2rdtdr, &
    182                 dt_gap, dt_particle, dt_particle_m, d_sum, e_a, e_int, &
     182                dt_gap, dt_particle, dt_particle_m, d_sum, epsilon, e_a, e_int,&
    183183                e_int_l, e_int_u, e_mean_int, e_s, err_ros, errmax, exp_arg, &
    184184                exp_term, fs_int, gg, g1, g2, g3, g4, integral, &
     
    639639!
    640640!--                   Determine dissipation class index of this gridbox
    641                       eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
    642                                     dissipation_classes ) + 1
    643                       IF ( hall_kernel .OR. diss(k,j,i) * 1.0E4 < 0.001 )  THEN
     641                      IF ( wang_kernel )  THEN
     642                         eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
     643                                       dissipation_classes ) + 1
     644                         epsilon = diss(k,j,i)
     645                      ELSE
     646                         epsilon = 0.0
     647                      ENDIF
     648                      IF ( hall_kernel  .OR.  epsilon * 1.0E4 < 0.001 )  THEN
    644649                         eclass = 0   ! Hall kernel is used
    645650                      ELSE
Note: See TracChangeset for help on using the changeset viewer.