OHATS04: Logbook

OHATS04: Chronological Logbook Entries


Return to Logbook Top Page
Entry Date Title Site Author #Graphics
1 Thu 29-Apr-2004Logbook creatednoneaster
2 Thu 29-Apr-2004System upmarshallmaclean
3 Thu 29-Apr-2004Sonic status in rainmarshallmaclean
4 Thu 29-Apr-2004Marshall flow distortion test marshallhorst
5 Thu 29-Apr-2004sonic ?? outputmarshallmaclean
6 Thu 29-Apr-2004CSAT wind tunnel testmarshallhorst
7 Thu 29-Apr-2004sonic didsmarshallhorst
8 Mon 03-May-2004Data system samplingmarshallmaclean2
9 Mon 03-May-2004Fixed duplicate channel number 203 in prep.configmarshallhorst
10 Tue 04-May-2004sonic tiltsmarshallhorst
11 Fri 07-May-2004Raised data-system boxmarshallmaclean
12 Mon 10-May-2004chemistry trailersmarshallhorst
13 Mon 10-May-2004Sonic azimuthsmarshallhorst
14 Mon 10-May-2004current resultsmarshallhorst
15 Wed 12-May-2004height error analysismarshallhorst
16 Fri 14-May-2004shut down Marshallmarshallhorst
17 Mon 17-May-2004Final resultsmarshallhorst

1: LOG, Site none, Thu 29-Apr-2004 11:19:51 MDT, Logbook created 
initial log file /home/isff/aster/projects/OHATS/logbook/tklog.log created by USER aster
2: LOG, Site marshall, Thu 29-Apr-2004 11:21:39 MDT, System up 
Sonics on goal-post, data system up, network operational
3: sonics, Site marshall, Thu 29-Apr-2004 11:22:36 MDT, Sonic status in rain 
3 sonics, channel 201,205,206 (all of which have firmware rev 3.0.a)
are having problems in the rain, reporting

	>Consec hopped exceeded

Here are the firmware rev levels for each channel:

4.0s: 200,204,207
4.0r: 203
3.0a: 201,202,205,206

4.0s is apparently the rain firmware upgrade.  These sonics also have nylon 
"water wicks" on each transducer.  However sonic on channel 207 is
missing one of the upper wicks.

Channel 205 is having the most difficulties, and has needed to be 
restarted several times.


To restart it use the "D" command from rserial -A ragwort 205

rserial -A ragwort 205


Acq sigs  60. Type x to abort....................................<
612a!!!!!!! 2 02554 e 02681 b a 00449 a 00467 c e 02203 f 01556 34?5?
 
>??
ET= 60 ts=i XD=d GN=858a TK=1 UP=5 FK=0 RN=1 IT=1 DR=102 rx=2 fx=038 BX=0 AH=1
AT=0 RS=1 BR=0 RI=1 GO=00000 HA=0 6X=3 3X=2 PD=2 SD=0 ?d sa=1
WM=o ar=0 ZZ=0 DC=1  ELo=010 010 011 ELb=010 010 010 TNo=0o=-2-2-2 C0b=-2-2-2 RC=0 tlo=8 8 9 tlb=8 8 8 DTR=01740 CA=1 TD=  duty=087     AQ= 60 AC=1 CD=0 SR=1 UX=0 MX=0 DTU=02320 DTC=01160 RD=o ss=1 XP=2 RF=018 DS=007 SN0367 01nov99 HF=005 JC=3 CB=3 MD=5 DF=05000 RNA=1 rev 3.0a cs=37637 &=1 os=
>Consec hopped exceeded.
 
Acq sigs  60. Type x to abort....................................<
012a e 01448 e 03158 b a 00443 a 00470 c f 02500 f 01159 3456
 
>Consec hopped exceeded.
 
Acq sigs  60. Type x to abort....................................<
012a 6 03069 0 03946 b a 00304 a 00492 c e 02129 a 01656 3456
 
>Consec hopped exceeded.
 
Acq sigs  60. Type x to abort....................................<
012a 1 01978 a 01907 b b 00450 a 00422 c f 02304 b 01553 3456
 
>Consec hopped exceeded.
 
Acq sigs  60. Type x to abort....................................<
012a f 02723 9 02779 b b 00410 a 00371 c e 02178 f 01639 3456
 
>Consec hopped exceeded.
 
