Changeset 1337 for palm


Ignore:
Timestamp:
Mar 25, 2014 3:11:48 PM (11 years ago)
Author:
heinze
Message:

Bugfix: REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1333 r1337  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Bugfix: REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    299299          IF ( .NOT. outflow_l  .OR.  i > nxl )  THEN
    300300
    301              tend(:,j,i) = 0.0
     301             tend(:,j,i) = 0.0_wp
    302302             IF ( timestep_scheme(1:5) == 'runge' )  THEN
    303303                IF ( ws_scheme_mom )  THEN
     
    354354                         intermediate_timestep_count_max )  THEN
    355355                   DO  k = nzb_u_inner(j,i)+1, nzt
    356                       tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tu_m(k,j,i)
     356                      tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i)
    357357                   ENDDO
    358358                ENDIF
     
    365365          IF ( .NOT. outflow_s  .OR.  j > nys )  THEN
    366366
    367              tend(:,j,i) = 0.0
     367             tend(:,j,i) = 0.0_wp
    368368             IF ( timestep_scheme(1:5) == 'runge' )  THEN
    369369                IF ( ws_scheme_mom )  THEN
     
    413413                         intermediate_timestep_count_max )  THEN
    414414                   DO  k = nzb_v_inner(j,i)+1, nzt
    415                       tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tv_m(k,j,i)
     415                      tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i)
    416416                   ENDDO
    417417                ENDIF
     
    422422!
    423423!--       Tendency terms for w-velocity component
    424           tend(:,j,i) = 0.0
     424          tend(:,j,i) = 0.0_wp
    425425          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    426426             IF ( ws_scheme_mom )  THEN
     
    471471                      intermediate_timestep_count_max )  THEN
    472472                DO  k = nzb_w_inner(j,i)+1, nzt-1
    473                    tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tw_m(k,j,i)
     473                   tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i)
    474474                ENDDO
    475475             ENDIF
     
    484484!
    485485!--          Tendency terms for potential temperature
    486              tend(:,j,i) = 0.0
     486             tend(:,j,i) = 0.0_wp
    487487             IF ( timestep_scheme(1:5) == 'runge' )  THEN
    488488                   IF ( ws_scheme_sca )  THEN
     
    516516!
    517517!--          Consideration of heat sources within the plant canopy
    518              IF ( plant_canopy  .AND.  cthf /= 0.0 )  THEN
     518             IF ( plant_canopy  .AND.  cthf /= 0.0_wp )  THEN
    519519                CALL plant_canopy_model( i, j, 4 )
    520520             ENDIF
     
    551551                         intermediate_timestep_count_max )  THEN
    552552                   DO  k = nzb_s_inner(j,i)+1, nzt
    553                       tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    554                                       5.3125 * tpt_m(k,j,i)
     553                      tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     554                                      5.3125_wp * tpt_m(k,j,i)
    555555                   ENDDO
    556556                ENDIF
     
    565565!
    566566!--          Tendency-terms for salinity
    567              tend(:,j,i) = 0.0
     567             tend(:,j,i) = 0.0_wp
    568568             IF ( timestep_scheme(1:5) == 'runge' ) &
    569569             THEN
     
    588588                                        - tsc(5) * rdf_sc(k) *                 &
    589589                                          ( sa(k,j,i) - sa_init(k) )
    590                 IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
     590                IF ( sa_p(k,j,i) < 0.0_wp )  sa_p(k,j,i) = 0.1_wp * sa(k,j,i)
    591591             ENDDO
    592592
     
    601601                         intermediate_timestep_count_max )  THEN
    602602                   DO  k = nzb_s_inner(j,i)+1, nzt
    603                       tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    604                                       5.3125 * tsa_m(k,j,i)
     603                      tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     604                                      5.3125_wp * tsa_m(k,j,i)
    605605                   ENDDO
    606606                ENDIF
     
    620620!
    621621!--          Tendency-terms for total water content / scalar
    622              tend(:,j,i) = 0.0
     622             tend(:,j,i) = 0.0_wp
    623623             IF ( timestep_scheme(1:5) == 'runge' ) &
    624624             THEN
     
    666666                                      - tsc(5) * rdf_sc(k) *                   &
    667667                                        ( q(k,j,i) - q_init(k) )
    668                 IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
     668                IF ( q_p(k,j,i) < 0.0_wp )  q_p(k,j,i) = 0.1_wp * q(k,j,i)
    669669             ENDDO
    670670
     
    679679                         intermediate_timestep_count_max )  THEN
    680680                   DO  k = nzb_s_inner(j,i)+1, nzt
    681                       tq_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    682                                      5.3125 * tq_m(k,j,i)
     681                      tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     682                                     5.3125_wp * tq_m(k,j,i)
    683683                   ENDDO
    684684                ENDIF
     
    692692!
    693693!--             Calculate prognostic equation for rain water content
    694                 tend(:,j,i) = 0.0
     694                tend(:,j,i) = 0.0_wp
    695695                IF ( timestep_scheme(1:5) == 'runge' ) &
    696696                THEN
     
    718718                                                       tsc(3) * tqr_m(k,j,i) ) &
    719719                                           - tsc(5) * rdf_sc(k) * qr(k,j,i)
    720                    IF ( qr_p(k,j,i) < 0.0 )  qr_p(k,j,i) = 0.0
     720                   IF ( qr_p(k,j,i) < 0.0_wp )  qr_p(k,j,i) = 0.0_wp
    721721                ENDDO
    722722!
     
    730730                            intermediate_timestep_count_max )  THEN
    731731                      DO  k = nzb_s_inner(j,i)+1, nzt
    732                          tqr_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    733                                         5.3125 * tqr_m(k,j,i)
     732                         tqr_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     733                                        5.3125_wp * tqr_m(k,j,i)
    734734                      ENDDO
    735735                   ENDIF
     
    738738!
    739739!--             Calculate prognostic equation for rain drop concentration.
    740                 tend(:,j,i) = 0.0
     740                tend(:,j,i) = 0.0_wp
    741741                IF ( timestep_scheme(1:5) == 'runge' )  THEN
    742742                   IF ( ws_scheme_sca )  THEN
     
    764764                                                       tsc(3) * tnr_m(k,j,i) ) &
    765765                                           - tsc(5) * rdf_sc(k) * nr(k,j,i)
    766                    IF ( nr_p(k,j,i) < 0.0 )  nr_p(k,j,i) = 0.0
     766                   IF ( nr_p(k,j,i) < 0.0_wp )  nr_p(k,j,i) = 0.0_wp
    767767                ENDDO
    768768!
     
    776776                            intermediate_timestep_count_max )  THEN
    777777                      DO  k = nzb_s_inner(j,i)+1, nzt
    778                          tnr_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    779                                         5.3125 * tnr_m(k,j,i)
     778                         tnr_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     779                                        5.3125_wp * tnr_m(k,j,i)
    780780                      ENDDO
    781781                   ENDIF
     
    793793!
    794794!--          Tendency-terms for TKE
    795              tend(:,j,i) = 0.0
     795             tend(:,j,i) = 0.0_wp
    796796             IF ( timestep_scheme(1:5) == 'runge'  &
    797797                 .AND.  .NOT. use_upstream_for_tke )  THEN
     
    830830                e_p(k,j,i) = e(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +       &
    831831                                                  tsc(3) * te_m(k,j,i) )
    832                 IF ( e_p(k,j,i) < 0.0 )  e_p(k,j,i) = 0.1 * e(k,j,i)
     832                IF ( e_p(k,j,i) < 0.0_wp )  e_p(k,j,i) = 0.1_wp * e(k,j,i)
    833833             ENDDO
    834834
     
    843843                         intermediate_timestep_count_max )  THEN
    844844                   DO  k = nzb_s_inner(j,i)+1, nzt
    845                       te_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    846                                      5.3125 * te_m(k,j,i)
     845                      te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     846                                     5.3125_wp * te_m(k,j,i)
    847847                   ENDDO
    848848                ENDIF
     
    880880    CALL cpu_log( log_point(5), 'u-equation', 'start' )
    881881
    882     tend = 0.0
     882    tend = 0.0_wp
    883883    IF ( timestep_scheme(1:5) == 'runge' )  THEN
    884884       IF ( ws_scheme_mom )  THEN
     
    946946             DO  j = nys, nyn
    947947                DO  k = nzb_u_inner(j,i)+1, nzt
    948                    tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tu_m(k,j,i)
     948                   tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i)
    949949                ENDDO
    950950             ENDDO
     
    959959    CALL cpu_log( log_point(6), 'v-equation', 'start' )
    960960
    961     tend = 0.0
     961    tend = 0.0_wp
    962962    IF ( timestep_scheme(1:5) == 'runge' )  THEN
    963963       IF ( ws_scheme_mom )  THEN
     
    10221022             DO  j = nysv, nyn
    10231023                DO  k = nzb_v_inner(j,i)+1, nzt
    1024                    tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tv_m(k,j,i)
     1024                   tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i)
    10251025                ENDDO
    10261026             ENDDO
     
    10941094             DO  j = nys, nyn
    10951095                DO  k = nzb_w_inner(j,i)+1, nzt-1
    1096                    tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tw_m(k,j,i)
     1096                   tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i)
    10971097                ENDDO
    10981098             ENDDO
     
    11201120             sbt = 1.0
    11211121          ENDIF
    1122           tend = 0.0
     1122          tend = 0.0_wp
    11231123          CALL advec_s_bc( pt, 'pt' )
    11241124
     
    11281128!--    pt-tendency terms with no communication
    11291129       IF ( scalar_advec /= 'bc-scheme' )  THEN
    1130           tend = 0.0
     1130          tend = 0.0_wp
    11311131          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    11321132             IF ( ws_scheme_sca )  THEN
     
    11561156!
    11571157!--    Consideration of heat sources within the plant canopy
    1158        IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
     1158       IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN
    11591159          CALL plant_canopy_model( 4 )
    11601160       ENDIF
     
    12011201                DO  j = nys, nyn
    12021202                   DO  k = nzb_s_inner(j,i)+1, nzt
    1203                       tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    1204                                       5.3125 * tpt_m(k,j,i)
     1203                      tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     1204                                      5.3125_wp * tpt_m(k,j,i)
    12051205                   ENDDO
    12061206                ENDDO
     
    12291229             sbt = 1.0
    12301230          ENDIF
    1231           tend = 0.0
     1231          tend = 0.0_wp
    12321232          CALL advec_s_bc( sa, 'sa' )
    12331233
     
    12621262                                        - tsc(5) * rdf_sc(k) *                 &
    12631263                                          ( sa(k,j,i) - sa_init(k) )
    1264                 IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
     1264                IF ( sa_p(k,j,i) < 0.0_wp )  sa_p(k,j,i) = 0.1_wp * sa(k,j,i)
    12651265             ENDDO
    12661266          ENDDO
     
    12831283                DO  j = nys, nyn
    12841284                   DO  k = nzb_s_inner(j,i)+1, nzt
    1285                       tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + &
    1286                                       5.3125 * tsa_m(k,j,i)
     1285                      tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + &
     1286                                      5.3125_wp * tsa_m(k,j,i)
    12871287                   ENDDO
    12881288                ENDDO
     
    13151315!
    13161316!--          Bott-Chlond scheme always uses Euler time step. Thus:
    1317              sbt = 1.0
    1318           ENDIF
    1319           tend = 0.0
     1317             sbt = 1.0_wp
     1318          ENDIF
     1319          tend = 0.0_wp
    13201320          CALL advec_s_bc( q, 'q' )
    13211321
     
    13251325!--    Scalar/q-tendency terms with no communication
    13261326       IF ( scalar_advec /= 'bc-scheme' )  THEN
    1327           tend = 0.0
     1327          tend = 0.0_wp
    13281328          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    13291329             IF ( ws_scheme_sca )  THEN
     
    13711371                                      - tsc(5) * rdf_sc(k) *                   &
    13721372                                        ( q(k,j,i) - q_init(k) )
    1373                 IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
     1373                IF ( q_p(k,j,i) < 0.0_wp )  q_p(k,j,i) = 0.1_wp * q(k,j,i)
    13741374             ENDDO
    13751375          ENDDO
     
    13921392                DO  j = nys, nyn
    13931393                   DO  k = nzb_s_inner(j,i)+1, nzt
    1394                       tq_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tq_m(k,j,i)
     1394                      tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tq_m(k,j,i)
    13951395                   ENDDO
    13961396                ENDDO
     
    14171417!
    14181418!--             Bott-Chlond scheme always uses Euler time step. Thus:
    1419                 sbt = 1.0
    1420              ENDIF
    1421              tend = 0.0
     1419                sbt = 1.0_wp
     1420             ENDIF
     1421             tend = 0.0_wp
    14221422             CALL advec_s_bc( e, 'e' )
    14231423
     
    14291429       IF ( scalar_advec /= 'bc-scheme'  .OR.  use_upstream_for_tke )  THEN
    14301430          IF ( use_upstream_for_tke )  THEN
    1431              tend = 0.0
     1431             tend = 0.0_wp
    14321432             CALL advec_s_up( e )
    14331433          ELSE
    1434              tend = 0.0
     1434             tend = 0.0_wp
    14351435             IF ( timestep_scheme(1:5) == 'runge' )  THEN
    14361436                IF ( ws_scheme_sca )  THEN
     
    14721472                e_p(k,j,i) = e(k,j,i) + dt_3d * ( sbt * tend(k,j,i) +          &
    14731473                                                  tsc(3) * te_m(k,j,i) )
    1474                 IF ( e_p(k,j,i) < 0.0 )  e_p(k,j,i) = 0.1 * e(k,j,i)
     1474                IF ( e_p(k,j,i) < 0.0_wp )  e_p(k,j,i) = 0.1_wp * e(k,j,i)
    14751475             ENDDO
    14761476          ENDDO
     
    14931493                DO  j = nys, nyn
    14941494                   DO  k = nzb_s_inner(j,i)+1, nzt
    1495                       te_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * te_m(k,j,i)
     1495                      te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * te_m(k,j,i)
    14961496                   ENDDO
    14971497                ENDDO
     
    15451545          CALL advec_u_ws_acc
    15461546       ELSE
    1547           tend = 0.0    ! to be removed later??
     1547          tend = 0.0_wp   ! to be removed later??
    15481548          CALL advec_u_pw
    15491549       ENDIF
     
    15971597                   tu_m(k,j,i) = tend(k,j,i)
    15981598                ELSEIF ( runge_step == 2 )  THEN
    1599                    tu_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tu_m(k,j,i)
     1599                   tu_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tu_m(k,j,i)
    16001600                ENDIF
    16011601             ENDIF
     
    16151615          CALL advec_v_ws_acc
    16161616       ELSE
    1617           tend = 0.0    ! to be removed later??
     1617          tend = 0.0_wp    ! to be removed later??
    16181618          CALL advec_v_pw
    16191619       END IF
     
    16641664                   tv_m(k,j,i) = tend(k,j,i)
    16651665                ELSEIF ( runge_step == 2 )  THEN
    1666                    tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tv_m(k,j,i)
     1666                   tv_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tv_m(k,j,i)
    16671667                ENDIF
    16681668             ENDIF
     
    16821682          CALL advec_w_ws_acc
    16831683       ELSE
    1684           tend = 0.0    ! to be removed later??
     1684          tend = 0.0_wp    ! to be removed later??
    16851685          CALL advec_w_pw
    16861686       ENDIF
     
    17271727                   tw_m(k,j,i) = tend(k,j,i)
    17281728                ELSEIF ( runge_step == 2 )  THEN
    1729                    tw_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tw_m(k,j,i)
     1729                   tw_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tw_m(k,j,i)
    17301730                ENDIF
    17311731             ENDIF
     
    17541754             sbt = 1.0
    17551755          ENDIF
    1756           tend = 0.0
     1756          tend = 0.0_wp
    17571757          CALL advec_s_bc( pt, 'pt' )
    17581758
     
    17621762!--    pt-tendency terms with no communication
    17631763       IF ( scalar_advec /= 'bc-scheme' )  THEN
    1764           tend = 0.0
     1764          tend = 0.0_wp
    17651765          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    17661766             IF ( ws_scheme_sca )  THEN
    17671767                CALL advec_s_ws_acc( pt, 'pt' )
    17681768             ELSE
    1769                 tend = 0.0    ! to be removed later??
     1769                tend = 0.0_wp    ! to be removed later??
    17701770                CALL advec_s_pw( pt )
    17711771             ENDIF
     
    17911791!
    17921792!--    Consideration of heat sources within the plant canopy
    1793        IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
     1793       IF ( plant_canopy .AND. ( cthf /= 0.0_wp ) ) THEN
    17941794          CALL plant_canopy_model( 4 )
    17951795       ENDIF
     
    18271827                      tpt_m(k,j,i) = tend(k,j,i)
    18281828                   ELSEIF ( runge_step == 2 )  THEN
    1829                       tpt_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tpt_m(k,j,i)
     1829                      tpt_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tpt_m(k,j,i)
    18301830                   ENDIF
    18311831                ENDIF
     
    18531853!
    18541854!--          Bott-Chlond scheme always uses Euler time step. Thus:
    1855              sbt = 1.0
    1856           ENDIF
    1857           tend = 0.0
     1855             sbt = 1.0_wp
     1856          ENDIF
     1857          tend = 0.0_wp
    18581858          CALL advec_s_bc( sa, 'sa' )
    18591859
     
    18631863!--    sa-tendency terms with no communication
    18641864       IF ( scalar_advec /= 'bc-scheme' )  THEN
    1865           tend = 0.0
     1865          tend = 0.0_wp
    18661866          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    18671867             IF ( ws_scheme_sca )  THEN
     
    18881888                                        - tsc(5) * rdf_sc(k) *                 &
    18891889                                          ( sa(k,j,i) - sa_init(k) )
    1890                 IF ( sa_p(k,j,i) < 0.0 )  sa_p(k,j,i) = 0.1 * sa(k,j,i)
     1890                IF ( sa_p(k,j,i) < 0.0_wp )  sa_p(k,j,i) = 0.1_wp * sa(k,j,i)
    18911891!
    18921892!--             Tendencies for the next Runge-Kutta step
     
    18941894                   tsa_m(k,j,i) = tend(k,j,i)
    18951895                ELSEIF ( runge_step == 2 )  THEN
    1896                    tsa_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tsa_m(k,j,i)
     1896                   tsa_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tsa_m(k,j,i)
    18971897                ENDIF
    18981898             ENDDO
     
    19241924!
    19251925!--          Bott-Chlond scheme always uses Euler time step. Thus:
    1926              sbt = 1.0
    1927           ENDIF
    1928           tend = 0.0
     1926             sbt = 1.0_wp
     1927          ENDIF
     1928          tend = 0.0_wp
    19291929          CALL advec_s_bc( q, 'q' )
    19301930
     
    19341934!--    Scalar/q-tendency terms with no communication
    19351935       IF ( scalar_advec /= 'bc-scheme' )  THEN
    1936           tend = 0.0
     1936          tend = 0.0_wp
    19371937          IF ( timestep_scheme(1:5) == 'runge' )  THEN
    19381938             IF ( ws_scheme_sca )  THEN
     
    19801980                                      - tsc(5) * rdf_sc(k) *                   &
    19811981                                        ( q(k,j,i) - q_init(k) )
    1982                 IF ( q_p(k,j,i) < 0.0 )  q_p(k,j,i) = 0.1 * q(k,j,i)
     1982                IF ( q_p(k,j,i) < 0.0_wp )  q_p(k,j,i) = 0.1_wp * q(k,j,i)
    19831983!
    19841984!--             Tendencies for the next Runge-Kutta step
     
    19861986                   tq_m(k,j,i) = tend(k,j,i)
    19871987                ELSEIF ( runge_step == 2 )  THEN
    1988                    tq_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tq_m(k,j,i)
     1988                   tq_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * tq_m(k,j,i)
    19891989                ENDIF
    19901990             ENDDO
     
    20102010!
    20112011!--             Bott-Chlond scheme always uses Euler time step. Thus:
    2012                 sbt = 1.0
    2013              ENDIF
    2014              tend = 0.0
     2012                sbt = 1.0_wp
     2013             ENDIF
     2014             tend = 0.0_wp
    20152015             CALL advec_s_bc( e, 'e' )
    20162016
     
    20222022       IF ( scalar_advec /= 'bc-scheme'  .OR.  use_upstream_for_tke )  THEN
    20232023          IF ( use_upstream_for_tke )  THEN
    2024              tend = 0.0
     2024             tend = 0.0_wp
    20252025             CALL advec_s_up( e )
    20262026          ELSE
     
    20292029                   CALL advec_s_ws_acc( e, 'e' )
    20302030                ELSE
    2031                    tend = 0.0    ! to be removed later??
     2031                   tend = 0.0_wp    ! to be removed later??
    20322032                   CALL advec_s_pw( e )
    20332033                ENDIF
    20342034             ELSE
    2035                 tend = 0.0    ! to be removed later??
     2035                tend = 0.0_wp    ! to be removed later??
    20362036                CALL advec_s_up( e )
    20372037             ENDIF
     
    20712071                   e_p(k,j,i) = e(k,j,i) + dt_3d * ( sbt * tend(k,j,i) +          &
    20722072                                                     tsc(3) * te_m(k,j,i) )
    2073                    IF ( e_p(k,j,i) < 0.0 )  e_p(k,j,i) = 0.1 * e(k,j,i)
     2073                   IF ( e_p(k,j,i) < 0.0_wp )  e_p(k,j,i) = 0.1_wp * e(k,j,i)
    20742074!
    20752075!--                Tendencies for the next Runge-Kutta step
     
    20772077                      te_m(k,j,i) = tend(k,j,i)
    20782078                   ELSEIF ( runge_step == 2 )  THEN
    2079                       te_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * te_m(k,j,i)
     2079                      te_m(k,j,i) = -9.5625_wp * tend(k,j,i) + 5.3125_wp * te_m(k,j,i)
    20802080                   ENDIF
    20812081                ENDIF
Note: See TracChangeset for help on using the changeset viewer.