Ignore:
Timestamp:
Aug 30, 2019 1:25:00 PM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_data_output_mod.f90

    r4182 r4205  
    2525! -----------------
    2626! $Id$
     27! - Correct x,y-coordinates of vertical surfaces in netcdf output
     28! - Change definition of azimuth angle, reference is north 0 degree
     29! - zenith angle is always defined, also for vertical surfaces where it is
     30!   90 degree, while azimuth angle is only defined for vertical surfaces, not
     31!   for horizontal ones
     32!
     33! 4182 2019-08-22 15:20:23Z scharf
    2734! Corrected "Former revisions" section
    2835!
     
    10681075!
    10691076!--       For vertical surfaces, zenith angles are not defined (fill value).
    1070 !--       Azimuth angle: eastward (0), westward (180), northward (270),
    1071 !--       southward (90).
     1077!--       Azimuth angle: northward (0), eastward (90), southward (180),
     1078!--       westward (270).
     1079!--       Note, for vertical surfaces, zenith angles are 90.0_wp.
    10721080          DO  l = 0, 3
    10731081             IF ( l == 0 )  THEN
    1074                 az    = 270.0_wp
     1082                az    = 0.0_wp
    10751083                off_x = 0.5_wp
    10761084                off_y = 0.0_wp
    10771085             ELSEIF ( l == 1 )  THEN
     1086                az    = 180.0_wp
     1087                off_x = 0.5_wp
     1088                off_y = 1.0_wp
     1089             ELSEIF ( l == 2 )  THEN
    10781090                az    = 90.0_wp
    1079                 off_x = 0.5_wp
    1080                 off_y = 0.0_wp
    1081              ELSEIF ( l == 2 )  THEN
    1082                 az    = 0.0_wp
    10831091                off_x = 0.0_wp
    10841092                off_y = 0.5_wp
    10851093             ELSEIF ( l == 3 )  THEN
    1086                 az    = 180.0_wp
    1087                 off_x = 0.0_wp
     1094                az    = 270.0_wp
     1095                off_x = 1.0_wp
    10881096                off_y = 0.5_wp
    10891097             ENDIF
     
    10951103                surfaces%zs(mm)      = zu(surf_def_v(l)%k(m))
    10961104                surfaces%azimuth(mm) = az
    1097                 surfaces%zenith(mm)  = surfaces%fillvalue
     1105                surfaces%zenith(mm)  = 90.0_wp
    10981106                i                    = i + 1
    10991107                mm                   = mm + 1
     
    11051113                surfaces%zs(mm)      = zu(surf_lsm_v(l)%k(m))
    11061114                surfaces%azimuth(mm) = az
    1107                 surfaces%zenith(mm)  = surfaces%fillvalue
     1115                surfaces%zenith(mm)  = 90.0_wp
    11081116                i                    = i + 1
    11091117                mm                   = mm + 1
     
    11151123                surfaces%zs(mm)      = zu(surf_usm_v(l)%k(m))
    11161124                surfaces%azimuth(mm) = az
    1117                 surfaces%zenith(mm)  = surfaces%fillvalue
     1125                surfaces%zenith(mm)  = 90.0_wp
    11181126                i                    = i + 1
    11191127                mm                   = mm + 1
Note: See TracChangeset for help on using the changeset viewer.