Acq sigs  60. Type x to abort....................................<
012a!!!!! e 02673 a 02742 b b 00460 a 00415 c! f456
4: array, Site marshall, Thu 29-Apr-2004 11:45:07 MDT, Marshall flow distortion test 
Analysis of HATS data showed that the spatial resolution of
finite-difference estimates of horizontal and vertical gradients can
be matched by making the vertical spacing of the two horizontal arrays
approximately equal to the sonic spacing within those arrays.  For OHATS,
WHOI engineers working on the array support structure (Mike Purcell and
Megan Carroll) have chosen the sonic spacing to be 58 cm.

The closest vertical spacing used for HATS was 1 m.  Comparison of wind
speed, TKE, and u* between the horizontal arrays and a nearby reference
tower do not show any evidence of flow distortion.  However it is possible
that the closer vertical spacing for the OHATS arrary may start to show
some evidence of flow distortion.  Peter Sullivan, John Wyngaard, and
I have discussed this issue and will address it from two approaches.
Peter found a potential flow solution for two parallel cylinders,
surrogates for the two horizontal masts, in Morse and Feshbach, Volume 2.
I have made preliminary steps to applying the technique of Wyngaards's
1981 investigation of flow distortion of a single cylinder to the two
cylinder problem, using Mathematica.  In the meantime, however, we have
also erected a subset of the planned OHATS array at Marshall to directly
measure flow distortion.

The Marshall array, or OHATS ops0, consists of two horizontal ASTER
masts at nominal heights of 4.75 m and 5.27 m.  (Because of interference
between the cross-bracing of the vertical ASTER masts and the attachment
hardware, we could not separate the horizontal masts by the desired 58 cm
and had to settle for 52 cm.)  The horizontal masts are oriented E-W so
that winds from the south are perpendicular to them.  About 6 m to the W
we erected a reference tower.  Nine sonic yokes were placed on each of
the two horizontal masts, with a spacing of 58 cm, and CSAT electronics
were mounted on the back of the masts in the intervening spaces.  Sonics
were mounted at the center three locations on both horizontal masts and
at the same two heights on the reference tower.

The sonics are designated as, e.g. u3, l2, etc, with u and l referring
to the upper and lower arrays and the ordinal number increasing from west
to east.  Thus u1 and l1 are the reference sonics and u3 and l3 are the
center sonics in the horizontal arrays.  The specific sonics are:

Loc	S/N	Channel
l1	855	200
l2	246	201
l3	247	202
l4	800	203
u1	856	204
u2	367	205
u3	369	206
u4	853	207

Sonics 853, 855, and 856 (u4, l1, u1) have the new CSI modification to
enhance their ability to operate in rain.  However, one of the nylon
rain socks is missing on an upper transducer of SN 853.
5: sonics, Site marshall, Thu 29-Apr-2004 11:59:23 MDT, sonic ?? output 
Output from "??" command on each sonic during ops0 flow distortion
test at Marshall.  Configuration as of April 29, 2004.
 
channel 200:
>??
SN0855 04mar04 rev 4.0s &=1 AC=1 AF=050 AH=1 AO=00300 ar=0 AQ= 60 BR=0 BX=0
 CF=1 C0o= 0 0 0 C0b= 0 0 0 CA=1 CD=0 cs=26307 DC=4 dl=015 DM=c DR=03465
 duty=087 DT=16240 ET= 60 FA=00050 FL=007 FX=038 GN=744a GO=00000 HA=0
 HG=01560 HH=02700 KT=0 LG=00832 LH=00100 MA=-020 MS=-010 MX=0 ND=1 NI=2
 ns=00203 OR=1 os=0 PD=2 RA=00020 RC=0 RF=00900 RH=015 RI=1 RS=1 RX=002 SD=0
 SL=035 SR=1 ss=1 T0123=1000 TD=a TF=02600 02600 02600 TK=1 TO= 0 0 0 TP=t
 ts=i UX=0 WM=o WR=006 XD=d xp=2 XX=00875 ZZ=0
  
 
 
