The data in this data set are available via DODS
at
http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
and through the DODS html form mechanism:
http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/.html
Below is the Ferret script (and following is the output) that produced this file from the original hourly data set (which is also available via DODS and on this Live Access Server at http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc).
! Description: Create daily filtered OASIS data file from hourly gridded data ! Automatically generated by hourly2daily.pl on Thu Sep 20 15:04:58 2001. ! ! All data sets have U_UNCORRECTED_HR variable. Use it to get axis time limits ! Filter to remove tidal and inertial signals ! Regrid from 1 hour to 1 day to reduce # file size and number of points ! ! ! Leslie / Fred's plweight.m for 33 hour filter ! DEFINE AXIS/X=-33:33:1 XAX1 DEFINE GRID/X=XAX1 G1 FILE/VAR=WEIGHTS33/GRID=G1 "./33hrfilt" SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33 LIST WEIGHTS33 ! ! Original hourly data file ! USE "http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc" DEFINE SYMBOL TBEG="`U_UNCORRECTED_HR,return=tstart`" DEFINE SYMBOL TEND="`U_UNCORRECTED_HR,return=tend`" DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY SPAWN "/bin/rm -f dailyM1.nc" LET W33=XSEQUENCE(WEIGHTS33[D=1]) LET RELATIVE_HUMIDITY_FILT = CONVOLVEL(RELATIVE_HUMIDITY_HR,W33) LET RELATIVE_HUMIDITY_DAY = RELATIVE_HUMIDITY_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Relative Humidity"/UNITS="%" RELATIVE_HUMIDITY_DAY STAT RELATIVE_HUMIDITY_DAY SAVE/APPEND/FILE="dailyM1.nc" RELATIVE_HUMIDITY_DAY CANCEL VAR RELATIVE_HUMIDITY_DAY RELATIVE_HUMIDITY_FILT LET AIR_TEMPERATURE_FILT = CONVOLVEL(AIR_TEMPERATURE_HR,W33) LET AIR_TEMPERATURE_DAY = AIR_TEMPERATURE_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily air_temperature"/UNITS="Celsius" AIR_TEMPERATURE_DAY STAT AIR_TEMPERATURE_DAY SAVE/APPEND/FILE="dailyM1.nc" AIR_TEMPERATURE_DAY CANCEL VAR AIR_TEMPERATURE_DAY AIR_TEMPERATURE_FILT LET WIND_DIR_FILT = CONVOLVEL(WIND_DIR_HR,W33) LET WIND_DIR_DAY = WIND_DIR_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Wind Direction Magnetic (from)"/UNITS="Degrees" WIND_DIR_DAY STAT WIND_DIR_DAY SAVE/APPEND/FILE="dailyM1.nc" WIND_DIR_DAY CANCEL VAR WIND_DIR_DAY WIND_DIR_FILT LET REL_WIND_DIR_FILT = CONVOLVEL(REL_WIND_DIR_HR,W33) LET REL_WIND_DIR_DAY = REL_WIND_DIR_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Relative wind direction"/UNITS="Degrees" REL_WIND_DIR_DAY STAT REL_WIND_DIR_DAY SAVE/APPEND/FILE="dailyM1.nc" REL_WIND_DIR_DAY CANCEL VAR REL_WIND_DIR_DAY REL_WIND_DIR_FILT LET V_UNCORRECTED_FILT = CONVOLVEL(V_UNCORRECTED_HR,W33) LET V_UNCORRECTED_DAY = V_UNCORRECTED_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily North-South Velocity Towards(+N/-S) True no mooring motion correction"/UNITS="cm/s" V_UNCORRECTED_DAY STAT V_UNCORRECTED_DAY SAVE/APPEND/FILE="dailyM1.nc" V_UNCORRECTED_DAY CANCEL VAR V_UNCORRECTED_DAY V_UNCORRECTED_FILT LET E1_FILT = CONVOLVEL(E1_HR,W33) LET E1_DAY = E1_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Beam1 Echo Intensity data per beam"/UNITS="Counts" E1_DAY STAT E1_DAY SAVE/APPEND/FILE="dailyM1.nc" E1_DAY CANCEL VAR E1_DAY E1_FILT LET V_FILT = CONVOLVEL(V_HR,W33) LET V_DAY = V_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily North-South Velocity Towards(+N/-S) True Mooring Motion corrected"/UNITS="cm/s" V_DAY STAT V_DAY SAVE/APPEND/FILE="dailyM1.nc" V_DAY CANCEL VAR V_DAY V_FILT LET E3_FILT = CONVOLVEL(E3_HR,W33) LET E3_DAY = E3_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Beam3 Echo Intensity data per beam"/UNITS="Counts" E3_DAY STAT E3_DAY SAVE/APPEND/FILE="dailyM1.nc" E3_DAY CANCEL VAR E3_DAY E3_FILT LET PRESSURE_FILT = CONVOLVEL(PRESSURE_HR,W33) LET PRESSURE_DAY = PRESSURE_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Pressure"/UNITS="db" PRESSURE_DAY STAT PRESSURE_DAY SAVE/APPEND/FILE="dailyM1.nc" PRESSURE_DAY CANCEL VAR PRESSURE_DAY PRESSURE_FILT LET WIND_SPEED_FILT = CONVOLVEL(WIND_SPEED_HR,W33) LET WIND_SPEED_DAY = WIND_SPEED_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Wind Speed"/UNITS="m/s" WIND_SPEED_DAY STAT WIND_SPEED_DAY SAVE/APPEND/FILE="dailyM1.nc" WIND_SPEED_DAY CANCEL VAR WIND_SPEED_DAY WIND_SPEED_FILT LET AIR_PRESS_FILT = CONVOLVEL(AIR_PRESS_HR,W33) LET AIR_PRESS_DAY = AIR_PRESS_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Air Pressure"/UNITS="mbar" AIR_PRESS_DAY STAT AIR_PRESS_DAY SAVE/APPEND/FILE="dailyM1.nc" AIR_PRESS_DAY CANCEL VAR AIR_PRESS_DAY AIR_PRESS_FILT LET SALINITY_FILT = CONVOLVEL(SALINITY_HR,W33) LET SALINITY_DAY = SALINITY_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Salinity"/UNITS="none" SALINITY_DAY STAT SALINITY_DAY SAVE/APPEND/FILE="dailyM1.nc" SALINITY_DAY CANCEL VAR SALINITY_DAY SALINITY_FILT LET TEMPERATURE_FILT = CONVOLVEL(TEMPERATURE_HR,W33) LET TEMPERATURE_DAY = TEMPERATURE_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Temperature"/UNITS="Celsius" TEMPERATURE_DAY STAT TEMPERATURE_DAY SAVE/APPEND/FILE="dailyM1.nc" TEMPERATURE_DAY CANCEL VAR TEMPERATURE_DAY TEMPERATURE_FILT LET BUOY_HEADING_FILT = CONVOLVEL(BUOY_HEADING_HR,W33) LET BUOY_HEADING_DAY = BUOY_HEADING_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Buoy heading Magnetic"/UNITS="Degrees" BUOY_HEADING_DAY STAT BUOY_HEADING_DAY SAVE/APPEND/FILE="dailyM1.nc" BUOY_HEADING_DAY CANCEL VAR BUOY_HEADING_DAY BUOY_HEADING_FILT LET U_UNCORRECTED_FILT = CONVOLVEL(U_UNCORRECTED_HR,W33) LET U_UNCORRECTED_DAY = U_UNCORRECTED_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily East-West Velocity Towards(+E/-W) True no mooring motion correction"/UNITS="cm/s" U_UNCORRECTED_DAY STAT U_UNCORRECTED_DAY SAVE/APPEND/FILE="dailyM1.nc" U_UNCORRECTED_DAY CANCEL VAR U_UNCORRECTED_DAY U_UNCORRECTED_FILT LET U_FILT = CONVOLVEL(U_HR,W33) LET U_DAY = U_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily East-West Velocity Towards(+E/-W) True Mooring Motion corrected"/UNITS="cm/s" U_DAY STAT U_DAY SAVE/APPEND/FILE="dailyM1.nc" U_DAY CANCEL VAR U_DAY U_FILT LET E2_FILT = CONVOLVEL(E2_HR,W33) LET E2_DAY = E2_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Beam2 Echo Intensity data per beam"/UNITS="Counts" E2_DAY STAT E2_DAY SAVE/APPEND/FILE="dailyM1.nc" E2_DAY CANCEL VAR E2_DAY E2_FILT LET E4_FILT = CONVOLVEL(E4_HR,W33) LET E4_DAY = E4_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Beam4 Echo Intensity data per beam"/UNITS="Counts" E4_DAY STAT E4_DAY SAVE/APPEND/FILE="dailyM1.nc" E4_DAY CANCEL VAR E4_DAY E4_FILT LET CONDUCTIVITY_FILT = CONVOLVEL(CONDUCTIVITY_HR,W33) LET CONDUCTIVITY_DAY = CONDUCTIVITY_FILT[GT=TIME_DAY] SET VAR/TITLE="Daily Conductivity"/UNITS="Siemens/m" CONDUCTIVITY_DAY STAT CONDUCTIVITY_DAY SAVE/APPEND/FILE="dailyM1.nc" CONDUCTIVITY_DAY CANCEL VAR CONDUCTIVITY_DAY CONDUCTIVITY_FILT QUIT
NOAA/PMEL TMAP
Program FERRET
Version 5.22 - 07/27/00
20-Sep-01 15:04
yes? ! Description: Create daily filtered OASIS data file from hourly gridded d
ata
yes? ! Automatically generated by hourly2daily.pl on Thu Sep 20 15:04:58 2001.
yes? !
yes? ! All data sets have U_UNCORRECTED_HR variable. Use it to get axis time l
imits
yes? ! Filter to remove tidal and inertial signals
yes? ! Regrid from 1 hour to 1 day to reduce # file size and number of points
yes? !
yes? !
yes? ! Leslie / Fred's plweight.m for 33 hour filter
yes? !
yes? DEFINE AXIS/X=-33:33:1 XAX1
yes? DEFINE GRID/X=XAX1 G1
yes? FILE/VAR=WEIGHTS33/GRID=G1 "./33hrfilt"
yes? SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33
yes? LIST WEIGHTS33
33 Hour Filter Weights
DATA SET: ./33hrfilt
-33 / 1: -0.00026
-32 / 2: -0.00112
-31 / 3: -0.00210
-30 / 4: -0.00315
-29 / 5: -0.00426
-28 / 6: -0.00536
-27 / 7: -0.00641
-26 / 8: -0.00735
-25 / 9: -0.00811
-24 / 10: -0.00864
-23 / 11: -0.00886
-22 / 12: -0.00872
-21 / 13: -0.00816
-20 / 14: -0.00714
-19 / 15: -0.00561
-18 / 16: -0.00356
-17 / 17: -0.00098
-16 / 18: 0.00213
-15 / 19: 0.00573
-14 / 20: 0.00979
-13 / 21: 0.01424
-12 / 22: 0.01901
-11 / 23: 0.02400
-10 / 24: 0.02910
-9 / 25: 0.03422
-8 / 26: 0.03922
-7 / 27: 0.04399
-6 / 28: 0.04841
-5 / 29: 0.05237
-4 / 30: 0.05576
-3 / 31: 0.05850
-2 / 32: 0.06051
-1 / 33: 0.06174
0 / 34: 0.06215
1 / 35: 0.06174
2 / 36: 0.06051
3 / 37: 0.05850
4 / 38: 0.05576
5 / 39: 0.05237
6 / 40: 0.04841
7 / 41: 0.04399
8 / 42: 0.03922
9 / 43: 0.03422
10 / 44: 0.02910
11 / 45: 0.02400
12 / 46: 0.01901
13 / 47: 0.01424
14 / 48: 0.00979
15 / 49: 0.00573
16 / 50: 0.00213
17 / 51: -0.00098
18 / 52: -0.00356
19 / 53: -0.00561
20 / 54: -0.00714
21 / 55: -0.00816
22 / 56: -0.00872
23 / 57: -0.00886
24 / 58: -0.00864
25 / 59: -0.00811
26 / 60: -0.00735
27 / 61: -0.00641
28 / 62: -0.00536
29 / 63: -0.00426
30 / 64: -0.00315
31 / 65: -0.00210
32 / 66: -0.00112
33 / 67: -0.00026
yes?
yes? !
yes? ! Original hourly data file
yes? !
yes? USE "http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hour
lyM1.nc"
*** NOTE: Evenly spaced axis has edges definition: LONGITUDE - ignored
*** NOTE: Evenly spaced axis has edges definition: LATITUDE - ignored
*** NOTE: Evenly spaced axis has edges definition: LONGITUDE_HR - ignored
*** NOTE: Evenly spaced axis has edges definition: LATITUDE_HR - ignored
*** NOTE: Evenly spaced axis has edges definition: LONGITUDE_MET_HR - ignored
*** NOTE: Evenly spaced axis has edges definition: LATITUDE_MET_HR - ignored
*** NOTE: Evenly spaced axis has edges definition: LONGITUDE_TS_HR - ignored
*** NOTE: Evenly spaced axis has edges definition: LATITUDE_TS_HR - ignored
yes? DEFINE SYMBOL TBEG="`U_UNCORRECTED_HR,return=tstart`"
!-> DEFINE SYMBOL TBEG="14-JAN-1992 00:00"
yes? DEFINE SYMBOL TEND="`U_UNCORRECTED_HR,return=tend`"
!-> DEFINE SYMBOL TEND="06-OCT-2000 23:00"
yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY
!-> DEFINE AXIS/T="14-JAN-1992 00:00":"06-OCT-2000 23:00":1/unit=day TIME_DAY
yes? SPAWN "/bin/rm -f dailyM1.nc"
yes? LET W33=XSEQUENCE(WEIGHTS33[D=1])
yes?
yes?
yes?
yes? LET RELATIVE_HUMIDITY_FILT = CONVOLVEL(RELATIVE_HUMIDITY_HR,W33)
yes? LET RELATIVE_HUMIDITY_DAY = RELATIVE_HUMIDITY_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Relative Humidity"/UNITS="%" RELATIVE_HUMIDITY_DAY
yes? STAT RELATIVE_HUMIDITY_DAY
Daily Relative Humidity
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2786
Minimum value: 46.749
Maximum value: 100.74
Mean value: 86.108 (unweighted average)
Standard deviation: 9.3259
yes? SAVE/APPEND/FILE="dailyM1.nc" RELATIVE_HUMIDITY_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR RELATIVE_HUMIDITY_DAY RELATIVE_HUMIDITY_FILT
yes?
yes?
yes? LET AIR_TEMPERATURE_FILT = CONVOLVEL(AIR_TEMPERATURE_HR,W33)
yes? LET AIR_TEMPERATURE_DAY = AIR_TEMPERATURE_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily air_temperature"/UNITS="Celsius" AIR_TEMPERATURE_DAY
yes? STAT AIR_TEMPERATURE_DAY
Daily air_temperature
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2608
Minimum value: 5.3359
Maximum value: 16.43
Mean value: 12.09 (unweighted average)
Standard deviation: 1.7578
yes? SAVE/APPEND/FILE="dailyM1.nc" AIR_TEMPERATURE_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR AIR_TEMPERATURE_DAY AIR_TEMPERATURE_FILT
yes?
yes?
yes? LET WIND_DIR_FILT = CONVOLVEL(WIND_DIR_HR,W33)
yes? LET WIND_DIR_DAY = WIND_DIR_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Wind Direction Magnetic (from)"/UNITS="Degrees" WIND_
DIR_DAY
yes? STAT WIND_DIR_DAY
Daily Wind Direction Magnetic (from)
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2449
Minimum value: 35.467
Maximum value: 320.04
Mean value: 229.33 (unweighted average)
Standard deviation: 58.424
yes? SAVE/APPEND/FILE="dailyM1.nc" WIND_DIR_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR WIND_DIR_DAY WIND_DIR_FILT
yes?
yes?
yes? LET REL_WIND_DIR_FILT = CONVOLVEL(REL_WIND_DIR_HR,W33)
yes? LET REL_WIND_DIR_DAY = REL_WIND_DIR_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Relative wind direction"/UNITS="Degrees" REL_WIND_DIR
_DAY
yes? STAT REL_WIND_DIR_DAY
Daily Relative wind direction
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2495
Minimum value: -2.187
Maximum value: 262.14
Mean value: 166.7 (unweighted average)
Standard deviation: 77.073
yes? SAVE/APPEND/FILE="dailyM1.nc" REL_WIND_DIR_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR REL_WIND_DIR_DAY REL_WIND_DIR_FILT
yes?
yes?
yes? LET V_UNCORRECTED_FILT = CONVOLVEL(V_UNCORRECTED_HR,W33)
yes? LET V_UNCORRECTED_DAY = V_UNCORRECTED_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily North-South Velocity Towards(+N/-S) True no mooring m
otion correction"/UNITS="cm/s" V_UNCORRECTED_DAY
yes? STAT V_UNCORRECTED_DAY
Daily North-South Velocity Towards(+N/-S) True no mooring motion correction
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 199547
Minimum value: -40.744
Maximum value: 45.43
Mean value: 0.78338 (unweighted average)
Standard deviation: 6.768
yes? SAVE/APPEND/FILE="dailyM1.nc" V_UNCORRECTED_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR V_UNCORRECTED_DAY V_UNCORRECTED_FILT
yes?
yes?
yes? LET E1_FILT = CONVOLVEL(E1_HR,W33)
yes? LET E1_DAY = E1_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Beam1 Echo Intensity data per beam"/UNITS="Counts" E1
_DAY
yes? STAT E1_DAY
Daily Beam1 Echo Intensity data per beam
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 186479
Minimum value: 31.996
Maximum value: 227.95
Mean value: 114.96 (unweighted average)
Standard deviation: 36.588
yes? SAVE/APPEND/FILE="dailyM1.nc" E1_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR E1_DAY E1_FILT
yes?
yes?
yes? LET V_FILT = CONVOLVEL(V_HR,W33)
yes? LET V_DAY = V_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily North-South Velocity Towards(+N/-S) True Mooring Moti
on corrected"/UNITS="cm/s" V_DAY
yes? STAT V_DAY
Daily North-South Velocity Towards(+N/-S) True Mooring Motion corrected
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 241153
Minimum value: -40.811
Maximum value: 45.486
Mean value: 1.3297 (unweighted average)
Standard deviation: 7.1524
yes? SAVE/APPEND/FILE="dailyM1.nc" V_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR V_DAY V_FILT
yes?
yes?
yes? LET E3_FILT = CONVOLVEL(E3_HR,W33)
yes? LET E3_DAY = E3_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Beam3 Echo Intensity data per beam"/UNITS="Counts" E3
_DAY
yes? STAT E3_DAY
Daily Beam3 Echo Intensity data per beam
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 186479
Minimum value: 33.811
Maximum value: 224.43
Mean value: 120.32 (unweighted average)
Standard deviation: 37.88
yes? SAVE/APPEND/FILE="dailyM1.nc" E3_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR E3_DAY E3_FILT
yes?
yes?
yes? LET PRESSURE_FILT = CONVOLVEL(PRESSURE_HR,W33)
yes? LET PRESSURE_DAY = PRESSURE_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Pressure"/UNITS="db" PRESSURE_DAY
yes? STAT PRESSURE_DAY
Daily Pressure
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -3.5 to 375
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 38280 (1*1*12*3190)
# flagged as bad data: 30171
Minimum value: 0.49877
Maximum value: 303.41
Mean value: 109.31 (unweighted average)
Standard deviation: 95.701
yes? SAVE/APPEND/FILE="dailyM1.nc" PRESSURE_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR PRESSURE_DAY PRESSURE_FILT
yes?
yes?
yes? LET WIND_SPEED_FILT = CONVOLVEL(WIND_SPEED_HR,W33)
yes? LET WIND_SPEED_DAY = WIND_SPEED_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Wind Speed"/UNITS="m/s" WIND_SPEED_DAY
yes? STAT WIND_SPEED_DAY
Daily Wind Speed
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2449
Minimum value: -0.9382
Maximum value: 15.43
Mean value: 5.1227 (unweighted average)
Standard deviation: 2.3114
yes? SAVE/APPEND/FILE="dailyM1.nc" WIND_SPEED_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR WIND_SPEED_DAY WIND_SPEED_FILT
yes?
yes?
yes? LET AIR_PRESS_FILT = CONVOLVEL(AIR_PRESS_HR,W33)
yes? LET AIR_PRESS_DAY = AIR_PRESS_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Air Pressure"/UNITS="mbar" AIR_PRESS_DAY
yes? STAT AIR_PRESS_DAY
Daily Air Pressure
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2740
Minimum value: 999.06
Maximum value: 1031.2
Mean value: 1018.1 (unweighted average)
Standard deviation: 4.8047
yes? SAVE/APPEND/FILE="dailyM1.nc" AIR_PRESS_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR AIR_PRESS_DAY AIR_PRESS_FILT
yes?
yes?
yes? LET SALINITY_FILT = CONVOLVEL(SALINITY_HR,W33)
yes? LET SALINITY_DAY = SALINITY_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Salinity"/UNITS="none" SALINITY_DAY
yes? STAT SALINITY_DAY
Daily Salinity
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -3.5 to 375
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 38280 (1*1*12*3190)
# flagged as bad data: 30495
Minimum value: 31.938
Maximum value: 37.123
Mean value: 33.824 (unweighted average)
Standard deviation: 0.26317
yes? SAVE/APPEND/FILE="dailyM1.nc" SALINITY_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR SALINITY_DAY SALINITY_FILT
yes?
yes?
yes? LET TEMPERATURE_FILT = CONVOLVEL(TEMPERATURE_HR,W33)
yes? LET TEMPERATURE_DAY = TEMPERATURE_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Temperature"/UNITS="Celsius" TEMPERATURE_DAY
yes? STAT TEMPERATURE_DAY
Daily Temperature
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -3.5 to 375
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 38280 (1*1*12*3190)
# flagged as bad data: 30352
Minimum value: 6.6932
Maximum value: 16.27
Mean value: 9.8923 (unweighted average)
Standard deviation: 1.7755
yes? SAVE/APPEND/FILE="dailyM1.nc" TEMPERATURE_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR TEMPERATURE_DAY TEMPERATURE_FILT
yes?
yes?
yes? LET BUOY_HEADING_FILT = CONVOLVEL(BUOY_HEADING_HR,W33)
yes? LET BUOY_HEADING_DAY = BUOY_HEADING_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Buoy heading Magnetic"/UNITS="Degrees" BUOY_HEADING_D
AY
yes? STAT BUOY_HEADING_DAY
Daily Buoy heading Magnetic
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): 3.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 3190 (1*1*1*3190)
# flagged as bad data: 2449
Minimum value: 28.672
Maximum value: 367.47
Mean value: 144.8 (unweighted average)
Standard deviation: 65.33
yes? SAVE/APPEND/FILE="dailyM1.nc" BUOY_HEADING_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR BUOY_HEADING_DAY BUOY_HEADING_FILT
yes?
yes?
yes? LET U_UNCORRECTED_FILT = CONVOLVEL(U_UNCORRECTED_HR,W33)
yes? LET U_UNCORRECTED_DAY = U_UNCORRECTED_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily East-West Velocity Towards(+E/-W) True no mooring mot
ion correction"/UNITS="cm/s" U_UNCORRECTED_DAY
yes? STAT U_UNCORRECTED_DAY
** unknown netCDF error code: -1
yes? SAVE/APPEND/FILE="dailyM1.nc" U_UNCORRECTED_DAY
** unknown netCDF error code: -1
yes? CANCEL VAR U_UNCORRECTED_DAY U_UNCORRECTED_FILT
yes?
yes?
yes? LET U_FILT = CONVOLVEL(U_HR,W33)
yes? LET U_DAY = U_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily East-West Velocity Towards(+E/-W) True Mooring Motion
corrected"/UNITS="cm/s" U_DAY
yes? STAT U_DAY
Daily East-West Velocity Towards(+E/-W) True Mooring Motion corrected
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 241153
Minimum value: -25.466
Maximum value: 32.722
Mean value: 0.72556 (unweighted average)
Standard deviation: 5.0842
yes? SAVE/APPEND/FILE="dailyM1.nc" U_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR U_DAY U_FILT
yes?
yes?
yes? LET E2_FILT = CONVOLVEL(E2_HR,W33)
yes? LET E2_DAY = E2_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Beam2 Echo Intensity data per beam"/UNITS="Counts" E2
_DAY
yes? STAT E2_DAY
Daily Beam2 Echo Intensity data per beam
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 186479
Minimum value: 31.896
Maximum value: 230.5
Mean value: 117.32 (unweighted average)
Standard deviation: 37.357
yes? SAVE/APPEND/FILE="dailyM1.nc" E2_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR E2_DAY E2_FILT
yes?
yes?
yes? LET E4_FILT = CONVOLVEL(E4_HR,W33)
yes? LET E4_DAY = E4_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Beam4 Echo Intensity data per beam"/UNITS="Counts" E4
_DAY
yes? STAT E4_DAY
Daily Beam4 Echo Intensity data per beam
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -2.5 to 502.5
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 322190 (1*1*101*3190)
# flagged as bad data: 186479
Minimum value: 22.466
Maximum value: 226.14
Mean value: 111.62 (unweighted average)
Standard deviation: 37.809
yes? SAVE/APPEND/FILE="dailyM1.nc" E4_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR E4_DAY E4_FILT
yes?
yes?
yes? LET CONDUCTIVITY_FILT = CONVOLVEL(CONDUCTIVITY_HR,W33)
yes? LET CONDUCTIVITY_DAY = CONDUCTIVITY_FILT[GT=TIME_DAY]
yes? SET VAR/TITLE="Daily Conductivity"/UNITS="Siemens/m" CONDUCTIVITY_DAY
yes? STAT CONDUCTIVITY_DAY
Daily Conductivity
LONGITUDE: 122W
LATITUDE: 36.8N
DEPTH (m): -3.5 to 375
TIME: 13-JAN-1992 12:00 to 07-OCT-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM1.nc
Total # of data points: 38280 (1*1*12*3190)
# flagged as bad data: 30291
Minimum value: 2.7764
Maximum value: 4.2387
Mean value: 3.6876 (unweighted average)
Standard deviation: 0.14001
yes? SAVE/APPEND/FILE="dailyM1.nc" CONDUCTIVITY_DAY
LISTing to file dailyM1.nc
yes? CANCEL VAR CONDUCTIVITY_DAY CONDUCTIVITY_FILT
yes? QUIT
A 33 hour low-pass filter is applied to the hourly time series. The plot below was created with these Ferret commands (33hrfilt is a local file containing the filter weights, the weight values are shown in the output above)
DEFINE AXIS/T=-33:33:1/UNITS=hour TAX1 DEFINE GRID/t=TAX1 g2 FILE/VAR=WEIGHTS33T/GRID=G2 "./33hrfilt" SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33T PLOT WEIGHTS33T FRAME/FILE=weights33.gif
The plot below shows an example of the high diurnal and tidal frequency energies contained in mooring time series data. The red line is the fourier transform of the above filter. Its sharp cutoff at 33 hours helps preserve low frequency energies and virtually eliminate the high frequency energies which could be aliased to lower frequencies if we did not apply this filter before sub-sampling to one day intervals. Below are the Ferret commands which produced this plot
USE "http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc" ! Define the time series at a point in space. ! Set the FFT, using explicit time specification. ! Plot the amplitude spectrum vs frequency. LET FFT_uwndtmp = U_UNCORRECTED_HR[Z=30] ! 30 meters depth LET FFT_uwndtim = MISSING(FFT_uwndtmp, 0.0) ! ffta needs data, assign to 0 LET FFT_uwndfft = ffta(FFT_uwndtim) SET VARIABLE/TITLE="U Velocity Energy Spectrum"/units=Amplitude FFT_uwndfft ! For easier interpretation, invert the frequency axis and plot the ! spectrum vs period: hours/cycle ! Get the frequency increment used in the FFT. LET FFT_nf = `FFT_uwndfft,return=lend` LET FFT_nyquist = 0.5 LET FFT_freq1 = FFT_nyquist/ FFT_nf ! Define a frequency axis. DEFINE AXIS/T=`FFT_freq1`:`FFT_nyquist`:`FFT_freq1` faxis DEFINE GRID/T=faxis gfftfreq LET a = T[g=gfftfreq] ! Define the period from the frequency axis. LET per = 1./a ! Plot as a "Y VS X" plot, showing the first 72 hours where the most energy is. ! The PPL ccommands clean up the plot appearance. SET WIN/SIZE=0.5 PLOT/VS/LINE/HLIMITS=0:72:3/VLIMITS=0:2/TITLE="Amplitude Spectrum"/SET_UP per[l=1:`FFT_nf`], FFT_uwndfft PPL XFOR (I2) PPL XLAB Period, hours/cycle PPL YLAB Energy of M4 ADCP U Velocity Component at 30m go unlabel 4 go unlabel 5 go unlabel 6 go unlabel 8 PPL PLOT ! Plot the filter's transform DEFINE AXIS/T=-33:33:1/UNITS=hour TAX1 DEFINE GRID/t=TAX1 g2 FILE/VAR=WEIGHTS33T/GRID=G2 "./33hrfilt" SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33T LET FFT_filt = ffta(WEIGHTS33T) LET FFT_filt_nf = `FFT_filt,return=lend` LET FFT_filt_freq1 = FFT_nyquist/ FFT_filt_nf ! Define a frequency axis. DEFINE AXIS/T=`FFT_filt_freq1`:`FFT_nyquist`:`FFT_filt_freq1` filtaxis DEFINE GRID/T=filtaxis gfiltfftfreq LET af = T[g=gfiltfftfreq] ! Define the period from the frequency axis. LET perfilt = 1./af ! Plot the filter over the transform of the data LET W33NORM = FFT_filt*31.43 SET VAR/TITLE="33 Hour Filter Transform" W33NORM PLOT/OVER/LINE/VS perfilt[l=1:`FFT_filt_nf`], W33NORM FRAME/FILE=w33_freq.gif
To use FFTA on data from other moorings which have multiple deployments, and hence gaps in the time axis you need to first regrid the data onto a regular time axis. Here are the Ferret command to do that:
set mem/size=20 def sym tst="`U_UNCORRECTED_HR,return=tstart`" def sym ten="`U_UNCORRECTED_HR,return=tend`" def axis/t=():():1/units=hour treg let u_reg = U_UNCORRECTED_HR[gt=treg] ! Can now pass u_reg to FFTA, or save it to a file...