Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (7 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

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

    r1321 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    121121                   ALLOCATE( e_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    122122                ENDIF
    123                 e_av = 0.0
     123                e_av = 0.0_wp
    124124
    125125             CASE ( 'lpt' )
     
    127127                   ALLOCATE( lpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    128128                ENDIF
    129                 lpt_av = 0.0
     129                lpt_av = 0.0_wp
    130130
    131131             CASE ( 'lwp*' )
     
    133133                   ALLOCATE( lwp_av(nysg:nyng,nxlg:nxrg) )
    134134                ENDIF
    135                 lwp_av = 0.0
     135                lwp_av = 0.0_wp
    136136
    137137             CASE ( 'nr' )
     
    139139                   ALLOCATE( nr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    140140                ENDIF
    141                 nr_av = 0.0
     141                nr_av = 0.0_wp
    142142
    143143             CASE ( 'p' )
     
    145145                   ALLOCATE( p_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    146146                ENDIF
    147                 p_av = 0.0
     147                p_av = 0.0_wp
    148148
    149149             CASE ( 'pc' )
     
    151151                   ALLOCATE( pc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    152152                ENDIF
    153                 pc_av = 0.0
     153                pc_av = 0.0_wp
    154154
    155155             CASE ( 'pr' )
     
    157157                   ALLOCATE( pr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    158158                ENDIF
    159                 pr_av = 0.0
     159                pr_av = 0.0_wp
    160160
    161161             CASE ( 'prr' )
     
    163163                   ALLOCATE( prr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    164164                ENDIF
    165                 prr_av = 0.0
     165                prr_av = 0.0_wp
    166166
    167167             CASE ( 'prr*' )
     
    169169                   ALLOCATE( precipitation_rate_av(nysg:nyng,nxlg:nxrg) )
    170170                ENDIF
    171                 precipitation_rate_av = 0.0
     171                precipitation_rate_av = 0.0_wp
    172172
    173173             CASE ( 'pt' )
     
    175175                   ALLOCATE( pt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    176176                ENDIF
    177                 pt_av = 0.0
     177                pt_av = 0.0_wp
    178178
    179179             CASE ( 'q' )
     
    181181                   ALLOCATE( q_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    182182                ENDIF
    183                 q_av = 0.0
     183                q_av = 0.0_wp
    184184
    185185             CASE ( 'qc' )
     
    187187                   ALLOCATE( qc_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    188188                ENDIF
    189                 qc_av = 0.0
     189                qc_av = 0.0_wp
    190190
    191191             CASE ( 'ql' )
     
    193193                   ALLOCATE( ql_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    194194                ENDIF
    195                 ql_av = 0.0
     195                ql_av = 0.0_wp
    196196
    197197             CASE ( 'ql_c' )
     
    199199                   ALLOCATE( ql_c_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    200200                ENDIF
    201                 ql_c_av = 0.0
     201                ql_c_av = 0.0_wp
    202202
    203203             CASE ( 'ql_v' )
     
    205205                   ALLOCATE( ql_v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    206206                ENDIF
    207                 ql_v_av = 0.0
     207                ql_v_av = 0.0_wp
    208208
    209209             CASE ( 'ql_vp' )
     
    211211                   ALLOCATE( ql_vp_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    212212                ENDIF
    213                 ql_vp_av = 0.0
     213                ql_vp_av = 0.0_wp
    214214
    215215             CASE ( 'qr' )
     
    217217                   ALLOCATE( qr_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    218218                ENDIF
    219                 qr_av = 0.0
     219                qr_av = 0.0_wp
    220220
    221221             CASE ( 'qsws*' )
     
    223223                   ALLOCATE( qsws_av(nysg:nyng,nxlg:nxrg) )
    224224                ENDIF
    225                 qsws_av = 0.0
     225                qsws_av = 0.0_wp
    226226
    227227             CASE ( 'qv' )
     
    229229                   ALLOCATE( qv_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    230230                ENDIF
    231                 qv_av = 0.0
     231                qv_av = 0.0_wp
    232232
    233233             CASE ( 'rho' )
     
    235235                   ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    236236                ENDIF
    237                 rho_av = 0.0
     237                rho_av = 0.0_wp
    238238
    239239             CASE ( 's' )
     
    241241                   ALLOCATE( s_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    242242                ENDIF
    243                 s_av = 0.0
     243                s_av = 0.0_wp
    244244
    245245             CASE ( 'sa' )
     
    247247                   ALLOCATE( sa_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    248248                ENDIF
    249                 sa_av = 0.0
     249                sa_av = 0.0_wp
    250250
    251251             CASE ( 'shf*' )
     
    253253                   ALLOCATE( shf_av(nysg:nyng,nxlg:nxrg) )
    254254                ENDIF
    255                 shf_av = 0.0
     255                shf_av = 0.0_wp
    256256
    257257             CASE ( 't*' )
     
    259259                   ALLOCATE( ts_av(nysg:nyng,nxlg:nxrg) )
    260260                ENDIF
    261                 ts_av = 0.0
     261                ts_av = 0.0_wp
    262262
    263263             CASE ( 'u' )
     
    265265                   ALLOCATE( u_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    266266                ENDIF
    267                 u_av = 0.0
     267                u_av = 0.0_wp
    268268
    269269             CASE ( 'u*' )
     
    271271                   ALLOCATE( us_av(nysg:nyng,nxlg:nxrg) )
    272272                ENDIF
    273                 us_av = 0.0
     273                us_av = 0.0_wp
    274274
    275275             CASE ( 'v' )
     
    277277                   ALLOCATE( v_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    278278                ENDIF
    279                 v_av = 0.0
     279                v_av = 0.0_wp
    280280
    281281             CASE ( 'vpt' )
     
    283283                   ALLOCATE( vpt_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    284284                ENDIF
    285                 vpt_av = 0.0
     285                vpt_av = 0.0_wp
    286286
    287287             CASE ( 'w' )
     
    289289                   ALLOCATE( w_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    290290                ENDIF
    291                 w_av = 0.0
     291                w_av = 0.0_wp
    292292
    293293             CASE ( 'z0*' )
     
    295295                   ALLOCATE( z0_av(nysg:nyng,nxlg:nxrg) )
    296296                ENDIF
    297                 z0_av = 0.0
     297                z0_av = 0.0_wp
    298298
    299299             CASE ( 'z0h*' )
     
    301301                   ALLOCATE( z0h_av(nysg:nyng,nxlg:nxrg) )
    302302                ENDIF
    303                 z0h_av = 0.0
     303                z0h_av = 0.0_wp
    304304
    305305             CASE DEFAULT
     
    380380                   DO  k = nzb, nzt+1
    381381                      psi = prt_start_index(k,j,i)
    382                       s_r3 = 0.0
    383                       s_r4 = 0.0
     382                      s_r3 = 0.0_wp
     383                      s_r4 = 0.0_wp
    384384                      DO  n = psi, psi+prt_count(k,j,i)-1
    385385                         s_r3 = s_r3 + particles(n)%radius**3 * &
     
    388388                                       particles(n)%weight_factor
    389389                      ENDDO
    390                       IF ( s_r3 /= 0.0 )  THEN
     390                      IF ( s_r3 /= 0.0_wp )  THEN
    391391                         mean_r = s_r4 / s_r3
    392392                      ELSE
    393                          mean_r = 0.0
     393                         mean_r = 0.0_wp
    394394                      ENDIF
    395395                      pr_av(k,j,i) = pr_av(k,j,i) + mean_r
Note: See TracChangeset for help on using the changeset viewer.