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

    r1323 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    8989!-- Default: eddy radius rc, eddy strength z,
    9090!--          position of eddy centre: ic, jc, kc1, kc2
    91     rc  =  4.0 * dx
     91    rc  =  4.0_wp * dx
    9292    ic  =  ( nx+1 ) / 2
    9393    jc  =  ic
     
    111111    DO  i = nxl, nxr
    112112       DO  j = nys, nyn
    113           x = ( i - ic - 0.5 ) * dx
    114           y = ( j - jc ) * dy
     113          x = ( i - ic - 0.5_wp ) * dx
     114          y = ( j - jc          ) * dy
    115115          radius = SQRT( x**2 + y**2 )
    116           IF ( radius <= 2.0 * rc )  THEN
    117              betrag = radius / ( 2.0 * rc ) * 0.08_wp
    118           ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    119              betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    120           ELSE
    121              betrag = 0.0
    122           ENDIF
    123 
    124           IF ( x == 0.0 )  THEN
    125              IF ( y > 0.0 )  THEN
     116          IF ( radius <= 2.0_wp * rc )  THEN
     117             betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
     118          ELSEIF ( radius > 2.0_wp * rc  .AND.  radius < 8.0_wp * rc )  THEN
     119             betrag = 0.08_wp * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
     120          ELSE
     121             betrag = 0.0_wp
     122          ENDIF
     123
     124          IF ( x == 0.0_wp )  THEN
     125             IF ( y > 0.0_wp )  THEN
    126126                alpha = pi / 2.0_wp
    127              ELSEIF ( y < 0.0 )  THEN
    128                 alpha = 3.0 * pi / 2.0_wp
    129              ENDIF
    130           ELSE
    131              IF ( x < 0.0 )  THEN
     127             ELSEIF ( y < 0.0_wp )  THEN
     128                alpha = 3.0_wp * pi / 2.0_wp
     129             ENDIF
     130          ELSE
     131             IF ( x < 0.0_wp )  THEN
    132132                alpha = ATAN( y / x ) + pi
    133133             ELSE
    134                 IF ( y < 0.0 )  THEN
     134                IF ( y < 0.0_wp )  THEN
    135135                   alpha = ATAN( y / x ) + 2.0_wp * pi
    136136                ELSE
     
    152152    DO  i = nxl, nxr
    153153       DO  j = nys, nyn
    154           x = ( i - ic ) * dx
    155           y = ( j - jc - 0.5) * dy
     154          x = ( i - ic          ) * dx
     155          y = ( j - jc - 0.5_wp ) * dy
    156156          radius = SQRT( x**2 + y**2 )
    157           IF ( radius <= 2.0 * rc )  THEN
     157          IF ( radius <= 2.0_wp * rc )  THEN
    158158             betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
    159           ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    160              betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    161           ELSE
    162              betrag = 0.0
    163           ENDIF
    164 
    165           IF ( x == 0.0 )  THEN
    166              IF ( y > 0.0 )  THEN
     159          ELSEIF ( radius > 2.0_wp * rc  .AND.  radius < 8.0_wp * rc )  THEN
     160             betrag = 0.08_wp * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
     161          ELSE
     162             betrag = 0.0_wp
     163          ENDIF
     164
     165          IF ( x == 0.0_wp )  THEN
     166             IF ( y > 0.0_wp )  THEN
    167167                alpha = pi / 2.0_wp
    168              ELSEIF ( y < 0.0 )  THEN
    169                 alpha = 3.0 * pi / 2.0_wp
    170              ENDIF
    171           ELSE
    172              IF ( x < 0.0 )  THEN
     168             ELSEIF ( y < 0.0_wp )  THEN
     169                alpha = 3.0_wp * pi / 2.0_wp
     170             ENDIF
     171          ELSE
     172             IF ( x < 0.0_wp )  THEN
    173173                alpha = ATAN( y / x ) + pi
    174174             ELSE
    175                 IF ( y < 0.0 )  THEN
     175                IF ( y < 0.0_wp )  THEN
    176176                   alpha = ATAN( y / x ) + 2.0_wp * pi
    177177                ELSE
Note: See TracChangeset for help on using the changeset viewer.