Changeset 2281 for palm


Ignore:
Timestamp:
Jun 13, 2017 11:34:50 AM (7 years ago)
Author:
suehring
Message:

Clean-up unnecessary index access to surface type

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r2258 r2281  
    56715671!--          For simplicity anterpolate within buildings and under elevated
    56725672!--          terrain too
    5673              DO  kk = kcb, kct
     5673             DO  kk = kcb, kct-1
    56745674!
    56755675!--             ijfc and kfc are precomputed in pmci_init_anterp_tophat
     
    56905690
    56915691             ENDDO
     5692             fc(kct,jj,ii) = fc(kct,jj,ii)
     5693
    56925694          ENDDO
    56935695       ENDDO
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r2233 r2281  
    2525! -----------------
    2626! $Id$
     27! Clean-up unnecessary index access to surface type
     28!
     29! 2233 2017-05-30 18:08:54Z suehring
    2730!
    2831! 2232 2017-05-30 17:47:52Z suehring
     
    12001203             i   = surf%i(m)           
    12011204             j   = surf%j(m)
    1202              k   = surf%k(m)
    12031205
    12041206             z_mo = surf%z_mo(m)
     
    13191321       ELSEIF ( TRIM( most_method ) == 'lookup' )  THEN
    13201322
    1321           !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo, li ) FIRSTPRIVATE( li_bnd ) LASTPRIVATE( li_bnd )
     1323          !$OMP PARALLEL DO PRIVATE( i, j, z_mo, li ) FIRSTPRIVATE( li_bnd ) LASTPRIVATE( li_bnd )
    13221324          DO  m = 1, surf%ns
    13231325
    13241326             i   = surf%i(m)           
    13251327             j   = surf%j(m)
    1326              k   = surf%k(m)
    13271328!
    13281329!--          If the bulk Richardson number is outside the range of the lookup
     
    14611462!--       Compute u* at upward-facing surfaces
    14621463          IF ( .NOT. downward )  THEN
    1463              !$OMP PARALLEL  DO PRIVATE( k, z_mo )
     1464             !$OMP PARALLEL  DO PRIVATE( z_mo )
    14641465             DO  m = 1, surf%ns
    14651466
    1466                 k    = surf%k(m)
    14671467                z_mo = surf%z_mo(m)
    14681468!
     
    14781478!--       any stability.
    14791479          ELSE
    1480              !$OMP PARALLEL  DO PRIVATE( k, z_mo )
     1480             !$OMP PARALLEL  DO PRIVATE( z_mo )
    14811481             DO  m = 1, surf%ns
    14821482
    1483                 k    = surf%k(m)
    14841483                z_mo = surf%z_mo(m)
    14851484!
     
    16111610!--    land-surface model, in order to compute aerodynamical resistance.
    16121611       IF ( surf_vertical )  THEN
    1613           !$OMP PARALLEL DO PRIVATE( i, j, k )
     1612          !$OMP PARALLEL DO PRIVATE( i, j )
    16141613          DO  m = 1, surf%ns
    16151614
    16161615             i          =  surf%i(m)           
    16171616             j          =  surf%j(m)
    1618              k          =  surf%k(m)
    16191617             surf%ts(m) = -surf%shf(m) / ( surf%us(m) + 1E-30_wp )
    16201618!
     
    17131711!--       Compute q* at vertical surfaces
    17141712          IF ( surf_vertical )  THEN
    1715              !$OMP PARALLEL DO PRIVATE( i, j, k )
     1713             !$OMP PARALLEL DO PRIVATE( i, j )
    17161714             DO  m = 1, surf%ns 
    17171715
    17181716                i          =  surf%i(m)           
    17191717                j          =  surf%j(m)
    1720                 k          =  surf%k(m)
    17211718                surf%qs(m) = -surf%qsws(m) / ( surf%us(m) + 1E-30_wp )
    17221719
     
    17331730!
    17341731!--          For a given scalar flux in the surface layer
    1735              !$OMP PARALLEL DO PRIVATE( i, j, k )
     1732             !$OMP PARALLEL DO PRIVATE( i, j )
    17361733             DO  m = 1, surf%ns 
    17371734                i   = surf%i(m)           
    17381735                j   = surf%j(m)
    1739                 k   = surf%k(m)
    17401736                surf%ss(m) = -surf%ssws(m) / ( surf%us(m) + 1E-30_wp )
    17411737             ENDDO
     
    17441740!--       At vertical surfaces
    17451741          IF ( surf_vertical )  THEN
    1746              !$OMP PARALLEL DO PRIVATE( i, j, k )
     1742             !$OMP PARALLEL DO PRIVATE( i, j )
    17471743             DO  m = 1, surf%ns 
    17481744                i          =  surf%i(m)           
    17491745                j          =  surf%j(m)
    1750                 k          =  surf%k(m)
    17511746                surf%ss(m) = -surf%ssws(m) / ( surf%us(m) + 1E-30_wp )
    17521747             ENDDO
     
    19241919          IF (  .NOT.  constant_scalarflux  .AND.  passive_scalar  .AND.       &
    19251920                .NOT.  downward )  THEN
    1926              !$OMP PARALLEL DO PRIVATE( i, j, k )
     1921             !$OMP PARALLEL DO PRIVATE( i, j )
    19271922             DO  m = 1, surf%ns   
    19281923
    19291924                i    = surf%i(m)           
    19301925                j    = surf%j(m)
    1931                 k    = surf%k(m)
    19321926                surf%ssws(m) = -surf%ss(m) * surf%us(m)
    19331927
     
    19381932          IF ( cloud_physics  .AND.  microphysics_seifert  .AND.              &
    19391933               .NOT. downward)  THEN
    1940              !$OMP PARALLEL DO PRIVATE( i, j, k )
     1934             !$OMP PARALLEL DO PRIVATE( i, j )
    19411935             DO  m = 1, surf%ns   
    19421936
    19431937                i    = surf%i(m)           
    19441938                j    = surf%j(m)
    1945                 k    = surf%k(m)
    19461939
    19471940                surf%qrsws(m) = -surf%qrs(m) * surf%us(m)
     
    20502043             ENDIF
    20512044
    2052              !$OMP PARALLEL DO PRIVATE( i, j, k, dum, z_mo )
     2045             !$OMP PARALLEL DO PRIVATE( i, j, dum, z_mo )
    20532046             DO  m = 1, surf%ns 
    20542047                i = surf%i(m)           
    20552048                j = surf%j(m)
    2056                 k = surf%k(m)
    20572049
    20582050                z_mo = surf%z_mo(m)
Note: See TracChangeset for help on using the changeset viewer.