Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (7 years ago)
Author:
raasch
Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:
1 edited

Legend:

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

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    9898          DO  j = dist_nys(dist_range), dist_nyn(dist_range)
    9999             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
    100                 randomnumber = 3.0 * disturbance_amplitude *                   &
    101                                ( random_function( iran ) - 0.5 )
     100                randomnumber = 3.0_wp * disturbance_amplitude *                &
     101                               ( random_function( iran ) - 0.5_wp )
    102102                IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.         &
    103103                     nyn >= j )                                                &
     
    113113             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
    114114#if defined( __nec )
    115                 randomnumber = 3.0 * disturbance_amplitude *                   &
    116                                ( RANDOM( 0 ) - 0.5 )
     115                randomnumber = 3.0_wp * disturbance_amplitude *                &
     116                               ( RANDOM( 0 ) - 0.5_wp )
    117117#else
    118118                CALL RANDOM_NUMBER( randomnumber )
    119                 randomnumber = 3.0 * disturbance_amplitude *                   &
    120                                 ( randomnumber - 0.5 )
     119                randomnumber = 3.0_wp * disturbance_amplitude *                &
     120                                ( randomnumber - 0.5_wp )
    121121#endif
    122122                IF ( nxl <= i .AND. nxr >= i .AND. nys <= j .AND. nyn >= j )   &
     
    145145             dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1)                  &
    146146                            + dist1(k,j+1,i) + dist1(k+1,j,i)                  &
    147                             ) / 12.0
     147                            ) / 12.0_wp
    148148          ENDDO
    149149          DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
    150150              dist2(k,j,i) = dist2(k,j,i) + ( dist1(k,j-1,i) + dist1(k-1,j,i)  &
    151                             + 6.0 * dist1(k,j,i)                               &
    152                             ) / 12.0
     151                            + 6.0_wp * dist1(k,j,i)                            &
     152                            ) / 12.0_wp
    153153          ENDDO
    154154        ENDDO
     
    165165             dist1(k,j,i) = ( dist2(k,j,i-1) + dist2(k,j,i+1) + dist2(k,j-1,i) &
    166166                            + dist2(k,j+1,i) + dist2(k+1,j,i) + dist2(k-1,j,i) &
    167                             + 6.0 * dist2(k,j,i)                               &
    168                             ) / 12.0
     167                            + 6.0_wp * dist2(k,j,i)                            &
     168                            ) / 12.0_wp
    169169          ENDDO
    170170       ENDDO
Note: See TracChangeset for help on using the changeset viewer.