The data in this data set are available via DODS/OPeNDAP
at
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM4.nc
and through the DODS html form mechanism:
http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/dailyFilteredM4.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/hourlyFilteredM4.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/hourlyM4.nc).
For information on the individual deployments and instruments that contributed to this data set please see the M4 hourly data set documentation.
Summary plots of all the variables in this data set
! Description: Create daily filtered OASIS data file from hourly gridded data ! Automatically generated by ./hourly2dailyfilter.pl on Mon Jun 13 07:54:13 2005. ! ! All data sets have AIR_TEMPERATURE_HR variable; ! use it to get axis time limits. ! Filter to remove tidal and inertial signals. ! Regrid from 1 hour to 1 day to reduce # file size and number of points ! Also save an intermediate hourly filtered file for finer scale use. ! SET MEMORY/SIZE=30 SET MODE DESPERATE ! ! Leslie / Fred's plweight.m for 33 hour filter ! DEFINE AXIS/X=-33:33:1 XAX1 DEFINE GRID/X=XAX1 G1 FILE/VAR=WEIGHTS33/GRID=G1 "./33hrfilt" SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33 LIST WEIGHTS33 ! ! Original hourly data file ! !USE "http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.nc" USE "hourlyM4.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 hourlyFilteredM4.nc" SPAWN "/bin/rm -f dailyFilteredM4.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="hourlyFilteredM4.nc" AIR_PRESS_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" AIR_TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" CONDUCTIVITY_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM1_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM2_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM3_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM4_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" GPS_LATITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" GPS_LONGITUDE_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" PRESSURE_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" RELATIVE_HUMIDITY_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" SALINITY_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" TEMPERATURE_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" U_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" V_COMPONENT_UNCORR_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" WIND_U_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM4.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="hourlyFilteredM4.nc" WIND_V_COMPONENT_FIL SAVE/APPEND/FILE="dailyFilteredM4.nc" WIND_V_COMPONENT_DAY CANCEL VAR WIND_V_COMPONENT_DAY CANCEL VAR WIND_V_COMPONENT_FIL QUIT
NOAA/PMEL TMAP FERRET v5.80 Solaris 5.8 - 01/03/05 13-Jun-05 07: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 Mon Jun 13 07:54:13 2005. yes? ! yes? ! All data sets have AIR_TEMPERATURE_HR variable; yes? ! use it to get axis time limits. yes? ! Filter to remove tidal and inertial signals. yes? ! Regrid from 1 hour to 1 day to reduce # file size and number of points yes? ! Also save an intermediate hourly filtered file for finer scale use. yes? ! yes? SET MEMORY/SIZE=30 Cached data cleared from memory cancel mode journal sp rm -f ferret.jnl yes? SET MODE DESPERATE yes? yes? ! yes? ! Leslie / Fred's plweight.m for 33 hour filter yes? ! yes? DEFINE AXIS/X=-33:33:1 XAX1 yes? DEFINE GRID/X=XAX1 G1 yes? FILE/VAR=WEIGHTS33/GRID=G1 "./33hrfilt" yes? SET VAR/TITLE="33 Hour Filter Weights" WEIGHTS33 yes? LIST WEIGHTS33 VARIABLE : 33 Hour Filter Weights FILENAME : 33hrfilt SUBSET : 67 points (X) -33 / 1: -0.00026 -32 / 2: -0.00112 -31 / 3: -0.00210 -30 / 4: -0.00315 -29 / 5: -0.00426 -28 / 6: -0.00536 -27 / 7: -0.00641 -26 / 8: -0.00735 -25 / 9: -0.00811 -24 / 10: -0.00864 -23 / 11: -0.00886 -22 / 12: -0.00872 -21 / 13: -0.00816 -20 / 14: -0.00714 -19 / 15: -0.00561 -18 / 16: -0.00356 -17 / 17: -0.00098 -16 / 18: 0.00213 -15 / 19: 0.00573 -14 / 20: 0.00979 -13 / 21: 0.01424 -12 / 22: 0.01901 -11 / 23: 0.02400 -10 / 24: 0.02910 -9 / 25: 0.03422 -8 / 26: 0.03922 -7 / 27: 0.04399 -6 / 28: 0.04841 -5 / 29: 0.05237 -4 / 30: 0.05576 -3 / 31: 0.05850 -2 / 32: 0.06051 -1 / 33: 0.06174 0 / 34: 0.06215 1 / 35: 0.06174 2 / 36: 0.06051 3 / 37: 0.05850 4 / 38: 0.05576 5 / 39: 0.05237 6 / 40: 0.04841 7 / 41: 0.04399 8 / 42: 0.03922 9 / 43: 0.03422 10 / 44: 0.02910 11 / 45: 0.02400 12 / 46: 0.01901 13 / 47: 0.01424 14 / 48: 0.00979 15 / 49: 0.00573 16 / 50: 0.00213 17 / 51: -0.00098 18 / 52: -0.00356 19 / 53: -0.00561 20 / 54: -0.00714 21 / 55: -0.00816 22 / 56: -0.00872 23 / 57: -0.00886 24 / 58: -0.00864 25 / 59: -0.00811 26 / 60: -0.00735 27 / 61: -0.00641 28 / 62: -0.00536 29 / 63: -0.00426 30 / 64: -0.00315 31 / 65: -0.00210 32 / 66: -0.00112 33 / 67: -0.00026 yes? yes? ! yes? ! Original hourly data file yes? ! yes? !USE "http://dods.mbari.org/cgi-bin/nph-nc/data/OASISdata/netcdf/hourlyM4.n c" yes? USE "hourlyM4.nc" yes? DEFINE SYMBOL TBEG="`AIR_TEMPERATURE_HR,return=tstart`" !- DEFINE SYMBOL TBEG="24-AUG-1999 00:00" yes? DEFINE SYMBOL TEND="`AIR_TEMPERATURE_HR,return=tend`" !- DEFINE SYMBOL TEND="08-SEP-2000 00:00" yes? DEFINE AXIS/T=($TBEG):($TEND):1/unit=day TIME_DAY !- DEFINE AXIS/T="24-AUG-1999 00:00":"08-SEP-2000 00:00":1/unit=day TIME_DAY yes? SPAWN "/bin/rm -f hourlyFilteredM4.nc" yes? SPAWN "/bin/rm -f dailyFilteredM4.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.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 960 Minimum value: 999.69 Maximum value: 1030.8 Mean value: 1017.8 (unweighted average) Standard deviation: 4.4489 yes? STAT AIR_PRESS_DAY Filtered Air Pressure LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 38 Minimum value: 999.79 Maximum value: 1030.7 Mean value: 1017.9 (unweighted average) Standard deviation: 4.4746 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" AIR_PRESS_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" AIR_PRESS_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 1094 Minimum value: 9.1288 Maximum value: 16.563 Mean value: 12.649 (unweighted average) Standard deviation: 1.2616 yes? STAT AIR_TEMPERATURE_DAY Filtered air_temperature LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 42 Minimum value: 9.1288 Maximum value: 16.472 Mean value: 12.661 (unweighted average) Standard deviation: 1.2632 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" AIR_TEMPERATURE_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" AIR_TEMPERATURE_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 109740 (1*1*12*9145) # flagged as bad data: 2870 Minimum value: 3.3512 Maximum value: 4.2816 Mean value: 3.6911 (unweighted average) Standard deviation: 0.1906 yes? STAT CONDUCTIVITY_DAY Filtered Conductivity LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 4584 (1*1*12*382) # flagged as bad data: 132 Minimum value: 3.3534 Maximum value: 4.2781 Mean value: 3.6912 (unweighted average) Standard deviation: 0.19076 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" CONDUCTIVITY_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" CONDUCTIVITY_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 768772 Minimum value: 40.438 Maximum value: 165.69 Mean value: 79.063 (unweighted average) Standard deviation: 33.797 yes? STAT ECHO_INTENSITY_BEAM1_DAY Filtered beam1 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 32128 Minimum value: 40.485 Maximum value: 162.42 Mean value: 78.971 (unweighted average) Standard deviation: 33.776 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM1_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" ECHO_INTENSITY_BEAM1_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 768772 Minimum value: 40.019 Maximum value: 176.56 Mean value: 80.869 (unweighted average) Standard deviation: 36.291 yes? STAT ECHO_INTENSITY_BEAM2_DAY Filtered beam2 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 32128 Minimum value: 40.043 Maximum value: 172.97 Mean value: 80.771 (unweighted average) Standard deviation: 36.265 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM2_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" ECHO_INTENSITY_BEAM2_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 768772 Minimum value: 42.979 Maximum value: 175.3 Mean value: 83.804 (unweighted average) Standard deviation: 35.623 yes? STAT ECHO_INTENSITY_BEAM3_DAY Filtered beam3 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 32128 Minimum value: 43.049 Maximum value: 171.76 Mean value: 83.707 (unweighted average) Standard deviation: 35.601 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM3_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" ECHO_INTENSITY_BEAM3_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 768772 Minimum value: 41.891 Maximum value: 172.32 Mean value: 81.853 (unweighted average) Standard deviation: 35.471 yes? STAT ECHO_INTENSITY_BEAM4_DAY Filtered beam4 Echo Intensity data per beam LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 32128 Minimum value: 41.906 Maximum value: 170.71 Mean value: 81.758 (unweighted average) Standard deviation: 35.448 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" ECHO_INTENSITY_BEAM4_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" ECHO_INTENSITY_BEAM4_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (decibars): 5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 9145 yes? STAT GPS_LATITUDE_DAY Filtered Latitude LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (decibars): 5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 382 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" GPS_LATITUDE_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" GPS_LATITUDE_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (decibars): 5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 9145 yes? STAT GPS_LONGITUDE_DAY Filtered Longitude LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (decibars): 5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 382 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" GPS_LONGITUDE_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" GPS_LONGITUDE_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 109740 (1*1*12*9145) # flagged as bad data: 1187 Minimum value: 0.55218 Maximum value: 348.35 Mean value: 130.23 (unweighted average) Standard deviation: 111.57 yes? STAT PRESSURE_DAY Filtered Pressure LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 4584 (1*1*12*382) # flagged as bad data: 63 Minimum value: 0.55233 Maximum value: 348.14 Mean value: 130.22 (unweighted average) Standard deviation: 111.58 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" PRESSURE_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" PRESSURE_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 1027 Minimum value: 64.094 Maximum value: 102.4 Mean value: 89.127 (unweighted average) Standard deviation: 7.2298 yes? STAT RELATIVE_HUMIDITY_DAY Filtered Relative Humidity LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 41 Minimum value: 65.92 Maximum value: 101.34 Mean value: 89.056 (unweighted average) Standard deviation: 7.2476 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" RELATIVE_HUMIDITY_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" RELATIVE_HUMIDITY_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 109740 (1*1*12*9145) # flagged as bad data: 4820 Minimum value: 32.963 Maximum value: 34.224 Mean value: 33.77 (unweighted average) Standard deviation: 0.31354 yes? STAT SALINITY_DAY Filtered Salinity LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 4584 (1*1*12*382) # flagged as bad data: 220 Minimum value: 32.997 Maximum value: 34.223 Mean value: 33.771 (unweighted average) Standard deviation: 0.3135 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" SALINITY_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" SALINITY_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 109740 (1*1*12*9145) # flagged as bad data: 2889 Minimum value: 5.8758 Maximum value: 17.074 Mean value: 9.9876 (unweighted average) Standard deviation: 2.4175 yes? STAT TEMPERATURE_DAY Filtered Temperature LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 to 375 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 4584 (1*1*12*382) # flagged as bad data: 132 Minimum value: 5.9065 Maximum value: 17.049 Mean value: 9.9892 (unweighted average) Standard deviation: 2.4187 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" TEMPERATURE_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" TEMPERATURE_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 820593 Minimum value: -38.496 Maximum value: 36.515 Mean value: -2.9808 (unweighted average) Standard deviation: 11.005 yes? STAT U_COMPONENT_UNCORR_DAY Filtered East-West Velocity Towards(+E/-W) True LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 34296 Minimum value: -37.28 Maximum value: 36.498 Mean value: -2.9397 (unweighted average) Standard deviation: 11.016 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" U_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" U_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 914500 (1*1*100*9145) # flagged as bad data: 820593 Minimum value: -33.932 Maximum value: 25.247 Mean value: -1.1826 (unweighted average) Standard deviation: 9.059 yes? STAT V_COMPONENT_UNCORR_DAY Filtered North-South Velocity Towards(+N/-S) True LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): 2.5 to 502.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 38200 (1*1*100*382) # flagged as bad data: 34296 Minimum value: -33.932 Maximum value: 24.764 Mean value: -1.1996 (unweighted average) Standard deviation: 9.0871 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" V_COMPONENT_UNCORR_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" V_COMPONENT_UNCORR_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 97 Minimum value: -8.8039 Maximum value: 13.481 Mean value: 3.6278 (unweighted average) Standard deviation: 3.6333 yes? STAT WIND_U_COMPONENT_DAY Filtered U Component of Wind Velocity LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 5 Minimum value: -7.4487 Maximum value: 12.948 Mean value: 3.6441 (unweighted average) Standard deviation: 3.6592 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" WIND_U_COMPONENT_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" WIND_U_COMPONENT_DAY LISTing to file dailyFilteredM4.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.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 23:30 to 08-SEP-2000 00:30 DATA SET: ./hourlyM4.nc Total # of data points: 9145 (1*1*1*9145) # flagged as bad data: 97 Minimum value: -10.942 Maximum value: 10.745 Mean value: -2.3403 (unweighted average) Standard deviation: 4.3523 yes? STAT WIND_V_COMPONENT_DAY Filtered V Component of Wind Velocity LONGITUDE: 122.4W LATITUDE: 36.2N DEPTH (m): -3.5 TIME: 23-AUG-1999 12:00 to 08-SEP-2000 12:00 DATA SET: ./hourlyM4.nc Total # of data points: 382 (1*1*1*382) # flagged as bad data: 5 Minimum value: -10.942 Maximum value: 10.116 Mean value: -2.332 (unweighted average) Standard deviation: 4.3457 yes? SAVE/APPEND/FILE="hourlyFilteredM4.nc" WIND_V_COMPONENT_FIL LISTing to file hourlyFilteredM4.nc yes? SAVE/APPEND/FILE="dailyFilteredM4.nc" WIND_V_COMPONENT_DAY LISTing to file dailyFilteredM4.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...