channel 201
>??
ET= 60 ts=i XD=d GN=848a TK=1 UP=5 FK=0 RN=1 IT=1 DR=102 rx=2 fx=038 BX=0 AH=1
AT=0 RS=1 BR=0 RI=1 GO=00000 HA=0 6X=3 3X=2 PD=2 SD=0 ?d sa=1
WM=o ar=2 ZZ=0 DC=1  ELo=010 010 010 ELb=010 010 010 TNo=471 d TNb=576 JD= 005
C0o=-2-2-2 C0b=-2-2-2 RC=0 tlo=8 8 8 tlb=8 8 8 DTR=01740 CA=1 TD=  duty=086
AQ= 60 AC=1 CD=0 SR=1 UX=0 MX=0 DTU=02320 DTC=01160 RD=o ss=1 XP=2 RF=018 DS=007 SN0246 26mar98 HF=005 JC=3 CB=3 MD=5 DF=05000 RNA=1 rev 3.0a cs=36550 &=1 os=
>Consec hopped exceeded.
 
channel 202
 
ET= 60 ts=i XD=d GN=858a TK=1 UP=5 FK=0 RN=1 IT=1 DR=102 rx=2 fx=038 BX=0 AH=1
AT=0 RS=1 BR=0 RI=1 GO=00000 HA=0 6X=3 3X=2 PD=2 SD=0 ?d sa=1
WM=o ar=1 ZZ=0 DC=1  ELo=010 010 011 ELb=010 010 010 TNo=181 d TNb=670 JD= 005
C0o=-2-2-2 C0b=-2-2-2 RC=0 tlo=8 8 9 tlb=8 8 8 DTR=01740 CA=1 TD=  duty=087
AQ= 60 AC=1 CD=0 SR=1 UX=0 MX=0 DTU=02320 DTC=01160 RD=o ss=1 XP=2 RF=018 DS=007 SN0247 26mar98 HF=005 JC=3 CB=3 MD=5 DF=05000 RNA=1 rev 3.0a cs=04332 &=1 os=
>Consec hopped exceeded.

channel 203
>??
SN0800 19sep03 rev 4.0r &=1 AC=1 AF=050 AH=1 AO=00300 ar=0 AQ= 60 BR=0 BX=0
 CF=1 C0o= 0 0 0 C0b= 0 0 0 CA=1 CD=0 cs=59636 DC=4 dl=015 DM=c duty=087
 DT=16240 ET= 60 FA=00050 FL=007 FX=038 GN=866a GO=00000 HA=0 HG=01560
 HH=02700 KT=0 LG=00832 LH=00100 MA=-020 MS=-010 MX=0 ND=1 NI=2 ns=00203 OR=1
 os=0 PD=2 RA=00020 RC=0 RF=00900 RH=015 RI=1 RS=1 RX=002 SD=0 SL=035 SR=1
 ss=1 T0123=1000 TD=a TF=02600 02600 02600 TK=1 TO= 0 0 0 TP=t ts=i UX=0 WM=o
 XD=d xp=2 XX=00875 ZZ=0
                                                                                
channel 204:
>??
SN0856 04mar04 rev 4.0s &=1 AC=1 AF=050 AH=1 AO=00300 ar=0 AQ= 60 BR=0 BX=0
 CF=1 C0o= 0 0 0 C0b= 0 0 0 CA=1 CD=0 cs=47929 DC=4 dl=015 DM=c DR=03465
 duty=087 DT=16240 ET= 60 FA=00050 FL=007 FX=038 GN=555a GO=00000 HA=0
 HG=01560 HH=02700 KT=0 LG=00832 LH=00100 MA=-020 MS=-010 MX=0 ND=1 NI=2
 ns=00203 OR=1 os=0 PD=2 RA=00020 RC=0 RF=00900 RH=015 RI=1 RS=1 RX=002 SD=0
 SL=035 SR=1 ss=1 T0123=1000 TD=a TF=02600 02600 02600 TK=1 TO= 0 0 0 TP=t
 ts=i UX=0 WM=o WR=006 XD=d xp=2 XX=00875 ZZ=0
                                                                                
channel 205:
>??
ET= 60 ts=i XD=d GN=868a TK=1 UP=5 FK=0 RN=1 IT=1 DR=102 rx=2 fx=038 BX=0 AH=1
AT=0 RS=1 BR=0 RI=1 GO=00000 HA=0 6X=3 3X=2 PD=2 SD=0 ?d sa=1
WM=o ar=0 ZZ=0 DC=1  ELo=010 010 011 ELb=010 010 010 TNo=bf7 d TNb=75f JD= 007
C0o=-2-2-2 C0b=-2-2-2 RC=0 tlo=8 8 9 tlb=8 8 8 DTR=01740 CA=1 TD=  duty=089
AQ= 60 AC=1 CD=0 SR=1 UX=0 MX=0 DTU=02320 DTC=01160 RD=o ss=1 XP=2 RF=018 DS=007 SN0367 01nov99 HF=005 JC=3 CB=3 MD=5 DF=05000 RNA=1 rev 3.0a cs=37637 &=1 os=
>Consec hopped exceeded.

