LAST CHANGED : 2022/11/15 23:49:10 CRUISE NAME(S) : P02E (rr2205) CRUISE DATES : 2022/06/13 to 2022/07/16 SHIP NAME : Roger Revelle PORTS : Honolulu, HI to San Diego, CA CHIEF SCIENTIST : Andreas Thurnherr / Sebastien Bigorre (Co-Chief) DATABASE NAME : a_rr DATA FILES : rr2022_164_19782.raw to rr2022_196_50400.raw STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] NOTE: No automated time-dependent heading correction exists check heading correction [ x ] calibration [ x ] edited [ x ] re-check heading correction [ x ] check editing [ x ] figures [ x ] INSTRUMENT : os150 ACQUISITION : PROGRAM : uhdas PROCESSING : python LOGGING : PARAMETERS : BT : bottom track mode (on or off) SI : sampling interval or averaging period for ensemble (sec) NB : number of bins BL : bin length (m) TD : transducer depth (m) BK : blanking length (m) HO : heading offset applied by DAS (deg) HB : heading bias (deg) CRPH : compensation for roll-pitch-heading, 1:on, 0:off) yy/mm/dd hh:mm:ss BT SI NB BL TD BK HO HB CRPH 2022/06/20 18:51:28 off 180 60 8 6 8 41.78 0.00 0001 HEADING : PRIMARY : heading from gyro CORRECTION : does not exist? POSITIONS : gps positions from seapath CALIBRATION : (check original processing parameters) additional rotation 0 final transducer angle is: (41.78) - (-0.28) = 42.06 applied scale factor 1 additional scale factor 1.005 PROCESSOR : Antonio Goulart --- final processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 150 ## (determined from "sonar"): instname = os150 ## (determined from "sonar"): pingtype = nb badbeam None beam_order None beamangle 30 configtype python cruisename RR2205 datatype uhdas dbname a_rr ens_len 180 fixfile a_rr.gps pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os150nb txy_file a_rr.agt xducer_dx -1 xducer_dy 11 yearbase 2022 ================================================== PROCESSING STEPS: ================================================== ----------------------------------------------------------------------------- # 1. Check heading correction ----------------------------------------------------------------------------- # Processing commands here forward start in the os150nb directory. Check cruise track $ plot_nav.py nav/a_rr.gps Check that we have a good heading correction for all of the ADCP data: $ figview.py cal/rotate/ens* There should be no poor fixes in the heading correction graphs. These are represented by red crosses, good fixes are represented by green dots. Gaps in the heading correction graphs are okay, no data was collected during these times. If there are red crosses on the heading correction graphs then they need to be patched with patch_hcorr.py. OBS: I had to go check the files on moli, on the cal/rotate directory for each of the two parts of the cruise. No red crosses spotted. ----------------------------------------------------------------------------- # 2. ADCP calibration, Part 1 ----------------------------------------------------------------------------- Calibrate the dataset in sum by checking the calibration values and applying a fix if necessary. When looking at a water track or bottom track calibration, we want to see statistics inside the following values: median_tolerance amplitude 0.997 to 1.003 phase/angle -0.05 to 0.05 dx,dy round to -2 to 2 The xy calibration is an estimate of the horizontal separation between the physical location of the ADCP and primary GPS device in meters. When looking at an xy calibration estimate we want the dx and dy values to be as close to zero as possible, as long as the signal is between 1000 and 5000. Less than 1000 means there is very little data, more than 5000 means there is too much change. If values are out of range for either calibration we make a bulk correction to the entire dataset. This should reduce the number of outliers to be edited out. If we make any changes to x or y they should be integers. ------------ $ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_rr.agt calculation done at 2022/11/15 23:49:09 xducer_dx = 0.499127 xducer_dy = -1.785002 signal = 1320.600288 ------------ $ catwt **watertrack** ------------ Number of edited points: 160 out of 161 amp = 1.0053 + 0.0000 (t - 183.1) phase = -0.27 + -0.0010 (t - 183.1) median mean std amplitude 1.0050 1.0053 0.0030 phase -0.2870 -0.2708 0.2031 ------------ dx, dy and phase are within good limits. Amplitude and phase need adjustments. amplitude: 1.005 phase: -0.28 We will first check the calibration figures to make sure that the watertrack estimates make sense: $ figview.py cal/watertrk/*png The calibration estimates are reasonable. We can go ahead and apply the calibration: $ quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.005 --rotate_angle -0.28 --auto Checking again: $ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_rr.agt calculation done at 2022/12/07 21:46:08 xducer_dx = 0.467742 xducer_dy = -1.772205 signal = 1320.600288 ------------ $ catwt **watertrack** ------------ Number of edited points: 160 out of 161 amp = 1.0004 + 0.0000 (t - 183.1) phase = 0.01 + -0.0009 (t - 183.1) median mean std amplitude 1.0000 1.0004 0.0033 phase -0.0070 0.0094 0.2042 ------------ Everything looks good for now, time to move to the next step. ----------------------------------------------------------------------------- # 3. Flagging Questionable Data ----------------------------------------------------------------------------- Flagging of biased/questionable data or artifacts in the dataset. We used dataviewer.py to look for problems with the data and flag the data as questionable. $ dataviewer.py -e There are no other sonars for comparison. Now we check the calibrations again to make sure that there were no large shifts since we made edits. $ catxy **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_rr.agt calculation done at 2022/12/07 21:46:08 xducer_dx = 0.467742 xducer_dy = -1.772205 signal = 1320.600288 ------------ $ catwt **watertrack** ------------ Number of edited points: 160 out of 161 amp = 1.0004 + 0.0000 (t - 183.1) phase = 0.01 + -0.0009 (t - 183.1) median mean std amplitude 1.0000 1.0004 0.0033 phase -0.0070 0.0094 0.2042 ------------ All calibration values are within reasonable limits. ----------------------------------------------------------------------------- # 4. Re-check all figures ----------------------------------------------------------------------------- Check all figures again to make sure that any problems were addressed and no new problems have appeared after making changes to the dataset. $ figview.py All of the figures look fine. ----------------------------------------------------------------------------- # 5. Check edited, calibrated dataset against original dataset ----------------------------------------------------------------------------- Compare the edited, calibrated dataset against the original dataset to make sure all problems have been dealt with and no new problems have appeared. $ dataviewer.py -c . ../os150nb.orig ----------------------------------------------------------------------------- # 6. Make plots and files ----------------------------------------------------------------------------- Create the figures and data files needed to finish processing and submit. Make the plots needed for web viewing, matlab files (legacy), and netCDF files, then check that the netCDF files are readable. Plots should be in 3 to 5 day chunks for a cruise longer than one week, or divided by geographic features where it makes sense. $ quick_web.py --interactive $ figview.py webpy $ quick_adcp.py --steps2rerun matfiles --auto $ adcp_nc.py adcpdb contour/os150nb RR2205 os150nb --ship_name Roger Revelle $ ncdump contour/os150nb.nc -h