LAST CHANGED : 2018/06/15 16:23:51 CRUISE NAME(S) : FK180528 CRUISE DATES : 2018/05/29 to 2018/06/15 SHIP NAME : Falkor PORTS : San Diego, CA to San Diego, CA CHIEF SCIENTIST : Dr. Joao Borges de Sousa DATABASE NAME : a_fk DATA FILES : fk2018_148_75537.raw to fk2018_165_57600.raw STATUS : to do done ------ ----------- averaged [ x ] loaded [ x ] NOTE: heading correction instrument exists NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) check heading correction [ x ] calibration [ x ] edited [ x ] re-check heading correction [ x ] check editing [ x ] figures [ x ] INSTRUMENT : wh300 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 2018/05/29 21:00:58 off 120 70 2 7 2 44.63 0.00 0001 2018/05/31 16:44:36 off 120 70 2 7 2 44.63 0.00 0001 2018/06/02 18:44:36 off 120 70 2 7 2 44.63 0.00 0001 2018/06/03 02:36:42 off 120 70 2 7 2 44.63 0.00 0001 2018/06/14 18:36:42 off 120 70 2 7 2 44.63 0.00 0001 HEADING : PRIMARY : heading from gyro CORRECTION : heading correction from seapath NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) POSITIONS : gps positions from cnav CALIBRATION : (check original processing parameters) 44.63 additional rotation 0.05 final transducer angle is: (original transducer angle) - (rotate_angle) 44.63 - (0.05 ) = 44.58 original scale factor 1 additional scale factor 1.005 COMMENTS : no gaps in heading correction PROCESSOR : U. Pinho --- processing parameters ---------- ## (determined from "sonar"): model = wh ## (determined from "sonar"): frequency = 300 ## (determined from "sonar"): instname = wh300 ## (determined from "sonar"): pingtype = bb beamangle 20 cruisename FK180528 datatype uhdas dbname a_fk ens_len 120 fixfile a_fk.gps frequency 300 hcorr_inst seapath instname wh300 model wh pingtype bb proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar wh300 txy_file a_fk.agt xducer_dx 2 xducer_dy 18 yearbase 2018 (1) Get some info about the data set: uhdas_info.py --overview --logfile uhdas_info.txt /home/data/falkor/FK180528/ADCP/ ================================================================== ======= start of report for /home/data/falkor/FK180528/ADCP/ ===== ================================================================== ------raw------- raw: cnav 201 files (fk2018_148_45010 - fk2018_165_59105) raw: gyro 201 files (fk2018_148_45010 - fk2018_165_59105) raw: posmv 201 files (fk2018_148_45010 - fk2018_165_59105) raw: seapath 201 files (fk2018_148_45010 - fk2018_165_59105) raw: soundspeed 201 files (fk2018_148_45010 - fk2018_165_59105) raw: tsg 201 files (fk2018_148_45010 - fk2018_165_59105) adcp: os75 .raw 196 files (fk2018_148_45010 - fk2018_165_57600) adcp: os75 .raw.log 196 files (fk2018_148_45010 - fk2018_165_57600) adcp: os75 .raw.log.bin 196 files (fk2018_148_45010 - fk2018_165_57600) adcp: wh300 .raw 199 files (fk2018_148_75537 - fk2018_165_57600) adcp: wh300 .raw.log 199 files (fk2018_148_75537 - fk2018_165_57600) adcp: wh300 .raw.log.bin 199 files (fk2018_148_75537 - fk2018_165_57600) -----rbin------- rbin: cnav:gps 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: gyro:hdg 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: posmv:gps 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: posmv:pmv 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: seapath:gps_sea 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: seapath:sea 201 files (fk2018_148_45010 - fk2018_165_59105) rbin: soundspeed:spd 201 files (fk2018_148_45010 - fk2018_165_59105) -----gbin------- gbin: os75 (gyro,posmv,cnav,time,soundspeed,seapath) gbin: wh300 (gyro,posmv,cnav,time,soundspeed,seapath) ------ database time ranges -------- proc: os75bb 148.524 - 148.528 (2018/05/29 to 2018/05/29) proc: os75nb 148.524 - 165.680 (2018/05/29 to 2018/06/15) proc: wh300 148.876 - 165.682 (2018/05/29 to 2018/06/15) ---------- end of report for /home/data/falkor/FK180528/ADCP/-------------- (1a) Sonars: os75nb and wh300 We see: - positions cnav and posmv - headings from gyro heading correction from seapath - posmv also available. - 201 files for all devices. (1b) From cruise_info.txt: HEADING : PRIMARY : heading from gyro CORRECTION : heading correction from posmv NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) POSITIONS : gps positions from gpsnav (2) Cruise summary: cd FK180528/wh300 plot_nav.py nav/a_fk.gps Ocean Pacific 30N between 130-135W, south of Galapagos. #----------------- # Post-processing: (3) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png - no need to apply patch_hcorr_app.py. (b) check calibration: **watertrack** ------------ Number of edited points: 129 out of 141 amp = 1.0141 + 0.0003 (t - 156.4) phase = -0.01 + 0.0079 (t - 156.4) median mean std amplitude 1.0120 1.0141 0.0124 phase -0.0670 -0.0102 0.7521 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_fk.agt calculation done at 2018/06/15 16:23:50 xducer_dx = -1.216145 xducer_dy = 0.455353 signal = 2395.380567 ------------ quick_adcp.py --steps2rerun rotate:apply_edit:navsteps:calib --rotate_amplitude 1.01 --rotate_angle -0.05 --auto **watertrack** ------------ Number of edited points: 139 out of 143 amp = 0.9989 + 0.0004 (t - 156.6) phase = 0.11 + 0.0043 (t - 156.6) median mean std amplitude 0.9980 0.9989 0.0075 phase 0.0580 0.1057 0.3387 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_fk.agt calculation done at 2018/09/17 23:48:45 xducer_dx = -0.919556 xducer_dy = 0.455556 signal = 2395.380570 ------------ (c) look at the data: dataviewer.py (5) edit "gautoedit" cd edit gautoedit.py -n5 # apply editing: cd .. quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto # check editing -- looks OK dataviewer.py **watertrack** ------------ Number of edited points: 139 out of 143 amp = 0.9989 + 0.0004 (t - 156.6) phase = 0.11 + 0.0045 (t - 156.6) median mean std amplitude 0.9990 0.9989 0.0075 phase 0.0580 0.1064 0.3371 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from a_fk.agt calculation done at 2018/09/18 00:36:40 xducer_dx = -0.917553 xducer_dy = 0.453218 signal = 2395.380570 ------------ ======================================================= After comparing with OS75, it appears an additional 0.995 (amp) and 0.1 (phase) are needed on the wh300 to make it agree with the os75nb. Redoing the watertrack calibraiton with a deeper reference layer is more consistent with the OS75 calibration reference layer and also yields the right kind of additional correction to make for the two to agree. - using a deeper reference layer: quick_adcp.py --steps2rerun navsteps:calib --rl_startbin 20 --rl_endbin 40 --auto ** watertrack ** Number of edited points: 143 out of 146 amp = 0.9958 + 0.0002 (t - 156.6) phase = 0.12 + 0.0090 (t - 156.6) median mean std amplitude 0.9950 0.9958 0.0063 phase 0.1040 0.1170 0.3561 - apply it: quick_adcp.py --steps2rerun navsteps:calib:rotate --rotate_amplitude .995 --rotate_angle 0.1 --auto **watertrack** ------------ Number of edited points: 143 out of 146 amp = 1.0008 + 0.0002 (t - 156.6) phase = 0.02 + 0.0090 (t - 156.6) median mean std amplitude 1.0000 1.0008 0.0064 phase 0.0090 0.0167 0.3557 ------------ Then dataviewer comparison looks much better. ======================================================= figview.py cal/watertrk/wtcal1.png (7) make web plots quick_web.py --interactive - view with a browser, look at webpy/index.html (8) Extract data: quick_adcp.py --steps2rerun matfiles --auto Extract NetCDF files as well. adcp_nc.py adcpdb contour/wh300_short FK180528 wh300 adcp_nc.py --long adcpdb contour/wh300_long FK180528 wh300 To check netcdf file: ncdump -h contour/wh300_short.nc