Ignore:
Timestamp:
Oct 4, 2017 8:26:59 AM (4 years ago)
Author:
raasch
Message:

upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny; no output if redundant ghost layer data to NetCDF files

File:
1 edited

Legend:

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

    r2318 r2512  
    2525! -----------------
    2626! $Id$
     27! upper bounds of 3d output changed from nx+1,ny+1 to nx,ny
     28! no output of ghost layer data
     29!
     30! 2318 2017-07-20 17:27:44Z suehring
    2731! Get topography top index via Function call
    2832!
     
    356360    LOGICAL      ::  found !<
    357361
    358     REAL(sp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb:nz_do3d) ::  local_pf !<
     362    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nz_do3d) ::  local_pf !<
    359363
    360364
     
    365369
    366370      CASE ( 'pcm_heatrate' )
    367          CALL exchange_horiz( pc_heating_rate, nbgp )
    368371         IF ( av == 0 )  THEN
    369             DO  i = nxlg, nxrg
    370                DO  j = nysg, nyng
     372            DO  i = nxl, nxr
     373               DO  j = nys, nyn
    371374                  DO  k = nzb_s_inner(j,i), nz_do3d
    372375                     local_pf(i,j,k) = pc_heating_rate(k-nzb_s_inner(j,i),j,i)
     
    379382      CASE ( 'pcm_lad' )
    380383         IF ( av == 0 )  THEN
    381             DO  i = nxlg, nxrg
    382                DO  j = nysg, nyng
     384            DO  i = nxl, nxr
     385               DO  j = nys, nyn
    383386                  DO  k = nzb_s_inner(j,i), nz_do3d
    384387                     local_pf(i,j,k) = lad_s(k-nzb_s_inner(j,i),j,i)
Note: See TracChangeset for help on using the changeset viewer.