The data in this data set are available via DODS/OPeNDAP
at
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM2.nc
and through the DODS html form mechanism:
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM2.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/hourlyFilteredM2.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/hourlyM2.nc).
For information on the individual deployments and instruments that contributed to this data set please see the M2 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 Thu Aug 27 09:01:38 2009. ! ! 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/hourlyM2.nc" USE "hourlyM2.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 hourlyFilteredM2.nc" SPAWN "/bin/rm -f dailyFilteredM2.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="hourlyFilteredM2.nc" AIR_PRESS_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" AIR_TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" CONDUCTIVITY_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM1_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM2_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM3_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM4_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" GPS_LATITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" GPS_LONGITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" PRESSURE_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" RELATIVE_HUMIDITY_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" SALINITY_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" U_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" V_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" WIND_U_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM2.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="hourlyFilteredM2.nc" WIND_V_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM2.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 27-Aug-09 09:01 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 Thu Aug 27 09:01:38 2009. 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/hourlyM2.n c" yes? USE "hourlyM2.nc" yes? DEFINE SYMBOL TBEG="`AIR_TEMPERATURE_HR,return=tstart`" !- DEFINE SYMBOL TBEG="24-SEP-1992 05:00" yes? DEFINE SYMBOL TEND="`AIR_TEMPERATURE_HR,return=tend`" !- DEFINE SYMBOL TEND="29-APR-2009 19:00" yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY !- DEFINE AXIS/T="24-SEP-1992 05:00":"29-APR-2009 19:00":1/unit=day TIME_DAY yes? SPAWN "/bin/rm -f hourlyFilteredM2.nc" yes? SPAWN "/bin/rm -f dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 66047 Minimum value: 995.23 Maximum value: 1033.8 Mean value: 1016.6 (unweighted average) Standard deviation: 4.5777 yes? STAT AIR_PRESS_DAY Filtered Air Pressure LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 2747 Minimum value: 999.34 Maximum value: 1033.8 Mean value: 1016.6 (unweighted average) Standard deviation: 4.5637 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" AIR_PRESS_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" AIR_PRESS_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 41865 Minimum value: 6.507 Maximum value: 18.089 Mean value: 12.592 (unweighted average) Standard deviation: 1.615 yes? STAT AIR_TEMPERATURE_DAY Filtered air_temperature LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 1736 Minimum value: 6.5388 Maximum value: 17.891 Mean value: 12.598 (unweighted average) Standard deviation: 1.6126 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" AIR_TEMPERATURE_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" AIR_TEMPERATURE_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 1745748 (1*1*12*145479) # flagged as bad data: 707927 Minimum value: 3.4003 Maximum value: 4.5134 Mean value: 3.7253 (unweighted average) Standard deviation: 0.16675 yes? STAT CONDUCTIVITY_DAY Filtered Conductivity LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 72756 (1*1*12*6063) # flagged as bad data: 29473 Minimum value: 3.4011 Maximum value: 4.5082 Mean value: 3.7253 (unweighted average) Standard deviation: 0.16684 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" CONDUCTIVITY_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" CONDUCTIVITY_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 7583128 Minimum value: 14.753 Maximum value: 267.59 Mean value: 123.52 (unweighted average) Standard deviation: 35.999 yes? STAT ECHO_INTENSITY_BEAM1_DAY Filtered beam1 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 315779 Minimum value: 16.226 Maximum value: 257.8 Mean value: 123.52 (unweighted average) Standard deviation: 36.206 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM1_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" ECHO_INTENSITY_BEAM1_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 7583128 Minimum value: 29.14 Maximum value: 223.66 Mean value: 124.6 (unweighted average) Standard deviation: 34.767 yes? STAT ECHO_INTENSITY_BEAM2_DAY Filtered beam2 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 315779 Minimum value: 30.7 Maximum value: 222.49 Mean value: 124.6 (unweighted average) Standard deviation: 34.934 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM2_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" ECHO_INTENSITY_BEAM2_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 7583128 Minimum value: 14.26 Maximum value: 236.62 Mean value: 124.72 (unweighted average) Standard deviation: 35.692 yes? STAT ECHO_INTENSITY_BEAM3_DAY Filtered beam3 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 315779 Minimum value: 15.052 Maximum value: 233.97 Mean value: 124.72 (unweighted average) Standard deviation: 35.88 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM3_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" ECHO_INTENSITY_BEAM3_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 7583128 Minimum value: 17.714 Maximum value: 223.38 Mean value: 123.17 (unweighted average) Standard deviation: 35.502 yes? STAT ECHO_INTENSITY_BEAM4_DAY Filtered beam4 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 315779 Minimum value: 17.73 Maximum value: 222.05 Mean value: 123.17 (unweighted average) Standard deviation: 35.675 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" ECHO_INTENSITY_BEAM4_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" ECHO_INTENSITY_BEAM4_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (decibars): 5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 79830 Minimum value: 36.669 Maximum value: 36.711 Mean value: 36.694 (unweighted average) Standard deviation: 0.0091309 yes? STAT GPS_LATITUDE_DAY Filtered Latitude LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (decibars): 5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 3322 Minimum value: 36.67 Maximum value: 36.711 Mean value: 36.694 (unweighted average) Standard deviation: 0.0091562 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" GPS_LATITUDE_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" GPS_LATITUDE_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (decibars): 5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 79830 Minimum value: -122.42 Maximum value: -122.38 Mean value: -122.4 (unweighted average) Standard deviation: 0.0083436 yes? STAT GPS_LONGITUDE_DAY Filtered Longitude LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (decibars): 5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 3322 Minimum value: -122.42 Maximum value: -122.38 Mean value: -122.4 (unweighted average) Standard deviation: 0.0083086 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" GPS_LONGITUDE_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" GPS_LONGITUDE_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 1745748 (1*1*12*145479) # flagged as bad data: 262167 Minimum value: -76.936 Maximum value: 1736.4 Mean value: 117.13 (unweighted average) Standard deviation: 102.04 yes? STAT PRESSURE_DAY Filtered Pressure LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 72756 (1*1*12*6063) # flagged as bad data: 10879 Minimum value: -73.273 Maximum value: 1694.5 Mean value: 117.12 (unweighted average) Standard deviation: 102.59 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" PRESSURE_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" PRESSURE_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 64686 Minimum value: 37.343 Maximum value: 102.51 Mean value: 87.984 (unweighted average) Standard deviation: 8.4938 yes? STAT RELATIVE_HUMIDITY_DAY Filtered Relative Humidity LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 2687 Minimum value: 39.843 Maximum value: 101.21 Mean value: 88.012 (unweighted average) Standard deviation: 8.4064 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" RELATIVE_HUMIDITY_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" RELATIVE_HUMIDITY_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 1745748 (1*1*12*145479) # flagged as bad data: 738150 Minimum value: 29.85 Maximum value: 34.266 Mean value: 33.713 (unweighted average) Standard deviation: 0.36507 yes? STAT SALINITY_DAY Filtered Salinity LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 72756 (1*1*12*6063) # flagged as bad data: 30734 Minimum value: 30.255 Maximum value: 34.259 Mean value: 33.713 (unweighted average) Standard deviation: 0.36505 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" SALINITY_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" SALINITY_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 1745748 (1*1*12*145479) # flagged as bad data: 334413 Minimum value: 5.997 Maximum value: 19.028 Mean value: 10.429 (unweighted average) Standard deviation: 2.1861 yes? STAT TEMPERATURE_DAY Filtered Temperature LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 to 375 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 72756 (1*1*12*6063) # flagged as bad data: 13910 Minimum value: 6.0059 Maximum value: 18.989 Mean value: 10.428 (unweighted average) Standard deviation: 2.1868 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" TEMPERATURE_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" TEMPERATURE_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 8788940 Minimum value: -55.913 Maximum value: 53.628 Mean value: -0.27527 (unweighted average) Standard deviation: 9.3088 yes? STAT U_COMPONENT_UNCORR_DAY Filtered East-West Velocity Towards(+E/-W) True LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 366021 Minimum value: -54.425 Maximum value: 50.835 Mean value: -0.2503 (unweighted average) Standard deviation: 9.3608 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" U_COMPONENT_UNCORR_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" U_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 14547900 (1*1*100*145479) # flagged as bad data: 8788940 Minimum value: -56.726 Maximum value: 62.17 Mean value: 0.90493 (unweighted average) Standard deviation: 12.102 yes? STAT V_COMPONENT_UNCORR_DAY Filtered North-South Velocity Towards(+N/-S) True LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): 2.5 to 502.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 606300 (1*1*100*6063) # flagged as bad data: 366021 Minimum value: -55.309 Maximum value: 62.17 Mean value: 0.91232 (unweighted average) Standard deviation: 12.185 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" V_COMPONENT_UNCORR_FIL **ERROR: insufficient memory: 14547900 words were requested. LISTing to file hourlyFilteredM2.nc *** NOTE: You can use SET MEMORY/SIZE=xxx to increase memory. *** NOTE: The "Memory use" section of the FERRET Users Guide has further tips. yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" V_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 36820 Minimum value: -15.812 Maximum value: 12.741 Mean value: 2.1464 (unweighted average) Standard deviation: 2.8587 yes? STAT WIND_U_COMPONENT_DAY Filtered U Component of Wind Velocity LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 1526 Minimum value: -13.715 Maximum value: 12.564 Mean value: 2.1818 (unweighted average) Standard deviation: 2.8195 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" WIND_U_COMPONENT_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" WIND_U_COMPONENT_DAY LISTing to file dailyFilteredM2.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: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 24-SEP-1992 04:30 to 29-APR-2009 19:30 DATA SET: ./hourlyM2.nc Total # of data points: 145479 (1*1*1*145479) # flagged as bad data: 36820 Minimum value: -15.912 Maximum value: 15.39 Mean value: -2.9673 (unweighted average) Standard deviation: 4.4908 yes? STAT WIND_V_COMPONENT_DAY Filtered V Component of Wind Velocity LONGITUDE: 122.4W LATITUDE: 36.7N DEPTH (m): -3.5 TIME: 23-SEP-1992 17:00 to 30-APR-2009 17:00 DATA SET: ./hourlyM2.nc Total # of data points: 6063 (1*1*1*6063) # flagged as bad data: 1526 Minimum value: -15.912 Maximum value: 15.077 Mean value: -2.9934 (unweighted average) Standard deviation: 4.5177 yes? SAVE/APPEND/FILE="hourlyFilteredM2.nc" WIND_V_COMPONENT_FIL LISTing to file hourlyFilteredM2.nc yes? SAVE/APPEND/FILE="dailyFilteredM2.nc" WIND_V_COMPONENT_DAY LISTing to file dailyFilteredM2.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...