Ignore:
Timestamp:
Mar 20, 2014 4:38:49 PM (10 years ago)
Author:
raasch
Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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
Note: See TracChangeset for help on using the changeset viewer.