channel 206:
>??
ET= 60 ts=i XD=d GN=888a TK=1 UP=5 FK=0 RN=1 IT=1 DR=102 rx=2 fx=038 BX=0 AH=1
AT=0 RS=1 BR=0 RI=1 GO=00000 HA=0 6X=3 3X=2 PD=2 SD=0 ?d sa=1
WM=o ar=0 ZZ=0 DC=1  ELo=010 010 010 ELb=010 011 010 TNo=670 d TNb=a64 JD= 007
C0o=-2-2-2 C0b=-2-2-2 RC=0 tlo=8 8 8 tlb=8 9 8 DTR=01740 CA=1 TD=  duty=087
AQ= 60 AC=1 CD=0 SR=1 UX=0 MX=0 DTU=02320 DTC=01160 RD=o ss=1 XP=2 RF=018 DS=007 SN0369 04nov99 HF=005 JC=3 CB=3 MD=5 DF=05000 RNA=1 rev 3.0a cs=26112 &=1 os=
>Consec hopped exceeded.
                                                                                
channel 207:
>??
SN0853 04mar04 rev 4.0s &=1 AC=1 AF=050 AH=1 AO=00300 ar=0 AQ= 60 BR=0 BX=0
 CF=1 C0o= 0 0 0 C0b= 0 0 0 CA=1 CD=0 cs=43838 DC=4 dl=015 DM=c DR=03465
 duty=087 DT=16240 ET= 60 FA=00050 FL=007 FX=038 GN=567a GO=00000 HA=0
 HG=01560 HH=02700 KT=0 LG=00832 LH=00100 MA=-020 MS=-010 MX=0 ND=1 NI=2
 ns=00203 OR=1 os=0 PD=2 RA=00020 RC=0 RF=00900 RH=015 RI=1 RS=1 RX=002 SD=0
 SL=035 SR=1 ss=1 T0123=1000 TD=a TF=02600 02600 02600 TK=1 TO= 0 0 0 TP=t
 ts=i UX=0 WM=o WR=006 XD=d xp=2 XX=00875 ZZ=0
6: sonics, Site marshall, Thu 29-Apr-2004 14:01:33 MDT, CSAT wind tunnel test 
Summary of CSAT wind tunnel tests for OHATS, 4/26-27/04
Reran CSAT853.  Original run had smallest slope (0.955) and found the
rain sock on one of the lower transducers to be extended beyond the
transducer face by about half its length.

l1fit to all data, 1.9 m/s < speed < 19.6 m/s (11 pts):

246: sonic speed =  4.3 cm/s + 0.976 pitot speed (rms = 4.1 cm/s)
247: sonic speed =  1.5 cm/s + 0.979 pitot speed (rms = 2.9 cm/s)
364: sonic speed =  5.8 cm/s + 0.970 pitot speed (rms = 1.7 cm/s)
367: sonic speed =  1.4 cm/s + 0.979 pitot speed (rms = 3.7 cm/s)
368: sonic speed =  7.2 cm/s + 0.982 pitot speed (rms = 1.9 cm/s)
369: sonic speed =  0.6 cm/s + 0.974 pitot speed (rms = 2.5 cm/s)
370: sonic speed =  4.2 cm/s + 0.980 pitot speed (rms = 5.1 cm/s)
373: sonic speed =  8.1 cm/s + 0.976 pitot speed (rms = 2.8 cm/s)
800: sonic speed =  4.6 cm/s + 0.968 pitot speed (rms = 2.9 cm/s)
853: sonic speed =  6.4 cm/s + 0.961 pitot speed (rms = 2.9 cm/s)
855: sonic speed =  2.7 cm/s + 0.973 pitot speed (rms = 3.2 cm/s)
856: sonic speed = -0.2 cm/s + 0.972 pitot speed (rms = 3.3 cm/s)

l1fit to  1.9 m/s < speed < 10 m/s (5 pts):

