Changeset 4205 for palm/trunk/SOURCE/surface_data_output_mod.f90

Timestamp:

Aug 30, 2019 1:25:00 PM (4 years ago)

Author:

suehring

Message:

plant canopy: Missing working precision + bugfix in calculation of wind speed; surface data output: Correct x,y-coordinates of vertical surfaces in netcdf output; Change definition of azimuth angle, reference is north 0 degree

File:

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

palm/trunk/SOURCE/surface_data_output_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - Correct x,y-coordinates of vertical surfaces in netcdf output
+! - Change definition of azimuth angle, reference is north 0 degree
+! - zenith angle is always defined, also for vertical surfaces where it is 
+!   90 degree, while azimuth angle is only defined for vertical surfaces, not
+!   for horizontal ones
+! 
+! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
 ! 
@@ -1068 +1075 @@
 !
 !--       For vertical surfaces, zenith angles are not defined (fill value).
-!--       Azimuth angle: eastward (0), westward (180), northward (270),
-!--       southward (90).
+!--       Azimuth angle: northward (0), eastward (90), southward (180), 
+!--       westward (270). 
+!--       Note, for vertical surfaces, zenith angles are 90.0_wp. 
           DO  l = 0, 3
              IF ( l == 0 )  THEN
-                az    = 270.0_wp
+                az    = 0.0_wp
                 off_x = 0.5_wp
                 off_y = 0.0_wp
              ELSEIF ( l == 1 )  THEN
+                az    = 180.0_wp
+                off_x = 0.5_wp
+                off_y = 1.0_wp
+             ELSEIF ( l == 2 )  THEN
                 az    = 90.0_wp
-                off_x = 0.5_wp
-                off_y = 0.0_wp
-             ELSEIF ( l == 2 )  THEN
-                az    = 0.0_wp
                 off_x = 0.0_wp
                 off_y = 0.5_wp
              ELSEIF ( l == 3 )  THEN
-                az    = 180.0_wp
-                off_x = 0.0_wp
+                az    = 270.0_wp
+                off_x = 1.0_wp
                 off_y = 0.5_wp
              ENDIF
@@ -1095 +1103 @@
                 surfaces%zs(mm)      = zu(surf_def_v(l)%k(m))
                 surfaces%azimuth(mm) = az
-                surfaces%zenith(mm)  = surfaces%fillvalue
+                surfaces%zenith(mm)  = 90.0_wp
                 i                    = i + 1
                 mm                   = mm + 1
@@ -1105 +1113 @@
                 surfaces%zs(mm)      = zu(surf_lsm_v(l)%k(m))
                 surfaces%azimuth(mm) = az
-                surfaces%zenith(mm)  = surfaces%fillvalue
+                surfaces%zenith(mm)  = 90.0_wp
                 i                    = i + 1
                 mm                   = mm + 1
@@ -1115 +1123 @@
                 surfaces%zs(mm)      = zu(surf_usm_v(l)%k(m))
                 surfaces%azimuth(mm) = az
-                surfaces%zenith(mm)  = surfaces%fillvalue
+                surfaces%zenith(mm)  = 90.0_wp
                 i                    = i + 1
                 mm                   = mm + 1