Changeset 2512 for palm/trunk/SOURCE/user_data_output_2d.f90

Timestamp:

Oct 4, 2017 8:26:59 AM (7 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:

• palm/trunk/SOURCE/user_data_output_2d.f90 (modified) (6 diffs)

palm/trunk/SOURCE/user_data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! ghost layer points removed from output array local_pf
+! 
+! 2233 2017-05-30 18:08:54Z suehring
 !
 ! 2232 2017-05-30 17:47:52Z suehring
@@ -85 +88 @@
     LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
 
-    REAL(wp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb:nzt+1) ::  local_pf !< 
+    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !< 
 
 
@@ -99 +102 @@
 !       CASE ( 'u2_xy', 'u2_xz', 'u2_yz' )
 !          IF ( av == 0 )  THEN
-!             DO  i = nxlg, nxrg
-!                DO  j = nysg, nyng
+!             DO  i = nxl, nxr
+!                DO  j = nys, nyn
 !                   DO  k = nzb_do, nzt_do
 !                      local_pf(i,j,k) = u2(k,j,i)
@@ -107 +110 @@
 !             ENDDO
 !          ELSE
-!             DO  i = nxlg, nxrg
-!                DO  j = nysg, nyng
+!             DO  i = nxl, nxr
+!                DO  j = nys, nyn
 !                   DO  k = nzb_do, nzt_do
 !                      local_pf(i,j,k) = u2_av(k,j,i)
@@ -118 +121 @@
 !          grid = 'zu'
 !
-!--    In case two-dimensional surface variables are outputted, the user
+!--    In case two-dimensional surface variables are output, the user
 !--    has to access related surface-type. Uncomment and extend following lines
 !--    appropriately (example output of vertical surface momentum flux of u-
@@ -160 +163 @@
 
  END SUBROUTINE user_data_output_2d
-