246: sonic speed = -4.8 cm/s + 0.991 pitot speed (rms = 1.3 cm/s)
247: sonic speed = -3.8 cm/s + 0.989 pitot speed (rms = 1.1 cm/s)
364: sonic speed =  2.8 cm/s + 0.976 pitot speed (rms = 1.1 cm/s)
367: sonic speed = -2.5 cm/s + 0.986 pitot speed (rms = 1.5 cm/s)
368: sonic speed =  4.0 cm/s + 0.988 pitot speed (rms = 0.6 cm/s)
369: sonic speed = -3.5 cm/s + 0.980 pitot speed (rms = 0.8 cm/s)
370: sonic speed = -1.5 cm/s + 0.988 pitot speed (rms = 1.6 cm/s)
373: sonic speed = -0.2 cm/s + 0.988 pitot speed (rms = 1.2 cm/s)
800: sonic speed = -1.4 cm/s + 0.978 pitot speed (rms = 1.0 cm/s)
853: sonic speed =  2.0 cm/s + 0.968 pitot speed (rms = 1.0 cm/s)
855: sonic speed =  1.2 cm/s + 0.975 pitot speed (rms = 0.9 cm/s)
856: sonic speed = -2.6 cm/s + 0.976 pitot speed (rms = 0.5 cm/s)
7: LOG, Site marshall, Thu 29-Apr-2004 14:54:02 MDT, sonic dids 
I changed the dids from l{1:4}.4.5m and u{1:4}.5m to the actual heights,
l{1:4}.4.75m and u{1:4}.5.27m.  Changes made in ops0/prep.config and 
ops0/covar.config.  Gordon said he would restart covars.
8: data_system, Site marshall, Mon 03-May-2004 08:57:50 MDT, Data system sampling 
The attached plots show the sampling jitter and drift in a
1 hour 15 minute series from the u.l1.4.75m sonic.

One can see the 9600 baud inter-character spacing in the vertical
spacing between the groups of points:    

	10bits/byte / 9600 bits/sec = .001 secs/byte.

The jitter is about +-0.005 sec. This is about 1/3rd the sampling
period, which isn't great.

The first plot is from u.l1.4.75m (CSAT SN 855), and the
second is u.u1.5.27m (CSAT SN 856).


The system clock was kept with NTP.  Here are the ntp log messages
around that time.


May  3 12:43:29 localhost ntpd[55]: offset -0.001755 sec freq 7.762 ppm error 0.001982 poll 10
May  3 13:43:29 localhost ntpd[55]: offset -0.001508 sec freq 7.755 ppm error 0.002179 poll 10
May  3 14:43:29 localhost ntpd[55]: offset -0.007275 sec freq 7.754 ppm error 0.025945 poll 10
9: LOG, Site marshall, Mon 03-May-2004 11:28:33 MDT, Fixed duplicate channel number 203 in prep.config 
In prep.config, two successive channels were designated 203, when they
should have been 202, 203.  I fixed this at 11:30 MDT and Gordon
restarted the processes.
10: sonics, Site marshall, Tue 04-May-2004 10:11:39 MDT, sonic tilts 
Marshall sonic tilts determined on data from 5/1 - 5/4 09:10 MDT
Assumed w.off=0
rm.azm=45 deg for reference sonics l1 and u1
rm.azm=120 deg for sonics in horizontal array

sonic	w.off	lean	leanaz
l1	0	2.4	-8.1
l2	0	0.1	123.8
l3	0	1.0	-37.3
l4	0	1.9	21.9
u1	0	1.7	-0.3
u2	0	0.5	131.3
u3	0	0.1	91.5
u4	0	0.2	-61.8

Gordon applied these tilts and I verified them by 
rerunning plot.tilt as requested:

Tom,

I've re-run the covars with the new tilts.

Could you check them with plot.tilt?

The tilt files are in $ASTER/projects/OHATS/cal_files
if you want to check/change them.

Gordon

5/13

Marshall sonic tilts determined on data from 5/1 - 5/12 15:00 MDT
Assumed w.off=0
rm.azm=90 deg for reference sonics l1 and u1
rm.azm=120 deg for sonics in horizontal array

sonic	w.off	lean	leanaz	
l1	0	2.6	0.9
l2	0	0.2	-97.8
l3	0	1.1	-47.9
l4	0	1.7	18.2
u1	0	1.5	-7.8
u2	0	0.3	123.3	
u3	0	0.1	160.8
u4	0	0.3	-84.7

:sagebrush2:cal_files: grep apr csat3*

