Changeset 4205 for palm/trunk/SOURCE/surface_data_output_mod.f90
- Timestamp:
- Aug 30, 2019 1:25:00 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_data_output_mod.f90
r4182 r4205 25 25 ! ----------------- 26 26 ! $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 27 34 ! Corrected "Former revisions" section 28 35 ! … … 1068 1075 ! 1069 1076 !-- 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. 1072 1080 DO l = 0, 3 1073 1081 IF ( l == 0 ) THEN 1074 az = 270.0_wp1082 az = 0.0_wp 1075 1083 off_x = 0.5_wp 1076 1084 off_y = 0.0_wp 1077 1085 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 1078 1090 az = 90.0_wp 1079 off_x = 0.5_wp1080 off_y = 0.0_wp1081 ELSEIF ( l == 2 ) THEN1082 az = 0.0_wp1083 1091 off_x = 0.0_wp 1084 1092 off_y = 0.5_wp 1085 1093 ELSEIF ( l == 3 ) THEN 1086 az = 180.0_wp1087 off_x = 0.0_wp1094 az = 270.0_wp 1095 off_x = 1.0_wp 1088 1096 off_y = 0.5_wp 1089 1097 ENDIF … … 1095 1103 surfaces%zs(mm) = zu(surf_def_v(l)%k(m)) 1096 1104 surfaces%azimuth(mm) = az 1097 surfaces%zenith(mm) = surfaces%fillvalue1105 surfaces%zenith(mm) = 90.0_wp 1098 1106 i = i + 1 1099 1107 mm = mm + 1 … … 1105 1113 surfaces%zs(mm) = zu(surf_lsm_v(l)%k(m)) 1106 1114 surfaces%azimuth(mm) = az 1107 surfaces%zenith(mm) = surfaces%fillvalue1115 surfaces%zenith(mm) = 90.0_wp 1108 1116 i = i + 1 1109 1117 mm = mm + 1 … … 1115 1123 surfaces%zs(mm) = zu(surf_usm_v(l)%k(m)) 1116 1124 surfaces%azimuth(mm) = az 1117 surfaces%zenith(mm) = surfaces%fillvalue1125 surfaces%zenith(mm) = 90.0_wp 1118 1126 i = i + 1 1119 1127 mm = mm + 1
Note: See TracChangeset
for help on using the changeset viewer.