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

    r2504 r2512  
    2525! -----------------
    2626! $Id$
     27! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny
     28! no output of ghost layer data
     29!
     30! 2504 2017-09-27 10:36:13Z maronga
    2731! Support roots and water under pavement. Added several pavement types.
    2832!
     
    38863890    LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
    38873891
    3888     REAL(wp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb:nzt+1) ::  local_pf !<
     3892    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !<
    38893893
    38903894
     
    39033907             ENDDO
    39043908          ELSE
    3905              DO  i = nxlg, nxrg
    3906                 DO  j = nysg, nyng
     3909             DO  i = nxl, nxr
     3910                DO  j = nys, nyn
    39073911                   local_pf(i,j,nzb+1) = c_liq_av(j,i)
    39083912                ENDDO
     
    39213925             ENDDO
    39223926          ELSE
    3923              DO  i = nxlg, nxrg
    3924                 DO  j = nysg, nyng
     3927             DO  i = nxl, nxr
     3928                DO  j = nys, nyn
    39253929                   local_pf(i,j,nzb+1) = c_soil_av(j,i)
    39263930                ENDDO
     
    39393943             ENDDO
    39403944          ELSE
    3941              DO  i = nxlg, nxrg
    3942                 DO  j = nysg, nyng
     3945             DO  i = nxl, nxr
     3946                DO  j = nys, nyn
    39433947                   local_pf(i,j,nzb+1) = c_veg_av(j,i)
    39443948                ENDDO
     
    39573961             ENDDO
    39583962          ELSE
    3959              DO  i = nxlg, nxrg
    3960                 DO  j = nysg, nyng
     3963             DO  i = nxl, nxr
     3964                DO  j = nys, nyn
    39613965                   local_pf(i,j,nzb+1) = ghf_av(j,i)
    39623966                ENDDO
     
    39753979             ENDDO
    39763980          ELSE
    3977              DO  i = nxlg, nxrg
    3978                 DO  j = nysg, nyng
     3981             DO  i = nxl, nxr
     3982                DO  j = nys, nyn
    39793983                   local_pf(i,j,nzb+1) = lai_av(j,i)
    39803984                ENDDO
     
    39933997             ENDDO
    39943998          ELSE
    3995              DO  i = nxlg, nxrg
    3996                 DO  j = nysg, nyng
     3999             DO  i = nxl, nxr
     4000                DO  j = nys, nyn
    39974001                   local_pf(i,j,nzb+1) = m_liq_av(j,i)
    39984002                ENDDO
     
    40134017             ENDDO
    40144018          ELSE
    4015              DO  i = nxlg, nxrg
    4016                 DO  j = nysg, nyng
     4019             DO  i = nxl, nxr
     4020                DO  j = nys, nyn
    40174021                   DO k = nzb_soil, nzt_soil
    40184022                      local_pf(i,j,k) = m_soil_av(k,j,i)
     
    40354039             ENDDO
    40364040          ELSE
    4037              DO  i = nxlg, nxrg
    4038                 DO  j = nysg, nyng
     4041             DO  i = nxl, nxr
     4042                DO  j = nys, nyn
    40394043                   local_pf(i,j,nzb+1) =  qsws_liq_av(j,i)
    40404044                ENDDO
     
    40534057             ENDDO
    40544058          ELSE
    4055              DO  i = nxlg, nxrg
    4056                 DO  j = nysg, nyng
     4059             DO  i = nxl, nxr
     4060                DO  j = nys, nyn
    40574061                   local_pf(i,j,nzb+1) =  qsws_soil_av(j,i)
    40584062                ENDDO
     
    40714075             ENDDO
    40724076          ELSE
    4073              DO  i = nxlg, nxrg
    4074                 DO  j = nysg, nyng
     4077             DO  i = nxl, nxr
     4078                DO  j = nys, nyn
    40754079                   local_pf(i,j,nzb+1) =  qsws_veg_av(j,i)
    40764080                ENDDO
     
    40904094             ENDDO
    40914095          ELSE
    4092              DO  i = nxlg, nxrg
    4093                 DO  j = nysg, nyng
     4096             DO  i = nxl, nxr
     4097                DO  j = nys, nyn
    40944098                   local_pf(i,j,nzb+1) = r_a_av(j,i)
    40954099                ENDDO
     
    41084112             ENDDO
    41094113          ELSE
    4110              DO  i = nxlg, nxrg
    4111                 DO  j = nysg, nyng
     4114             DO  i = nxl, nxr
     4115                DO  j = nys, nyn
    41124116                   local_pf(i,j,nzb+1) = r_s_av(j,i)
    41134117                ENDDO
     
    41284132             ENDDO
    41294133          ELSE
    4130              DO  i = nxlg, nxrg
    4131                 DO  j = nysg, nyng
     4134             DO  i = nxl, nxr
     4135                DO  j = nys, nyn
    41324136                   DO k = nzb_soil, nzt_soil
    41334137                      local_pf(i,j,k) = t_soil_av(k,j,i)
     
    41774181    LOGICAL      ::  found !<
    41784182
    4179     REAL(sp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb_soil:nzt_soil) ::  local_pf !<
     4183    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_soil:nzt_soil) ::  local_pf !<
    41804184
    41814185
Note: See TracChangeset for help on using the changeset viewer.