Ignore:
Timestamp:
Mar 25, 2014 7:45:13 PM (10 years ago)
Author:
kanani
Message:

REAL constants defined as wp-kind

File:
1 edited

Legend:

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

    r1321 r1340  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants defined as wp-kind
    2323!
    2424! Former revisions:
     
    100100                                              v(k,j+1,i)      +  &
    101101                                              v(k,j+1,i-1) )     &
    102                                             / 4.0 )**2        +  &
     102                                            / 4.0_wp )**2     +  &
    103103                                          ( ( w(k-1,j,i-1)    +  &
    104104                                              w(k-1,j,i)      +  &
    105105                                              w(k,j,i-1)      +  &
    106106                                              w(k,j,i) )         &
    107                                             / 4.0 )**2 )      *  &
     107                                            / 4.0_wp )**2 )   *  &
    108108                                    u(k,j,i)
    109109                   ENDDO
     
    123123                                              u(k,j,i)        +  &
    124124                                              u(k,j,i+1) )       &
    125                                             / 4.0 )**2        +  &
     125                                            / 4.0_wp )**2     +  &
    126126                                              v(k,j,i)**2     +  &
    127127                                          ( ( w(k-1,j-1,i)    +  &
     
    129129                                              w(k,j-1,i)      +  &
    130130                                              w(k,j,i) )         &
    131                                             / 4.0 )**2 )      *  &
     131                                            / 4.0_wp )**2 )   *  &
    132132                                    v(k,j,i) 
    133133                   ENDDO
     
    147147                                              u(k+1,j,i)      +  &
    148148                                              u(k+1,j,i+1) )     &
    149                                             / 4.0 )**2        +  &
     149                                            / 4.0_wp )**2     +  &
    150150                                          ( ( v(k,j,i)        +  &
    151151                                              v(k,j+1,i)      +  &
    152152                                              v(k+1,j,i)      +  &
    153153                                              v(k+1,j+1,i) )     &
    154                                             / 4.0 )**2        +  &
     154                                            / 4.0_wp )**2     +  &
    155155                                              w(k,j,i)**2 )   *  &
    156156                                    w(k,j,i)
     
    183183                                    SQRT( ( ( u(k,j,i)        +       &
    184184                                              u(k,j,i+1) )            &
    185                                             / 2.0 )**2        +       &
     185                                            / 2.0_wp )**2     +       &
    186186                                          ( ( v(k,j,i)        +       &
    187187                                              v(k,j+1,i) )            &
    188                                             / 2.0 )**2        +       &
     188                                            / 2.0_wp )**2     +       &
    189189                                          ( ( w(k-1,j,i)      +       &
    190190                                              w(k,j,i) )              &
    191                                             / 2.0 )**2 )      *       &
     191                                            / 2.0_wp )**2 )   *       &
    192192                                    ( q(k,j,i) - sls(k,j,i) )
    193193                   ENDDO
     
    202202                   DO  k = nzb_s_inner(j,i)+1, pch_index
    203203                      tend(k,j,i) = tend(k,j,i) -                     &
    204                                     2.0 * cdc(k,j,i) * lad_s(k,j,i) * &
     204                                    2.0_wp * cdc(k,j,i) * lad_s(k,j,i) * &
    205205                                    SQRT( ( ( u(k,j,i)              + &
    206206                                              u(k,j,i+1) )            &
    207                                             / 2.0 )**2              + &
     207                                            / 2.0_wp )**2           + &
    208208                                          ( ( v(k,j,i)              + &
    209209                                              v(k,j+1,i) )            &
    210                                             / 2.0 )**2              + &
     210                                            / 2.0_wp )**2           + &
    211211                                          ( ( w(k,j,i)              + &
    212212                                              w(k+1,j,i) )            &
    213                                             / 2.0 )**2 )            * &
     213                                            / 2.0_wp )**2 )         * &
    214214                                    e(k,j,i)
    215215                   ENDDO
     
    267267                                        v(k,j+1,i)  +        &
    268268                                        v(k,j+1,i-1) )       &
    269                                       / 4.0 )**2    +        &
     269                                      / 4.0_wp )**2 +        &
    270270                                    ( ( w(k-1,j,i-1) +       &
    271271                                        w(k-1,j,i)   +       &
    272272                                        w(k,j,i-1)   +       &
    273273                                        w(k,j,i) )           &
    274                                       / 4.0 )**2 ) *         &
     274                                      / 4.0_wp )**2 ) *      &
    275275                              u(k,j,i)
    276276          ENDDO
     
    286286                                        u(k,j,i)     +       &
    287287                                        u(k,j,i+1) )         &
    288                                       / 4.0 )**2     +       &
     288                                      / 4.0_wp )**2  +       &
    289289                                        v(k,j,i)**2  +       &
    290290                                    ( ( w(k-1,j-1,i) +       &
     
    292292                                        w(k,j-1,i)   +       &
    293293                                        w(k,j,i) )           &
    294                                       / 4.0 )**2 ) *         &
     294                                      / 4.0_wp )**2 ) *      &
    295295                              v(k,j,i)
    296296          ENDDO
     
    306306                                        u(k+1,j,i)  +        &
    307307                                        u(k+1,j,i+1) )       &
    308                                       / 4.0 )**2    +        &
     308                                      / 4.0_wp )**2 +        &
    309309                                    ( ( v(k,j,i)    +        &
    310310                                        v(k,j+1,i)  +        &
    311311                                        v(k+1,j,i)  +        &
    312312                                        v(k+1,j+1,i) )       &
    313                                       / 4.0 )**2    +        &
     313                                      / 4.0_wp )**2 +        &
    314314                                        w(k,j,i)**2 ) *      &
    315315                              w(k,j,i)
     
    336336                              SQRT( ( ( u(k,j,i)        +       &
    337337                                        u(k,j,i+1) )            &
    338                                       / 2.0 )**2        +       &
     338                                      / 2.0_wp )**2     +       &
    339339                                    ( ( v(k,j,i)        +       &
    340340                                        v(k,j+1,i) )            &
    341                                       / 2.0 )**2        +       &
     341                                      / 2.0_wp )**2     +       &
    342342                                    ( ( w(k-1,j,i)      +       &
    343343                                        w(k,j,i) )              &
    344                                       / 2.0 )**2 )      *       &
     344                                      / 2.0_wp )**2 )   *       &
    345345                              ( q(k,j,i) - sls(k,j,i) )
    346346             ENDDO   
     
    351351          DO  k = nzb_s_inner(j,i)+1, pch_index   
    352352             tend(k,j,i) = tend(k,j,i) -                        &
    353                               2.0 * cdc(k,j,i) * lad_s(k,j,i) * &
     353                              2.0_wp * cdc(k,j,i) * lad_s(k,j,i) * &
    354354                              SQRT( ( ( u(k,j,i)           +    &
    355355                                        u(k,j,i+1) )            &
    356                                       / 2.0 )**2           +    & 
     356                                      / 2.0_wp )**2        +    & 
    357357                                    ( ( v(k,j,i)           +    &
    358358                                        v(k,j+1,i) )            &
    359                                       / 2.0 )**2           +    &
     359                                      / 2.0_wp )**2        +    &
    360360                                    ( ( w(k,j,i)           +    &
    361361                                        w(k+1,j,i) )            &
    362                                       / 2.0 )**2 )         *    &
     362                                      / 2.0_wp )**2 )      *    &
    363363                              e(k,j,i)
    364364          ENDDO
Note: See TracChangeset for help on using the changeset viewer.