Changeset 2512 for palm/trunk/SOURCE/urban_surface_mod.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/urban_surface_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! upper bounds of 3d output changed from nx+1,ny+1 to nx,ny
+! no output of ghost layer data
+! 
+! 2350 2017-08-15 11:48:26Z kanani
 ! Bugfix and error message for nopointer version. 
 ! Additional "! defined(__nopointer)" as workaround to enable compilation of 
@@ -2424 +2428 @@
         INTEGER(iwp), INTENT(IN)       ::  nzt_do    !< vertical upper limit of the data output (usually nz_do3d)
         LOGICAL, INTENT(OUT)           ::  found     !< 
-        REAL(sp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb_do:nzt_do) ::  local_pf   !< sp - it has to correspond to module data_output_3d
-        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg)     ::  temp_pf    !< temp array for urban surface output procedure
+        REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< sp - it has to correspond to module data_output_3d
+        REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf    !< temp array for urban surface output procedure
         
         CHARACTER (len=varnamelength)                          :: var, surfid
@@ -2900 +2904 @@
               
         END SELECT
-        
-!--     fill out array local_pf which is subsequently treated by data_output_3d
-        CALL exchange_horiz( temp_pf, nbgp )
-!
-!--  To Do: why reversed loop order
-        DO j = nysg,nyng
-            DO i = nxlg,nxrg
-                DO k = nzb_do, nzt_do
+!
+!--     Rearrange dimensions for NetCDF output
+        DO  j = nys, nyn
+            DO  i = nxl, nxr
+                DO  k = nzb_do, nzt_do
                     local_pf(i,j,k) = temp_pf(k,j,i)
                 ENDDO