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

palm/trunk/SOURCE/urban_surface_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Enable 3D data output also with 64-bit precision
+!
+! 4750 2020-10-16 14:27:48Z suehring
 ! - bugfix in openmp directive
 ! - make t_green_h and t_green_v public (required in indoor model)
@@ -2455 +2458 @@
     LOGICAL, INTENT(OUT)           ::  found     !<
 
-    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !< sp - it has to correspond to module data_output_3d
-    REAL(sp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr)     ::  temp_pf   !< temp array for urban surface output procedure
+    REAL(wp), 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
 
     found = .TRUE.
@@ -2523 +2526 @@
                 j = surf_usm_h(l)%j(m)
                 k = surf_usm_h(l)%k(m)
-                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = sp) )
+                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = wp) )
              ENDDO
           ELSE
@@ -2530 +2533 @@
                 j = surf_usm_v(l)%j(m)
                 k = surf_usm_v(l)%k(m)
-                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = sp) + 1.0_sp )
+                temp_pf(0,j,i) = MAX( temp_pf(0,j,i), REAL( k, KIND = wp) + 1.0_sp )
              ENDDO
           ENDIF