csat3.l1:2004 apr 29 00:00 	 0.00	 0.00	 0.00	2.6	0.9	269.1	0 
csat3.l2:2004 apr 29 00:00 	 0.00	 0.00	 0.00	0.2	-97.8	268.9	0 
csat3.l3:2004 apr 29 00:00 	 0.00	 0.00	 0.00	1.1	-47.9	268.9	0 
csat3.l4:2004 apr 29 00:00 	 0.00	 0.00	 0.00	1.7	18.2	269.4	0 
csat3.u1:2004 apr 29 00:00 	 0.00	 0.00	 0.00	1.5	-7.8	268.6	0 
csat3.u2:2004 apr 29 00:00 	 0.00	 0.00	 0.00	0.3	123.3	269.7	0 
csat3.u3:2004 apr 29 00:00 	 0.00	 0.00	 0.00	0.1	160.8	269.2	0 
csat3.u4:2004 apr 29 00:00 	 0.00	 0.00	 0.00	0.3	-84.7	268.6	0 

5/17

Checked the tilt correction by running plot.tilt on the corrected data
(without the azimuth correction!).  Two sonics (l2.4.75m and u2.5.27m)
had b coefficients of 0 0 -0.001; the rest had coefficients of 0 0 0.
11: data_system, Site marshall, Fri 07-May-2004 15:14:23 MDT, Raised data-system box 
Kurt had noticed water in the bottom of the data-system
weather box after the rains at the end of April.

So today at around 3:00pm Tom raised the data-system box
onto the tower.  This required disconnecting all the
data cables. We've checked with rserial to channel 200,204
that the cables are connected as before.


He also add some extra cable between the charge monitor and
the data-system - previously it was stretched tight.
12: array, Site marshall, Mon 10-May-2004 12:13:40 MDT, chemistry trailers 
There are 2 chemistry trailers to the north of the sonic array (which is
oriented for winds from the south).  The height of the trailers is 125"
and the distance to the horizontal sonic array is 80 paces/strides or
about 240'.  I also used the range finder function on the data scope
to estimate the distance in units of trailer height.  The distance was
roughly 24 (+/- 1) trailer heights, consistant with the paced distance.

Thus the flow distortion due to the trailers is about (1/24)^2 = 0.002.
13: sonics, Site marshall, Mon 10-May-2004 12:20:26 MDT, Sonic azimuths 
I used the data scope to measure the orientation of the sonics.

5/7:

Calibrated the data scope at Marshall.  This took three tries, the last
two using the monopod to help hold it steady as I rotated in 45 degree
increments.  The final try was graded "good 9".

5/8
Entered declinatation 10.2 E

5/10:

Measured azimuths for all 8 sonics and repeated the process 4 times:

Sonic					mean	median	std dev

u1	358.8	358.0	358.3	359.8	358.7	358.6	0.8
u2	000.1	359.2	358.7	000.7	359.7	359.7	0.9
u3	359.0	359.3	358.8	000.3	359.4	359.2	0.7
u4	359.1	358.8	358.4	358.3	358.7	358.6	0.4

l1	359.1	358.4	359.0	359.5	359.0	359.1	0.5
l2	359.0	358.2	358.8	359.2	358.8	358.9	0.4
l3	359.0	358.5	358.8	000.5	359.2	358.9	0.9
l4	359.5	359.3	358.7	000.9	359.6	359.4	0.9

Also measured a single azimuth 10 times.  In numerical order, the
samples were:

000.3
000.6
000.7
000.9
001.0
001.2
001.3
001.3
001.3
001.3

The median is 1.1.  Deleting the 0.3, the mean is 1.07, the standard
deviation is 0.3, and the range is +0.23, -0.47.
14: data analysis, Site marshall, Mon 10-May-2004 14:08:27 MDT, current results 
Marshall results to date, 5/1/2004 8:45 - 5/10/2004 12:45

linear fits of data from horizontal array sonics to reference sonics

wind speed (5 min averages):
spd.u2.5.27m = 0.980 spd.u1.5.27m + 0.02 m/s	nmr = 0.056
spd.u3.5.27m = 0.974 spd.u1.5.27m + 0.03 m/s	nmr = 0.059
spd.u4.5.27m = 0.984 spd.u1.5.27m + 0.03 m/s	nmr = 0.057

spd.l2.4.75m = 0.996 spd.l1.4.75m - 0.00 m/s	nmr = 0.063
spd.l3.4.75m = 0.986 spd.l1.4.75m + 0.01 m/s	nmr = 0.059
spd.l4.4.75m = 0.988 spd.l1.4.75m + 0.00 m/s	nmr = 0.059

