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

    r1321 r1353  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    119119
    120120
    121        ptnudge = 0.0; qnudge = 0.0; tnudge = 0.0; unudge = 0.0
    122        vnudge = 0.0; wnudge = 0.0; timenudge = 0.0
     121       ptnudge = 0.0_wp; qnudge = 0.0_wp; tnudge = 0.0_wp; unudge = 0.0_wp
     122       vnudge = 0.0_wp; wnudge = 0.0_wp; timenudge = 0.0_wp
    123123
    124124       t = 0
     
    193193             fac = ( highheight - zu(k) ) / ( highheight - lowheight )
    194194
    195              tnudge(k,t)  = fac * lowtnudge + ( 1.0 - fac ) * hightnudge
    196              unudge(k,t)  = fac * lowunudge + ( 1.0 - fac ) * highunudge
    197              vnudge(k,t)  = fac * lowvnudge + ( 1.0 - fac ) * highvnudge
    198              wnudge(k,t)  = fac * lowwnudge + ( 1.0 - fac ) * highwnudge
    199              ptnudge(k,t) = fac * lowptnudge + ( 1.0 - fac ) * highptnudge
    200              qnudge(k,t)  = fac * lowqnudge + ( 1.0 - fac ) * highqnudge
     195             tnudge(k,t)  = fac * lowtnudge  + ( 1.0_wp - fac ) * hightnudge
     196             unudge(k,t)  = fac * lowunudge  + ( 1.0_wp - fac ) * highunudge
     197             vnudge(k,t)  = fac * lowvnudge  + ( 1.0_wp - fac ) * highvnudge
     198             wnudge(k,t)  = fac * lowwnudge  + ( 1.0_wp - fac ) * highwnudge
     199             ptnudge(k,t) = fac * lowptnudge + ( 1.0_wp - fac ) * highptnudge
     200             qnudge(k,t)  = fac * lowqnudge  + ( 1.0_wp - fac ) * highqnudge
    201201          ENDDO
    202202
     
    208208!--    Prevent nudging if nudging profiles exhibt too small values
    209209!--    not used so far
    210        lptnudge  = ANY( ABS( ptnudge ) > 1e-8 )
    211        lqnudge   = ANY( ABS( qnudge ) > 1e-8 )
    212        lunudge   = ANY( ABS( unudge ) > 1e-8 )
    213        lvnudge   = ANY( ABS( vnudge ) > 1e-8 )
    214        lwnudge   = ANY( ABS( wnudge ) > 1e-8 )
     210       lptnudge  = ANY( ABS( ptnudge ) > 1.0e-8_wp )
     211       lqnudge   = ANY( ABS( qnudge )  > 1.0e-8_wp )
     212       lunudge   = ANY( ABS( unudge )  > 1.0e-8_wp )
     213       lvnudge   = ANY( ABS( vnudge )  > 1.0e-8_wp )
     214       lwnudge   = ANY( ABS( wnudge )  > 1.0e-8_wp )
    215215
    216216    END SUBROUTINE init_nudge
Note: See TracChangeset for help on using the changeset viewer.