Ignore:
Timestamp:
Oct 28, 2015 5:14:10 PM (6 years ago)
Author:
raasch
Message:

FORTRAN an OpenMP errors removed
misplaced cpp-directive fixed
small E- and F-FORMAT changes to avoid informative compiler messages about insufficient field width

File:
1 edited

Legend:

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

    • Property svn:keywords set to Id
    r1692 r1697  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! FORTRAN and OpenMP errors removed
    2222!
    2323! Former revisions:
     
    2525! $Id$
    2626!
    27 !
     27! 1696 2015-10-27 10:03:34Z maronga
    2828! Modularized and completely re-written version of prandtl_fluxes.f90. In the
    2929! course of the re-writing two additional methods have been implemented. See
     
    9090!>    that this method cannot be used in case of roughness heterogeneity
    9191!>
    92 !> @todo limiting of uv_total might be required is some cases
     92!> @todo limiting of uv_total might be required in some cases
    9393!> @todo (re)move large_scale_forcing actions
    9494!> @todo check/optimize OpenMP and OpenACC directives
     
    387387                       l        !< look index
    388388
    389        REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: rib                
     389       REAL(wp), DIMENSION(nysg:nyng,nxlg:nxrg) :: rib !< Bulk Richardson number
    390390                       
    391391       REAL(wp)     :: f,      & !< Function for Newton iteration: f = Ri - [...]/[...]^2 = 0
     
    394394                       ol_m,   & !< Previous value of L for Newton iteration
    395395                       ol_old, & !< Previous time step value of L
    396                        ol_u,   & !< Upper bound of L for Newton iteration
    397                        rib       !< Bulk Richardson number
     396                       ol_u      !< Upper bound of L for Newton iteration
    398397
    399398!
     
    686685       IMPLICIT NONE
    687686
     687       !$OMP PARALLEL DO PRIVATE( k, z_mo )
     688       !$acc kernels loop
    688689       DO  i = nxlg, nxrg
    689690          DO  j = nysg, nyng
     
    694695!
    695696!--          Compute u* at the scalars' grid points
    696              !$OMP PARALLEL DO PRIVATE( a, b, k, uv_total, z_mo )
    697              !$acc kernels loop
    698697             us(j,i) = kappa * uv_total(j,i) / ( LOG( z_mo / z0(j,i) )         &
    699698                                          - psi_m( z_mo / ol(j,i) )            &
     
    764763          ENDIF
    765764
    766           !$OMP PARALLEL DO PRIVATE( a, b, k, z_mo )
     765          !$OMP PARALLEL DO PRIVATE( k, z_mo )
    767766          !$acc kernels loop
    768767          DO  i = nxlg, nxrg
     
    817816             ENDIF
    818817
    819              !$OMP PARALLEL DO PRIVATE( a, b, k, z_mo )
     818             !$OMP PARALLEL DO PRIVATE( e_s, k, z_mo )
    820819             !$acc kernels loop independent
    821820             DO  i = nxlg, nxrg
     
    859858       IF ( cloud_physics .AND. icloud_scheme == 0 .AND. precipitation )  THEN
    860859
    861           !$OMP PARALLEL DO PRIVATE( a, b, k, z_mo )
     860          !$OMP PARALLEL DO PRIVATE( k, z_mo )
    862861          !$acc kernels loop independent
    863862          DO  i = nxlg, nxrg
Note: See TracChangeset for help on using the changeset viewer.