Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (8 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/tridia_solver.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    170170                   IF ( i >= 0  .AND.  i <= nnxh )  THEN
    171171                      ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
    172                                                   REAL( nx+1 ) ) ) / ( dx * dx ) + &
     172                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
    173173                                2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    174                                                   REAL( ny+1 ) ) ) / ( dy * dy )
     174                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
    175175                   ELSE
    176176                      ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
    177                                                   REAL( nx+1 ) ) ) / ( dx * dx ) + &
     177                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
    178178                                2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    179                                                   REAL( ny+1 ) ) ) / ( dy * dy )
     179                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
    180180                   ENDIF
    181181                ELSE
    182182                   IF ( i >= 0  .AND.  i <= nnxh )  THEN
    183183                      ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
    184                                                   REAL( nx+1 ) ) ) / ( dx * dx ) + &
     184                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
    185185                                2.0 * ( 1.0 - COS( ( 2.0 * pi * ( ny+1-j ) ) / &
    186                                                   REAL( ny+1 ) ) ) / ( dy * dy )
     186                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
    187187                   ELSE
    188188                      ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
    189                                                   REAL( nx+1 ) ) ) / ( dx * dx ) + &
     189                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
    190190                                2.0 * ( 1.0 - COS( ( 2.0 * pi * ( ny+1-j ) ) / &
    191                                                   REAL( ny+1 ) ) ) / ( dy * dy )
     191                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
    192192                   ENDIF
    193193                ENDIF
     
    558558             IF ( i >= 0 .AND. i <= nnxh ) THEN
    559559                l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
    560                                          REAL( nx+1 ) ) ) * ddx2 + &
     560                                         REAL( nx+1, KIND=wp ) ) ) * ddx2 + &
    561561                       2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    562                                          REAL( ny+1 ) ) ) * ddy2
     562                                         REAL( ny+1, KIND=wp ) ) ) * ddy2
    563563             ELSE
    564564                l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
    565                                          REAL( nx+1 ) ) ) * ddx2 + &
     565                                         REAL( nx+1, KIND=wp ) ) ) * ddx2 + &
    566566                       2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
    567                                          REAL( ny+1 ) ) ) * ddy2
     567                                         REAL( ny+1, KIND=wp ) ) ) * ddy2
    568568             ENDIF
    569569          ENDDO
Note: See TracChangeset for help on using the changeset viewer.