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/surface_coupler.f90

    r1325 r1353  
    2020! Current revisions:
    2121! ------------------
     22! REAL constants provided with KIND-attribute
    2223!
    2324! Former revisions:
     
    206207!
    207208!--       Send heat flux at bottom surface to the ocean
    208           total_2d_a = 0.0
    209           total_2d   = 0.0
     209          total_2d_a = 0.0_wp
     210          total_2d   = 0.0_wp
    210211          total_2d(nys:nyn,nxl:nxr) = shf(nys:nyn,nxl:nxr)
    211212
     
    216217!--       Send humidity flux at bottom surface to the ocean
    217218          IF ( humidity )  THEN
    218              total_2d_a = 0.0
    219              total_2d   = 0.0
     219             total_2d_a = 0.0_wp
     220             total_2d   = 0.0_wp
    220221             total_2d(nys:nyn,nxl:nxr) = qsws(nys:nyn,nxl:nxr)
    221222
     
    236237!
    237238!--       Send momentum flux (u) at bottom surface to the ocean
    238           total_2d_a = 0.0
    239           total_2d   = 0.0
     239          total_2d_a = 0.0_wp
     240          total_2d   = 0.0_wp
    240241          total_2d(nys:nyn,nxl:nxr) = usws(nys:nyn,nxl:nxr)
    241242          CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, &
     
    244245!
    245246!--       Send momentum flux (v) at bottom surface to the ocean
    246           total_2d_a = 0.0
    247           total_2d   = 0.0
     247          total_2d_a = 0.0_wp
     248          total_2d   = 0.0_wp
    248249          total_2d(nys:nyn,nxl:nxr) = vsws(nys:nyn,nxl:nxr)
    249250          CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, &
     
    338339!
    339340!--       Send surface temperature to atmosphere
    340           total_2d_o = 0.0
    341           total_2d   = 0.0
     341          total_2d_o = 0.0_wp
     342          total_2d   = 0.0_wp
    342343          total_2d(nys:nyn,nxl:nxr) = pt(nzt,nys:nyn,nxl:nxr)
    343344
     
    367368!
    368369!--       Send u to atmosphere
    369           total_2d_o = 0.0 
    370           total_2d   = 0.0
     370          total_2d_o = 0.0_wp 
     371          total_2d   = 0.0_wp
    371372          total_2d(nys:nyn,nxl:nxr) = u(nzt,nys:nyn,nxl:nxr)
    372373          CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, &
     
    375376!
    376377!--       Send v to atmosphere
    377           total_2d_o = 0.0
    378           total_2d   = 0.0
     378          total_2d_o = 0.0_wp
     379          total_2d   = 0.0_wp
    379380          total_2d(nys:nyn,nxl:nxr) = v(nzt,nys:nyn,nxl:nxr)
    380381          CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, &
     
    396397!--       following Steinhorn (1991), JPO 21, pp. 1681-1683:
    397398!--       S'w' = -S * evaporation / ( rho_water * ( 1 - S ) )
    398           saswst = -1.0 * sa(nzt,:,:) * qswst_remote /  &
    399                     ( rho(nzt,:,:) * ( 1.0 - sa(nzt,:,:) ) )
     399          saswst = -1.0_wp * sa(nzt,:,:) * qswst_remote /  &
     400                    ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) ) )
    400401       ENDIF
    401402
     
    474475       dny2 = 2 * ( dny / 2 )
    475476
    476        total_2d_a = 0.0
     477       total_2d_a = 0.0_wp
    477478!
    478479!--    Interpolation from ocean-grid-layer to atmosphere-grid-layer
Note: See TracChangeset for help on using the changeset viewer.