turbulent kinetic energy (5 min averages):
tke.u2.5.27m = 1.00  tke.u1.5.27m - 0.001 m2/s2	nmr = 0.073
tke.u3.5.27m = 0.983 tke.u1.5.27m + 0.001 m2/s2	nmr = 0.075
tke.u4.5.27m = 0.999 tke.u1.5.27m - 0.001 m2/s2	nmr = 0.076

tke.l2.4.75m = 1.01 tke.l1.4.75m - 0.002 m2/s2	nmr = 0.075
tke.l3.4.75m = 1.01 tke.l1.4.75m - 0.002 m2/s2	nmr = 0.077
tke.l4.4.75m = 1.01 tke.l1.4.75m - 0.002 m2/s2	nmr = 0.078

friction velocity (30 min averages):
u*.u2.5.27m = 1.03 u*.u1.5.27m + 0.01 m/s	nmr = 0.231
u*.u3.5.27m = 1.02 u*.u1.5.27m + 0.01 m/s	nmr = 0.224
u*.u4.5.27m = 1.02 u*.u1.5.27m + 0.01 m/s	nmr = 0.225

u*.l2.4.75m = 1.02  u*.l1.4.75m + 0.001 m/s	nmr = 0.221
u*.l3.4.75m = 0.965 u*.l1.4.75m + 0.005 m/s	nmr = 0.300
u*.l4.4.75m = 0.962 u*.l1.4.75m + 0.00k m/s	nmr = 0.229

nmr = mean absolute deviation from linear fit/
                          standard deviation of data


For comparison, here are the same data for the fourth HATS configuration;
sonic spacing = 50 cm for d array, 62.5 cm for the s array, vertical
spacing = 1 m.

Speed (30 min averages)			nmr
d1.4m spd = 0.994 spd.0.4m + 0.02 m/s 	0.045
d2.4m spd = 0.980 spd.0.4m + 0.05 m/s	0.057
d3.4m spd = 0.995 spd.0.4m - 0.00 m/s	0.038	
d4.4m spd = 0.991 spd.0.4m - 0.00 m/s	0.044
d5.4m spd = 0.993 spd.0.4m + 0.08 m/s	0.050
d6.4m spd = 0.998 spd.0.4m + 0.03 m/s	0.051
d7.4m spd = 0.992 spd.0.4m + 0.03 m/s	0.037
d8.4m spd = 0.984 spd.0.4m + 0.01 m/s	0.055
d9.4m spd = 0.985 spd.0.4m + 0.04 m/s	0.053

s1.5m spd = 0.982 spd.0.5m + 0.05 m/s	0.053
s2.5m spd = 0.989 spd.0.5m + 0.08 m/s	0.040
s3.5m spd = 0.995 spd.0.5m + 0.06 m/s	0.048
s4.5m spd = 0.996 spd.0.5m + 0.05 m/s	0.039
s5.5m spd = 0.982 spd.0.5m + 0.08 m/s	0.047

Turbulent kinetic emergy (30 min averages):
d1.4m tke = 0.972 tke.0.4m + 0.020 m2/s2	0.072
d2.4m tke = 0.967 tke.0.4m + 0.047 m2/s2	0.071
d3.4m tke = 0.973 tke.0.4m - 0.044 m2/s2	0.071
d4.4m tke = 0.972 tke.0.4m - 0.003 m2/s2	0.069
d5.4m tke = 0.984 tke.0.4m + 0.079 m2/s2	0.067
d6.4m tke = 0.979 tke.0.4m + 0.028 m2/s2	0.065
d7.4m tke = 0.971 tke.0.4m + 0.033 m2/s2	0.063
d8.4m tke = 0.972 tke.0.4m + 0.013 m2/s2	0.061
d9.4m tke = 0.973 tke.0.4m + 0.039 m2/s2	0.059

s1.5m tke = 0.974 tke.0.5m + 2.5e-4 m2/s2	0.067
s2.5m tke = 0.986 tke.0.5m + 1.5e-4 m2/s2	0.062
s3.5m tke = 1.000 tke.0.5m + 0.8e-4 m2/s2	0.061
s4.5m tke = 1.000 tke.0.5m + 0.6e-4 m2/s2	0.060
s5.5m tke = 0.989 tke.0.5m - 0.2e-4 m2/s2	0.057

