The data in this data set are available via DODS
at
http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.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/hourlyM4.nc).
! Description: Create daily filtered OASIS data file from hourly gridded data ! Automatically generated by hourly2daily.pl on Thu Sep 20 09:53:44 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/hourlyM4.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 dailyM4.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="dailyM4.nc" RELATIVE_HUMIDITY_DAY 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="dailyM4.nc" AIR_TEMPERATURE_DAY 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="dailyM4.nc" WIND_DIR_DAY 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="dailyM4.nc" REL_WIND_DIR_DAY 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"/UNITS="cm/s" V_UNCORRECTED_DAY STAT V_UNCORRECTED_DAY SAVE/APPEND/FILE="dailyM4.nc" V_UNCORRECTED_DAY 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="dailyM4.nc" E1_DAY 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="dailyM4.nc" E3_DAY 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="dailyM4.nc" PRESSURE_DAY 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="dailyM4.nc" WIND_SPEED_DAY 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="dailyM4.nc" AIR_PRESS_DAY 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="dailyM4.nc" SALINITY_DAY 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="dailyM4.nc" TEMPERATURE_DAY 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="dailyM4.nc" BUOY_HEADING_DAY 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"/UNITS="cm/s" U_UNCORRECTED_DAY STAT U_UNCORRECTED_DAY SAVE/APPEND/FILE="dailyM4.nc" U_UNCORRECTED_DAY 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="dailyM4.nc" E2_DAY 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="dailyM4.nc" E4_DAY 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="dailyM4.nc" CONDUCTIVITY_DAY QUIT
NOAA/PMEL TMAP
Program FERRET
Version 5.22 - 07/27/00
20-Sep-01 09:53
yes? ! Description: Create daily filtered OASIS data file from hourly gridded d
ata
yes? ! Automatically generated by hourly2daily.pl on Thu Sep 20 09:53:44 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
lyM4.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="24-AUG-1999 00:00"
yes? DEFINE SYMBOL TEND="`U_UNCORRECTED_HR,return=tend`"
!-> DEFINE SYMBOL TEND="08-SEP-2000 00:00"
yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY
!-> DEFINE AXIS/T="24-AUG-1999 00:00":"08-SEP-2000 00:00":1/unit=day TIME_DAY
yes? SPAWN "/bin/rm -f dailyM4.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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 14
Minimum value: 65.96
Maximum value: 101.35
Mean value: 89.166 (unweighted average)
Standard deviation: 7.2087
yes? SAVE/APPEND/FILE="dailyM4.nc" RELATIVE_HUMIDITY_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 7
Minimum value: 9.139
Maximum value: 16.486
Mean value: 12.738 (unweighted average)
Standard deviation: 1.2874
yes? SAVE/APPEND/FILE="dailyM4.nc" AIR_TEMPERATURE_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 5
Minimum value: 88.154
Maximum value: 334.22
Mean value: 267.22 (unweighted average)
Standard deviation: 61.172
yes? SAVE/APPEND/FILE="dailyM4.nc" WIND_DIR_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 5
Minimum value: 56.693
Maximum value: 265.92
Mean value: 160.8 (unweighted average)
Standard deviation: 35.801
yes? SAVE/APPEND/FILE="dailyM4.nc" REL_WIND_DIR_DAY
LISTing to file dailyM4.nc
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"/UNITS="cm/s
" V_UNCORRECTED_DAY
yes? STAT V_UNCORRECTED_DAY
Daily North-South Velocity Towards(+N/-S) True
LONGITUDE: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 34587
Minimum value: -33.932
Maximum value: 24.787
Mean value: -1.0932 (unweighted average)
Standard deviation: 9.047
yes? SAVE/APPEND/FILE="dailyM4.nc" V_UNCORRECTED_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 32418
Minimum value: 40.123
Maximum value: 162.42
Mean value: 78.782 (unweighted average)
Standard deviation: 33.712
yes? SAVE/APPEND/FILE="dailyM4.nc" E1_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 32418
Minimum value: 43.015
Maximum value: 171.76
Mean value: 83.51 (unweighted average)
Standard deviation: 35.535
yes? SAVE/APPEND/FILE="dailyM4.nc" E3_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -3.5 to 375
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 4584 (1*1*12*382)
# flagged as bad data: 63
Minimum value: 0.55233
Maximum value: 348.14
Mean value: 130.22 (unweighted average)
Standard deviation: 111.58
yes? SAVE/APPEND/FILE="dailyM4.nc" PRESSURE_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 5
Minimum value: 1.4037
Maximum value: 14.379
Mean value: 6.5565 (unweighted average)
Standard deviation: 3.1076
yes? SAVE/APPEND/FILE="dailyM4.nc" WIND_SPEED_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 5
Minimum value: 999.79
Maximum value: 1030.7
Mean value: 1017.8 (unweighted average)
Standard deviation: 4.4727
yes? SAVE/APPEND/FILE="dailyM4.nc" AIR_PRESS_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -3.5 to 375
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 4584 (1*1*12*382)
# flagged as bad data: 220
Minimum value: 32.997
Maximum value: 34.223
Mean value: 33.771 (unweighted average)
Standard deviation: 0.3135
yes? SAVE/APPEND/FILE="dailyM4.nc" SALINITY_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -3.5 to 375
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 4584 (1*1*12*382)
# flagged as bad data: 132
Minimum value: 5.9065
Maximum value: 17.049
Mean value: 9.9892 (unweighted average)
Standard deviation: 2.4187
yes? SAVE/APPEND/FILE="dailyM4.nc" TEMPERATURE_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): 3.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 382 (1*1*1*382)
# flagged as bad data: 5
Minimum value: 49.577
Maximum value: 312.67
Mean value: 182.14 (unweighted average)
Standard deviation: 44.931
yes? SAVE/APPEND/FILE="dailyM4.nc" BUOY_HEADING_DAY
LISTing to file dailyM4.nc
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"/UNITS="cm/s"
U_UNCORRECTED_DAY
yes? STAT U_UNCORRECTED_DAY
Daily East-West Velocity Towards(+E/-W) True
LONGITUDE: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 34587
Minimum value: -37.28
Maximum value: 36.498
Mean value: -2.9224 (unweighted average)
Standard deviation: 10.948
yes? SAVE/APPEND/FILE="dailyM4.nc" U_UNCORRECTED_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 32418
Minimum value: 40.037
Maximum value: 172.97
Mean value: 80.571 (unweighted average)
Standard deviation: 36.192
yes? SAVE/APPEND/FILE="dailyM4.nc" E2_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -2.5 to 502.5
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 38582 (1*1*101*382)
# flagged as bad data: 32418
Minimum value: 41.909
Maximum value: 170.71
Mean value: 81.572 (unweighted average)
Standard deviation: 35.379
yes? SAVE/APPEND/FILE="dailyM4.nc" E4_DAY
LISTing to file dailyM4.nc
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: 122.4W
LATITUDE: 36.2N
DEPTH (m): -3.5 to 375
TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00
DATA SET: http://dods.shore.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc
Total # of data points: 4584 (1*1*12*382)
# flagged as bad data: 132
Minimum value: 3.3534
Maximum value: 4.2781
Mean value: 3.6912 (unweighted average)
Standard deviation: 0.19076
yes? SAVE/APPEND/FILE="dailyM4.nc" CONDUCTIVITY_DAY
LISTing to file dailyM4.nc
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...