Changeset 1322 for palm/trunk


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

Location:
palm/trunk/SOURCE
Files:
43 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r1321 r1322  
    269269diffusivities.o: modules.o mod_kinds.o
    270270disturb_field.o: modules.o cpulog.o mod_kinds.o random_function.o
    271 disturb_heatflux.o: modules.o mod_kinds.o
     271disturb_heatflux.o: modules.o cpulog.o mod_kinds.o
    272272eqn_state_seawater.o: modules.o mod_kinds.o
    273273exchange_horiz.o: modules.o cpulog.o mod_kinds.o
  • palm/trunk/SOURCE/advec_ws.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    195195!
    196196!--    Set the appropriate factors for scalar and momentum advection.
    197        adv_sca_5 = 1./60.
    198        adv_sca_3 = 1./12.
    199        adv_sca_1 = 1./2.
    200        adv_mom_5 = 1./120.
    201        adv_mom_3 = 1./24.
    202        adv_mom_1 = 1./4.
     197       adv_sca_5 = 1.0_wp /  60.0_wp
     198       adv_sca_3 = 1.0_wp /  12.0_wp
     199       adv_sca_1 = 1.0_wp /   2.0_wp
     200       adv_mom_5 = 1.0_wp / 120.0_wp
     201       adv_mom_3 = 1.0_wp /  24.0_wp
     202       adv_mom_1 = 1.0_wp /   4.0_wp
    203203!         
    204204!--    Arrays needed for statical evaluation of fluxes.
     
    11741174                              + ( flux_r(k) *                                 &
    11751175                                ( u_comp(k) - 2.0 * hom(k,1,1,0) )            &
    1176                               / ( u_comp(k) - gu + 1.0E-20      )             &
     1176                              / ( u_comp(k) - gu + 1.0E-20_wp   )             &
    11771177                              +   diss_r(k) *                                 &
    11781178                                  ABS( u_comp(k) - 2.0 * hom(k,1,1,0) )       &
    1179                               / ( ABS( u_comp(k) - gu ) + 1.0E-20 ) )         &
     1179                              / ( ABS( u_comp(k) - gu ) + 1.0E-20_wp ) )      &
    11801180                              *   weight_substep(intermediate_timestep_count)
    11811181!
     
    12781278                              + ( flux_r(k) *                                 &
    12791279                                ( u_comp(k) - 2.0 * hom(k,1,1,0) )            &
    1280                               / ( u_comp(k) - gu + 1.0E-20      )             &
     1280                              / ( u_comp(k) - gu + 1.0E-20_wp   )             &
    12811281                              +   diss_r(k) *                                 &
    12821282                                  ABS( u_comp(k) - 2.0 * hom(k,1,1,0) )       &
    1283                               / ( ABS( u_comp(k) - gu ) + 1.0E-20 ) )         &
     1283                              / ( ABS( u_comp(k) - gu ) + 1.0E-20_wp ) )      &
    12841284                              *   weight_substep(intermediate_timestep_count)
    12851285!
     
    16311631             + ( flux_n(k)                                                    &
    16321632             * ( v_comp(k) - 2.0 * hom(k,1,2,0) )                             &
    1633              / ( v_comp(k) - gv + 1.0E-20 )                                   &
     1633             / ( v_comp(k) - gv + 1.0E-20_wp )                                &
    16341634             +   diss_n(k)                                                    &
    16351635             *   ABS( v_comp(k) - 2.0 * hom(k,1,2,0) )                        &
    1636              / ( ABS( v_comp(k) - gv ) +1.0E-20 ) )                           &
     1636             / ( ABS( v_comp(k) - gv ) +1.0E-20_wp ) )                        &
    16371637             *   weight_substep(intermediate_timestep_count)
    16381638!
     
    17401740             + ( flux_n(k)                                                    &
    17411741             * ( v_comp(k) - 2.0 * hom(k,1,2,0) )                             &
    1742              / ( v_comp(k) - gv + 1.0E-20 )                                   &
     1742             / ( v_comp(k) - gv + 1.0E-20_wp )                                &
    17431743             +   diss_n(k)                                                    &
    17441744             *   ABS( v_comp(k) - 2.0 * hom(k,1,2,0) )                        &
    1745              / ( ABS( v_comp(k) - gv ) +1.0E-20 ) )                           &
     1745             / ( ABS( v_comp(k) - gv ) +1.0E-20_wp ) )                        &
    17461746             *   weight_substep(intermediate_timestep_count)
    17471747!
     
    33123312                              + ( flux_r(k) *                                 &
    33133313                                ( u_comp(k) - 2.0 * hom(k,1,1,0) )            &
    3314                               / ( u_comp(k) - gu + 1.0E-20      )             &
     3314                              / ( u_comp(k) - gu + 1.0E-20_wp    )            &
    33153315                              +   diss_r(k) *                                 &
    33163316                                  ABS( u_comp(k) - 2.0 * hom(k,1,1,0) )       &
    3317                               / ( ABS( u_comp(k) - gu ) + 1.0E-20 ) )         &
     3317                              / ( ABS( u_comp(k) - gu ) + 1.0E-20_wp ) )      &
    33183318                              *   weight_substep(intermediate_timestep_count)
    33193319!
     
    34183418                              + ( flux_r(k) *                                 &
    34193419                                ( u_comp(k) - 2.0 * hom(k,1,1,0) )            &
    3420                               / ( u_comp(k) - gu + 1.0E-20      )             &
     3420                              / ( u_comp(k) - gu + 1.0E-20_wp   )             &
    34213421                              +   diss_r(k) *                                 &
    34223422                                  ABS( u_comp(k) - 2.0 * hom(k,1,1,0) )       &
    3423                               / ( ABS( u_comp(k) - gu ) + 1.0E-20 ) )         &
     3423                              / ( ABS( u_comp(k) - gu ) + 1.0E-20_wp ) )      &
    34243424                              *   weight_substep(intermediate_timestep_count)
    34253425!
     
    37493749!                              + ( flux_r    *                                 &
    37503750!                                ( u_comp    - 2.0 * hom(k,1,1,0) )            &
    3751 !                              / ( u_comp    - gu + 1.0E-20      )             &
     3751!                              / ( u_comp    - gu + 1.0E-20_wp   )             &
    37523752!                              +   diss_r    *                                 &
    37533753!                                  ABS( u_comp    - 2.0 * hom(k,1,1,0) )       &
    3754 !                              / ( ABS( u_comp    - gu ) + 1.0E-20 ) )         &
     3754!                              / ( ABS( u_comp    - gu ) + 1.0E-20_wp ) )      &
    37553755!                              *   weight_substep(intermediate_timestep_count)
    37563756!
     
    41084108                      + ( flux_n(k)                                           &
    41094109                      * ( v_comp(k) - 2.0 * hom(k,1,2,0) )                    &
    4110                       / ( v_comp(k) - gv + 1.0E-20 )                          &
     4110                      / ( v_comp(k) - gv + 1.0E-20_wp )                       &
    41114111                      +   diss_n(k)                                           &
    41124112                      *   ABS( v_comp(k) - 2.0 * hom(k,1,2,0) )               &
    4113                       / ( ABS( v_comp(k) - gv ) +1.0E-20 ) )                  &
     4113                      / ( ABS( v_comp(k) - gv ) +1.0E-20_wp ) )               &
    41144114                      *   weight_substep(intermediate_timestep_count)
    41154115!
     
    42214221                         + ( flux_n(k)                                        &
    42224222                         * ( v_comp(k) - 2.0 * hom(k,1,2,0) )                 &
    4223                          / ( v_comp(k) - gv + 1.0E-20 )                       &
     4223                         / ( v_comp(k) - gv + 1.0E-20_wp )                    &
    42244224                         +   diss_n(k)                                        &
    42254225                         *   ABS( v_comp(k) - 2.0 * hom(k,1,2,0) )            &
    4226                          / ( ABS( v_comp(k) - gv ) +1.0E-20 ) )               &
     4226                         / ( ABS( v_comp(k) - gv ) +1.0E-20_wp ) )            &
    42274227                         *   weight_substep(intermediate_timestep_count)
    42284228!
     
    45544554!                      + ( flux_n                                           &
    45554555!                      * ( v_comp - 2.0 * hom(k,1,2,0) )                    &
    4556 !                      / ( v_comp - gv + 1.0E-20 )                          &
     4556!                      / ( v_comp - gv + 1.0E-20_wp )                       &
    45574557!                      +   diss_n                                           &
    45584558!                      *   ABS( v_comp - 2.0 * hom(k,1,2,0) )               &
    4559 !                      / ( ABS( v_comp - gv ) +1.0E-20 ) )                  &
     4559!                      / ( ABS( v_comp - gv ) +1.0E-20_wp ) )               &
    45604560!                      *   weight_substep(intermediate_timestep_count)
    45614561!
  • palm/trunk/SOURCE/average_3d_data.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    9999                DO  j = nysg, nyng
    100100                   DO  k = nzb, nzt+1
    101                       e_av(k,j,i) = e_av(k,j,i) / REAL( average_count_3d )
     101                      e_av(k,j,i) = e_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    102102                   ENDDO
    103103                ENDDO
     
    107107             DO  i = nxlg, nxrg
    108108                DO  j = nysg, nyng
    109                    qsws_av(j,i) = qsws_av(j,i) / REAL( average_count_3d )
     109                   qsws_av(j,i) = qsws_av(j,i) / REAL( average_count_3d, KIND=wp )
    110110                ENDDO
    111111             ENDDO
     
    115115                DO  j = nysg, nyng
    116116                   DO  k = nzb, nzt+1
    117                       lpt_av(k,j,i) = lpt_av(k,j,i) / REAL( average_count_3d )
     117                      lpt_av(k,j,i) = lpt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    118118                   ENDDO
    119119                ENDDO
     
    123123             DO  i = nxlg, nxrg
    124124                DO  j = nysg, nyng
    125                    lwp_av(j,i) = lwp_av(j,i) / REAL( average_count_3d )
     125                   lwp_av(j,i) = lwp_av(j,i) / REAL( average_count_3d, KIND=wp )
    126126                ENDDO
    127127             ENDDO
     
    131131                DO  j = nysg, nyng
    132132                   DO  k = nzb, nzt+1
    133                       nr_av(k,j,i) = nr_av(k,j,i) / REAL( average_count_3d )
     133                      nr_av(k,j,i) = nr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    134134                   ENDDO
    135135                ENDDO
     
    140140                DO  j = nysg, nyng
    141141                   DO  k = nzb, nzt+1
    142                       p_av(k,j,i) = p_av(k,j,i) / REAL( average_count_3d )
     142                      p_av(k,j,i) = p_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    143143                   ENDDO
    144144                ENDDO
     
    149149                DO  j = nys, nyn
    150150                   DO  k = nzb, nzt+1
    151                       pc_av(k,j,i) = pc_av(k,j,i) / REAL( average_count_3d )
     151                      pc_av(k,j,i) = pc_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    152152                   ENDDO
    153153                ENDDO
     
    158158                DO  j = nys, nyn
    159159                   DO  k = nzb, nzt+1
    160                       pr_av(k,j,i) = pr_av(k,j,i) / REAL( average_count_3d )
     160                      pr_av(k,j,i) = pr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    161161                   ENDDO
    162162                ENDDO
     
    167167                DO  j = nysg, nyng
    168168                   precipitation_rate_av(j,i) = precipitation_rate_av(j,i) /   &
    169                                                 REAL( average_count_3d )
     169                                                REAL( average_count_3d, KIND=wp )
    170170                ENDDO
    171171             ENDDO
     
    175175                DO  j = nysg, nyng
    176176                   DO  k = nzb, nzt+1
    177                       pt_av(k,j,i) = pt_av(k,j,i) / REAL( average_count_3d )
     177                      pt_av(k,j,i) = pt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    178178                   ENDDO
    179179                ENDDO
     
    184184                DO  j = nysg, nyng
    185185                   DO  k = nzb, nzt+1
    186                       q_av(k,j,i) = q_av(k,j,i) / REAL( average_count_3d )
     186                      q_av(k,j,i) = q_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    187187                   ENDDO
    188188                ENDDO
     
    193193                DO  j = nysg, nyng
    194194                   DO  k = nzb, nzt+1
    195                       qc_av(k,j,i) = qc_av(k,j,i) / REAL( average_count_3d )
     195                      qc_av(k,j,i) = qc_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    196196                   ENDDO
    197197                ENDDO
     
    202202                DO  j = nysg, nyng
    203203                   DO  k = nzb, nzt+1
    204                       ql_av(k,j,i) = ql_av(k,j,i) / REAL( average_count_3d )
     204                      ql_av(k,j,i) = ql_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    205205                   ENDDO
    206206                ENDDO
     
    211211                DO  j = nysg, nyng
    212212                   DO  k = nzb, nzt+1
    213                       ql_c_av(k,j,i) = ql_c_av(k,j,i) / REAL( average_count_3d )
     213                      ql_c_av(k,j,i) = ql_c_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    214214                   ENDDO
    215215                ENDDO
     
    220220                DO  j = nysg, nyng
    221221                   DO  k = nzb, nzt+1
    222                       ql_v_av(k,j,i) = ql_v_av(k,j,i) / REAL( average_count_3d )
     222                      ql_v_av(k,j,i) = ql_v_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    223223                   ENDDO
    224224                ENDDO
     
    230230                   DO  k = nzb, nzt+1
    231231                      ql_vp_av(k,j,i) = ql_vp_av(k,j,i) /                      &
    232                                         REAL( average_count_3d )
     232                                        REAL( average_count_3d, KIND=wp )
    233233                   ENDDO
    234234                ENDDO
     
    239239                DO  j = nysg, nyng
    240240                   DO  k = nzb, nzt+1
    241                       qr_av(k,j,i) = qr_av(k,j,i) / REAL( average_count_3d )
     241                      qr_av(k,j,i) = qr_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    242242                   ENDDO
    243243                ENDDO
     
    248248                DO  j = nysg, nyng
    249249                   DO  k = nzb, nzt+1
    250                       qv_av(k,j,i) = qv_av(k,j,i) / REAL( average_count_3d )
     250                      qv_av(k,j,i) = qv_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    251251                   ENDDO
    252252                ENDDO
     
    257257                DO  j = nysg, nyng
    258258                   DO  k = nzb, nzt+1
    259                       rho_av(k,j,i) = rho_av(k,j,i) / REAL( average_count_3d )
     259                      rho_av(k,j,i) = rho_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    260260                   ENDDO
    261261                ENDDO
     
    266266                DO  j = nysg, nyng
    267267                   DO  k = nzb, nzt+1
    268                       s_av(k,j,i) = s_av(k,j,i) / REAL( average_count_3d )
     268                      s_av(k,j,i) = s_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    269269                   ENDDO
    270270                ENDDO
     
    275275                DO  j = nysg, nyng
    276276                   DO  k = nzb, nzt+1
    277                       sa_av(k,j,i) = sa_av(k,j,i) / REAL( average_count_3d )
     277                      sa_av(k,j,i) = sa_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    278278                   ENDDO
    279279                ENDDO
     
    283283             DO  i = nxlg, nxrg
    284284                DO  j = nysg, nyng
    285                    shf_av(j,i) = shf_av(j,i) / REAL( average_count_3d )
     285                   shf_av(j,i) = shf_av(j,i) / REAL( average_count_3d, KIND=wp )
    286286                ENDDO
    287287             ENDDO
     
    290290             DO  i = nxlg, nxrg
    291291                DO  j = nysg, nyng
    292                    ts_av(j,i) = ts_av(j,i) / REAL( average_count_3d )
     292                   ts_av(j,i) = ts_av(j,i) / REAL( average_count_3d, KIND=wp )
    293293                ENDDO
    294294             ENDDO
     
    298298                DO  j = nysg, nyng
    299299                   DO  k = nzb, nzt+1
    300                       u_av(k,j,i) = u_av(k,j,i) / REAL( average_count_3d )
     300                      u_av(k,j,i) = u_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    301301                   ENDDO
    302302                ENDDO
     
    306306             DO  i = nxlg, nxrg
    307307                DO  j = nysg, nyng
    308                    us_av(j,i) = us_av(j,i) / REAL( average_count_3d )
     308                   us_av(j,i) = us_av(j,i) / REAL( average_count_3d, KIND=wp )
    309309                ENDDO
    310310             ENDDO
     
    314314                DO  j = nysg, nyng
    315315                   DO  k = nzb, nzt+1
    316                       v_av(k,j,i) = v_av(k,j,i) / REAL( average_count_3d )
     316                      v_av(k,j,i) = v_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    317317                   ENDDO
    318318                ENDDO
     
    323323                DO  j = nysg, nyng
    324324                   DO  k = nzb, nzt+1
    325                       vpt_av(k,j,i) = vpt_av(k,j,i) / REAL( average_count_3d )
     325                      vpt_av(k,j,i) = vpt_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    326326                   ENDDO
    327327                ENDDO
     
    332332                DO  j = nysg, nyng
    333333                   DO  k = nzb, nzt+1
    334                       w_av(k,j,i) = w_av(k,j,i) / REAL( average_count_3d )
     334                      w_av(k,j,i) = w_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    335335                   ENDDO
    336336                ENDDO
     
    340340             DO  i = nxlg, nxrg
    341341                DO  j = nysg, nyng
    342                    z0_av(j,i) = z0_av(j,i) / REAL( average_count_3d )
     342                   z0_av(j,i) = z0_av(j,i) / REAL( average_count_3d, KIND=wp )
    343343                ENDDO
    344344             ENDDO
     
    347347             DO  i = nxlg, nxrg
    348348                DO  j = nysg, nyng
    349                    z0h_av(j,i) = z0h_av(j,i) / REAL( average_count_3d )
     349                   z0h_av(j,i) = z0h_av(j,i) / REAL( average_count_3d, KIND=wp )
    350350                ENDDO
    351351             ENDDO
  • palm/trunk/SOURCE/calc_radiation.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! exponent 4.0 changed to integer
    2323!
    2424! Former revisions:
     
    127127
    128128             temperature     = pt(nzb,j,i) * t_d_pt(nzb) + l_d_cp * ql(nzb,j,i)
    129              blackbody_emission(nzb) = sigma * temperature**4.0
     129             blackbody_emission(nzb) = sigma * temperature**4
    130130
    131131             DO  k = nzb_2d(j,i)+1, nzt
     
    139139
    140140                temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
    141                 blackbody_emission(k) = sigma * temperature**4.0
     141                blackbody_emission(k) = sigma * temperature**4
    142142
    143143             ENDDO
     
    149149             temperature       = pt(nzt+1,j,i) * t_d_pt(nzt+1) + l_d_cp *      &
    150150                                 ql(nzt+1,j,i)
    151              blackbody_emission(nzt+1) = sigma * temperature**4.0
     151             blackbody_emission(nzt+1) = sigma * temperature**4
    152152
    153153!
     
    267267
    268268       temperature     = pt(nzb,j,i) * t_d_pt(nzb) + l_d_cp * ql(nzb,j,i)
    269        blackbody_emission(nzb) = sigma * temperature**4.0
     269       blackbody_emission(nzb) = sigma * temperature**4
    270270
    271271       DO  k = nzb_2d(j,i)+1, nzt
     
    277277
    278278          temperature     = pt(k,j,i) * t_d_pt(k) + l_d_cp * ql(k,j,i)
    279           blackbody_emission(k) = sigma * temperature**4.0
     279          blackbody_emission(k) = sigma * temperature**4
    280280
    281281       ENDDO
     
    286286       temperature       = pt(nzt+1,j,i) * t_d_pt(nzt+1) + l_d_cp *            &
    287287                           ql(nzt+1,j,i)
    288        blackbody_emission(nzt+1) = sigma * temperature**4.0
     288       blackbody_emission(nzt+1) = sigma * temperature**4
    289289
    290290!
  • palm/trunk/SOURCE/check_parameters.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! some REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    259259    IF ( coupling_mode /= 'uncoupled')  THEN
    260260
    261        IF ( dt_coupling == 9999999.9 )  THEN
     261       IF ( dt_coupling == 9999999.9_wp )  THEN
    262262          message_string = 'dt_coupling is not set but required for coup' // &
    263263                           'ling mode "' //  TRIM( coupling_mode ) // '"'
     
    920920!
    921921!--    In case of no given gradients for ug, choose a zero gradient
    922        IF ( ug_vertical_gradient_level(1) == -9999999.9 )  THEN
     922       IF ( ug_vertical_gradient_level(1) == -9999999.9_wp )  THEN
    923923          ug_vertical_gradient_level(1) = 0.0
    924924       ENDIF 
     
    984984!
    985985!--    In case of no given gradients for vg, choose a zero gradient
    986        IF ( vg_vertical_gradient_level(1) == -9999999.9 )  THEN
     986       IF ( vg_vertical_gradient_level(1) == -9999999.9_wp )  THEN
    987987          vg_vertical_gradient_level(1) = 0.0
    988988       ENDIF
     
    991991!--    Let the initial wind profiles be the calculated ug/vg profiles or
    992992!--    interpolate them from wind profile data (if given)
    993        IF ( u_profile(1) == 9999999.9  .AND.  v_profile(1) == 9999999.9 )  THEN
     993       IF ( u_profile(1) == 9999999.9_wp  .AND.  v_profile(1) == 9999999.9_wp )  THEN
    994994
    995995          u_init = ug
     
    10181018             ENDIF
    10191019
    1020              IF ( kk < 100 .AND. uv_heights(kk+1) /= 9999999.9 )  THEN
     1020             IF ( kk < 100 .AND. uv_heights(kk+1) /= 9999999.9_wp )  THEN
    10211021                u_init(k) = u_profile(kk) + ( zu(k) - uv_heights(kk) ) /       &
    10221022                                       ( uv_heights(kk+1) - uv_heights(kk) ) * &
     
    11001100!--    In case of no given temperature gradients, choose gradient of neutral
    11011101!--    stratification
    1102        IF ( pt_vertical_gradient_level(1) == -9999999.9 )  THEN
     1102       IF ( pt_vertical_gradient_level(1) == -9999999.9_wp )  THEN
    11031103          pt_vertical_gradient_level(1) = 0.0
    11041104       ENDIF
     
    12361236!--       In case of no given leaf area density gradients, choose a vanishing
    12371237!--       gradient
    1238           IF ( lad_vertical_gradient_level(1) == -9999999.9 ) THEN
     1238          IF ( lad_vertical_gradient_level(1) == -9999999.9_wp ) THEN
    12391239             lad_vertical_gradient_level(1) = 0.0
    12401240          ENDIF
     
    12471247!-- Initialize large scale subsidence if required
    12481248    If ( large_scale_subsidence )  THEN
    1249        IF ( subs_vertical_gradient_level(1) /= -9999999.9 .AND. &
     1249       IF ( subs_vertical_gradient_level(1) /= -9999999.9_wp .AND. &
    12501250                                     .NOT. large_scale_forcing )  THEN
    12511251          CALL init_w_subsidence
     
    12551255!--    are read in from file LSF_DATA
    12561256
    1257        IF ( subs_vertical_gradient_level(1) == -9999999.9 .AND. &
     1257       IF ( subs_vertical_gradient_level(1) == -9999999.9_wp .AND. &
    12581258                                     .NOT. large_scale_forcing )  THEN
    12591259          message_string = 'There is no default large scale vertical ' // &
     
    12641264       ENDIF
    12651265    ELSE
    1266         IF ( subs_vertical_gradient_level(1) /= -9999999.9 )  THEN
     1266        IF ( subs_vertical_gradient_level(1) /= -9999999.9_wp )  THEN
    12671267           message_string = 'Enable usage of large scale subsidence by ' // &
    12681268                            'setting large_scale_subsidence = .T..'
     
    12841284    ELSEIF ( reference_state == 'single_value' )  THEN
    12851285       use_single_reference_value = .TRUE.
    1286        IF ( pt_reference == 9999999.9 )  pt_reference = pt_surface
     1286       IF ( pt_reference == 9999999.9_wp )  pt_reference = pt_surface
    12871287       vpt_reference = pt_reference * ( 1.0 + 0.61 * q_surface )
    12881288    ELSE
     
    15281528    ENDIF
    15291529
    1530     IF ( surface_heatflux == 9999999.9  )  THEN
     1530    IF ( surface_heatflux == 9999999.9_wp  )  THEN
    15311531       constant_heatflux = .FALSE.
    15321532       IF ( large_scale_forcing )  THEN
     
    15481548    ENDIF
    15491549
    1550     IF ( top_heatflux     == 9999999.9 )  constant_top_heatflux = .FALSE.
     1550    IF ( top_heatflux     == 9999999.9_wp )  constant_top_heatflux = .FALSE.
    15511551
    15521552    IF ( neutral )  THEN
    15531553
    1554        IF ( surface_heatflux /= 0.0  .AND.  surface_heatflux /= 9999999.9 ) &
     1554       IF ( surface_heatflux /= 0.0  .AND.  surface_heatflux /= 9999999.9_wp ) &
    15551555       THEN
    15561556          message_string = 'heatflux must not be set for pure neutral flow'
     
    15581558       ENDIF
    15591559
    1560        IF ( top_heatflux /= 0.0  .AND.  top_heatflux /= 9999999.9 ) &
     1560       IF ( top_heatflux /= 0.0  .AND.  top_heatflux /= 9999999.9_wp ) &
    15611561       THEN
    15621562          message_string = 'heatflux must not be set for pure neutral flow'
     
    15661566    ENDIF
    15671567
    1568     IF ( top_momentumflux_u /= 9999999.9  .AND.  &
    1569          top_momentumflux_v /= 9999999.9 )  THEN
     1568    IF ( top_momentumflux_u /= 9999999.9_wp  .AND.  &
     1569         top_momentumflux_v /= 9999999.9_wp )  THEN
    15701570       constant_top_momentumflux = .TRUE.
    1571     ELSEIF (  .NOT. ( top_momentumflux_u == 9999999.9  .AND.  &
    1572            top_momentumflux_v == 9999999.9 ) )  THEN
     1571    ELSEIF (  .NOT. ( top_momentumflux_u == 9999999.9_wp  .AND.  &
     1572           top_momentumflux_v == 9999999.9_wp ) )  THEN
    15731573       message_string = 'both, top_momentumflux_u AND top_momentumflux_v ' // &
    15741574                        'must be set'
     
    16171617       ENDIF
    16181618
    1619        IF ( top_salinityflux == 9999999.9 )  constant_top_salinityflux = .FALSE.
    1620        IF ( ibc_sa_t == 1  .AND.   top_salinityflux == 9999999.9 )  THEN
     1619       IF ( top_salinityflux == 9999999.9_wp )  constant_top_salinityflux = .FALSE.
     1620       IF ( ibc_sa_t == 1  .AND.   top_salinityflux == 9999999.9_wp )  THEN
    16211621          message_string = 'boundary condition: bc_sa_t = "' // &
    16221622                           TRIM( bc_sa_t ) // '" requires to set ' // &
     
    16671667       ENDIF
    16681668
    1669        IF ( surface_waterflux == 9999999.9  )  THEN
     1669       IF ( surface_waterflux == 9999999.9_wp  )  THEN
    16701670          constant_waterflux = .FALSE.
    16711671          IF ( large_scale_forcing )  THEN
     
    18171817!-- Set the default intervals for data output, if necessary
    18181818!-- NOTE: dt_dosp has already been set in package_parin
    1819     IF ( dt_data_output /= 9999999.9 )  THEN
    1820        IF ( dt_dopr           == 9999999.9 )  dt_dopr           = dt_data_output
    1821        IF ( dt_dopts          == 9999999.9 )  dt_dopts          = dt_data_output
    1822        IF ( dt_do2d_xy        == 9999999.9 )  dt_do2d_xy        = dt_data_output
    1823        IF ( dt_do2d_xz        == 9999999.9 )  dt_do2d_xz        = dt_data_output
    1824        IF ( dt_do2d_yz        == 9999999.9 )  dt_do2d_yz        = dt_data_output
    1825        IF ( dt_do3d           == 9999999.9 )  dt_do3d           = dt_data_output
    1826        IF ( dt_data_output_av == 9999999.9 )  dt_data_output_av = dt_data_output
     1819    IF ( dt_data_output /= 9999999.9_wp )  THEN
     1820       IF ( dt_dopr           == 9999999.9_wp )  dt_dopr           = dt_data_output
     1821       IF ( dt_dopts          == 9999999.9_wp )  dt_dopts          = dt_data_output
     1822       IF ( dt_do2d_xy        == 9999999.9_wp )  dt_do2d_xy        = dt_data_output
     1823       IF ( dt_do2d_xz        == 9999999.9_wp )  dt_do2d_xz        = dt_data_output
     1824       IF ( dt_do2d_yz        == 9999999.9_wp )  dt_do2d_yz        = dt_data_output
     1825       IF ( dt_do3d           == 9999999.9_wp )  dt_do3d           = dt_data_output
     1826       IF ( dt_data_output_av == 9999999.9_wp )  dt_data_output_av = dt_data_output
    18271827       DO  mid = 1, max_masks
    1828           IF ( dt_domask(mid) == 9999999.9 )  dt_domask(mid)    = dt_data_output
     1828          IF ( dt_domask(mid) == 9999999.9_wp )  dt_domask(mid)    = dt_data_output
    18291829       ENDDO
    18301830    ENDIF
     
    18321832!
    18331833!-- Set the default skip time intervals for data output, if necessary
    1834     IF ( skip_time_dopr    == 9999999.9 ) &
     1834    IF ( skip_time_dopr    == 9999999.9_wp ) &
    18351835                                       skip_time_dopr    = skip_time_data_output
    1836     IF ( skip_time_dosp    == 9999999.9 ) &
     1836    IF ( skip_time_dosp    == 9999999.9_wp ) &
    18371837                                       skip_time_dosp    = skip_time_data_output
    1838     IF ( skip_time_do2d_xy == 9999999.9 ) &
     1838    IF ( skip_time_do2d_xy == 9999999.9_wp ) &
    18391839                                       skip_time_do2d_xy = skip_time_data_output
    1840     IF ( skip_time_do2d_xz == 9999999.9 ) &
     1840    IF ( skip_time_do2d_xz == 9999999.9_wp ) &
    18411841                                       skip_time_do2d_xz = skip_time_data_output
    1842     IF ( skip_time_do2d_yz == 9999999.9 ) &
     1842    IF ( skip_time_do2d_yz == 9999999.9_wp ) &
    18431843                                       skip_time_do2d_yz = skip_time_data_output
    1844     IF ( skip_time_do3d    == 9999999.9 ) &
     1844    IF ( skip_time_do3d    == 9999999.9_wp ) &
    18451845                                       skip_time_do3d    = skip_time_data_output
    1846     IF ( skip_time_data_output_av == 9999999.9 ) &
     1846    IF ( skip_time_data_output_av == 9999999.9_wp ) &
    18471847                                skip_time_data_output_av = skip_time_data_output
    18481848    DO  mid = 1, max_masks
    1849        IF ( skip_time_domask(mid) == 9999999.9 ) &
     1849       IF ( skip_time_domask(mid) == 9999999.9_wp ) &
    18501850                                skip_time_domask(mid)    = skip_time_data_output
    18511851    ENDDO
     
    18601860    ENDIF
    18611861
    1862     IF ( averaging_interval_pr == 9999999.9 )  THEN
     1862    IF ( averaging_interval_pr == 9999999.9_wp )  THEN
    18631863       averaging_interval_pr = averaging_interval
    18641864    ENDIF
     
    18701870    ENDIF
    18711871
    1872     IF ( averaging_interval_sp == 9999999.9 )  THEN
     1872    IF ( averaging_interval_sp == 9999999.9_wp )  THEN
    18731873       averaging_interval_sp = averaging_interval
    18741874    ENDIF
     
    18821882!
    18831883!-- Set the default interval for profiles entering the temporal average
    1884     IF ( dt_averaging_input_pr == 9999999.9 )  THEN
     1884    IF ( dt_averaging_input_pr == 9999999.9_wp )  THEN
    18851885       dt_averaging_input_pr = dt_averaging_input
    18861886    ENDIF
     
    18891889!-- Set the default interval for the output of timeseries to a reasonable
    18901890!-- value (tries to minimize the number of calls of flow_statistics)
    1891     IF ( dt_dots == 9999999.9 )  THEN
     1891    IF ( dt_dots == 9999999.9_wp )  THEN
    18921892       IF ( averaging_interval_pr == 0.0 )  THEN
    18931893          dt_dots = MIN( dt_run_control, dt_dopr )
     
    19161916!-- Set the default value for the integration interval of precipitation amount
    19171917    IF ( precipitation )  THEN
    1918        IF ( precipitation_amount_interval == 9999999.9 )  THEN
     1918       IF ( precipitation_amount_interval == 9999999.9_wp )  THEN
    19191919          precipitation_amount_interval = dt_do2d_xy
    19201920       ELSE
     
    32773277!
    32783278!-- Determine upper and lower hight level indices for random perturbations
    3279     IF ( disturbance_level_b == -9999999.9 )  THEN
     3279    IF ( disturbance_level_b == -9999999.9_wp )  THEN
    32803280       IF ( ocean ) THEN
    32813281          disturbance_level_b     = zu((nzt*2)/3)
     
    33023302    ENDIF
    33033303
    3304     IF ( disturbance_level_t == -9999999.9 )  THEN
     3304    IF ( disturbance_level_t == -9999999.9_wp )  THEN
    33053305       IF ( ocean )  THEN
    33063306          disturbance_level_t     = zu(nzt-3)
     
    34193419!-- In case of turbulent inflow calculate the index of the recycling plane
    34203420    IF ( turbulent_inflow )  THEN
    3421        IF ( recycling_width == 9999999.9 )  THEN
     3421       IF ( recycling_width == 9999999.9_wp )  THEN
    34223422!
    34233423!--       Set the default value for the width of the recycling domain
     
    34823482!-- Set time for the next user defined restart (time_restart is the
    34833483!-- internal parameter for steering restart events)
    3484     IF ( restart_time /= 9999999.9 )  THEN
     3484    IF ( restart_time /= 9999999.9_wp )  THEN
    34853485       IF ( restart_time > time_since_reference_point )  THEN
    34863486          time_restart = restart_time
     
    34903490!--    In case of a restart run, set internal parameter to default (no restart)
    34913491!--    if the NAMELIST-parameter restart_time is omitted
    3492        time_restart = 9999999.9
     3492       time_restart = 9999999.9_wp
    34933493    ENDIF
    34943494
  • palm/trunk/SOURCE/cpulog.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    165165#if defined( __lc ) || defined( __decalpha ) || defined( __nec )
    166166       CALL SYSTEM_CLOCK( count, count_rate )
    167        mtime = REAL( count ) / REAL( count_rate )
     167       mtime = REAL( count, KIND=wp ) / REAL( count_rate, KIND=wp )
    168168#elif defined( __ibm )
    169        mtime = IRTC( ) * 1E-9
     169       mtime = IRTC( ) * 1E-9_wp
    170170#else
    171171       message_string = 'no time measurement defined on this host'
     
    322322!--       Get total time in order to calculate CPU-time per gridpoint and timestep
    323323          IF ( nr_timesteps_this_run /= 0 )  THEN
    324              average_cputime = log_point(1)%sum / REAL( (nx+1) * (ny+1) * nz ) / &
    325                                REAL( nr_timesteps_this_run ) * 1E6  ! in micro-sec
     324             average_cputime = log_point(1)%sum / REAL( (nx+1) * (ny+1) * nz, KIND=wp ) / &
     325                               REAL( nr_timesteps_this_run, KIND=wp ) * 1E6_wp  ! in micro-sec
    326326          ELSE
    327327             average_cputime = -1.0
  • palm/trunk/SOURCE/data_output_profiles.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    102102    ELSE
    103103       IF ( average_count_pr > 0 )  THEN
    104           hom_sum = hom_sum / REAL( average_count_pr )
     104          hom_sum = hom_sum / REAL( average_count_pr, KIND=wp )
    105105       ELSE
    106106!
  • palm/trunk/SOURCE/data_output_spectra.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    113113       ENDIF
    114114       IF ( average_count_sp /= 1 )  THEN
    115           spectrum_x = spectrum_x / REAL( average_count_sp )
    116           spectrum_y = spectrum_y / REAL( average_count_sp )
     115          spectrum_x = spectrum_x / REAL( average_count_sp, KIND=wp )
     116          spectrum_y = spectrum_y / REAL( average_count_sp, KIND=wp )
    117117          average_count_sp = 0
    118118       ENDIF
  • palm/trunk/SOURCE/diffusion_w.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    132132!
    133133!--             Interpolate eddy diffusivities on staggered gridpoints
    134                 kmxp = 0.25 * &
     134                kmxp = 0.25_wp *                                              &
    135135                       ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    136                 kmxm = 0.25 * &
     136                kmxm = 0.25_wp *                                              &
    137137                       ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    138                 kmyp = 0.25 * &
     138                kmyp = 0.25_wp *                                              &
    139139                       ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    140                 kmym = 0.25 * &
     140                kmym = 0.25_wp *                                              &
    141141                       ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    142142
     
    152152                      &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)     &
    153153                      &   ) * ddy                                              &
    154                       & + 2.0 * (                                              &
     154                      & + 2.0_wp * (                                           &
    155155                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
    156156                      & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k)  &
     
    165165!
    166166!--                Interpolate eddy diffusivities on staggered gridpoints
    167                    kmxp = 0.25 * &
     167                   kmxp = 0.25_wp *                                            &
    168168                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    169                    kmxm = 0.25 * &
     169                   kmxm = 0.25_wp *                                            &
    170170                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    171                    kmyp = 0.25 * &
     171                   kmyp = 0.25_wp *                                            &
    172172                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    173                    kmym = 0.25 * &
     173                   kmym = 0.25_wp *                                            &
    174174                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    175175
     
    262262!
    263263!--                Interpolate eddy diffusivities on staggered gridpoints
    264                    kmxp = 0.25 * &
     264                   kmxp = 0.25_wp *                                            &
    265265                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    266                    kmxm = 0.25 * &
     266                   kmxm = 0.25_wp *                                            &
    267267                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    268                    kmyp = 0.25 * &
     268                   kmyp = 0.25_wp *                                            &
    269269                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    270                    kmym = 0.25 * &
     270                   kmym = 0.25_wp *                                            &
    271271                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    272272
     
    282282                         &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)    &
    283283                         &   ) * ddy                                             &
    284                          & + 2.0 * (                                             &
     284                         & + 2.0_wp * (                                          &
    285285                         &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
    286286                         & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k) &
     
    297297!
    298298!--                Interpolate eddy diffusivities on staggered gridpoints
    299                    kmxp = 0.25 * &
     299                   kmxp = 0.25_wp *                                            &
    300300                          ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    301                    kmxm = 0.25 * &
     301                   kmxm = 0.25_wp *                                            &
    302302                          ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    303                    kmyp = 0.25 * &
     303                   kmyp = 0.25_wp *                                            &
    304304                          ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    305                    kmym = 0.25 * &
     305                   kmym = 0.25_wp *                                            &
    306306                          ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    307307
     
    378378!
    379379!--       Interpolate eddy diffusivities on staggered gridpoints
    380           kmxp = 0.25 * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    381           kmxm = 0.25 * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    382           kmyp = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    383           kmym = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
     380          kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
     381          kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
     382          kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
     383          kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    384384
    385385          tend(k,j,i) = tend(k,j,i)                                            &
     
    394394                      &   - kmym * ( v(k+1,j,i) - v(k,j,i)   ) * ddzu(k+1)     &
    395395                      &   ) * ddy                                              &
    396                       & + 2.0 * (                                              &
     396                      & + 2.0_wp * (                                           &
    397397                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
    398398                      & - km(k,j,i)   * ( w(k,j,i)   - w(k-1,j,i) ) * ddzw(k)  &
     
    423423!
    424424!--          Interpolate eddy diffusivities on staggered gridpoints
    425              kmxp = 0.25 * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
    426              kmxm = 0.25 * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
    427              kmyp = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
    428              kmym = 0.25 * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
     425             kmxp = 0.25_wp * ( km(k,j,i)+km(k,j,i+1)+km(k+1,j,i)+km(k+1,j,i+1) )
     426             kmxm = 0.25_wp * ( km(k,j,i)+km(k,j,i-1)+km(k+1,j,i)+km(k+1,j,i-1) )
     427             kmyp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j+1,i)+km(k+1,j+1,i) )
     428             kmym = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )
    429429
    430430             tend(k,j,i) = tend(k,j,i)                                         &
  • palm/trunk/SOURCE/diffusivities.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    141141                IF ( dvar_dz > 0.0 ) THEN
    142142                   IF ( use_single_reference_value )  THEN
    143                       l_stable = 0.76 * sqrt_e /                               &
    144                                  SQRT( g / var_reference * dvar_dz ) + 1E-5
     143                      l_stable = 0.76_wp * sqrt_e /                            &
     144                                 SQRT( g / var_reference * dvar_dz ) + 1E-5_wp
    145145                   ELSE
    146                       l_stable = 0.76 * sqrt_e /                               &
    147                                  SQRT( g / var(k,j,i) * dvar_dz ) + 1E-5
     146                      l_stable = 0.76_wp * sqrt_e /                            &
     147                                 SQRT( g / var(k,j,i) * dvar_dz ) + 1E-5_wp
    148148                   ENDIF
    149149                ELSE
     
    162162      !
    163163      !--       Compute diffusion coefficients for momentum and heat
    164                 km(k,j,i) = 0.1 * l * sqrt_e
    165                 kh(k,j,i) = ( 1.0 + 2.0 * l / ll ) * km(k,j,i)
     164                km(k,j,i) = 0.1_wp * l * sqrt_e
     165                kh(k,j,i) = ( 1.0_wp + 2.0_wp * l / ll ) * km(k,j,i)
    166166
    167167             ENDIF
  • palm/trunk/SOURCE/disturb_field.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    9898          DO  j = dist_nys(dist_range), dist_nyn(dist_range)
    9999             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
    100                 randomnumber = 3.0 * disturbance_amplitude *                   &
    101                                ( random_function( iran ) - 0.5 )
     100                randomnumber = 3.0_wp * disturbance_amplitude *                &
     101                               ( random_function( iran ) - 0.5_wp )
    102102                IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.         &
    103103                     nyn >= j )                                                &
     
    113113             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
    114114#if defined( __nec )
    115                 randomnumber = 3.0 * disturbance_amplitude *                   &
    116                                ( RANDOM( 0 ) - 0.5 )
     115                randomnumber = 3.0_wp * disturbance_amplitude *                &
     116                               ( RANDOM( 0 ) - 0.5_wp )
    117117#else
    118118                CALL RANDOM_NUMBER( randomnumber )
    119                 randomnumber = 3.0 * disturbance_amplitude *                   &
    120                                 ( randomnumber - 0.5 )
     119                randomnumber = 3.0_wp * disturbance_amplitude *                &
     120                                ( randomnumber - 0.5_wp )
    121121#endif
    122122                IF ( nxl <= i .AND. nxr >= i .AND. nys <= j .AND. nyn >= j )   &
     
    145145             dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1)                  &
    146146                            + dist1(k,j+1,i) + dist1(k+1,j,i)                  &
    147                             ) / 12.0
     147                            ) / 12.0_wp
    148148          ENDDO
    149149          DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
    150150              dist2(k,j,i) = dist2(k,j,i) + ( dist1(k,j-1,i) + dist1(k-1,j,i)  &
    151                             + 6.0 * dist1(k,j,i)                               &
    152                             ) / 12.0
     151                            + 6.0_wp * dist1(k,j,i)                            &
     152                            ) / 12.0_wp
    153153          ENDDO
    154154        ENDDO
     
    165165             dist1(k,j,i) = ( dist2(k,j,i-1) + dist2(k,j,i+1) + dist2(k,j-1,i) &
    166166                            + dist2(k,j+1,i) + dist2(k+1,j,i) + dist2(k-1,j,i) &
    167                             + 6.0 * dist2(k,j,i)                               &
    168                             ) / 12.0
     167                            + 6.0_wp * dist2(k,j,i)                            &
     168                            ) / 12.0_wp
    169169          ENDDO
    170170       ENDDO
  • palm/trunk/SOURCE/disturb_heatflux.f90

    r1321 r1322  
    1919!
    2020! Current revisions:
    21 ! -----------------
    22 !
     21! ------------------
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    7979    DO  i = 0, nx
    8080       DO  j = 0, ny
    81           randomnumber = random_gauss( iran, 5.0 )
     81          randomnumber = random_gauss( iran, 5.0_wp )
    8282          IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.  nyn >= j )   &
    8383          THEN
  • palm/trunk/SOURCE/fft_xy.f90

    r1321 r1322  
    1616!
    1717! Copyright 1997-2014 Leibniz Universitaet Hannover
    18 !--------------------------------------------------------------------------------!
     18!------------------------------------------------------------------------------!
    1919!
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    413413
    414414                   DO  i = 0, (nx+1)/2
    415                       ar(i,j,k) = REAL( cwork(i) )
     415                      ar(i,j,k) = REAL( cwork(i), KIND=wp )
    416416                   ENDDO
    417417                   DO  i = 1, (nx+1)/2 - 1
     
    441441
    442442                   DO  i = 0, nx
    443                       ar(i,j,k) = REAL( cwork(i) )
     443                      ar(i,j,k) = REAL( cwork(i), KIND=wp )
    444444                   ENDDO
    445445
     
    519519
    520520                      DO  i = 0, (nx+1)/2
    521                          ar_2d(i,j) = REAL( x_out(i) ) / ( nx+1 )
     521                         ar_2d(i,j) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
    522522                      ENDDO
    523523                      DO  i = 1, (nx+1)/2 - 1
     
    528528
    529529                      DO  i = 0, (nx+1)/2
    530                          ar(i,j,k) = REAL( x_out(i) ) / ( nx+1 )
     530                         ar(i,j,k) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
    531531                      ENDDO
    532532                      DO  i = 1, (nx+1)/2 - 1
     
    689689
    690690                   DO  i = 0, (nx+1)/2
    691                       ar(i,j,k)      = REAL( ar_tmp(i,j,k) )  * dnx
     691                      ar(i,j,k)      = REAL( ar_tmp(i,j,k), KIND=wp )  * dnx
    692692                   ENDDO
    693693
     
    794794             CALL FFTN( cwork, ishape )
    795795             DO  i = 0, (nx+1)/2
    796                 ar(i) = REAL( cwork(i) )
     796                ar(i) = REAL( cwork(i), KIND=wp )
    797797             ENDDO
    798798             DO  i = 1, (nx+1)/2 - 1
     
    813813
    814814             DO  i = 0, nx
    815                 ar(i) = REAL( cwork(i) )
     815                ar(i) = REAL( cwork(i), KIND=wp )
    816816             ENDDO
    817817
     
    862862
    863863             DO  i = 0, (nx+1)/2
    864                 ar(i) = REAL( x_out(i) ) / ( nx+1 )
     864                ar(i) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
    865865             ENDDO
    866866             DO  i = 1, (nx+1)/2 - 1
     
    10561056
    10571057                   DO  j = 0, (ny+1)/2
    1058                       ar_tr(j,i,k) = REAL( cwork(j) )
     1058                      ar_tr(j,i,k) = REAL( cwork(j), KIND=wp )
    10591059                   ENDDO
    10601060                   DO  j = 1, (ny+1)/2 - 1
     
    10841084
    10851085                   DO  j = 0, ny
    1086                       ar(j,i,k) = REAL( cwork(j) )
     1086                      ar(j,i,k) = REAL( cwork(j), KIND=wp )
    10871087                   ENDDO
    10881088
     
    11601160
    11611161                   DO  j = 0, (ny+1)/2
    1162                       ar_tr(j,i,k) = REAL( y_out(j) ) / (ny+1)
     1162                      ar_tr(j,i,k) = REAL( y_out(j), KIND=wp ) / (ny+1)
    11631163                   ENDDO
    11641164                   DO  j = 1, (ny+1)/2 - 1
     
    13051305
    13061306                   DO  j = 0, (ny+1)/2
    1307                       ar(j,i,k)      = REAL( ar_tmp(j,i,k) )  * dny
     1307                      ar(j,i,k)      = REAL( ar_tmp(j,i,k), KIND=wp )  * dny
    13081308                   ENDDO
    13091309
     
    14121412
    14131413             DO  j = 0, (ny+1)/2
    1414                 ar(j) = REAL( cwork(j) )
     1414                ar(j) = REAL( cwork(j), KIND=wp )
    14151415             ENDDO
    14161416             DO  j = 1, (ny+1)/2 - 1
     
    14311431
    14321432             DO  j = 0, ny
    1433                 ar(j) = REAL( cwork(j) )
     1433                ar(j) = REAL( cwork(j), KIND=wp )
    14341434             ENDDO
    14351435
     
    14801480
    14811481             DO  j = 0, (ny+1)/2
    1482                 ar(j) = REAL( y_out(j) ) / (ny+1)
     1482                ar(j) = REAL( y_out(j), KIND=wp ) / (ny+1)
    14831483             ENDDO
    14841484             DO  j = 1, (ny+1)/2 - 1
     
    16741674             DO  k = 1, nz
    16751675                DO  i = 0, (nx+1)/2
    1676                    ar(i,k) = REAL( work(i+1,k) )
     1676                   ar(i,k) = REAL( work(i+1,k), KIND=wp )
    16771677                ENDDO
    16781678                DO  i = 1, (nx+1)/2 - 1
     
    18171817             DO  k = 1, nz
    18181818                DO  j = 0, (ny+1)/2
    1819                    ar(j,k) = REAL( work(j+1,k) )
     1819                   ar(j,k) = REAL( work(j+1,k), KIND=wp )
    18201820                ENDDO
    18211821                DO  j = 1, (ny+1)/2 - 1
  • palm/trunk/SOURCE/flow_statistics.f90

    r1321 r1322  
    2121! Current revisions:
    2222! -----------------
    23 !
     23! REAL constants defined as wp-kind
    2424!
    2525! Former revisions:
     
    11931193          dptdz(k) = ( hom(k,1,4,sr) - hom(k-1,1,4,sr) ) * ddzu(k)
    11941194       ENDDO
    1195        dptdz_threshold = 0.2 / 100.0
     1195       dptdz_threshold = 0.2_wp / 100.0_wp
    11961196
    11971197       IF ( ocean )  THEN
     
    28592859          dptdz(k) = ( hom(k,1,4,sr) - hom(k-1,1,4,sr) ) * ddzu(k)
    28602860       ENDDO
    2861        dptdz_threshold = 0.2 / 100.0
     2861       dptdz_threshold = 0.2_wp / 100.0_wp
    28622862
    28632863       IF ( ocean )  THEN
  • palm/trunk/SOURCE/header.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute,
     23! some REAL constants defined as wp-kind
    2324!
    2425! Former revisions:
     
    454455    WRITE ( io, 203 )  simulated_time_at_begin, end_time
    455456
    456     IF ( time_restart /= 9999999.9  .AND. &
     457    IF ( time_restart /= 9999999.9_wp  .AND. &
    457458         simulated_time_at_begin == simulated_time )  THEN
    458        IF ( dt_restart == 9999999.9 )  THEN
     459       IF ( dt_restart == 9999999.9_wp )  THEN
    459460          WRITE ( io, 204 )  ' Restart at:       ',time_restart
    460461       ELSE
     
    472473                                              simulated_time_at_begin )
    473474       ENDIF
    474        WRITE ( io, 206 )  simulated_time, log_point_s(10)%sum, &
    475                           log_point_s(10)%sum / REAL( i ),    &
     475       WRITE ( io, 206 )  simulated_time, log_point_s(10)%sum,      &
     476                          log_point_s(10)%sum / REAL( i, KIND=wp ), &
    476477                          cpuseconds_per_simulated_second
    477        IF ( time_restart /= 9999999.9  .AND.  time_restart < end_time )  THEN
    478           IF ( dt_restart == 9999999.9 )  THEN
     478       IF ( time_restart /= 9999999.9_wp  .AND.  time_restart < end_time )  THEN
     479          IF ( dt_restart == 9999999.9_wp )  THEN
    479480             WRITE ( io, 204 )  ' Next restart at:     ',time_restart
    480481          ELSE
     
    524525          bh  = INT( building_height / dz )
    525526
    526           IF ( building_wall_left == 9999999.9 )  THEN
     527          IF ( building_wall_left == 9999999.9_wp )  THEN
    527528             building_wall_left = ( nx + 1 - blx ) / 2 * dx
    528529          ENDIF
     
    530531          bxr = bxl + blx
    531532
    532           IF ( building_wall_south == 9999999.9 )  THEN
     533          IF ( building_wall_south == 9999999.9_wp )  THEN
    533534             building_wall_south = ( ny + 1 - bly ) / 2 * dy
    534535          ENDIF
     
    541542       CASE ( 'single_street_canyon' )
    542543          ch  = NINT( canyon_height / dz )
    543           IF ( canyon_width_x /= 9999999.9 )  THEN
     544          IF ( canyon_width_x /= 9999999.9_wp )  THEN
    544545!
    545546!--          Street canyon in y direction
    546547             cwx = NINT( canyon_width_x / dx )
    547              IF ( canyon_wall_left == 9999999.9 )  THEN
     548             IF ( canyon_wall_left == 9999999.9_wp )  THEN
    548549                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
    549550             ENDIF
     
    552553             WRITE ( io, 272 )  'y', canyon_height, ch, 'u', cxl, cxr
    553554
    554           ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
     555          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
    555556!
    556557!--          Street canyon in x direction
    557558             cwy = NINT( canyon_width_y / dy )
    558              IF ( canyon_wall_south == 9999999.9 )  THEN
     559             IF ( canyon_wall_south == 9999999.9_wp )  THEN
    559560                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
    560561             ENDIF
     
    11881189!
    11891190!-- Timeseries
    1190     IF ( dt_dots /= 9999999.9 )  THEN
     1191    IF ( dt_dots /= 9999999.9_wp )  THEN
    11911192       WRITE ( io, 340 )
    11921193
     
    12021203!
    12031204!-- Dvrp-output
    1204     IF ( dt_dvrp /= 9999999.9 )  THEN
     1205    IF ( dt_dvrp /= 9999999.9_wp )  THEN
    12051206       WRITE ( io, 360 )  dt_dvrp, TRIM( dvrp_output ), TRIM( dvrp_host ), &
    12061207                          TRIM( dvrp_username ), TRIM( dvrp_directory )
     
    12541255!
    12551256!-- Spectra output
    1256     IF ( dt_dosp /= 9999999.9 )  THEN
     1257    IF ( dt_dosp /= 9999999.9_wp )  THEN
    12571258       WRITE ( io, 370 )
    12581259
     
    13771378!
    13781379!-- Initial wind profiles
    1379     IF ( u_profile(1) /= 9999999.9 )  WRITE ( io, 427 )
     1380    IF ( u_profile(1) /= 9999999.9_wp )  WRITE ( io, 427 )
    13801381
    13811382!
     
    16271628          ENDIF
    16281629       ENDIF
    1629        IF ( dt_write_particle_data /= 9999999.9 )  THEN
     1630       IF ( dt_write_particle_data /= 9999999.9_wp )  THEN
    16301631          WRITE ( io, 485 )  dt_write_particle_data
    16311632          output_format = ' '
     
    16411642          WRITE ( io, 344 )  output_format
    16421643       ENDIF
    1643        IF ( dt_dopts /= 9999999.9 )  WRITE ( io, 494 )  dt_dopts
     1644       IF ( dt_dopts /= 9999999.9_wp )  WRITE ( io, 494 )  dt_dopts
    16441645       IF ( write_particle_statistics )  WRITE ( io, 486 )
    16451646
     
    16471648
    16481649       DO  i = 1, number_of_particle_groups
    1649           IF ( i == 1  .AND.  density_ratio(i) == 9999999.9 )  THEN
     1650          IF ( i == 1  .AND.  density_ratio(i) == 9999999.9_wp )  THEN
    16501651             WRITE ( io, 490 )  i, 0.0
    16511652             WRITE ( io, 492 )
  • palm/trunk/SOURCE/init_1d_model.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    777777          energy = energy + 0.5 * ( u1d(k)**2 + v1d(k)**2 )
    778778       ENDDO
    779        energy = energy / REAL( nzt - nzb + 1 )
     779       energy = energy / REAL( nzt - nzb + 1, KIND=wp )
    780780
    781781       uv_total = SQRT( u1d(nzb+1)**2 + v1d(nzb+1)**2 )
  • palm/trunk/SOURCE/init_3d_model.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
    23 !
     22! REAL constants defined as wp-kind
    2423! module interfaces removed
    2524!
     
    16041603    IF ( TRIM(timestep_scheme) == 'runge-kutta-3' )  THEN      ! for RK3-method
    16051604
    1606        weight_substep(1) = 1./6.
    1607        weight_substep(2) = 3./10.
    1608        weight_substep(3) = 8./15.
    1609 
    1610        weight_pres(1)    = 1./3.
    1611        weight_pres(2)    = 5./12.
    1612        weight_pres(3)    = 1./4.
     1605       weight_substep(1) = 1._wp/6._wp
     1606       weight_substep(2) = 3._wp/10._wp
     1607       weight_substep(3) = 8._wp/15._wp
     1608
     1609       weight_pres(1)    = 1._wp/3._wp
     1610       weight_pres(2)    = 5._wp/12._wp
     1611       weight_pres(3)    = 1._wp/4._wp
    16131612
    16141613    ELSEIF ( TRIM(timestep_scheme) == 'runge-kutta-2' )  THEN  ! for RK2-method
    16151614
    1616        weight_substep(1) = 1./2.
    1617        weight_substep(2) = 1./2.
     1615       weight_substep(1) = 1._wp/2._wp
     1616       weight_substep(2) = 1._wp/2._wp
    16181617         
    1619        weight_pres(1)    = 1./2.
    1620        weight_pres(2)    = 1./2.       
     1618       weight_pres(1)    = 1._wp/2._wp
     1619       weight_pres(2)    = 1._wp/2._wp       
    16211620
    16221621    ELSE                                     ! for Euler-method
     
    16711670          DO  k = dp_level_ind_b+1, nzt
    16721671             dp_smooth_factor(k) = 0.5 * ( 1.0 + SIN( pi * &
    1673                   ( REAL( k - dp_level_ind_b ) /  &
    1674                     REAL( nzt - dp_level_ind_b ) - 0.5 ) ) )
     1672                  ( REAL( k - dp_level_ind_b, KIND=wp ) /  &
     1673                    REAL( nzt - dp_level_ind_b, KIND=wp ) - 0.5 ) ) )
    16751674          ENDDO
    16761675       ENDIF
     
    16861685       DO  i = nxl, nxr
    16871686          IF ( ( i * dx ) < pt_damping_width )  THEN
    1688              ptdf_x(i) = pt_damping_factor * ( SIN( pi * 0.5 *        &
    1689                          REAL( pt_damping_width - i * dx ) / (        &
    1690                          REAL( pt_damping_width )            ) ) )**2
     1687             ptdf_x(i) = pt_damping_factor * ( SIN( pi * 0.5 *              &
     1688                         REAL( pt_damping_width - i * dx, KIND=wp ) / (                 &
     1689                         REAL( pt_damping_width, KIND=wp )            ) ) )**2
    16911690          ENDIF
    16921691       ENDDO
     
    16941693       DO  i = nxl, nxr
    16951694          IF ( ( i * dx ) > ( nx * dx - pt_damping_width ) )  THEN
    1696              ptdf_x(i) = pt_damping_factor *                                      &
    1697                          SIN( pi * 0.5 * ( ( i - nx ) * dx + pt_damping_width ) / &
    1698                                          REAL( pt_damping_width ) )**2
     1695             ptdf_x(i) = pt_damping_factor *                                   &
     1696                         SIN( pi * 0.5 *                                       &
     1697                              ( ( i - nx ) * dx + pt_damping_width ) /         &
     1698                              REAL( pt_damping_width, KIND=wp ) )**2
    16991699          ENDIF
    17001700       ENDDO
     
    17021702       DO  j = nys, nyn
    17031703          IF ( ( j * dy ) > ( ny * dy - pt_damping_width ) )  THEN
    1704              ptdf_y(j) = pt_damping_factor *                                      &
    1705                          SIN( pi * 0.5 * ( ( j - ny ) * dy + pt_damping_width ) / &
    1706                                          REAL( pt_damping_width ) )**2
     1704             ptdf_y(j) = pt_damping_factor *                                   &
     1705                         SIN( pi * 0.5 *                                       &
     1706                              ( ( j - ny ) * dy + pt_damping_width ) /         &
     1707                              REAL( pt_damping_width, KIND=wp ) )**2
    17071708          ENDIF
    17081709       ENDDO
     
    17101711       DO  j = nys, nyn
    17111712          IF ( ( j * dy ) < pt_damping_width )  THEN
    1712              ptdf_y(j) = pt_damping_factor *                             &
    1713                          SIN( pi * 0.5 * ( pt_damping_width - j * dy ) / &
    1714                                          REAL( pt_damping_width ) )**2
     1713             ptdf_y(j) = pt_damping_factor *                                   &
     1714                         SIN( pi * 0.5 *                                       &
     1715                             ( pt_damping_width - j * dy ) /                   &
     1716                             REAL( pt_damping_width, KIND=wp ) )**2
    17151717          ENDIF
    17161718       ENDDO
  • palm/trunk/SOURCE/init_advec.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    8585       ALLOCATE( aex(intervals), bex(intervals), dex(intervals), eex(intervals) )
    8686
    87        delt  = 1.0 / REAL( intervals )
     87       delt  = 1.0_wp / REAL( intervals, KIND=wp )
    8888       sterm = delt * 0.5
    8989
  • palm/trunk/SOURCE/init_cloud_physics.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    9696    l_d_rv = l_v / r_v
    9797
    98     schmidt_p_1d3 = schmidt**( 1.0 / 3.0 )
     98    schmidt_p_1d3 = schmidt**( 1.0_wp / 3.0_wp )
    9999
    100100    pirho_l  = pi * rho_l / 6.0
     
    125125       ENDIF
    126126       hyp(k)    = surface_pressure * 100.0 * &
    127                    ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0/0.286)
     127                   ( (t_surface - g/cp * zu(k)) / t_surface )**(1.0_wp/0.286_wp)
    128128       pt_d_t(k) = ( 100000.0 / hyp(k) )**0.286
    129129       t_d_pt(k) = 1.0 / pt_d_t(k)
  • palm/trunk/SOURCE/init_dvrp.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    710710!--       Define a default colourtable for particles
    711711          DO  i = 1, 11
    712              interval_values_dvrp_prt(1,i) = i - 1.0
    713              interval_values_dvrp_prt(2,i) = REAL( i )
    714              interval_h_dvrp_prt(:,i) = 270.0 - ( i - 1.0 ) * 9.0
     712             interval_values_dvrp_prt(1,i) = i - 1.0_wp
     713             interval_values_dvrp_prt(2,i) = REAL( i, KIND=wp )
     714             interval_h_dvrp_prt(:,i) = 270.0_wp - ( i - 1.0_wp ) * 9.0_wp
    715715          ENDDO
    716716
    717717          DO  i = 12, 22
    718              interval_values_dvrp_prt(1,i) = i - 1.0
    719              interval_values_dvrp_prt(2,i) = REAL( i )
    720              interval_h_dvrp_prt(:,i) = 70.0 - ( i - 12.0 ) * 9.5
     718             interval_values_dvrp_prt(1,i) = i - 1.0_wp
     719             interval_values_dvrp_prt(2,i) = REAL( i, KIND=wp )
     720             interval_h_dvrp_prt(:,i) = 70.0_wp - ( i - 12.0_wp ) * 9.5_wp
    721721          ENDDO
    722722
  • palm/trunk/SOURCE/init_grid.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    377377!-- Compute the grid-dependent mixing length.
    378378    DO  k = 1, nzt
    379        l_grid(k)  = ( dx * dy * dzw(k) )**0.33333333333333
     379       l_grid(k)  = ( dx * dy * dzw(k) )**0.33333333333333_wp
    380380    ENDDO
    381381
     
    516516          bh  = NINT( building_height / dz )
    517517
    518           IF ( building_wall_left == 9999999.9 )  THEN
     518          IF ( building_wall_left == 9999999.9_wp )  THEN
    519519             building_wall_left = ( nx + 1 - blx ) / 2 * dx
    520520          ENDIF
     
    522522          bxr = bxl + blx
    523523
    524           IF ( building_wall_south == 9999999.9 )  THEN
     524          IF ( building_wall_south == 9999999.9_wp )  THEN
    525525             building_wall_south = ( ny + 1 - bly ) / 2 * dy
    526526          ENDIF
     
    547547!--       Single quasi-2D street canyon of infinite length in x or y direction.
    548548!--       The canyon is centered in the other direction by default.
    549           IF ( canyon_width_x /= 9999999.9 )  THEN
     549          IF ( canyon_width_x /= 9999999.9_wp )  THEN
    550550!
    551551!--          Street canyon in y direction
    552552             cwx = NINT( canyon_width_x / dx )
    553              IF ( canyon_wall_left == 9999999.9 )  THEN
     553             IF ( canyon_wall_left == 9999999.9_wp )  THEN
    554554                canyon_wall_left = ( nx + 1 - cwx ) / 2 * dx
    555555             ENDIF
     
    557557             cxr = cxl + cwx
    558558
    559           ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
     559          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
    560560!
    561561!--          Street canyon in x direction
    562562             cwy = NINT( canyon_width_y / dy )
    563              IF ( canyon_wall_south == 9999999.9 )  THEN
     563             IF ( canyon_wall_south == 9999999.9_wp )  THEN
    564564                canyon_wall_south = ( ny + 1 - cwy ) / 2 * dy
    565565             ENDIF
     
    578578!
    579579!--       Street canyon size has to meet some requirements
    580           IF ( canyon_width_x /= 9999999.9 )  THEN
     580          IF ( canyon_width_x /= 9999999.9_wp )  THEN
    581581             IF ( ( cxl < 1 ) .OR. ( cxr > nx-1 ) .OR. ( cwx < 3 ) .OR.  &
    582582               ( ch < 3 ) )  THEN
     
    587587                CALL message( 'init_grid', 'PA0205', 1, 2, 0, 6, 0 )
    588588             ENDIF
    589           ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
     589          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
    590590             IF ( ( cys < 1 ) .OR. ( cyn > ny-1 ) .OR. ( cwy < 3 ) .OR.  &
    591591               ( ch < 3 ) )  THEN
     
    597597             ENDIF
    598598          ENDIF
    599           IF ( canyon_width_x /= 9999999.9 .AND. canyon_width_y /= 9999999.9 ) &
     599          IF ( canyon_width_x /= 9999999.9_wp .AND. canyon_width_y /= 9999999.9_wp ) &
    600600               THEN
    601601             message_string = 'inconsistent canyon parameters:' //     & 
     
    606606
    607607          nzb_local = ch
    608           IF ( canyon_width_x /= 9999999.9 )  THEN
     608          IF ( canyon_width_x /= 9999999.9_wp )  THEN
    609609             nzb_local(:,cxl+1:cxr-1) = 0
    610           ELSEIF ( canyon_width_y /= 9999999.9 )  THEN
     610          ELSEIF ( canyon_width_y /= 9999999.9_wp )  THEN
    611611             nzb_local(cys+1:cyn-1,:) = 0
    612612          ENDIF
  • palm/trunk/SOURCE/init_ocean.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    8788!-- and the reference density (used later in buoyancy term)
    8889!-- First step: Calculate pressure using reference density
    89     hyp(nzt+1) = surface_pressure * 100.0
     90    hyp(nzt+1) = surface_pressure * 100.0_wp
    9091
    9192    hyp(nzt)      = hyp(nzt+1) + rho_surface * g * 0.5 * dzu(nzt+1)
  • palm/trunk/SOURCE/init_pegrid.f90

    r1321 r1322  
    1919!
    2020! Current revisions:
    21 ! -----------------
     21! ------------------
     22! REAL functions provided with KIND-attribute
    2223!
    2324! Former revisions:
     
    174175!
    175176!--    Automatic determination of the topology
    176        numproc_sqr = SQRT( REAL( numprocs ) )
     177       numproc_sqr = SQRT( REAL( numprocs, KIND=wp ) )
    177178       pdims(1)    = MAX( numproc_sqr , 1 )
    178179       DO  WHILE ( MOD( numprocs , pdims(1) ) /= 0 )
  • palm/trunk/SOURCE/init_pt_anomaly.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    9394             radius = SQRT( x**2 + y**2 + z**2 )
    9495             IF ( radius <= rc )  THEN
    95                 betrag = 5.0 * EXP( -( radius * 0.001 / 2.0 )**2 )
     96                betrag = 5.0 * EXP( -( radius * 0.001_wp / 2.0_wp )**2 )
    9697             ELSE
    9798                betrag = 0.0
  • palm/trunk/SOURCE/init_rankine.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    111112          radius = SQRT( x**2 + y**2 )
    112113          IF ( radius <= 2.0 * rc )  THEN
    113              betrag = radius / ( 2.0 * rc ) * 0.08
     114             betrag = radius / ( 2.0 * rc ) * 0.08_wp
    114115          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    115              betrag = 0.08 * EXP( -( radius - 2.0 * rc ) / 2.0 )
     116             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    116117          ELSE
    117118             betrag = 0.0
     
    120121          IF ( x == 0.0 )  THEN
    121122             IF ( y > 0.0 )  THEN
    122                 alpha = pi / 2.0
     123                alpha = pi / 2.0_wp
    123124             ELSEIF ( y < 0.0 )  THEN
    124                 alpha = 3.0 * pi / 2.0
     125                alpha = 3.0 * pi / 2.0_wp
    125126             ENDIF
    126127          ELSE
     
    129130             ELSE
    130131                IF ( y < 0.0 )  THEN
    131                    alpha = ATAN( y / x ) + 2.0 * pi
     132                   alpha = ATAN( y / x ) + 2.0_wp * pi
    132133                ELSE
    133134                   alpha = ATAN( y / x )
     
    152153          radius = SQRT( x**2 + y**2 )
    153154          IF ( radius <= 2.0 * rc )  THEN
    154              betrag = radius / ( 2.0 * rc ) * 0.08
     155             betrag = radius / ( 2.0_wp * rc ) * 0.08_wp
    155156          ELSEIF ( radius > 2.0 * rc  .AND.  radius < 8.0 * rc )  THEN
    156              betrag = 0.08 * EXP( -( radius - 2.0 * rc ) / 2.0 )
     157             betrag = 0.08 * EXP( -( radius - 2.0_wp * rc ) / 2.0_wp )
    157158          ELSE
    158159             betrag = 0.0
     
    161162          IF ( x == 0.0 )  THEN
    162163             IF ( y > 0.0 )  THEN
    163                 alpha = pi / 2.0
     164                alpha = pi / 2.0_wp
    164165             ELSEIF ( y < 0.0 )  THEN
    165                 alpha = 3.0 * pi / 2.0
     166                alpha = 3.0 * pi / 2.0_wp
    166167             ENDIF
    167168          ELSE
     
    170171             ELSE
    171172                IF ( y < 0.0 )  THEN
    172                    alpha = ATAN( y / x ) + 2.0 * pi
     173                   alpha = ATAN( y / x ) + 2.0_wp * pi
    173174                ELSE
    174175                   alpha = ATAN( y / x )
  • palm/trunk/SOURCE/init_slope.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    109110!
    110111!--       Compute temperatures in the rotated coordinate system
    111           alpha    = alpha + alpha_surface / 180.0 * pi
     112          alpha    = alpha + alpha_surface / 180.0_wp * pi
    112113          pt_value = pt_surface + radius * SIN( alpha ) * &
    113                                   pt_vertical_gradient(1) / 100.0
     114                                  pt_vertical_gradient(1) / 100.0_wp
    114115          pt_slope_ref(k,i) = pt_value
    115116       ENDDO               
     
    120121!-- used for the cyclic boundary in x-direction
    121122    pt_slope_offset = (nx+1) * dx * sin_alpha_surface * &
    122                       pt_vertical_gradient(1) / 100.0
     123                      pt_vertical_gradient(1) / 100.0_wp
    123124
    124125
  • palm/trunk/SOURCE/local_tremain.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    6970#if defined( __ibm )
    7071
    71     actual_wallclock_time = IRTC( ) * 1E-9
     72    actual_wallclock_time = IRTC( ) * 1E-9_wp
    7273    remaining_time = maximum_cpu_time_allowed - &
    7374                     ( actual_wallclock_time - initial_wallclock_time )
     
    7677
    7778    CALL SYSTEM_CLOCK( count, count_rate )
    78     actual_wallclock_time = REAL( count ) / REAL( count_rate )
     79    actual_wallclock_time = REAL( count, KIND=wp ) / REAL( count_rate, KIND=wp )
    7980    remaining_time = maximum_cpu_time_allowed - &
    8081                     ( actual_wallclock_time - initial_wallclock_time )
  • palm/trunk/SOURCE/local_tremain_ini.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    6768#if defined( __ibm )
    6869
    69     initial_wallclock_time = IRTC( ) * 1E-9
     70    initial_wallclock_time = IRTC( ) * 1E-9_wp
    7071
    7172#elif defined( __lc )
    7273
    7374    CALL SYSTEM_CLOCK( count, count_rate )
    74     initial_wallclock_time = REAL( count ) / REAL( count_rate )
     75    initial_wallclock_time = REAL( count, KIND=wp ) / REAL( count_rate, KIND=wp )
    7576
    7677#else
  • palm/trunk/SOURCE/lpm_advec.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    186187!
    187188!--          Determine the sublayer. Further used as index.
    188              height_p     = ( particles(n)%z - z0_av_global )           &
    189                                   * REAL( number_of_sublayers )        &
     189             height_p     = ( particles(n)%z - z0_av_global )            &
     190                                  * REAL( number_of_sublayers, KIND=wp ) &
    190191                                  * d_z_p_z0
    191192!
     
    210211             us_int = 0.5 * ( us(j,i) + us(j,i-1) ) 
    211212
    212              u_int  = -usws(j,i) / ( us_int * kappa + 1E-10 )           &
     213             u_int  = -usws(j,i) / ( us_int * kappa + 1E-10_wp )           &
    213214                      * log_z_z0_int
    214215
     
    277278              us_int = 0.5 * ( us(j,i) + us(j-1,i) )   
    278279
    279               v_int  = -vsws(j,i) / ( us_int * kappa + 1E-10 )             &
     280              v_int  = -vsws(j,i) / ( us_int * kappa + 1E-10_wp )             &
    280281                       * log_z_z0_int
    281282
     
    10931094          ENDIF
    10941095
    1095           vv_int = ( 1.0 / 3.0 ) * ( aa + bb + cc )
    1096 
    1097           fs_int = ( 2.0 / 3.0 ) * e_mean_int / &
    1098                       ( vv_int + ( 2.0 / 3.0 ) * e_mean_int )
     1096          vv_int = ( 1.0_wp / 3.0_wp ) * ( aa + bb + cc )
     1097
     1098          fs_int = ( 2.0_wp / 3.0_wp ) * e_mean_int / &
     1099                      ( vv_int + ( 2.0_wp / 3.0_wp ) * e_mean_int )
    10991100
    11001101!
     
    11261127!--          from becoming unrealistically large.
    11271128             particles(n)%rvar1 = SQRT( 2.0 * sgs_wfu_part * e_int ) * &
    1128                                         ( random_gauss( iran_part, 5.0 ) - 1.0 )
     1129                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
    11291130             particles(n)%rvar2 = SQRT( 2.0 * sgs_wfv_part * e_int ) * &
    1130                                         ( random_gauss( iran_part, 5.0 ) - 1.0 )
     1131                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
    11311132             particles(n)%rvar3 = SQRT( 2.0 * sgs_wfw_part * e_int ) * &
    1132                                         ( random_gauss( iran_part, 5.0 ) - 1.0 )
     1133                                        ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp )
    11331134
    11341135          ELSE
     
    11681169                                  ) * dt_particle / 2.0          +             &
    11691170                                  SQRT( fs_int * c_0 * diss_int ) *            &
    1170                                   ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
     1171                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
    11711172                                  SQRT( dt_particle )
    11721173
     
    11771178                                    particles(n)%rvar2 /                       &
    11781179                                    ( 2.0 * sgs_wfv_part * e_int ) + de_dy_int &
    1179                                   ) * dt_particle / 2.0          +             &
     1180                                  ) * dt_particle / 2.0_wp          +             &
    11801181                                  SQRT( fs_int * c_0 * diss_int ) *            &
    1181                                   ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
     1182                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
    11821183                                  SQRT( dt_particle )
    11831184
     
    11881189                                    particles(n)%rvar3 /                       &
    11891190                                    ( 2.0 * sgs_wfw_part * e_int ) + de_dz_int &
    1190                                   ) * dt_particle / 2.0          +             &
     1191                                  ) * dt_particle / 2.0_wp         +             &
    11911192                                  SQRT( fs_int * c_0 * diss_int ) *            &
    1192                                   ( random_gauss( iran_part, 5.0 ) - 1.0 ) *   &
     1193                                  ( random_gauss( iran_part, 5.0_wp ) - 1.0_wp ) *   &
    11931194                                  SQRT( dt_particle )
    11941195
     
    12521253                                ( v_int - particles(n)%speed_y )**2 +   &
    12531254                                ( w_int - particles(n)%speed_z )**2 ) / &
    1254                                          molecular_viscosity )**0.687   &
     1255                                         molecular_viscosity )**0.687_wp   &
    12551256                        )
    12561257             exp_term = EXP( -exp_arg * dt_particle )
  • palm/trunk/SOURCE/lpm_collision_kernels.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
     22! REAL constants defined as wp_kind
    2223!
    2324! Former revisions:
     
    167168!--       Calculate the radius class bounds with logarithmic distances
    168169!--       in the interval [1.0E-6, 2.0E-4] m
    169           rclass_lbound = LOG( 1.0E-6 )
    170           rclass_ubound = LOG( 2.0E-4 )
     170          rclass_lbound = LOG( 1.0E-6_wp )
     171          rclass_ubound = LOG( 2.0E-4_wp )
    171172          radclass(1)   = 1.0E-6
    172173          DO  i = 2, radius_classes
     
    182183!--       Set the class bounds for dissipation in interval [0.0, 0.1] m**2/s**3
    183184          DO  i = 1, dissipation_classes
    184              epsclass(i) = 0.1 * REAL( i ) / dissipation_classes
     185             epsclass(i) = 0.1 * REAL( i, KIND=wp ) / dissipation_classes
    185186!             IF ( myid == 0 )  THEN
    186187!                PRINT*, 'i=', i, ' eps = ', epsclass(i)
     
    197198
    198199             epsilon = epsclass(k)
    199              urms    = 2.02 * ( epsilon / 0.04 )**( 1.0 / 3.0 )
     200             urms    = 2.02 * ( epsilon / 0.04_wp )**( 1.0_wp / 3.0_wp )
    200201
    201202             CALL turbsd
     
    313314          epsilon = 0.0
    314315       ENDIF
    315        urms    = 2.02 * ( epsilon / 0.04 )**( 0.33333333333 )
     316       urms    = 2.02 * ( epsilon / 0.04_wp )**( 0.33333333333_wp )
    316317
    317318       IF ( wang_kernel  .AND.  epsilon > 1.0E-7 )  THEN
     
    431432       ENDIF
    432433
    433        lambda    = urms * SQRT( 15.0 * molecular_viscosity / epsilon )    ! in m
    434        lambda_re = urms**2 * SQRT( 15.0 / epsilon / molecular_viscosity )
     434       lambda    = urms * SQRT( 15.0_wp * molecular_viscosity / epsilon ) ! in m
     435       lambda_re = urms**2 * SQRT( 15.0_wp / epsilon / molecular_viscosity )
    435436       tl        = urms**2 / epsilon                       ! in s
    436437       lf        = 0.5 * urms**3 / epsilon                 ! in m
    437438       tauk      = SQRT( molecular_viscosity / epsilon )                  ! in s
    438        eta       = ( molecular_viscosity**3 / epsilon )**0.25             ! in m
     439       eta       = ( molecular_viscosity**3 / epsilon )**0.25_wp          ! in m
    439440       vk        = eta / tauk
    440441
    441442       ao = ( 11.0 + 7.0 * lambda_re ) / ( 205.0 + lambda_re )
    442        tt = SQRT( 2.0 * lambda_re / ( SQRT( 15.0 ) * ao ) ) * tauk   ! in s
     443       tt = SQRT( 2.0 * lambda_re / ( SQRT( 15.0_wp ) * ao ) ) * tauk   ! in s
    443444
    444445       CALL fallg    ! gives winf in m/s
     
    452453!--    Calculate wr (from Aayala 2008b, page 38f)
    453454       z   = tt / tl
    454        be  = SQRT( 2.0 ) * lambda / lf
     455       be  = SQRT( 2.0_wp ) * lambda / lf
    455456       bbb = SQRT( 1.0 - 2.0 * be**2 )
    456457       d1  = ( 1.0 + bbb ) / ( 2.0 * bbb )
     
    504505             wrtur2xy = vrms1xy**2 + vrms2xy**2 - 2.0 * v1v2xy   ! in m**2/s**2
    505506             IF ( wrtur2xy < 0.0 )  wrtur2xy = 0.0
    506              wrgrav2  = pi / 8.0 * ( winf(j) - winf(i) )**2
    507              wrfin    = SQRT( ( 2.0 / pi ) * ( wrtur2xy + wrgrav2) )  ! in m/s
     507             wrgrav2  = pi / 8.0_wp * ( winf(j) - winf(i) )**2
     508             wrfin    = SQRT( ( 2.0_wp / pi ) * ( wrtur2xy + wrgrav2) ) ! in m/s
    508509
    509510!
     
    518519                   5.3406 * sst
    519520             IF ( xx < 0.0 )  xx = 0.0
    520              yy = 0.1886 * EXP( 20.306 / lambda_re )
     521             yy = 0.1886 * EXP( 20.306_wp / lambda_re )
    521522
    522523             c1_gr  =  xx / ( g / vk * tauk )**yy
    523524
    524              ao_gr  = ao + ( pi / 8.0) * ( g / vk * tauk )**2
     525             ao_gr  = ao + ( pi / 8.0_wp) * ( g / vk * tauk )**2
    525526             fao_gr = 20.115 * SQRT( ao_gr / lambda_re )
    526527             rc     = SQRT( fao_gr * ABS( st(j) - st(i) ) ) * eta   ! in cm
     
    652653          stb   = 32.0 * rho_a * ( rho_l - rho_a) * g / (3.0 * eta * eta)
    653654          phy   = sigma**3 * rho_a**2 / ( eta**4 * g * ( rho_l - rho_a ) )
    654           py    = phy**( 1.0 / 6.0 )
     655          py    = phy**( 1.0_wp / 6.0_wp )
    655656
    656657       ENDIF
     
    685686             ENDIF
    686687
    687              x = LOG( 16.0 * bond * py / 3.0 )
     688             x = LOG( 16.0 * bond * py / 3.0_wp )
    688689             y = 0.0
    689690
     
    695696
    696697             IF ( radclass(j) > 0.0035 )  THEN
    697                 winf(j) = xrey * eta / ( 2.0 * rho_a * 0.0035 )
     698                winf(j) = xrey * eta / ( 2.0 * rho_a * 0.0035_wp )
    698699             ELSE
    699700                winf(j) = xrey * eta / ( 2.0 * rho_a * radclass(j) )
     
    799800       ALLOCATE( ira(1:radius_classes) )
    800801       DO  j = 1, radius_classes
    801           particle_radius = radclass(j) * 1.0E6
     802          particle_radius = radclass(j) * 1.0E6_wp
    802803          DO  k = 1, 15
    803804             IF ( particle_radius < r0(k) )  THEN
     
    822823             IF ( ir < 16 )  THEN
    823824                IF ( ir >= 2 )  THEN
    824                    pp = ( ( radclass(j) * 1.0E06 ) - r0(ir-1) ) / &
     825                   pp = ( ( radclass(j) * 1.0E06_wp ) - r0(ir-1) ) / &
    825826                        ( r0(ir) - r0(ir-1) )
    826827                   qq = ( rq- rat(iq-1) ) / ( rat(iq) - rat(iq-1) )
     
    930931
    931932       DO  j = 1, radius_classes
    932           particle_radius = radclass(j) * 1.0E6
     933          particle_radius = radclass(j) * 1.0E6_wp
    933934          DO  k = 1, 7
    934935             IF ( particle_radius < r0(k) )  THEN
     
    10391040       REAL(wp), DIMENSION(1:9,1:19), SAVE ::  ef = 0.0          !:
    10401041
    1041        mean_rm = mean_r * 1.0E06
    1042        rm      = r      * 1.0E06
     1042       mean_rm = mean_r * 1.0E06_wp
     1043       rm      = r      * 1.0E06_wp
    10431044
    10441045       IF ( first )  THEN
  • palm/trunk/SOURCE/lpm_droplet_collision.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp_kind
    2323!
    2424! Former revisions:
     
    207207!--                Determine dissipation class index of this gridbox
    208208                   IF ( wang_kernel )  THEN
    209                       eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
     209                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
    210210                                    dissipation_classes ) + 1
    211211                      epsilon = diss(k,j,i)
     
    213213                      epsilon = 0.0
    214214                   ENDIF
    215                    IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
     215                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
    216216                      eclass = 0   ! Hall kernel is used
    217217                   ELSE
     
    270270!--                   Change of the current droplet radius
    271271                      rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
    272                                        sum3 )**0.33333333333333
     272                                       sum3 )**0.33333333333333_wp
    273273
    274274                      IF ( weight(n-is+1) < 0.0 )  THEN
     
    347347!--                   Change of the current droplet radius
    348348                      rad(n-is+1) = ( (r3 + dt_3d * ddV * (sum1 - sum2 * r3) )/&
    349                                        sum3 )**0.33333333333333
     349                                       sum3 )**0.33333333333333_wp
    350350
    351351                      IF ( weight(n-is+1) < 0.0 )  THEN
     
    534534!
    535535!--                   Change in radius due to collision
    536                       delta_r = effective_coll_efficiency / 3.0             &
     536                      delta_r = effective_coll_efficiency / 3.0_wp          &
    537537                                * pi * sl_r3 * ddx * ddy / dz               &
    538538                                * SQRT( ( u_int - particles(n)%speed_x )**2 &
     
    551551                      IF ( delta_v >= sl_r3  .AND.  sl_r3 > 0.0 )  THEN
    552552
    553                          delta_r = ( ( sl_r3/particles(n)%weight_factor )   &
    554                                      + particles(n)%radius**3 )**( 1./3. )  &
     553                         delta_r = ( ( sl_r3/particles(n)%weight_factor )           &
     554                                     + particles(n)%radius**3 )**( 1.0_wp/3.0_wp )  &
    555555                                     - particles(n)%radius
    556556
     
    612612
    613613                   IF ( wang_kernel )  THEN
    614                       eclass = INT( diss(k,j,i) * 1.0E4 / 1000.0 * &
     614                      eclass = INT( diss(k,j,i) * 1.0E4_wp / 1000.0_wp * &
    615615                                    dissipation_classes ) + 1
    616616                      epsilon = diss(k,j,i)
     
    618618                      epsilon = 0.0
    619619                   ENDIF
    620                    IF ( hall_kernel .OR. epsilon * 1.0E4 < 0.001 )  THEN
     620                   IF ( hall_kernel .OR. epsilon * 1.0E4_wp < 0.001 )  THEN
    621621                      eclass = 0   ! Hall kernel is used
    622622                   ELSE
     
    631631
    632632                   particles(psi)%radius = ( particles(psi)%radius**3 / &
    633                                              sum3 )**0.33333333333333
     633                                             sum3 )**0.33333333333333_wp
    634634                   particles(psi)%weight_factor = particles(psi)%weight_factor &
    635635                                                  * sum3
     
    652652
    653653                   particles(psi)%radius = ( particles(psi)%radius**3 / &
    654                                              sum3 )**0.33333333333333
     654                                             sum3 )**0.33333333333333_wp
    655655                   particles(psi)%weight_factor = particles(psi)%weight_factor &
    656656                                                  * sum3
  • palm/trunk/SOURCE/lpm_droplet_condensation.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    153153 
    154154!-- Parameters for Rosenbrock method
    155     REAL(wp), PARAMETER ::  a21 = 2.0          !:
    156     REAL(wp), PARAMETER ::  a31 = 48.0/25.0    !:
    157     REAL(wp), PARAMETER ::  a32 = 6.0/25.0     !:
    158     REAL(wp), PARAMETER ::  b1 = 19.0/9.0      !:
    159     REAL(wp), PARAMETER ::  b2 = 0.5           !:
    160     REAL(wp), PARAMETER ::  b3 = 25.0/108.0    !:
    161     REAL(wp), PARAMETER ::  b4 = 125.0/108.0   !:
    162     REAL(wp), PARAMETER ::  c21 = -8.0         !:
    163     REAL(wp), PARAMETER ::  c31 = 372.0/25.0   !:
    164     REAL(wp), PARAMETER ::  c32 = 12.0/5.0     !:
    165     REAL(wp), PARAMETER ::  c41 = -112.0/125.0 !:
    166     REAL(wp), PARAMETER ::  c42 = -54.0/125.0  !:
    167     REAL(wp), PARAMETER ::  c43 = -2.0/5.0     !:
    168     REAL(wp), PARAMETER ::  errcon = 0.1296    !:
    169     REAL(wp), PARAMETER ::  e1 = 17.0/54.0     !:
    170     REAL(wp), PARAMETER ::  e2 = 7.0/36.0      !:
    171     REAL(wp), PARAMETER ::  e3 = 0.0           !:
    172     REAL(wp), PARAMETER ::  e4 = 125.0/108.0   !:
    173     REAL(wp), PARAMETER ::  gam = 0.5          !:
    174     REAL(wp), PARAMETER ::  grow = 1.5         !:
    175     REAL(wp), PARAMETER ::  pgrow = -0.25      !:
    176     REAL(wp), PARAMETER ::  pshrnk = -1.0/3.0  !:
    177     REAL(wp), PARAMETER ::  shrnk = 0.5        !:
     155    REAL(wp), PARAMETER ::  a21 = 2.0             !:
     156    REAL(wp), PARAMETER ::  a31 = 48.0/25.0_wp    !:
     157    REAL(wp), PARAMETER ::  a32 = 6.0/25.0_wp     !:
     158    REAL(wp), PARAMETER ::  b1 = 19.0/9.0_wp      !:
     159    REAL(wp), PARAMETER ::  b2 = 0.5              !:
     160    REAL(wp), PARAMETER ::  b3 = 25.0/108.0_wp    !:
     161    REAL(wp), PARAMETER ::  b4 = 125.0/108.0_wp   !:
     162    REAL(wp), PARAMETER ::  c21 = -8.0            !:
     163    REAL(wp), PARAMETER ::  c31 = 372.0/25.0_wp   !:
     164    REAL(wp), PARAMETER ::  c32 = 12.0/5.0_wp     !:
     165    REAL(wp), PARAMETER ::  c41 = -112.0/125.0_wp !:
     166    REAL(wp), PARAMETER ::  c42 = -54.0/125.0_wp  !:
     167    REAL(wp), PARAMETER ::  c43 = -2.0/5.0_wp     !:
     168    REAL(wp), PARAMETER ::  errcon = 0.1296       !:
     169    REAL(wp), PARAMETER ::  e1 = 17.0/54.0_wp     !:
     170    REAL(wp), PARAMETER ::  e2 = 7.0/36.0_wp      !:
     171    REAL(wp), PARAMETER ::  e3 = 0.0              !:
     172    REAL(wp), PARAMETER ::  e4 = 125.0/108.0_wp   !:
     173    REAL(wp), PARAMETER ::  gam = 0.5             !:
     174    REAL(wp), PARAMETER ::  grow = 1.5            !:
     175    REAL(wp), PARAMETER ::  pgrow = -0.25         !:
     176    REAL(wp), PARAMETER ::  pshrnk = -1.0/3.0_wp  !:
     177    REAL(wp), PARAMETER ::  shrnk = 0.5           !:
    178178
    179179
  • palm/trunk/SOURCE/lpm_init.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    252252!
    253253!--    Calculate vertical depth of the sublayers
    254        height_int  = ( z_p - z0_av_global ) / REAL( number_of_sublayers )
     254       height_int  = ( z_p - z0_av_global ) / REAL( number_of_sublayers, KIND=wp )
    255255!
    256256!--    Precalculate LOG(z/z0)
  • palm/trunk/SOURCE/microphysics.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    522522!
    523523!--             Weight averaged radius of cloud droplets:
    524                 rc = 0.5 * ( xc * dpirho_l )**( 1.0 / 3.0 )
     524                rc = 0.5 * ( xc * dpirho_l )**( 1.0 / 3.0_wp )
    525525
    526526                alpha_cc = ( a_1 + a_2 * nu_c ) / ( 1.0 + a_3 * nu_c )
     
    529529!
    530530!--             Mixing length (neglecting distance to ground and stratification)
    531                 l_mix = ( dx * dy * dzu(k) )**( 1.0 / 3.0 )
     531                l_mix = ( dx * dy * dzu(k) )**( 1.0 / 3.0_wp )
    532532!
    533533!--             Limit dissipation rate according to Seifert, Nuijens and
     
    536536!
    537537!--             Compute Taylor-microscale Reynolds number:
    538                 re_lambda = 6.0 / 11.0 * ( l_mix / c_const )**( 2.0 / 3.0 ) *  &
    539                             SQRT( 15.0 / kin_vis_air ) * epsilon**( 1.0 / 6.0 )
     538                re_lambda = 6.0 / 11.0 * ( l_mix / c_const )**( 2.0 / 3.0_wp ) *  &
     539                            SQRT( 15.0 / kin_vis_air ) * epsilon**( 1.0 / 6.0_wp )
    540540!
    541541!--             The factor of 1.0E4 is needed to convert the dissipation rate
     
    664664!
    665665!--          Weight averaged diameter of rain drops:
    666              dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0 )
     666             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0_wp )
    667667!
    668668!--          Collisional breakup rate (Seifert, 2008):
     
    760760!
    761761!--          Weight averaged diameter of rain drops:
    762              dr = ( xr * dpirho_l )**( 1.0 / 3.0 )
     762             dr = ( xr * dpirho_l )**( 1.0 / 3.0_wp )
    763763!
    764764!--          Compute ventilation factor and intercept parameter
     
    772772!--             Slope parameter of gamma distribution (Seifert, 2008):
    773773                lambda_r = ( ( mu_r + 3.0 ) * ( mu_r + 2.0 ) *       &
    774                              ( mu_r + 1.0 ) )**( 1.0 / 3.0 ) / dr
     774                             ( mu_r + 1.0 ) )**( 1.0 / 3.0_wp ) / dr
    775775
    776776                mu_r_2   = mu_r + 2.0
     
    849849!
    850850!--    Sedimentation of cloud droplets (Heus et al., 2010):
    851        sed_qc_const = k_st * ( 3.0 / ( 4.0 * pi * rho_l ))**( 2.0 / 3.0 ) *   &
     851       sed_qc_const = k_st * ( 3.0 / ( 4.0 * pi * rho_l ))**( 2.0 / 3.0_wp ) *   &
    852852                     EXP( 5.0 * LOG( sigma_gc )**2 )
    853853
     
    856856       DO  k = nzt, nzb_s_inner(j,i)+1, -1
    857857          IF ( qc_1d(k) > eps_sb )  THEN
    858              sed_qc(k) = sed_qc_const * nc_1d(k)**( -2.0 / 3.0 ) * &
    859                         ( qc_1d(k) * hyrho(k) )**( 5.0 / 3.0 )
     858             sed_qc(k) = sed_qc_const * nc_1d(k)**( -2.0 / 3.0_wp ) * &
     859                        ( qc_1d(k) * hyrho(k) )**( 5.0 / 3.0_wp )
    860860          ELSE
    861861             sed_qc(k) = 0.0
     
    940940!
    941941!--          Weight averaged diameter of rain drops:
    942              dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0 )
     942             dr = ( hyrho(k) * qr_1d(k) / nr_1d(k) * dpirho_l )**( 1.0 / 3.0_wp )
    943943!
    944944!--          Shape parameter of gamma distribution (Milbrandt and Yau, 2005;
     
    948948!--          Slope parameter of gamma distribution (Seifert, 2008):
    949949             lambda_r = ( ( mu_r + 3.0 ) * ( mu_r + 2.0 ) *          &
    950                         ( mu_r + 1.0 ) )**( 1.0 / 3.0 ) / dr
     950                        ( mu_r + 1.0 ) )**( 1.0 / 3.0_wp ) / dr
    951951
    952952             w_nr(k) = MAX( 0.1, MIN( 20.0, a_term - b_term * ( 1.0 +          &
  • palm/trunk/SOURCE/netcdf.f90

    r1321 r1322  
    2323! Current revisions:
    2424! ------------------
    25 !
     25! Forgotten ONLY-attribute added to USE-statements
    2626!
    2727! Former revisions:
     
    49244924!------------------------------------------------------------------------------!
    49254925
    4926     USE control_parameters
     4926    USE control_parameters,                                                    &
     4927        ONLY:  netcdf_data_format
     4928
    49274929    USE netcdf
     4930
    49284931    USE netcdf_control
     4932
    49294933    USE pegrid
    49304934
     
    49624966!------------------------------------------------------------------------------!
    49634967
    4964     USE control_parameters
     4968    USE control_parameters,                                                    &
     4969        ONLY:  message_string
     4970
    49654971    USE netcdf
     4972
    49664973    USE netcdf_control
     4974
    49674975    USE pegrid
    49684976
  • palm/trunk/SOURCE/poismg.f90

    r1321 r1322  
    2323! Current revisions:
    2424! -----------------
    25 !
     25! REAL constants defined as wp-kind
    2626!
    2727! Former revisions:
     
    505505                                        ( r(k,j,i) - r(k-1,j+1,i+1) )
    506506
    507              f_mg(kc,jc,ic) = 1.0 / 64.0 * (                            &
     507             f_mg(kc,jc,ic) = 1.0 / 64.0_wp * (                         &
    508508                              8.0 * r(k,j,i)                            &
    509509                            + 4.0 * ( rkjim   + rkjip   +               &
     
    521521                            +       ( r(k+1,j-1,i-1) + r(k+1,j+1,i-1) + &
    522522                                      r(k+1,j-1,i+1) + r(k+1,j+1,i+1) ) &
    523                                            )
    524 
    525 !             f_mg(kc,jc,ic) = 1.0 / 64.0 * (                            &
     523                                              )
     524
     525!             f_mg(kc,jc,ic) = 1.0 / 64.0_wp * (                         &
    526526!                              8.0 * r(k,j,i)                            &
    527527!                            + 4.0 * ( r(k,j,i-1)     + r(k,j,i+1)     + &
     
    539539!                            +       ( r(k+1,j-1,i-1) + r(k+1,j+1,i-1) + &
    540540!                                      r(k+1,j-1,i+1) + r(k+1,j+1,i+1) ) &
    541 !                                           )
     541!                                              )
    542542          ENDDO
    543543       ENDDO
  • palm/trunk/SOURCE/read_var_list.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    906906!-- Calculate the temporal average of vertical profiles, if neccessary
    907907    IF ( average_count_pr /= 0 )  THEN
    908        hom_sum = hom_sum / REAL( average_count_pr )
     908       hom_sum = hom_sum / REAL( average_count_pr, KIND=wp )
    909909    ENDIF
    910910
  • palm/trunk/SOURCE/surface_coupler.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    386386!--       * latent heat of vaporization in m2/s2, or 540 cal/g, or 40.65 kJ/mol
    387387!--       /(rho_atm(=1.0)*c_p)
    388           tswst = tswst + qswst_remote * 2.2626108E6 / 1005.0
     388          tswst = tswst + qswst_remote * 2.2626108E6_wp / 1005.0_wp
    389389!
    390390!--        ...and convert it to a salinity flux at the sea surface (top)
     
    399399!--    (constants are the specific heat capacities for air and water)
    400400!--    now tswst is the ocean top heat flux
    401        tswst = tswst / rho(nzt,:,:) * 1005.0 / 4218.0
     401       tswst = tswst / rho(nzt,:,:) * 1005.0_wp / 4218.0_wp
    402402
    403403!
  • palm/trunk/SOURCE/temperton_fft.f90

    r1321 r1322  
    44! Current revisions:
    55! -----------------
    6 !
     6! REAL constants defined as wp-kind
    77!
    88! Former revisions:
     
    11741174    je = jd + 2*m*inc2
    11751175    z = 1.0_wp/REAL(n)
    1176     zsin45 = z*SQRT(0.5)
     1176    zsin45 = z*SQRT(0.5_wp)
    11771177
    11781178    DO l = 1, la
     
    1628162880  CONTINUE
    16291629    ibase = 0
    1630     sin45 = SQRT(0.5)
     1630    sin45 = SQRT(0.5_wp)
    16311631    DO l = 1, la
    16321632       i = ibase
  • palm/trunk/SOURCE/timestep.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    156156             v_gtrans = v_gtrans + hom(k,1,2,0)
    157157          ENDDO
    158           u_gtrans = u_gtrans / REAL( nzt - nzb )
    159           v_gtrans = v_gtrans / REAL( nzt - nzb )
     158          u_gtrans = u_gtrans / REAL( nzt - nzb, KIND=wp )
     159          v_gtrans = v_gtrans / REAL( nzt - nzb, KIND=wp )
    160160       ELSE
    161161!
     
    173173          ENDDO
    174174          !$acc end parallel
    175           uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
    176           uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb) )
     175          uv_gtrans_l(1) = u_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
     176          uv_gtrans_l(2) = v_gtrans_l / REAL( (nxr-nxl+1)*(nyn-nys+1)*(nzt-nzb), KIND=wp )
    177177#if defined( __parallel )
    178178          IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    179179          CALL MPI_ALLREDUCE( uv_gtrans_l, uv_gtrans, 2, MPI_REAL, MPI_SUM, &
    180180                              comm2d, ierr )
    181           u_gtrans = uv_gtrans(1) / REAL( numprocs )
    182           v_gtrans = uv_gtrans(2) / REAL( numprocs )
     181          u_gtrans = uv_gtrans(1) / REAL( numprocs, KIND=wp )
     182          v_gtrans = uv_gtrans(2) / REAL( numprocs, KIND=wp )
    183183#else
    184184          u_gtrans = uv_gtrans_l(1)
  • palm/trunk/SOURCE/timestep_scheme_steering.f90

    r1321 r1322  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    4848        ONLY:  intermediate_timestep_count, timestep_scheme, tsc
    4949
     50    USE kinds
     51
    5052    IMPLICIT NONE
    5153
     
    5759       IF ( timestep_scheme == 'runge-kutta-2' )  THEN
    5860          IF ( intermediate_timestep_count == 1 )  THEN
    59              tsc(1:5) = (/ 1.0, 1.0,  0.0, 0.0, 0.0 /)
     61             tsc(1:5) = (/ 1.0_wp, 1.0_wp,  0.0_wp, 0.0_wp, 0.0_wp /)
    6062          ELSE
    61              tsc(1:5) = (/ 1.0, 0.5, -0.5, 0.0, 1.0 /)
     63             tsc(1:5) = (/ 1.0_wp, 0.5_wp, -0.5_wp, 0.0_wp, 1.0_wp /)
    6264          ENDIF
    6365       ELSE
    6466          IF ( intermediate_timestep_count == 1 )  THEN
    65              tsc(1:5) = (/ 1.0,   1.0/3.0,        0.0, 0.0, 0.0 /)
     67             tsc(1:5) = (/ 1.0_wp,  1.0_wp /  3.0_wp,           0.0_wp, 0.0_wp, 0.0_wp /)
    6668          ELSEIF ( intermediate_timestep_count == 2 )  THEN
    67              tsc(1:5) = (/ 1.0, 15.0/16.0, -25.0/48.0, 0.0, 0.0 /)
     69             tsc(1:5) = (/ 1.0_wp, 15.0_wp / 16.0_wp, -25.0_wp/48.0_wp, 0.0_wp, 0.0_wp /)
    6870          ELSE
    69              tsc(1:5) = (/ 1.0,  8.0/15.0,   1.0/15.0, 0.0, 1.0 /)
     71             tsc(1:5) = (/ 1.0_wp,  8.0_wp / 15.0_wp,   1.0_wp/15.0_wp, 0.0_wp, 1.0_wp /)
    7072          ENDIF         
    7173       ENDIF
     
    7476!
    7577!--    Euler scheme
    76        tsc(1:5) = (/ 1.0, 1.0, 0.0, 0.0, 1.0 /)
     78       tsc(1:5) = (/ 1.0_wp, 1.0_wp, 0.0_wp, 0.0_wp, 1.0_wp /)
    7779
    7880    ENDIF
  • 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
  • palm/trunk/SOURCE/user_3d_data_averaging.f90

    r1321 r1322  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL functions provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    118118!                DO  j = nysg, nyng
    119119!                   DO  k = nzb, nzt+1
    120 !                      u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d )
     120!                      u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    121121!                   ENDDO
    122122!                ENDDO
Note: See TracChangeset for help on using the changeset viewer.