Changeset 4768 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Nov 2, 2020 7:11:23 PM (4 years ago)

Author:

suehring

Message:

Enable 3D data output also with 64-bit precision

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! Enable 3D data output also with 64-bit precision
+! 
+! 4559 2020-06-11 08:51:48Z raasch
 ! file re-formatted to follow the PALM coding standard
 !
@@ -170 +173 @@
     REAL(wp)     ::  s_r3      !< sum( particle-radius**3 )
 
-    REAL(sp), DIMENSION(:,:,:), ALLOCATABLE ::  local_pf  !< output array
+    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  local_pf  !< output array
 
     REAL(wp), DIMENSION(:,:,:), POINTER ::  to_be_resorted  !< pointer to array which shall be
@@ -261 +264 @@
 !
 !--    Before each output, set array local_pf to fill value
-       local_pf = fill_value
+       local_pf = REAL( fill_value, KIND=wp )
 !
 !--    Set masking flag for topography for not resorted arrays
@@ -674 +677 @@
 !                count = (/ nxr-nxl+2, nyn-nys+2, nzt_do-nzb_do+1, 1 /) )
 !          ELSE
-             nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                          &
-                                     local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                      &
-                                     start = (/ nxl+1, nys+1, nzb_do+1, do3d_time_count(av) /),    &
-                                     count = (/ nxr-nxl+1, nyn-nys+1, nzt_do-nzb_do+1, 1 /)        &
-                                   )
+!
+!--       Call for parallel runs
+          nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                             &
+                                   local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                        &
+                                   start = (/ nxl+1, nys+1, nzb_do+1, do3d_time_count(av) /),      &
+                                   count = (/ nxr-nxl+1, nyn-nys+1, nzt_do-nzb_do+1, 1 /)          &
+                                 )
 !          ENDIF
           CALL netcdf_handle_error( 'data_output_3d', 386 )
@@ -685 +690 @@
 #else
 #if defined( __netcdf )
+!
+!--    Call for non-parallel runs
        nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_do3d(av,ivar),                                &
                                local_pf(nxl:nxr,nys:nyn,nzb_do:nzt_do),                            &