Friction velocity (30 min averages):
d1.4m u* = 1.014 u*.0.4m + 0.001 m/s 	0.156
d2.4m u* = 1.003 u*.0.4m - 0.001 m/s	0.183
d3.4m u* = 1.000 u*.0.4m + 0.000 m/s	0.169	
d4.4m u* = 0.993 u*.0.4m - 0.000 m/s	0.173
d5.4m u* = 0.996 u*.0.4m + 0.001 m/s	0.179
d6.4m u* = 0.990 u*.0.4m + 0.001 m/s	0.180
d7.4m u* = 0.993 u*.0.4m + 0.001 m/s	0.161
d8.4m u* = 0.989 u*.0.4m + 0.002 m/s	0.170
d9.4m u* = 1.002 u*.0.4m + 0.001 m/s	0.171

s1.5m u* = 0.990 u*.0.5m + 0.001 m/s	0.223
s2.5m u* = 1.004 u*.0.5m + 0.002 m/s	0.216
s3.5m u* = 0.988 u*.0.5m + 0.002 m/s	0.220
s4.5m u* = 1.021 u*.0.5m + 0.002 m/s	0.215
s5.5m u* = 1.015 u*.0.5m + 0.001 m/s	0.210
15: data analysis, Site marshall, Wed 12-May-2004 09:23:22 MDT, height error analysis 
What is error in wind speed due to a height difference between the
reference sonic and the horizontal array?

U = (u*/k)*ln(z/zo)

dU/dz = u*/kz = U/(z*ln(z/zo))

dU/U = dz/(z*ln(z/zo))
    
     = (1 cm/500 cm)/ln(500/3) = 0.04% per cm height error
                               
                               = 0.1% per inch height error

Thus a 1% error in U requires a 10 inch height error
16: LOG, Site marshall, Fri 14-May-2004 12:21:40 MDT, shut down Marshall 
We shut down data collection at Marshall around noon today (Friday, May 14).
We plan to tear down the sensors and data system this afternoon.
17: data analysis, Site marshall, Mon 17-May-2004 09:46:51 MDT, Final results 
Final Marshall results, 5/1/2004 8:45 - 5/12/2004 13:45

linear fits of data from horizontal array sonics to reference sonics

wind speed (942 5-min averages):
spd.u2.5.27m = 0.987 spd.u1.5.27m + 0.00 m/s	nmr = 0.042
spd.u3.5.27m = 0.983 spd.u1.5.27m + 0.01 m/s	nmr = 0.044
spd.u4.5.27m = 0.988 spd.u1.5.27m + 0.03 m/s	nmr = 0.043

spd.l2.4.75m = 0.994 spd.l1.4.75m + 0.00 m/s	nmr = 0.045
spd.l3.4.75m = 0.992 spd.l1.4.75m - 0.01 m/s	nmr = 0.044
spd.l4.4.75m = 0.984 spd.l1.4.75m + 0.01 m/s	nmr = 0.044

turbulent kinetic energy (942 5-min averages; TKE < 5 m^2/s^2):
tke.u2.5.27m = 0.984 tke.u1.5.27m + 0.001 m2/s2	nmr = 0.069
tke.u3.5.27m = 0.978 tke.u1.5.27m + 0.002 m2/s2	nmr = 0.070
tke.u4.5.27m = 0.993 tke.u1.5.27m + 0.001 m2/s2	nmr = 0.071

tke.l2.4.75m = 0.994 tke.l1.4.75m + 0.000 m2/s2	nmr = 0.072
tke.l3.4.75m = 0.995 tke.l1.4.75m + 0.000 m2/s2	nmr = 0.073
tke.l4.4.75m = 0.984 tke.l1.4.75m + 0.000 m2/s2	nmr = 0.074

friction velocity (359 15 min averages):
u*.u2.5.27m = 0.999 u*.u1.5.27m + 0.01 m/s	nmr = 0.202
u*.u3.5.27m = 0.995 u*.u1.5.27m + 0.01 m/s	nmr = 0.209
u*.u4.5.27m = 0.965 u*.u1.5.27m + 0.01 m/s	nmr = 0.215

u*.l2.4.75m = 1.01  u*.l1.4.75m + 0.01 m/s	nmr = 0.204
u*.l3.4.75m = 0.986 u*.l1.4.75m + 0.01 m/s	nmr = 0.222
u*.l4.4.75m = 0.960 u*.l1.4.75m + 0.01 m/s	nmr = 0.208

nmr = mean absolute deviation from linear fit/
                          standard deviation of data