Changeset 1562 for palm/trunk/SOURCE


Ignore:
Timestamp:
Mar 9, 2015 1:10:37 PM (9 years ago)
Author:
heinze
Message:

Bugfix: REAL constants provided with KIND-attribute in call of intrinsic function like MAX and MIN

File:
1 edited

Legend:

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

    r1561 r1562  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Bugfix: REAL constants provided with KIND-attribute in call of
     23! intrinsic function like MAX and MIN
    2324!
    2425! Former revisions:
     
    757758!--          Calculate ratio of upwind gradients. Note, Min/Max is just to
    758759!--          avoid if statements.
    759              rl     = ( MAX( 0.0, u(k,j,i) - u_gtrans ) *                     &
     760             rl     = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) *                  &
    760761                           ABS( ( sk(k,j,i-1) - sk(k,j,i-2)            ) /    &
    761762                                ( sk(k,j,i)   - sk(k,j,i-1) + 1E-20_wp )      &
    762763                              ) +                                             &
    763                         MIN( 0.0, u(k,j,i) - u_gtrans ) *                     &
     764                        MIN( 0.0_wp, u(k,j,i) - u_gtrans ) *                  &
    764765                           ABS( ( sk(k,j,i)   - sk(k,j,i+1)            ) /    &
    765766                                ( sk(k,j,i-1) - sk(k,j,i)   + 1E-20_wp )      &
     
    767768                      ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp )
    768769
    769              rr     = ( MAX( 0.0, u(k,j,i+1) - u_gtrans ) *                   &
     770             rr     = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) *                &
    770771                           ABS( ( sk(k,j,i)   - sk(k,j,i-1)            ) /    &
    771772                                ( sk(k,j,i+1) - sk(k,j,i)   + 1E-20_wp )      &
    772773                              ) +                                             &
    773                         MIN( 0.0, u(k,j,i+1) - u_gtrans ) *                   &
     774                        MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) *                &
    774775                           ABS( ( sk(k,j,i+1) - sk(k,j,i+2)            ) /    &
    775776                                ( sk(k,j,i)   - sk(k,j,i+1) + 1E-20_wp )      &
     
    777778                      ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp )
    778779
    779              rs     = ( MAX( 0.0, v(k,j,i) - v_gtrans ) *                     &
     780             rs     = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) *                  &
    780781                           ABS( ( sk(k,j-1,i) - sk(k,j-2,i)            ) /    &
    781782                                ( sk(k,j,i)   - sk(k,j-1,i) + 1E-20_wp )      &
    782783                              ) +                                             &
    783                         MIN( 0.0, v(k,j,i) - v_gtrans ) *                     &
     784                        MIN( 0.0_wp, v(k,j,i) - v_gtrans ) *                  &
    784785                           ABS( ( sk(k,j,i)   - sk(k,j+1,i)            ) /    &
    785786                                ( sk(k,j-1,i) - sk(k,j,i)   + 1E-20_wp )      &
     
    787788                      ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp )
    788789
    789              rn     = ( MAX( 0.0, v(k,j+1,i) - v_gtrans ) *                   &
     790             rn     = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) *                &
    790791                           ABS( ( sk(k,j,i)   - sk(k,j-1,i)            ) /    &
    791792                                ( sk(k,j+1,i) - sk(k,j,i)   + 1E-20_wp )      &
    792793                              ) +                                             &
    793                        MIN( 0.0, v(k,j+1,i) - v_gtrans ) *                    &
     794                       MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) *                 &
    794795                           ABS( ( sk(k,j+1,i) - sk(k,j+2,i)            ) /    &
    795796                                ( sk(k,j,i)   - sk(k,j+1,i) + 1E-20_wp )      &
     
    803804             k_mmm  = k - 3 * ibit8
    804805
    805              rd     = ( MAX( 0.0, w(k-1,j,i) ) *                              &
     806             rd     = ( MAX( 0.0_wp, w(k-1,j,i) ) *                           &
    806807                           ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i)           ) /  &
    807808                                ( sk(k-1,j,i)  - sk(k_mm,j,i) + 1E-20_wp )    &
    808809                              ) +                                             &
    809                         MIN( 0.0, w(k-1,j,i) ) *                              &
     810                        MIN( 0.0_wp, w(k-1,j,i) ) *                           &
    810811                           ABS( ( sk(k-1,j,i) - sk(k,j,i)            ) /      &
    811812                                ( sk(k_mm,j,i) - sk(k-1,j,i)   + 1E-20_wp )   &
     
    813814                      ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp )
    814815 
    815              rt     = ( MAX( 0.0, w(k,j,i) ) *                                &
     816             rt     = ( MAX( 0.0_wp, w(k,j,i) ) *                             &
    816817                           ABS( ( sk(k,j,i)   - sk(k-1,j,i)            ) /    &
    817818                                ( sk(k+1,j,i) - sk(k,j,i)   + 1E-20_wp )      &
    818819                              ) +                                             &
    819                         MIN( 0.0, w(k,j,i) ) *                                &
     820                        MIN( 0.0_wp, w(k,j,i) ) *                             &
    820821                           ABS( ( sk(k+1,j,i) - sk(k_pp,j,i)           ) /    &
    821822                                ( sk(k,j,i)   - sk(k+1,j,i) + 1E-20_wp )      &
     
    986987!--          Calculate ratio of upwind gradients. Note, Min/Max is just to
    987988!--          avoid if statements.
    988              rl     = ( MAX( 0.0, u(k,j,i) - u_gtrans ) *                     &
     989             rl     = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) *                  &
    989990                           ABS( ( sk(k,j,i-1) - sk(k,j,i-2)            ) /    &
    990991                                ( sk(k,j,i)   - sk(k,j,i-1) + 1E-20_wp )      &
    991992                              ) +                                             &
    992                         MIN( 0.0, u(k,j,i) - u_gtrans ) *                     &
     993                        MIN( 0.0_wp, u(k,j,i) - u_gtrans ) *                  &
    993994                           ABS( ( sk(k,j,i)   - sk(k,j,i+1)            ) /    &
    994995                                ( sk(k,j,i-1) - sk(k,j,i)   + 1E-20_wp )      &
     
    996997                      ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp )
    997998
    998              rr     = ( MAX( 0.0, u(k,j,i+1) - u_gtrans ) *                   &
     999             rr     = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) *                &
    9991000                           ABS( ( sk(k,j,i)   - sk(k,j,i-1)            ) /    &
    10001001                                ( sk(k,j,i+1) - sk(k,j,i)   + 1E-20_wp )      &
    10011002                              ) +                                             &
    1002                         MIN( 0.0, u(k,j,i+1) - u_gtrans ) *                   &
     1003                        MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) *                &
    10031004                           ABS( ( sk(k,j,i+1) - sk(k,j,i+2)            ) /    &
    10041005                                ( sk(k,j,i)   - sk(k,j,i+1) + 1E-20_wp )      &
     
    10061007                      ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp )
    10071008
    1008              rs     = ( MAX( 0.0, v(k,j,i) - v_gtrans ) *                     &
     1009             rs     = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) *                  &
    10091010                           ABS( ( sk(k,j-1,i) - sk(k,j-2,i)            ) /    &
    10101011                                ( sk(k,j,i)   - sk(k,j-1,i) + 1E-20_wp )      &
    10111012                              ) +                                             &
    1012                         MIN( 0.0, v(k,j,i) - v_gtrans ) *                     &
     1013                        MIN( 0.0_wp, v(k,j,i) - v_gtrans ) *                  &
    10131014                           ABS( ( sk(k,j,i)   - sk(k,j+1,i)            ) /    &
    10141015                                ( sk(k,j-1,i) - sk(k,j,i)   + 1E-20_wp )      &
     
    10161017                      ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp )
    10171018
    1018              rn     = ( MAX( 0.0, v(k,j+1,i) - v_gtrans ) *                   &
     1019             rn     = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) *                &
    10191020                           ABS( ( sk(k,j,i)   - sk(k,j-1,i)            ) /    &
    10201021                                ( sk(k,j+1,i) - sk(k,j,i)   + 1E-20_wp )      &
    10211022                              ) +                                             &
    1022                        MIN( 0.0, v(k,j+1,i) - v_gtrans ) *                    &
     1023                       MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) *                 &
    10231024                           ABS( ( sk(k,j+1,i) - sk(k,j+2,i)            ) /    &
    10241025                                ( sk(k,j,i)   - sk(k,j+1,i) + 1E-20_wp )      &
     
    10321033             k_mmm  = k - 3 * ibit8
    10331034
    1034              rd     = ( MAX( 0.0, w(k-1,j,i) ) *                              &
     1035             rd     = ( MAX( 0.0_wp, w(k-1,j,i) ) *                           &
    10351036                           ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i)           ) /  &
    10361037                                ( sk(k-1,j,i)  - sk(k_mm,j,i) + 1E-20_wp )    &
    10371038                              ) +                                             &
    1038                         MIN( 0.0, w(k-1,j,i) ) *                              &
     1039                        MIN( 0.0_wp, w(k-1,j,i) ) *                           &
    10391040                           ABS( ( sk(k-1,j,i) - sk(k,j,i)            ) /      &
    10401041                                ( sk(k_mm,j,i) - sk(k-1,j,i)   + 1E-20_wp )   &
     
    10421043                      ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp )
    10431044 
    1044              rt     = ( MAX( 0.0, w(k,j,i) ) *                                &
     1045             rt     = ( MAX( 0.0_wp, w(k,j,i) ) *                             &
    10451046                           ABS( ( sk(k,j,i)   - sk(k-1,j,i)            ) /    &
    10461047                                ( sk(k+1,j,i) - sk(k,j,i)   + 1E-20_wp )      &
    10471048                              ) +                                             &
    1048                         MIN( 0.0, w(k,j,i) ) *                                &
     1049                        MIN( 0.0_wp, w(k,j,i) ) *                             &
    10491050                           ABS( ( sk(k+1,j,i) - sk(k_pp,j,i)           ) /    &
    10501051                                ( sk(k,j,i)   - sk(k+1,j,i) + 1E-20_wp )      &
     
    28932894!--                Calculate ratio of upwind gradients. Note, Min/Max is just
    28942895!--                to avoid if statements.
    2895                    rl     = ( MAX( 0.0, u(k,j,i) - u_gtrans ) *               &
     2896                   rl     = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    28962897                               ABS( ( sk(k,j,i-1) - sk(k,j,i-2)            ) /&
    28972898                                    ( sk(k,j,i)   - sk(k,j,i-1) + 1E-20_wp )  &
    28982899                                  ) +                                         &
    2899                               MIN( 0.0, u(k,j,i) - u_gtrans ) *               &
     2900                              MIN( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    29002901                               ABS( ( sk(k,j,i)   - sk(k,j,i+1)            ) /&
    29012902                                    ( sk(k,j,i-1) - sk(k,j,i)   + 1E-20_wp )  &
     
    29032904                            ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp )
    29042905
    2905                    rr     = ( MAX( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     2906                   rr     = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    29062907                               ABS( ( sk(k,j,i)   - sk(k,j,i-1)            ) /&
    29072908                                    ( sk(k,j,i+1) - sk(k,j,i)   + 1E-20_wp )  &
    29082909                                  ) +                                         &
    2909                               MIN( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     2910                              MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    29102911                               ABS( ( sk(k,j,i+1) - sk(k,j,i+2)            ) /&
    29112912                                    ( sk(k,j,i)   - sk(k,j,i+1) + 1E-20_wp )  &
     
    29132914                            ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp )
    29142915
    2915                    rs     = ( MAX( 0.0, v(k,j,i) - v_gtrans ) *               &
     2916                   rs     = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    29162917                               ABS( ( sk(k,j-1,i) - sk(k,j-2,i)            ) /&
    29172918                                    ( sk(k,j,i)   - sk(k,j-1,i) + 1E-20_wp )  &
    29182919                                  ) +                                         &
    2919                               MIN( 0.0, v(k,j,i) - v_gtrans ) *               &
     2920                              MIN( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    29202921                               ABS( ( sk(k,j,i)   - sk(k,j+1,i)            ) /&
    29212922                                    ( sk(k,j-1,i) - sk(k,j,i)   + 1E-20_wp )  &
     
    29232924                            ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp )
    29242925
    2925                    rn     = ( MAX( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     2926                   rn     = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    29262927                               ABS( ( sk(k,j,i)   - sk(k,j-1,i)            ) /&
    29272928                                    ( sk(k,j+1,i) - sk(k,j,i)   + 1E-20_wp )  &
    29282929                                  ) +                                         &
    2929                               MIN( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     2930                              MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    29302931                               ABS( ( sk(k,j+1,i) - sk(k,j+2,i)            ) /&
    29312932                                    ( sk(k,j,i)   - sk(k,j+1,i) + 1E-20_wp )  &
     
    29392940                   k_mmm  = k - 3 * ibit8
    29402941
    2941                    rd     = ( MAX( 0.0, w(k-1,j,i) ) *                        &
     2942                   rd     = ( MAX( 0.0_wp, w(k-1,j,i) ) *                     &
    29422943                            ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i)           ) / &
    29432944                                 ( sk(k-1,j,i)  - sk(k_mm,j,i) + 1E-20_wp )   &
    29442945                               ) +                                            &
    2945                               MIN( 0.0, w(k-1,j,i) ) *                        &
     2946                              MIN( 0.0_wp, w(k-1,j,i) ) *                     &
    29462947                            ABS( ( sk(k-1,j,i) - sk(k,j,i)            ) /     &
    29472948                                 ( sk(k_mm,j,i) - sk(k-1,j,i)   + 1E-20_wp )  &
     
    29492950                            ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp )
    29502951 
    2951                    rt     = ( MAX( 0.0, w(k,j,i) ) *                          &
     2952                   rt     = ( MAX( 0.0_wp, w(k,j,i) ) *                       &
    29522953                            ABS( ( sk(k,j,i)   - sk(k-1,j,i)            ) /   &
    29532954                                 ( sk(k+1,j,i) - sk(k,j,i)   + 1E-20_wp )     &
    29542955                               ) +                                            &
    2955                               MIN( 0.0, w(k,j,i) ) *                          &
     2956                              MIN( 0.0_wp, w(k,j,i) ) *                       &
    29562957                            ABS( ( sk(k+1,j,i) - sk(k_pp,j,i)           ) /   &
    29572958                                 ( sk(k,j,i)   - sk(k+1,j,i) + 1E-20_wp )     &
     
    31183119!--                Calculate ratio of upwind gradients. Note, Min/Max is just
    31193120!--                to avoid if statements.
    3120                    rl     = ( MAX( 0.0, u(k,j,i) - u_gtrans ) *               &
     3121                   rl     = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    31213122                               ABS( ( sk(k,j,i-1) - sk(k,j,i-2)            ) /&
    31223123                                    ( sk(k,j,i)   - sk(k,j,i-1) + 1E-20_wp )  &
    31233124                                  ) +                                         &
    3124                               MIN( 0.0, u(k,j,i) - u_gtrans ) *               &
     3125                              MIN( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    31253126                               ABS( ( sk(k,j,i)   - sk(k,j,i+1)            ) /&
    31263127                                    ( sk(k,j,i-1) - sk(k,j,i)   + 1E-20_wp )  &
     
    31283129                            ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp )
    31293130
    3130                    rr     = ( MAX( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     3131                   rr     = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    31313132                               ABS( ( sk(k,j,i)   - sk(k,j,i-1)            ) /&
    31323133                                    ( sk(k,j,i+1) - sk(k,j,i)   + 1E-20_wp )  &
    31333134                                  ) +                                         &
    3134                               MIN( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     3135                              MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    31353136                               ABS( ( sk(k,j,i+1) - sk(k,j,i+2)            ) /&
    31363137                                    ( sk(k,j,i)   - sk(k,j,i+1) + 1E-20_wp )  &
     
    31383139                            ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp )
    31393140
    3140                    rs     = ( MAX( 0.0, v(k,j,i) - v_gtrans ) *               &
     3141                   rs     = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    31413142                               ABS( ( sk(k,j-1,i) - sk(k,j-2,i)            ) /&
    31423143                                    ( sk(k,j,i)   - sk(k,j-1,i) + 1E-20_wp )  &
    31433144                                  ) +                                         &
    3144                               MIN( 0.0, v(k,j,i) - v_gtrans ) *               &
     3145                              MIN( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    31453146                               ABS( ( sk(k,j,i)   - sk(k,j+1,i)            ) /&
    31463147                                    ( sk(k,j-1,i) - sk(k,j,i)   + 1E-20_wp )  &
     
    31483149                            ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp )
    31493150
    3150                    rn     = ( MAX( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     3151                   rn     = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    31513152                               ABS( ( sk(k,j,i)   - sk(k,j-1,i)            ) /&
    31523153                                    ( sk(k,j+1,i) - sk(k,j,i)   + 1E-20_wp )  &
    31533154                                  ) +                                         &
    3154                               MIN( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     3155                              MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    31553156                               ABS( ( sk(k,j+1,i) - sk(k,j+2,i)            ) /&
    31563157                                    ( sk(k,j,i)   - sk(k,j+1,i) + 1E-20_wp )  &
     
    31643165                   k_mmm  = k - 3 * ibit8
    31653166
    3166                    rd     = ( MAX( 0.0, w(k-1,j,i) ) *                        &
     3167                   rd     = ( MAX( 0.0_wp, w(k-1,j,i) ) *                     &
    31673168                            ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i)           ) / &
    31683169                                 ( sk(k-1,j,i)  - sk(k_mm,j,i) + 1E-20_wp )   &
    31693170                               ) +                                            &
    3170                               MIN( 0.0, w(k-1,j,i) ) *                        &
     3171                              MIN( 0.0_wp, w(k-1,j,i) ) *                     &
    31713172                            ABS( ( sk(k-1,j,i) - sk(k,j,i)            ) /     &
    31723173                                 ( sk(k_mm,j,i) - sk(k-1,j,i)   + 1E-20_wp )  &
     
    31743175                            ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp )
    31753176 
    3176                    rt     = ( MAX( 0.0, w(k,j,i) ) *                          &
     3177                   rt     = ( MAX( 0.0_wp, w(k,j,i) ) *                       &
    31773178                            ABS( ( sk(k,j,i)   - sk(k-1,j,i)            ) /   &
    31783179                                 ( sk(k+1,j,i) - sk(k,j,i)   + 1E-20_wp )     &
    31793180                               ) +                                            &
    3180                               MIN( 0.0, w(k,j,i) ) *                          &
     3181                              MIN( 0.0_wp, w(k,j,i) ) *                       &
    31813182                            ABS( ( sk(k+1,j,i) - sk(k_pp,j,i)           ) /   &
    31823183                                 ( sk(k,j,i)   - sk(k+1,j,i) + 1E-20_wp )     &
     
    36313632!--                Calculate ratio of upwind gradients. Note, Min/Max is just
    36323633!--                to avoid if statements.
    3633                    rl     = ( MAX( 0.0, u(k,j,i) - u_gtrans ) *               &
     3634                   rl     = ( MAX( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    36343635                               ABS( ( sk(k,j,i-1) - sk(k,j,i-2)            ) /&
    36353636                                    ( sk(k,j,i)   - sk(k,j,i-1) + 1E-20_wp )  &
    36363637                                  ) +                                         &
    3637                               MIN( 0.0, u(k,j,i) - u_gtrans ) *               &
     3638                              MIN( 0.0_wp, u(k,j,i) - u_gtrans ) *            &
    36383639                               ABS( ( sk(k,j,i)   - sk(k,j,i+1)            ) /&
    36393640                                    ( sk(k,j,i-1) - sk(k,j,i)   + 1E-20_wp )  &
     
    36413642                            ) / ABS( u(k,j,i) - u_gtrans + 1E-20_wp )
    36423643
    3643                    rr     = ( MAX( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     3644                   rr     = ( MAX( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    36443645                               ABS( ( sk(k,j,i)   - sk(k,j,i-1)            ) /&
    36453646                                    ( sk(k,j,i+1) - sk(k,j,i)   + 1E-20_wp )  &
    36463647                                  ) +                                         &
    3647                               MIN( 0.0, u(k,j,i+1) - u_gtrans ) *             &
     3648                              MIN( 0.0_wp, u(k,j,i+1) - u_gtrans ) *          &
    36483649                               ABS( ( sk(k,j,i+1) - sk(k,j,i+2)            ) /&
    36493650                                    ( sk(k,j,i)   - sk(k,j,i+1) + 1E-20_wp )  &
     
    36513652                            ) / ABS( u(k,j,i+1) - u_gtrans + 1E-20_wp )
    36523653
    3653                    rs     = ( MAX( 0.0, v(k,j,i) - v_gtrans ) *               &
     3654                   rs     = ( MAX( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    36543655                               ABS( ( sk(k,j-1,i) - sk(k,j-2,i)            ) /&
    36553656                                    ( sk(k,j,i)   - sk(k,j-1,i) + 1E-20_wp )  &
    36563657                                  ) +                                         &
    3657                               MIN( 0.0, v(k,j,i) - v_gtrans ) *               &
     3658                              MIN( 0.0_wp, v(k,j,i) - v_gtrans ) *            &
    36583659                               ABS( ( sk(k,j,i)   - sk(k,j+1,i)            ) /&
    36593660                                    ( sk(k,j-1,i) - sk(k,j,i)   + 1E-20_wp )  &
     
    36613662                            ) / ABS( v(k,j,i) - v_gtrans + 1E-20_wp )
    36623663
    3663                    rn     = ( MAX( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     3664                   rn     = ( MAX( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    36643665                               ABS( ( sk(k,j,i)   - sk(k,j-1,i)            ) /&
    36653666                                    ( sk(k,j+1,i) - sk(k,j,i)   + 1E-20_wp )  &
    36663667                                  ) +                                         &
    3667                               MIN( 0.0, v(k,j+1,i) - v_gtrans ) *             &
     3668                              MIN( 0.0_wp, v(k,j+1,i) - v_gtrans ) *          &
    36683669                               ABS( ( sk(k,j+1,i) - sk(k,j+2,i)            ) /&
    36693670                                    ( sk(k,j,i)   - sk(k,j+1,i) + 1E-20_wp )  &
     
    36773678                   k_mmm  = k - 3 * ibit8
    36783679
    3679                    rd     = ( MAX( 0.0, w(k-1,j,i) ) *                        &
     3680                   rd     = ( MAX( 0.0_wp, w(k-1,j,i) ) *                     &
    36803681                            ABS( ( sk(k_mm,j,i) - sk(k_mmm,j,i)           ) / &
    36813682                                 ( sk(k-1,j,i)  - sk(k_mm,j,i) + 1E-20_wp )   &
    36823683                               ) +                                            &
    3683                               MIN( 0.0, w(k-1,j,i) ) *                        &
     3684                              MIN( 0.0_wp, w(k-1,j,i) ) *                     &
    36843685                            ABS( ( sk(k-1,j,i) - sk(k,j,i)            ) /     &
    36853686                                 ( sk(k_mm,j,i) - sk(k-1,j,i)   + 1E-20_wp )  &
     
    36873688                            ) * ibit8 / ABS( w(k-1,j,i) + 1E-20_wp )
    36883689 
    3689                    rt     = ( MAX( 0.0, w(k,j,i) ) *                          &
     3690                   rt     = ( MAX( 0.0_wp, w(k,j,i) ) *                       &
    36903691                            ABS( ( sk(k,j,i)   - sk(k-1,j,i)            ) /   &
    36913692                                 ( sk(k+1,j,i) - sk(k,j,i)   + 1E-20_wp )     &
    36923693                               ) +                                            &
    3693                               MIN( 0.0, w(k,j,i) ) *                          &
     3694                              MIN( 0.0_wp, w(k,j,i) ) *                       &
    36943695                            ABS( ( sk(k+1,j,i) - sk(k_pp,j,i)           ) /   &
    36953696                                 ( sk(k,j,i)   - sk(k+1,j,i) + 1E-20_wp )     &
Note: See TracChangeset for help on using the changeset viewer.