The data in this data set are available via DODS/OPeNDAP
at
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM3.nc
and through the DODS html form mechanism:
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM3.nc.html
In addition to this daily gridded data set the intermediate
hourly filtered data is available through DODS/OPeNDAP at the
url:
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyFilteredM3.nc
(The hourly filtered data set is not available through this LAS.)
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.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM3.nc).
For information on the individual deployments and instruments that contributed to this data set please see the M3 hourly data set documentation.
Summary plots of all the variables in this data set
! Description: Create daily filtered OASIS data file from hourly gridded data ! Automatically generated by ./hourly2dailyfilter.pl on Mon Jun 13 10:11:14 2005. ! ! All data sets have AIR_TEMPERATURE_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 ! Also save an intermediate hourly filtered file for finer scale use. ! SET MEMORY/SIZE=30 SET MODE DESPERATE ! ! 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.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM3.nc" USE "hourlyM3.nc" DEFINE SYMBOL TBEG="`AIR_TEMPERATURE_HR,return=tstart`" DEFINE SYMBOL TEND="`AIR_TEMPERATURE_HR,return=tend`" DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY SPAWN "/bin/rm -f hourlyFilteredM3.nc" SPAWN "/bin/rm -f dailyFilteredM3.nc" LET W33=XSEQUENCE(WEIGHTS33[D=1]) !------------------------------------------------------------------ ! Variable AIR_PRESS_HR: 1 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = AIR_PRESS_HR[L=@CIA] + AIR_PRESS_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET AIR_PRESS_FIL AIR_PRESS_HR[L=@FLN] * GAP_MASK ! LET AIR_PRESS_FIL = CONVOLVEL(AIR_PRESS_FIL,W33) LET AIR_PRESS_FIL = CONVOLVEL(AIR_PRESS_HR,W33) LET AIR_PRESS_DAY = AIR_PRESS_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Air Pressure"/UNITS="mbar" AIR_PRESS_DAY SET MEM/SIZE=30 STAT AIR_PRESS_FIL STAT AIR_PRESS_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" AIR_PRESS_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" AIR_PRESS_DAY CANCEL VAR AIR_PRESS_DAY CANCEL VAR AIR_PRESS_FIL !------------------------------------------------------------------ ! Variable AIR_TEMPERATURE_HR: 2 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = AIR_TEMPERATURE_HR[L=@CIA] + AIR_TEMPERATURE_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET AIR_TEMPERATURE_FIL AIR_TEMPERATURE_HR[L=@FLN] * GAP_MASK ! LET AIR_TEMPERATURE_FIL = CONVOLVEL(AIR_TEMPERATURE_FIL,W33) LET AIR_TEMPERATURE_FIL = CONVOLVEL(AIR_TEMPERATURE_HR,W33) LET AIR_TEMPERATURE_DAY = AIR_TEMPERATURE_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered air_temperature"/UNITS="Celsius" AIR_TEMPERATURE_DAY SET MEM/SIZE=30 STAT AIR_TEMPERATURE_FIL STAT AIR_TEMPERATURE_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" AIR_TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" AIR_TEMPERATURE_DAY CANCEL VAR AIR_TEMPERATURE_DAY CANCEL VAR AIR_TEMPERATURE_FIL !------------------------------------------------------------------ ! Variable CONDUCTIVITY_HR: 3 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = CONDUCTIVITY_HR[L=@CIA] + CONDUCTIVITY_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET CONDUCTIVITY_FIL CONDUCTIVITY_HR[L=@FLN] * GAP_MASK ! LET CONDUCTIVITY_FIL = CONVOLVEL(CONDUCTIVITY_FIL,W33) LET CONDUCTIVITY_FIL = CONVOLVEL(CONDUCTIVITY_HR,W33) LET CONDUCTIVITY_DAY = CONDUCTIVITY_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Conductivity"/UNITS="Siemens/m" CONDUCTIVITY_DAY SET MEM/SIZE=30 STAT CONDUCTIVITY_FIL STAT CONDUCTIVITY_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" CONDUCTIVITY_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" CONDUCTIVITY_DAY CANCEL VAR CONDUCTIVITY_DAY CANCEL VAR CONDUCTIVITY_FIL !------------------------------------------------------------------ ! Variable ECHO_INTENSITY_BEAM1_HR: 4 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = ECHO_INTENSITY_BEAM1_HR[L=@CIA] + ECHO_INTENSITY_BEAM1_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET ECHO_INTENSITY_BEAM1_FIL ECHO_INTENSITY_BEAM1_HR[L=@FLN] * GAP_MASK ! LET ECHO_INTENSITY_BEAM1_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM1_FIL,W33) LET ECHO_INTENSITY_BEAM1_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM1_HR,W33) LET ECHO_INTENSITY_BEAM1_DAY = ECHO_INTENSITY_BEAM1_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered beam1 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM1_DAY SET MEM/SIZE=30 STAT ECHO_INTENSITY_BEAM1_FIL STAT ECHO_INTENSITY_BEAM1_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM1_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM1_DAY CANCEL VAR ECHO_INTENSITY_BEAM1_DAY CANCEL VAR ECHO_INTENSITY_BEAM1_FIL !------------------------------------------------------------------ ! Variable ECHO_INTENSITY_BEAM2_HR: 5 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = ECHO_INTENSITY_BEAM2_HR[L=@CIA] + ECHO_INTENSITY_BEAM2_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET ECHO_INTENSITY_BEAM2_FIL ECHO_INTENSITY_BEAM2_HR[L=@FLN] * GAP_MASK ! LET ECHO_INTENSITY_BEAM2_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM2_FIL,W33) LET ECHO_INTENSITY_BEAM2_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM2_HR,W33) LET ECHO_INTENSITY_BEAM2_DAY = ECHO_INTENSITY_BEAM2_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered beam2 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM2_DAY SET MEM/SIZE=30 STAT ECHO_INTENSITY_BEAM2_FIL STAT ECHO_INTENSITY_BEAM2_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM2_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM2_DAY CANCEL VAR ECHO_INTENSITY_BEAM2_DAY CANCEL VAR ECHO_INTENSITY_BEAM2_FIL !------------------------------------------------------------------ ! Variable ECHO_INTENSITY_BEAM3_HR: 6 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = ECHO_INTENSITY_BEAM3_HR[L=@CIA] + ECHO_INTENSITY_BEAM3_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET ECHO_INTENSITY_BEAM3_FIL ECHO_INTENSITY_BEAM3_HR[L=@FLN] * GAP_MASK ! LET ECHO_INTENSITY_BEAM3_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM3_FIL,W33) LET ECHO_INTENSITY_BEAM3_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM3_HR,W33) LET ECHO_INTENSITY_BEAM3_DAY = ECHO_INTENSITY_BEAM3_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered beam3 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM3_DAY SET MEM/SIZE=30 STAT ECHO_INTENSITY_BEAM3_FIL STAT ECHO_INTENSITY_BEAM3_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM3_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM3_DAY CANCEL VAR ECHO_INTENSITY_BEAM3_DAY CANCEL VAR ECHO_INTENSITY_BEAM3_FIL !------------------------------------------------------------------ ! Variable ECHO_INTENSITY_BEAM4_HR: 7 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = ECHO_INTENSITY_BEAM4_HR[L=@CIA] + ECHO_INTENSITY_BEAM4_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET ECHO_INTENSITY_BEAM4_FIL ECHO_INTENSITY_BEAM4_HR[L=@FLN] * GAP_MASK ! LET ECHO_INTENSITY_BEAM4_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM4_FIL,W33) LET ECHO_INTENSITY_BEAM4_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM4_HR,W33) LET ECHO_INTENSITY_BEAM4_DAY = ECHO_INTENSITY_BEAM4_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered beam4 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM4_DAY SET MEM/SIZE=30 STAT ECHO_INTENSITY_BEAM4_FIL STAT ECHO_INTENSITY_BEAM4_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM4_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM4_DAY CANCEL VAR ECHO_INTENSITY_BEAM4_DAY CANCEL VAR ECHO_INTENSITY_BEAM4_FIL !------------------------------------------------------------------ ! Variable GPS_LATITUDE_HR: 8 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = GPS_LATITUDE_HR[L=@CIA] + GPS_LATITUDE_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET GPS_LATITUDE_FIL GPS_LATITUDE_HR[L=@FLN] * GAP_MASK ! LET GPS_LATITUDE_FIL = CONVOLVEL(GPS_LATITUDE_FIL,W33) LET GPS_LATITUDE_FIL = CONVOLVEL(GPS_LATITUDE_HR,W33) LET GPS_LATITUDE_DAY = GPS_LATITUDE_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Latitude"/UNITS="Degrees North" GPS_LATITUDE_DAY SET MEM/SIZE=30 STAT GPS_LATITUDE_FIL STAT GPS_LATITUDE_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" GPS_LATITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" GPS_LATITUDE_DAY CANCEL VAR GPS_LATITUDE_DAY CANCEL VAR GPS_LATITUDE_FIL !------------------------------------------------------------------ ! Variable GPS_LONGITUDE_HR: 9 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = GPS_LONGITUDE_HR[L=@CIA] + GPS_LONGITUDE_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET GPS_LONGITUDE_FIL GPS_LONGITUDE_HR[L=@FLN] * GAP_MASK ! LET GPS_LONGITUDE_FIL = CONVOLVEL(GPS_LONGITUDE_FIL,W33) LET GPS_LONGITUDE_FIL = CONVOLVEL(GPS_LONGITUDE_HR,W33) LET GPS_LONGITUDE_DAY = GPS_LONGITUDE_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Longitude"/UNITS="Degrees East" GPS_LONGITUDE_DAY SET MEM/SIZE=30 STAT GPS_LONGITUDE_FIL STAT GPS_LONGITUDE_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" GPS_LONGITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" GPS_LONGITUDE_DAY CANCEL VAR GPS_LONGITUDE_DAY CANCEL VAR GPS_LONGITUDE_FIL !------------------------------------------------------------------ ! Variable PRESSURE_HR: 10 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = PRESSURE_HR[L=@CIA] + PRESSURE_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET PRESSURE_FIL PRESSURE_HR[L=@FLN] * GAP_MASK ! LET PRESSURE_FIL = CONVOLVEL(PRESSURE_FIL,W33) LET PRESSURE_FIL = CONVOLVEL(PRESSURE_HR,W33) LET PRESSURE_DAY = PRESSURE_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Pressure"/UNITS="db" PRESSURE_DAY SET MEM/SIZE=30 STAT PRESSURE_FIL STAT PRESSURE_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" PRESSURE_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" PRESSURE_DAY CANCEL VAR PRESSURE_DAY CANCEL VAR PRESSURE_FIL !------------------------------------------------------------------ ! Variable RELATIVE_HUMIDITY_HR: 11 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = RELATIVE_HUMIDITY_HR[L=@CIA] + RELATIVE_HUMIDITY_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET RELATIVE_HUMIDITY_FIL RELATIVE_HUMIDITY_HR[L=@FLN] * GAP_MASK ! LET RELATIVE_HUMIDITY_FIL = CONVOLVEL(RELATIVE_HUMIDITY_FIL,W33) LET RELATIVE_HUMIDITY_FIL = CONVOLVEL(RELATIVE_HUMIDITY_HR,W33) LET RELATIVE_HUMIDITY_DAY = RELATIVE_HUMIDITY_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Relative Humidity"/UNITS="%" RELATIVE_HUMIDITY_DAY SET MEM/SIZE=30 STAT RELATIVE_HUMIDITY_FIL STAT RELATIVE_HUMIDITY_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" RELATIVE_HUMIDITY_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" RELATIVE_HUMIDITY_DAY CANCEL VAR RELATIVE_HUMIDITY_DAY CANCEL VAR RELATIVE_HUMIDITY_FIL !------------------------------------------------------------------ ! Variable SALINITY_HR: 12 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = SALINITY_HR[L=@CIA] + SALINITY_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET SALINITY_FIL SALINITY_HR[L=@FLN] * GAP_MASK ! LET SALINITY_FIL = CONVOLVEL(SALINITY_FIL,W33) LET SALINITY_FIL = CONVOLVEL(SALINITY_HR,W33) LET SALINITY_DAY = SALINITY_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Salinity"/UNITS="" SALINITY_DAY SET MEM/SIZE=30 STAT SALINITY_FIL STAT SALINITY_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" SALINITY_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" SALINITY_DAY CANCEL VAR SALINITY_DAY CANCEL VAR SALINITY_FIL !------------------------------------------------------------------ ! Variable TEMPERATURE_HR: 13 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = TEMPERATURE_HR[L=@CIA] + TEMPERATURE_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET TEMPERATURE_FIL TEMPERATURE_HR[L=@FLN] * GAP_MASK ! LET TEMPERATURE_FIL = CONVOLVEL(TEMPERATURE_FIL,W33) LET TEMPERATURE_FIL = CONVOLVEL(TEMPERATURE_HR,W33) LET TEMPERATURE_DAY = TEMPERATURE_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered Temperature"/UNITS="Celsius" TEMPERATURE_DAY SET MEM/SIZE=30 STAT TEMPERATURE_FIL STAT TEMPERATURE_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" TEMPERATURE_DAY CANCEL VAR TEMPERATURE_DAY CANCEL VAR TEMPERATURE_FIL !------------------------------------------------------------------ ! Variable U_COMPONENT_UNCORR_HR: 14 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = U_COMPONENT_UNCORR_HR[L=@CIA] + U_COMPONENT_UNCORR_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET U_COMPONENT_UNCORR_FIL U_COMPONENT_UNCORR_HR[L=@FLN] * GAP_MASK ! LET U_COMPONENT_UNCORR_FIL = CONVOLVEL(U_COMPONENT_UNCORR_FIL,W33) LET U_COMPONENT_UNCORR_FIL = CONVOLVEL(U_COMPONENT_UNCORR_HR,W33) LET U_COMPONENT_UNCORR_DAY = U_COMPONENT_UNCORR_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered East-West Velocity Towards(+E/-W) True"/UNITS="cm/s" U_COMPONENT_UNCORR_DAY SET MEM/SIZE=30 STAT U_COMPONENT_UNCORR_FIL STAT U_COMPONENT_UNCORR_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" U_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" U_COMPONENT_UNCORR_DAY CANCEL VAR U_COMPONENT_UNCORR_DAY CANCEL VAR U_COMPONENT_UNCORR_FIL !------------------------------------------------------------------ ! Variable V_COMPONENT_UNCORR_HR: 15 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = V_COMPONENT_UNCORR_HR[L=@CIA] + V_COMPONENT_UNCORR_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET V_COMPONENT_UNCORR_FIL V_COMPONENT_UNCORR_HR[L=@FLN] * GAP_MASK ! LET V_COMPONENT_UNCORR_FIL = CONVOLVEL(V_COMPONENT_UNCORR_FIL,W33) LET V_COMPONENT_UNCORR_FIL = CONVOLVEL(V_COMPONENT_UNCORR_HR,W33) LET V_COMPONENT_UNCORR_DAY = V_COMPONENT_UNCORR_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered North-South Velocity Towards(+N/-S) True"/UNITS="cm/s" V_COMPONENT_UNCORR_DAY SET MEM/SIZE=30 STAT V_COMPONENT_UNCORR_FIL STAT V_COMPONENT_UNCORR_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" V_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" V_COMPONENT_UNCORR_DAY CANCEL VAR V_COMPONENT_UNCORR_DAY CANCEL VAR V_COMPONENT_UNCORR_FIL !------------------------------------------------------------------ ! Variable WIND_U_COMPONENT_HR: 16 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = WIND_U_COMPONENT_HR[L=@CIA] + WIND_U_COMPONENT_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET WIND_U_COMPONENT_FIL WIND_U_COMPONENT_HR[L=@FLN] * GAP_MASK ! LET WIND_U_COMPONENT_FIL = CONVOLVEL(WIND_U_COMPONENT_FIL,W33) LET WIND_U_COMPONENT_FIL = CONVOLVEL(WIND_U_COMPONENT_HR,W33) LET WIND_U_COMPONENT_DAY = WIND_U_COMPONENT_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered U Component of Wind Velocity"/UNITS="m/s (eastward)" WIND_U_COMPONENT_DAY SET MEM/SIZE=30 STAT WIND_U_COMPONENT_FIL STAT WIND_U_COMPONENT_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" WIND_U_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" WIND_U_COMPONENT_DAY CANCEL VAR WIND_U_COMPONENT_DAY CANCEL VAR WIND_U_COMPONENT_FIL !------------------------------------------------------------------ ! Variable WIND_V_COMPONENT_HR: 17 of 17 ! ! These commands need too much memory for the whole deployment ! Save for when this code is converted to processing by deployment ! ! LET GAP_ALLOWED = 2 ! LET GAP_SIZE = WIND_V_COMPONENT_HR[L=@CIA] + WIND_V_COMPONENT_HR[L=@CIB] ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 ! LET WIND_V_COMPONENT_FIL WIND_V_COMPONENT_HR[L=@FLN] * GAP_MASK ! LET WIND_V_COMPONENT_FIL = CONVOLVEL(WIND_V_COMPONENT_FIL,W33) LET WIND_V_COMPONENT_FIL = CONVOLVEL(WIND_V_COMPONENT_HR,W33) LET WIND_V_COMPONENT_DAY = WIND_V_COMPONENT_FIL[GT=TIME_DAY] SET VAR/TITLE="Filtered V Component of Wind Velocity"/UNITS="m/s (northward)" WIND_V_COMPONENT_DAY SET MEM/SIZE=30 STAT WIND_V_COMPONENT_FIL STAT WIND_V_COMPONENT_DAY SAVE/APPEND/FILE="hourlyFilteredM3.nc" WIND_V_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM3.nc" WIND_V_COMPONENT_DAY CANCEL VAR WIND_V_COMPONENT_DAY CANCEL VAR WIND_V_COMPONENT_FIL QUIT
NOAA/PMEL TMAP FERRET v5.80 Solaris 5.8 - 01/03/05 13-Jun-05 10:11 cancel mode journal sp rm -f ferret.jnl yes? ! Description: Create daily filtered OASIS data file from hourly gridded da ta yes? ! Automatically generated by ./hourly2dailyfilter.pl on Mon Jun 13 10:11:14 2005. yes? ! yes? ! All data sets have AIR_TEMPERATURE_HR variable; yes? ! use it to get axis time limits. 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? ! Also save an intermediate hourly filtered file for finer scale use. yes? ! yes? SET MEMORY/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? SET MODE DESPERATE 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 VARIABLE : 33 Hour Filter Weights FILENAME : 33hrfilt SUBSET : 67 points (X) -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.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM3.n c" yes? USE "hourlyM3.nc" yes? DEFINE SYMBOL TBEG="`AIR_TEMPERATURE_HR,return=tstart`" !- DEFINE SYMBOL TBEG="24-MAR-1998 00:00" yes? DEFINE SYMBOL TEND="`AIR_TEMPERATURE_HR,return=tend`" !- DEFINE SYMBOL TEND="08-MAR-2000 00:00" yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY !- DEFINE AXIS/T="24-MAR-1998 00:00":"08-MAR-2000 00:00":1/unit=day TIME_DAY yes? SPAWN "/bin/rm -f hourlyFilteredM3.nc" yes? SPAWN "/bin/rm -f dailyFilteredM3.nc" yes? LET W33=XSEQUENCE(WEIGHTS33[D=1]) yes? yes? !------------------------------------------------------------------ yes? ! Variable AIR_PRESS_HR: 1 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = AIR_PRESS_HR[L=@CIA] + AIR_PRESS_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET AIR_PRESS_FIL AIR_PRESS_HR[L=@FLN] * GAP_MASK yes? ! LET AIR_PRESS_FIL = CONVOLVEL(AIR_PRESS_FIL,W33) yes? yes? LET AIR_PRESS_FIL = CONVOLVEL(AIR_PRESS_HR,W33) yes? LET AIR_PRESS_DAY = AIR_PRESS_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Air Pressure"/UNITS="mbar" AIR_PRESS_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT AIR_PRESS_FIL CONVOLVEL(AIR_PRESS_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 9340 Minimum value: 998.61 Maximum value: 1030.8 Mean value: 1016.8 (unweighted average) Standard deviation: 4.8351 yes? STAT AIR_PRESS_DAY Filtered Air Pressure LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 389 Minimum value: 998.64 Maximum value: 1030.8 Mean value: 1016.8 (unweighted average) Standard deviation: 4.8819 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" AIR_PRESS_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" AIR_PRESS_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR AIR_PRESS_DAY yes? CANCEL VAR AIR_PRESS_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable AIR_TEMPERATURE_HR: 2 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = AIR_TEMPERATURE_HR[L=@CIA] + AIR_TEMPERATURE_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET AIR_TEMPERATURE_FIL AIR_TEMPERATURE_HR[L=@FLN] * GAP_MASK yes? ! LET AIR_TEMPERATURE_FIL = CONVOLVEL(AIR_TEMPERATURE_FIL,W33) yes? yes? LET AIR_TEMPERATURE_FIL = CONVOLVEL(AIR_TEMPERATURE_HR,W33) yes? LET AIR_TEMPERATURE_DAY = AIR_TEMPERATURE_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered air_temperature"/UNITS="Celsius" AIR_TEMPERATURE_DA Y yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT AIR_TEMPERATURE_FIL CONVOLVEL(AIR_TEMPERATURE_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 9229 Minimum value: 8.4178 Maximum value: 16.862 Mean value: 13.152 (unweighted average) Standard deviation: 1.6697 yes? STAT AIR_TEMPERATURE_DAY Filtered air_temperature LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 389 Minimum value: 8.5768 Maximum value: 16.849 Mean value: 13.161 (unweighted average) Standard deviation: 1.6763 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" AIR_TEMPERATURE_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" AIR_TEMPERATURE_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR AIR_TEMPERATURE_DAY yes? CANCEL VAR AIR_TEMPERATURE_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable CONDUCTIVITY_HR: 3 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = CONDUCTIVITY_HR[L=@CIA] + CONDUCTIVITY_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET CONDUCTIVITY_FIL CONDUCTIVITY_HR[L=@FLN] * GAP_MASK yes? ! LET CONDUCTIVITY_FIL = CONVOLVEL(CONDUCTIVITY_FIL,W33) yes? yes? LET CONDUCTIVITY_FIL = CONVOLVEL(CONDUCTIVITY_HR,W33) yes? LET CONDUCTIVITY_DAY = CONDUCTIVITY_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Conductivity"/UNITS="Siemens/m" CONDUCTIVITY_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT CONDUCTIVITY_FIL CONVOLVEL(CONDUCTIVITY_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 205932 (1*1*12*17161) # flagged as bad data: 190146 Minimum value: -41.988 Maximum value: 7.9241 Mean value: 3.7396 (unweighted average) Standard deviation: 2.4242 yes? STAT CONDUCTIVITY_DAY Filtered Conductivity LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 8592 (1*1*12*716) # flagged as bad data: 7933 Minimum value: -24.632 Maximum value: 5.7579 Mean value: 3.7513 (unweighted average) Standard deviation: 2.0131 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" CONDUCTIVITY_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" CONDUCTIVITY_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR CONDUCTIVITY_DAY yes? CANCEL VAR CONDUCTIVITY_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable ECHO_INTENSITY_BEAM1_HR: 4 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = ECHO_INTENSITY_BEAM1_HR[L=@CIA] + ECHO_INTENSITY_BEAM1_HR[ L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET ECHO_INTENSITY_BEAM1_FIL ECHO_INTENSITY_BEAM1_HR[L=@FLN] * GAP_MASK yes? ! LET ECHO_INTENSITY_BEAM1_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM1_FIL,W33) yes? yes? LET ECHO_INTENSITY_BEAM1_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM1_HR,W33) yes? LET ECHO_INTENSITY_BEAM1_DAY = ECHO_INTENSITY_BEAM1_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered beam1 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM1_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT ECHO_INTENSITY_BEAM1_FIL CONVOLVEL(ECHO_INTENSITY_BEAM1_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1373078 Minimum value: 39.869 Maximum value: 173.78 Mean value: 79.056 (unweighted average) Standard deviation: 32.045 yes? STAT ECHO_INTENSITY_BEAM1_DAY Filtered beam1 Echo Intensity data per beam LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 57317 Minimum value: 39.967 Maximum value: 169.9 Mean value: 78.985 (unweighted average) Standard deviation: 32.031 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM1_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM1_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR ECHO_INTENSITY_BEAM1_DAY yes? CANCEL VAR ECHO_INTENSITY_BEAM1_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable ECHO_INTENSITY_BEAM2_HR: 5 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = ECHO_INTENSITY_BEAM2_HR[L=@CIA] + ECHO_INTENSITY_BEAM2_HR[ L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET ECHO_INTENSITY_BEAM2_FIL ECHO_INTENSITY_BEAM2_HR[L=@FLN] * GAP_MASK yes? ! LET ECHO_INTENSITY_BEAM2_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM2_FIL,W33) yes? yes? LET ECHO_INTENSITY_BEAM2_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM2_HR,W33) yes? LET ECHO_INTENSITY_BEAM2_DAY = ECHO_INTENSITY_BEAM2_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered beam2 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM2_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT ECHO_INTENSITY_BEAM2_FIL CONVOLVEL(ECHO_INTENSITY_BEAM2_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1373078 Minimum value: 41.97 Maximum value: 175.82 Mean value: 78.637 (unweighted average) Standard deviation: 32.81 yes? STAT ECHO_INTENSITY_BEAM2_DAY Filtered beam2 Echo Intensity data per beam LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 57317 Minimum value: 42.075 Maximum value: 173.48 Mean value: 78.566 (unweighted average) Standard deviation: 32.794 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM2_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM2_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR ECHO_INTENSITY_BEAM2_DAY yes? CANCEL VAR ECHO_INTENSITY_BEAM2_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable ECHO_INTENSITY_BEAM3_HR: 6 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = ECHO_INTENSITY_BEAM3_HR[L=@CIA] + ECHO_INTENSITY_BEAM3_HR[ L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET ECHO_INTENSITY_BEAM3_FIL ECHO_INTENSITY_BEAM3_HR[L=@FLN] * GAP_MASK yes? ! LET ECHO_INTENSITY_BEAM3_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM3_FIL,W33) yes? yes? LET ECHO_INTENSITY_BEAM3_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM3_HR,W33) yes? LET ECHO_INTENSITY_BEAM3_DAY = ECHO_INTENSITY_BEAM3_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered beam3 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM3_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT ECHO_INTENSITY_BEAM3_FIL CONVOLVEL(ECHO_INTENSITY_BEAM3_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1373078 Minimum value: 44.94 Maximum value: 177.47 Mean value: 80.74 (unweighted average) Standard deviation: 31.982 yes? STAT ECHO_INTENSITY_BEAM3_DAY Filtered beam3 Echo Intensity data per beam LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 57317 Minimum value: 45.003 Maximum value: 173.53 Mean value: 80.671 (unweighted average) Standard deviation: 31.968 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM3_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM3_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR ECHO_INTENSITY_BEAM3_DAY yes? CANCEL VAR ECHO_INTENSITY_BEAM3_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable ECHO_INTENSITY_BEAM4_HR: 7 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = ECHO_INTENSITY_BEAM4_HR[L=@CIA] + ECHO_INTENSITY_BEAM4_HR[ L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET ECHO_INTENSITY_BEAM4_FIL ECHO_INTENSITY_BEAM4_HR[L=@FLN] * GAP_MASK yes? ! LET ECHO_INTENSITY_BEAM4_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM4_FIL,W33) yes? yes? LET ECHO_INTENSITY_BEAM4_FIL = CONVOLVEL(ECHO_INTENSITY_BEAM4_HR,W33) yes? LET ECHO_INTENSITY_BEAM4_DAY = ECHO_INTENSITY_BEAM4_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered beam4 Echo Intensity data per beam"/UNITS="Counts" ECHO_INTENSITY_BEAM4_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT ECHO_INTENSITY_BEAM4_FIL CONVOLVEL(ECHO_INTENSITY_BEAM4_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1373078 Minimum value: 43.189 Maximum value: 181.68 Mean value: 81.606 (unweighted average) Standard deviation: 33.146 yes? STAT ECHO_INTENSITY_BEAM4_DAY Filtered beam4 Echo Intensity data per beam LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 57317 Minimum value: 43.189 Maximum value: 177.02 Mean value: 81.532 (unweighted average) Standard deviation: 33.133 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" ECHO_INTENSITY_BEAM4_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" ECHO_INTENSITY_BEAM4_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR ECHO_INTENSITY_BEAM4_DAY yes? CANCEL VAR ECHO_INTENSITY_BEAM4_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable GPS_LATITUDE_HR: 8 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = GPS_LATITUDE_HR[L=@CIA] + GPS_LATITUDE_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET GPS_LATITUDE_FIL GPS_LATITUDE_HR[L=@FLN] * GAP_MASK yes? ! LET GPS_LATITUDE_FIL = CONVOLVEL(GPS_LATITUDE_FIL,W33) yes? yes? LET GPS_LATITUDE_FIL = CONVOLVEL(GPS_LATITUDE_HR,W33) yes? LET GPS_LATITUDE_DAY = GPS_LATITUDE_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Latitude"/UNITS="Degrees North" GPS_LATITUDE_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT GPS_LATITUDE_FIL CONVOLVEL(GPS_LATITUDE_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (decibars): 5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 10586 Minimum value: 36.546 Maximum value: 36.588 Mean value: 36.562 (unweighted average) Standard deviation: 0.011813 yes? STAT GPS_LATITUDE_DAY Filtered Latitude LONGITUDE: 123W LATITUDE: 36.6N DEPTH (decibars): 5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 442 Minimum value: 36.546 Maximum value: 36.585 Mean value: 36.562 (unweighted average) Standard deviation: 0.011838 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" GPS_LATITUDE_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" GPS_LATITUDE_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR GPS_LATITUDE_DAY yes? CANCEL VAR GPS_LATITUDE_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable GPS_LONGITUDE_HR: 9 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = GPS_LONGITUDE_HR[L=@CIA] + GPS_LONGITUDE_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET GPS_LONGITUDE_FIL GPS_LONGITUDE_HR[L=@FLN] * GAP_MASK yes? ! LET GPS_LONGITUDE_FIL = CONVOLVEL(GPS_LONGITUDE_FIL,W33) yes? yes? LET GPS_LONGITUDE_FIL = CONVOLVEL(GPS_LONGITUDE_HR,W33) yes? LET GPS_LONGITUDE_DAY = GPS_LONGITUDE_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Longitude"/UNITS="Degrees East" GPS_LONGITUDE_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT GPS_LONGITUDE_FIL CONVOLVEL(GPS_LONGITUDE_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (decibars): 5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 10586 Minimum value: -122.98 Maximum value: -122.93 Mean value: -122.95 (unweighted average) Standard deviation: 0.015021 yes? STAT GPS_LONGITUDE_DAY Filtered Longitude LONGITUDE: 123W LATITUDE: 36.6N DEPTH (decibars): 5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 442 Minimum value: -122.98 Maximum value: -122.93 Mean value: -122.95 (unweighted average) Standard deviation: 0.015054 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" GPS_LONGITUDE_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" GPS_LONGITUDE_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR GPS_LONGITUDE_DAY yes? CANCEL VAR GPS_LONGITUDE_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable PRESSURE_HR: 10 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = PRESSURE_HR[L=@CIA] + PRESSURE_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET PRESSURE_FIL PRESSURE_HR[L=@FLN] * GAP_MASK yes? ! LET PRESSURE_FIL = CONVOLVEL(PRESSURE_FIL,W33) yes? yes? LET PRESSURE_FIL = CONVOLVEL(PRESSURE_HR,W33) yes? LET PRESSURE_DAY = PRESSURE_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Pressure"/UNITS="db" PRESSURE_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT PRESSURE_FIL CONVOLVEL(PRESSURE_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 205932 (1*1*12*17161) # flagged as bad data: 134612 Minimum value: 0.64742 Maximum value: 300 Mean value: 96.819 (unweighted average) Standard deviation: 97.796 yes? STAT PRESSURE_DAY Filtered Pressure LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 8592 (1*1*12*716) # flagged as bad data: 5624 Minimum value: 0.65556 Maximum value: 300 Mean value: 96.752 (unweighted average) Standard deviation: 97.831 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" PRESSURE_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" PRESSURE_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR PRESSURE_DAY yes? CANCEL VAR PRESSURE_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable RELATIVE_HUMIDITY_HR: 11 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = RELATIVE_HUMIDITY_HR[L=@CIA] + RELATIVE_HUMIDITY_HR[L=@CIB ] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET RELATIVE_HUMIDITY_FIL RELATIVE_HUMIDITY_HR[L=@FLN] * GAP_MASK yes? ! LET RELATIVE_HUMIDITY_FIL = CONVOLVEL(RELATIVE_HUMIDITY_FIL,W33) yes? yes? LET RELATIVE_HUMIDITY_FIL = CONVOLVEL(RELATIVE_HUMIDITY_HR,W33) yes? LET RELATIVE_HUMIDITY_DAY = RELATIVE_HUMIDITY_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Relative Humidity"/UNITS="%" RELATIVE_HUMIDITY_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT RELATIVE_HUMIDITY_FIL CONVOLVEL(RELATIVE_HUMIDITY_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 16176 Minimum value: 72.474 Maximum value: 100.18 Mean value: 94.371 (unweighted average) Standard deviation: 5.6885 yes? STAT RELATIVE_HUMIDITY_DAY Filtered Relative Humidity LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 675 Minimum value: 77.875 Maximum value: 100.18 Mean value: 94.313 (unweighted average) Standard deviation: 5.6449 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" RELATIVE_HUMIDITY_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" RELATIVE_HUMIDITY_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR RELATIVE_HUMIDITY_DAY yes? CANCEL VAR RELATIVE_HUMIDITY_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable SALINITY_HR: 12 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = SALINITY_HR[L=@CIA] + SALINITY_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET SALINITY_FIL SALINITY_HR[L=@FLN] * GAP_MASK yes? ! LET SALINITY_FIL = CONVOLVEL(SALINITY_FIL,W33) yes? yes? LET SALINITY_FIL = CONVOLVEL(SALINITY_HR,W33) yes? LET SALINITY_DAY = SALINITY_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Salinity"/UNITS="" SALINITY_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT SALINITY_FIL CONVOLVEL(SALINITY_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 205932 (1*1*12*17161) # flagged as bad data: 190590 Minimum value: 31.94 Maximum value: 34.025 Mean value: 33.289 (unweighted average) Standard deviation: 0.30064 yes? STAT SALINITY_DAY Filtered Salinity LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 8592 (1*1*12*716) # flagged as bad data: 7953 Minimum value: 32.009 Maximum value: 34.02 Mean value: 33.289 (unweighted average) Standard deviation: 0.30156 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" SALINITY_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" SALINITY_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR SALINITY_DAY yes? CANCEL VAR SALINITY_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable TEMPERATURE_HR: 13 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = TEMPERATURE_HR[L=@CIA] + TEMPERATURE_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET TEMPERATURE_FIL TEMPERATURE_HR[L=@FLN] * GAP_MASK yes? ! LET TEMPERATURE_FIL = CONVOLVEL(TEMPERATURE_FIL,W33) yes? yes? LET TEMPERATURE_FIL = CONVOLVEL(TEMPERATURE_HR,W33) yes? LET TEMPERATURE_DAY = TEMPERATURE_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered Temperature"/UNITS="Celsius" TEMPERATURE_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT TEMPERATURE_FIL CONVOLVEL(TEMPERATURE_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 205932 (1*1*12*17161) # flagged as bad data: 136334 Minimum value: 6.3830 Maximum value: 17.571 Mean value: 10.965 (unweighted average) Standard deviation: 2.6154 yes? STAT TEMPERATURE_DAY Filtered Temperature LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 to 375 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 8592 (1*1*12*716) # flagged as bad data: 5696 Minimum value: 6.3839 Maximum value: 17.475 Mean value: 10.969 (unweighted average) Standard deviation: 2.6187 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" TEMPERATURE_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" TEMPERATURE_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR TEMPERATURE_DAY yes? CANCEL VAR TEMPERATURE_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable U_COMPONENT_UNCORR_HR: 14 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = U_COMPONENT_UNCORR_HR[L=@CIA] + U_COMPONENT_UNCORR_HR[L=@C IB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET U_COMPONENT_UNCORR_FIL U_COMPONENT_UNCORR_HR[L=@FLN] * GAP_MASK yes? ! LET U_COMPONENT_UNCORR_FIL = CONVOLVEL(U_COMPONENT_UNCORR_FIL,W33) yes? yes? LET U_COMPONENT_UNCORR_FIL = CONVOLVEL(U_COMPONENT_UNCORR_HR,W33) yes? LET U_COMPONENT_UNCORR_DAY = U_COMPONENT_UNCORR_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered East-West Velocity Towards(+E/-W) True"/UNITS="cm/s " U_COMPONENT_UNCORR_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT U_COMPONENT_UNCORR_FIL CONVOLVEL(U_COMPONENT_UNCORR_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1503154 Minimum value: -44.367 Maximum value: 40.353 Mean value: -1.7562 (unweighted average) Standard deviation: 11.506 yes? STAT U_COMPONENT_UNCORR_DAY Filtered East-West Velocity Towards(+E/-W) True LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 62743 Minimum value: -42.893 Maximum value: 39.954 Mean value: -1.7598 (unweighted average) Standard deviation: 11.49 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" U_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" U_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR U_COMPONENT_UNCORR_DAY yes? CANCEL VAR U_COMPONENT_UNCORR_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable V_COMPONENT_UNCORR_HR: 15 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = V_COMPONENT_UNCORR_HR[L=@CIA] + V_COMPONENT_UNCORR_HR[L=@C IB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET V_COMPONENT_UNCORR_FIL V_COMPONENT_UNCORR_HR[L=@FLN] * GAP_MASK yes? ! LET V_COMPONENT_UNCORR_FIL = CONVOLVEL(V_COMPONENT_UNCORR_FIL,W33) yes? yes? LET V_COMPONENT_UNCORR_FIL = CONVOLVEL(V_COMPONENT_UNCORR_HR,W33) yes? LET V_COMPONENT_UNCORR_DAY = V_COMPONENT_UNCORR_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered North-South Velocity Towards(+N/-S) True"/UNITS="cm /s" V_COMPONENT_UNCORR_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT V_COMPONENT_UNCORR_FIL CONVOLVEL(V_COMPONENT_UNCORR_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 1716100 (1*1*100*17161) # flagged as bad data: 1503154 Minimum value: -55.421 Maximum value: 28.341 Mean value: -1.4823 (unweighted average) Standard deviation: 11.443 yes? STAT V_COMPONENT_UNCORR_DAY Filtered North-South Velocity Towards(+N/-S) True LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): 2.5 to 502.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 71600 (1*1*100*716) # flagged as bad data: 62743 Minimum value: -53.118 Maximum value: 26.925 Mean value: -1.4727 (unweighted average) Standard deviation: 11.438 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" V_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" V_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR V_COMPONENT_UNCORR_DAY yes? CANCEL VAR V_COMPONENT_UNCORR_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable WIND_U_COMPONENT_HR: 16 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = WIND_U_COMPONENT_HR[L=@CIA] + WIND_U_COMPONENT_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET WIND_U_COMPONENT_FIL WIND_U_COMPONENT_HR[L=@FLN] * GAP_MASK yes? ! LET WIND_U_COMPONENT_FIL = CONVOLVEL(WIND_U_COMPONENT_FIL,W33) yes? yes? LET WIND_U_COMPONENT_FIL = CONVOLVEL(WIND_U_COMPONENT_HR,W33) yes? LET WIND_U_COMPONENT_DAY = WIND_U_COMPONENT_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered U Component of Wind Velocity"/UNITS="m/s (eastward) " WIND_U_COMPONENT_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT WIND_U_COMPONENT_FIL CONVOLVEL(WIND_U_COMPONENT_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 2969 Minimum value: -10.199 Maximum value: 13.108 Mean value: 3.0511 (unweighted average) Standard deviation: 3.4158 yes? STAT WIND_U_COMPONENT_DAY Filtered U Component of Wind Velocity LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 126 Minimum value: -9.9948 Maximum value: 12.943 Mean value: 3.0768 (unweighted average) Standard deviation: 3.4505 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" WIND_U_COMPONENT_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" WIND_U_COMPONENT_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR WIND_U_COMPONENT_DAY yes? CANCEL VAR WIND_U_COMPONENT_FIL yes? yes? !------------------------------------------------------------------ yes? ! Variable WIND_V_COMPONENT_HR: 17 of 17 yes? ! yes? ! These commands need too much memory for the whole deployment yes? ! Save for when this code is converted to processing by deployment yes? ! yes? ! LET GAP_ALLOWED = 2 yes? ! LET GAP_SIZE = WIND_V_COMPONENT_HR[L=@CIA] + WIND_V_COMPONENT_HR[L=@CIB] yes? ! LET GAP_MASK = IF GAP_SIZE LE GAP_ALLOWED THEN 1 yes? ! LET WIND_V_COMPONENT_FIL WIND_V_COMPONENT_HR[L=@FLN] * GAP_MASK yes? ! LET WIND_V_COMPONENT_FIL = CONVOLVEL(WIND_V_COMPONENT_FIL,W33) yes? yes? LET WIND_V_COMPONENT_FIL = CONVOLVEL(WIND_V_COMPONENT_HR,W33) yes? LET WIND_V_COMPONENT_DAY = WIND_V_COMPONENT_FIL[GT=TIME_DAY] yes? SET VAR/TITLE="Filtered V Component of Wind Velocity"/UNITS="m/s (northward )" WIND_V_COMPONENT_DAY yes? SET MEM/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? STAT WIND_V_COMPONENT_FIL CONVOLVEL(WIND_V_COMPONENT_HR,W33) LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 23:30 to 08-MAR-2000 00:30 DATA SET: ./hourlyM3.nc Total # of data points: 17161 (1*1*1*17161) # flagged as bad data: 2969 Minimum value: -12.766 Maximum value: 10.6 Mean value: -2.8968 (unweighted average) Standard deviation: 4.4369 yes? STAT WIND_V_COMPONENT_DAY Filtered V Component of Wind Velocity LONGITUDE: 123W LATITUDE: 36.6N DEPTH (m): -3.5 TIME: 23-MAR-1998 12:00 to 08-MAR-2000 12:00 DATA SET: ./hourlyM3.nc Total # of data points: 716 (1*1*1*716) # flagged as bad data: 126 Minimum value: -12.764 Maximum value: 10.471 Mean value: -2.9184 (unweighted average) Standard deviation: 4.4225 yes? SAVE/APPEND/FILE="hourlyFilteredM3.nc" WIND_V_COMPONENT_FIL LISTing to file hourlyFilteredM3.nc yes? SAVE/APPEND/FILE="dailyFilteredM3.nc" WIND_V_COMPONENT_DAY LISTing to file dailyFilteredM3.nc yes? CANCEL VAR WIND_V_COMPONENT_DAY yes? CANCEL VAR WIND_V_COMPONENT_FIL yes? 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.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...