Changeset 1244

Timestamp:

Oct 31, 2013 8:16:56 AM (11 years ago)

Author:

raasch

Message:

bugfix for index bounds in case of 3d-parallel output

File:

• palm/trunk/SOURCE/data_output_3d.f90 (modified) (2 diffs)

palm/trunk/SOURCE/data_output_3d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! Bugfix for index bounds in case of 3d-parallel output
 !
 ! Former revisions:
@@ -560 +560 @@
                 IF ( nxr == nx  .AND.  nyn /= ny )  THEN
                    nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,if), &
-                                  local_pf(nxl:nxrg,nys:nyn,nzb:nz_do3d),    &
+                                     local_pf(nxl:nxr+1,nys:nyn,nzb:nz_do3d), &
                       start = (/ nxl+1, nys+1, nzb+1, do3d_time_count(av) /), &
-                      count = (/ nxr-nxl+1+nbgp, nyn-nys+1, nz_do3d-nzb+1, 1 /) )
+                      count = (/ nxr-nxl+2, nyn-nys+1, nz_do3d-nzb+1, 1 /) )
                 ELSEIF ( nxr /= nx  .AND.  nyn == ny )  THEN
                    nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,if), &
-                                  local_pf(nxl:nxr,nys:nyng,nzb:nz_do3d),    &
+                                     local_pf(nxl:nxr,nys:nyn+1,nzb:nz_do3d), &
                       start = (/ nxl+1, nys+1, nzb+1, do3d_time_count(av) /), &
-                      count = (/ nxr-nxl+1, nyn-nys+1+nbgp, nz_do3d-nzb+1, 1 /) )
+                      count = (/ nxr-nxl+1, nyn-nys+2, nz_do3d-nzb+1, 1 /) )
                 ELSEIF ( nxr == nx  .AND.  nyn == ny )  THEN
                    nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,if), &
-                                  local_pf(nxl:nxrg,nys:nyng,nzb:nz_do3d),  &
+                                   local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d), &
                       start = (/ nxl+1, nys+1, nzb+1, do3d_time_count(av) /), &
-                      count = (/ nxr-nxl+1+nbgp, nyn-nys+1+nbgp, nz_do3d-nzb+1, 1 /) )
+                      count = (/ nxr-nxl+2, nyn-nys+2, nz_do3d-nzb+1, 1 /) )
                 ELSE
                    nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,if), &
-                                  local_pf(nxl:nxr,nys:nyn,nzb:nz_do3d),      &
+                                       local_pf(nxl:nxr,nys:nyn,nzb:nz_do3d), &
                       start = (/ nxl+1, nys+1, nzb+1, do3d_time_count(av) /), &
                       count = (/ nxr-nxl+1, nyn-nys+1, nz_do3d-nzb+1, 1 /) )