The data in this data set are available via DODS/OPeNDAP
at
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM0.nc
and through the DODS html form mechanism:
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM0.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/hourlyFilteredM0.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/hourlyM0.nc).
For information on the individual deployments and instruments that contributed to this data set please see the M0 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 Tue Mar 13 11:54:47 2007. ! ! 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/hourlyM0.nc" USE "hourlyM0.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 hourlyFilteredM0.nc" SPAWN "/bin/rm -f dailyFilteredM0.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="hourlyFilteredM0.nc" AIR_PRESS_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" AIR_TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" CONDUCTIVITY_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM1_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM2_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM3_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM4_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" GPS_LATITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" GPS_LONGITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" PRESSURE_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" RELATIVE_HUMIDITY_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" SALINITY_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" U_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" V_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" WIND_U_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM0.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="hourlyFilteredM0.nc" WIND_V_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM0.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-Mar-07 11:54 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 Tue Mar 13 11:54:47 2007. 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/hourlyM0.n c" yes? USE "hourlyM0.nc" yes? DEFINE SYMBOL TBEG="`AIR_TEMPERATURE_HR,return=tstart`" !- DEFINE SYMBOL TBEG="03-JUN-2004 18:00" yes? DEFINE SYMBOL TEND="`AIR_TEMPERATURE_HR,return=tend`" !- DEFINE SYMBOL TEND="31-JUL-2006 14:00" yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY !- DEFINE AXIS/T="03-JUN-2004 18:00":"31-JUL-2006 14:00":1/unit=day TIME_DAY yes? SPAWN "/bin/rm -f hourlyFilteredM0.nc" yes? SPAWN "/bin/rm -f dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 3436 Minimum value: 996.09 Maximum value: 1031.3 Mean value: 1016.3 (unweighted average) Standard deviation: 4.5486 yes? STAT AIR_PRESS_DAY Filtered Air Pressure LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 148 Minimum value: 996.09 Maximum value: 1031 Mean value: 1016.3 (unweighted average) Standard deviation: 4.5788 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" AIR_PRESS_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" AIR_PRESS_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 3234 Minimum value: 5.5646 Maximum value: 16.509 Mean value: 11.812 (unweighted average) Standard deviation: 1.715 yes? STAT AIR_TEMPERATURE_DAY Filtered air_temperature LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 139 Minimum value: 5.8176 Maximum value: 16.502 Mean value: 11.818 (unweighted average) Standard deviation: 1.7266 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" AIR_TEMPERATURE_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" AIR_TEMPERATURE_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 94545 (1*1*5*18909) # flagged as bad data: 8789 Minimum value: 1.7809 Maximum value: 4.3726 Mean value: 3.8226 (unweighted average) Standard deviation: 0.15156 yes? STAT CONDUCTIVITY_DAY Filtered Conductivity LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 3945 (1*1*5*789) # flagged as bad data: 371 Minimum value: 1.8082 Maximum value: 4.3652 Mean value: 3.8225 (unweighted average) Standard deviation: 0.1519 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" CONDUCTIVITY_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" CONDUCTIVITY_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 39489 Minimum value: 55.166 Maximum value: 173.96 Mean value: 107.98 (unweighted average) Standard deviation: 22.86 yes? STAT ECHO_INTENSITY_BEAM1_DAY Filtered beam1 Echo Intensity data per beam LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1669 Minimum value: 55.621 Maximum value: 173.96 Mean value: 107.87 (unweighted average) Standard deviation: 22.804 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM1_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" ECHO_INTENSITY_BEAM1_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 39489 Minimum value: 52.369 Maximum value: 177.81 Mean value: 109.37 (unweighted average) Standard deviation: 23.848 yes? STAT ECHO_INTENSITY_BEAM2_DAY Filtered beam2 Echo Intensity data per beam LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1669 Minimum value: 52.793 Maximum value: 177.81 Mean value: 109.25 (unweighted average) Standard deviation: 23.783 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM2_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" ECHO_INTENSITY_BEAM2_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 39489 Minimum value: 55.92 Maximum value: 177.89 Mean value: 111.96 (unweighted average) Standard deviation: 23.079 yes? STAT ECHO_INTENSITY_BEAM3_DAY Filtered beam3 Echo Intensity data per beam LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1669 Minimum value: 55.92 Maximum value: 177.89 Mean value: 111.84 (unweighted average) Standard deviation: 23.02 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM3_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" ECHO_INTENSITY_BEAM3_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 39489 Minimum value: 52.684 Maximum value: 164.85 Mean value: 101.28 (unweighted average) Standard deviation: 22.663 yes? STAT ECHO_INTENSITY_BEAM4_DAY Filtered beam4 Echo Intensity data per beam LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1669 Minimum value: 52.71 Maximum value: 163.61 Mean value: 101.12 (unweighted average) Standard deviation: 22.597 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" ECHO_INTENSITY_BEAM4_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" ECHO_INTENSITY_BEAM4_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (decibars): 1 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 8730 Minimum value: 36.834 Maximum value: 36.835 Mean value: 36.835 (unweighted average) Standard deviation: 0.00012712 yes? STAT GPS_LATITUDE_DAY Filtered Latitude LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (decibars): 1 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 367 Minimum value: 36.834 Maximum value: 36.835 Mean value: 36.835 (unweighted average) Standard deviation: 0.0001248 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" GPS_LATITUDE_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" GPS_LATITUDE_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (decibars): 1 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 8730 Minimum value: -121.9 Maximum value: -121.9 Mean value: -121.9 (unweighted average) Standard deviation: 0.00026175 yes? STAT GPS_LONGITUDE_DAY Filtered Longitude LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (decibars): 1 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 367 Minimum value: -121.9 Maximum value: -121.9 Mean value: -121.9 (unweighted average) Standard deviation: 0.00026279 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" GPS_LONGITUDE_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" GPS_LONGITUDE_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 94545 (1*1*5*18909) # flagged as bad data: 6625 Minimum value: -0.21543 Maximum value: 54.996 Mean value: 25.09 (unweighted average) Standard deviation: 19.604 yes? STAT PRESSURE_DAY Filtered Pressure LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 3945 (1*1*5*789) # flagged as bad data: 283 Minimum value: -0.10174 Maximum value: 54.991 Mean value: 25.102 (unweighted average) Standard deviation: 19.604 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" PRESSURE_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" PRESSURE_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 3843 Minimum value: 34.111 Maximum value: 100.01 Mean value: 87.007 (unweighted average) Standard deviation: 8.95 yes? STAT RELATIVE_HUMIDITY_DAY Filtered Relative Humidity LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 164 Minimum value: 35.802 Maximum value: 99.835 Mean value: 86.913 (unweighted average) Standard deviation: 8.9668 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" RELATIVE_HUMIDITY_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" RELATIVE_HUMIDITY_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 94545 (1*1*5*18909) # flagged as bad data: 13365 Minimum value: 31.099 Maximum value: 34.033 Mean value: 33.49 (unweighted average) Standard deviation: 0.31149 yes? STAT SALINITY_DAY Filtered Salinity LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 3945 (1*1*5*789) # flagged as bad data: 561 Minimum value: 31.242 Maximum value: 34.033 Mean value: 33.489 (unweighted average) Standard deviation: 0.31298 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" SALINITY_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" SALINITY_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 94545 (1*1*5*18909) # flagged as bad data: 7362 Minimum value: 8.7235 Maximum value: 17.912 Mean value: 11.825 (unweighted average) Standard deviation: 1.6022 yes? STAT TEMPERATURE_DAY Filtered Temperature LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 to 62.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 3945 (1*1*5*789) # flagged as bad data: 315 Minimum value: 8.7601 Maximum value: 17.83 Mean value: 11.825 (unweighted average) Standard deviation: 1.6011 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" TEMPERATURE_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" TEMPERATURE_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 41072 Minimum value: -23.15 Maximum value: 32.284 Mean value: -0.56269 (unweighted average) Standard deviation: 4.4035 yes? STAT U_COMPONENT_UNCORR_DAY Filtered East-West Velocity Towards(+E/-W) True LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1738 Minimum value: -20.665 Maximum value: 29.061 Mean value: -0.58658 (unweighted average) Standard deviation: 4.4142 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" U_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" U_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 207999 (1*1*11*18909) # flagged as bad data: 41072 Minimum value: -23.989 Maximum value: 27.767 Mean value: 1.3642 (unweighted average) Standard deviation: 4.5984 yes? STAT V_COMPONENT_UNCORR_DAY Filtered North-South Velocity Towards(+N/-S) True LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): 2.5 to 57.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 8679 (1*1*11*789) # flagged as bad data: 1738 Minimum value: -21.145 Maximum value: 27.485 Mean value: 1.3814 (unweighted average) Standard deviation: 4.5912 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" V_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" V_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 1245 Minimum value: -11.385 Maximum value: 12.963 Mean value: 0.90428 (unweighted average) Standard deviation: 2.6751 yes? STAT WIND_U_COMPONENT_DAY Filtered U Component of Wind Velocity LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 53 Minimum value: -11.226 Maximum value: 11.934 Mean value: 0.92266 (unweighted average) Standard deviation: 2.7139 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" WIND_U_COMPONENT_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" WIND_U_COMPONENT_DAY LISTing to file dailyFilteredM0.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: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 17:30 to 31-JUL-2006 14:30 DATA SET: ./hourlyM0.nc Total # of data points: 18909 (1*1*1*18909) # flagged as bad data: 1245 Minimum value: -8.7828 Maximum value: 10.593 Mean value: 0.19672 (unweighted average) Standard deviation: 2.0425 yes? STAT WIND_V_COMPONENT_DAY Filtered V Component of Wind Velocity LONGITUDE: 121.9W LATITUDE: 36.8N DEPTH (m): -3.5 TIME: 03-JUN-2004 06:00 to 01-AUG-2006 06:00 DATA SET: ./hourlyM0.nc Total # of data points: 789 (1*1*1*789) # flagged as bad data: 53 Minimum value: -8.5149 Maximum value: 9.8627 Mean value: 0.20511 (unweighted average) Standard deviation: 2.0116 yes? SAVE/APPEND/FILE="hourlyFilteredM0.nc" WIND_V_COMPONENT_FIL LISTing to file hourlyFilteredM0.nc yes? SAVE/APPEND/FILE="dailyFilteredM0.nc" WIND_V_COMPONENT_DAY LISTing to file dailyFilteredM0.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...