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/timestep.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    156156             v_gtrans = v_gtrans + hom(k,1,2,0)
    157157          ENDDO
    158           u_gtrans = u_gtrans / REAL( nzt - nzb )
    159           v_gtrans = v_gtrans / REAL( nzt - nzb )
     158          u_gtrans = u_gtrans / REAL( nzt - nzb, KIND=wp )
     159          v_gtrans = v_gtrans / REAL( nzt - nzb, KIND=wp )
    160160       ELSE
    161161!
     
    173173          ENDDO
    174174          !$acc end parallel
    175           uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
    176           uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
     175          uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
     176          uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
    177177#if defined( __parallel )
    178178          IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    179179          CALL MPI_ALLREDUCE( uv_gtrans_l, uv_gtrans, 2, MPI_REAL, MPI_SUM, &
    180180                              comm2d, ierr )
    181           u_gtrans = uv_gtrans(1) / REAL( numprocs )
    182           v_gtrans = uv_gtrans(2) / REAL( numprocs )
     181          u_gtrans = uv_gtrans(1) / REAL( numprocs, KIND=wp )
     182          v_gtrans = uv_gtrans(2) / REAL( numprocs, KIND=wp )
    183183#else
    184184          u_gtrans = uv_gtrans_l(1)
Note: See TracChangeset for help on using the changeset viewer.