Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (7 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1329 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    111111              pts_value_l(0:number_of_particle_groups,dopts_num) )
    112112
    113     pts_value_l = 0.0
    114     pts_value_l(:,16) = 9999999.9    ! for calculation of minimum radius
     113    pts_value_l = 0.0_wp
     114    pts_value_l(:,16) = 9999999.9_wp    ! for calculation of minimum radius
    115115
    116116!
     
    135135          pts_value_l(0,11) = pts_value_l(0,11) + particles(n)%rvar3 ! mean sgsw
    136136       ENDIF
    137        IF ( particles(n)%speed_z > 0.0 )  THEN
    138           pts_value_l(0,12) = pts_value_l(0,12) + 1.0  ! # of upward moving prts
     137       IF ( particles(n)%speed_z > 0.0_wp )  THEN
     138          pts_value_l(0,12) = pts_value_l(0,12) + 1.0_wp  ! # of upward moving prts
    139139          pts_value_l(0,13) = pts_value_l(0,13) + &
    140140                                            particles(n)%speed_z ! mean w upw.
     
    170170             pts_value_l(j,11) = pts_value_l(j,11) + particles(n)%rvar3
    171171          ENDIF
    172           IF ( particles(n)%speed_z > 0.0 )  THEN
    173              pts_value_l(j,12) = pts_value_l(j,12) + 1.0
     172          IF ( particles(n)%speed_z > 0.0_wp )  THEN
     173             pts_value_l(j,12) = pts_value_l(j,12) + 1.0_wp
    174174             pts_value_l(j,13) = pts_value_l(j,13) + particles(n)%speed_z
    175175          ELSE
     
    179179          pts_value_l(j,16) = MIN( pts_value(j,16), particles(n)%radius )
    180180          pts_value_l(j,17) = MAX( pts_value(j,17), particles(n)%radius )
    181           pts_value_l(j,18) = pts_value_l(j,18) + 1.0
    182           pts_value_l(j,19) = pts_value_l(j,19) + 1.0
     181          pts_value_l(j,18) = pts_value_l(j,18) + 1.0_wp
     182          pts_value_l(j,19) = pts_value_l(j,19) + 1.0_wp
    183183
    184184       ENDIF
     
    221221    DO  j = 0, inum
    222222
    223        IF ( pts_value(j,1) > 0.0 )  THEN
     223       IF ( pts_value(j,1) > 0.0_wp )  THEN
    224224
    225225          pts_value(j,2:15) = pts_value(j,2:15) / pts_value(j,1)
    226           IF ( pts_value(j,12) > 0.0  .AND.  pts_value(j,12) < 1.0 )  THEN
     226          IF ( pts_value(j,12) > 0.0_wp  .AND.  pts_value(j,12) < 1.0_wp )  THEN
    227227             pts_value(j,13) = pts_value(j,13) / pts_value(j,12)
    228              pts_value(j,14) = pts_value(j,14) / ( 1.0 - pts_value(j,12) )
    229           ELSEIF ( pts_value(j,12) == 0.0 )  THEN
    230              pts_value(j,13) = -1.0
     228             pts_value(j,14) = pts_value(j,14) / ( 1.0_wp - pts_value(j,12) )
     229          ELSEIF ( pts_value(j,12) == 0.0_wp )  THEN
     230             pts_value(j,13) = -1.0_wp
    231231          ELSE
    232              pts_value(j,14) = -1.0
     232             pts_value(j,14) = -1.0_wp
    233233          ENDIF
    234234
     
    323323    DO  j = 0, inum
    324324
    325        IF ( pts_value(j,1) > 0.0 )  THEN
     325       IF ( pts_value(j,1) > 0.0_wp )  THEN
    326326          pts_value(j,20:28) = pts_value(j,20:28) / pts_value(j,1)
    327327       ENDIF
Note: See TracChangeset for help on using the changeset viewer.