LAST CHANGED : 2018/12/28 20:01:53 CRUISE NAME(S) : S04P_2018 (nbp1802) CRUISE DATES : 2018/03/10 to 2018/05/10 SHIP NAME : N.B. Palmer PORTS : Hobart,Tasmania - Punta Arenas, Chile CHIEF SCIENTIST : Alison Macdonald DATABASE NAME : aship DATA FILES : np2018_068_22866.raw to np2018_129_64800.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 : os75 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/03/10 06:26:07 off 300 60 16 7 8 45.00 0.00 0001 HEADING : PRIMARY : heading from gyro CORRECTION : heading correction from seapath1 NOTE: time-dependent heading corrections applied IN the ensembles (see cal/rotate/ens_hcorr.ang) POSITIONS : gps positions from gpsnav CALIBRATION : (check original processing parameters) additional rotation 0 final transducer angle is: 45.0 (original transducer angle) - (rotate_angle) 45.0 - 0 original scale factor 1 additional scale factor 1.002 GPS x sonar offset; adjustment applied: .single ping processing: (q_py.cnt) - original: dx = 2, dy = -3, - applied: --xducer_dx 5, # ('dy' was the same as original --xducer_dy -3) .post-processing: --xducer_dy 2 .final: --xducer_dx 5 --xducer_dx -1 COMMENTS : just a few gaps in heading correction some scattering layers PROCESSOR : U. Pinho --- processing parameters ---------- ## (determined from "sonar"): model = os ## (determined from "sonar"): frequency = 75 ## (determined from "sonar"): instname = os75 ## (determined from "sonar"): pingtype = nb badbeam None beamangle 30 configtype python cruisename nbp1802 datatype uhdas dbname aship ens_len 300 fixfile aship.gps hcorr_inst seapath1 pgmin 50 pingpref None proc_engine python ref_method refsm refuv_smoothwin 3 refuv_source nav sonar os75nb txy_file aship.agt xducer_dx 5 xducer_dy -3 yearbase 2018 (1) Get some info: uhdas_info.py --overview --logfile uhdas_info.txt /home/data/palmer/nbp1802 ================================================================== ======= start of report for /home/data/palmer/nbp1802 ===== ================================================================== ------raw------- raw: gpsnav 739 files (np2018_068_22866 - np2018_129_64800) raw: gyro 739 files (np2018_068_22866 - np2018_129_64800) raw: seapath1 739 files (np2018_068_22866 - np2018_129_64800) raw: seapath2 739 files (np2018_068_22866 - np2018_129_64800) raw: sndspd 1 files (np2018_068_22866 - np2018_068_22866) adcp: os38 .raw 739 files (np2018_068_22866 - np2018_129_64800) adcp: os38 .raw.log 739 files (np2018_068_22866 - np2018_129_64800) adcp: os38 .raw.log.bin 739 files (np2018_068_22866 - np2018_129_64800) adcp: os75 .raw 739 files (np2018_068_22866 - np2018_129_64800) adcp: os75 .raw.log 739 files (np2018_068_22866 - np2018_129_64800) adcp: os75 .raw.log.bin 739 files (np2018_068_22866 - np2018_129_64800) -----rbin------- rbin: gpsnav:gps 739 files (np2018_068_22866 - np2018_129_64800) rbin: gyro:hdg 739 files (np2018_068_22866 - np2018_129_64800) rbin: seapath1:gps_sea 739 files (np2018_068_22866 - np2018_129_64800) rbin: seapath1:sea 739 files (np2018_068_22866 - np2018_129_64800) rbin: seapath2:gps_sea 739 files (np2018_068_22866 - np2018_129_64800) rbin: seapath2:sea 739 files (np2018_068_22866 - np2018_129_64800) -----gbin------- gbin: os38 (gyro,seapath2,gpsnav,time,sndspd,seapath1) ------ database time ranges -------- proc: os38bb No database found proc: os38nb 68.268 - 129.796 (2018/03/10 to 2018/05/10) proc: os75bb No database found proc: os75nb 68.268 - 129.796 (2018/03/10 to 2018/05/10) ---------- end of report for /home/data/palmer/nbp1802-------------- # What we see: - GPS messages from: gpsnav,seapath1 and seapath2 - Heading messages : seapath1, seapath2 and gyro - All other feeds : 739 files - proc: os38nb 68.268 - 129.796 (2018/03/10 to 2018/05/10) - proc: os75nb 68.268 - 129.796 (2018/03/10 to 2018/05/10) cruise 61.52 days long. # ****** # Single ping adcp_database_maker.py was applied # sonar - gps offset applied: # original dx=2, dy=-3, applied --xducer_dx 5, # (dy remain the same --xducer_dy -3) ### POST-PROCESSING ----------------------------------------------------------- cd os75nb plot_nav.py nav/aship.gps Go-Ship cruise from Hobart, Tasmania, to Chile, along 67 degrees South and Antartica. (2) review the data (a) check calibration:heading correction device: figview.py cal/rotate/*png - Gaps in the heading correction patch_hcorr_app.py applied. (b) check calibration: **watertrack** ------------ Number of edited points: 212 out of 239 amp = 1.0026 + 0.0001 (t - 101.3) phase = 0.04 + -0.0023 (t - 101.3) median mean std amplitude 1.0020 1.0026 0.0083 phase 0.0565 0.0402 0.4222 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2018/12/28 22:35:45 xducer_dx = 0.787936 xducer_dy = 1.348335 signal = 6312.485595 ------------ # uvship was also applied: quick_adcp.py --steps2rerun navsteps:calib --refuv_source uvship --auto # reference layer: 30 to 50: quick_adcp.py --steps2rerun navsteps:calib --rl_startbin 30 --rl_endbin 50 --auto **watertrack** ------------ Number of edited points: 35 out of 60 amp = 1.0030 + -0.0000 (t - 110.8) phase = 0.15 + -0.0018 (t - 110.8) median mean std amplitude 1.0020 1.0030 0.0056 phase 0.1840 0.1455 0.3941 ------------ (3) edit: dataviewer.py -e quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto check calibration: **watertrack** ------------ Number of edited points: 199 out of 212 amp = 1.0012 + 0.0001 (t - 101.4) phase = 0.05 + -0.0024 (t - 101.4) median mean std amplitude 1.0020 1.0012 0.0075 phase 0.0400 0.0520 0.4020 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2018/12/29 00:52:29 xducer_dx = -0.275101 xducer_dy = 2.477279 signal = 6312.485595 ------------ quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_amplitude 1.002 --xducer_dy 2 --auto **watertrack** ------------ Number of edited points: 200 out of 212 amp = 0.9992 + 0.0001 (t - 101.5) phase = 0.05 + -0.0024 (t - 101.5) median mean std amplitude 1.0000 0.9992 0.0074 phase 0.0450 0.0509 0.4010 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2018/12/29 01:01:19 xducer_dx = -0.290542 xducer_dy = 1.282775 signal = 6312.485595 ------------ # check editing -- looks OK dataviewer.py (4) small editing again: dataviewer.py -e quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto check calibration: **watertrack** ------------ Number of edited points: 199 out of 210 amp = 0.9992 + 0.0001 (t - 101.4) phase = 0.05 + -0.0024 (t - 101.4) median mean std amplitude 1.0000 0.9992 0.0074 phase 0.0450 0.0479 0.3976 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2018/12/29 01:15:11 xducer_dx = -0.280875 xducer_dy = 1.274689 signal = 6312.485595 ------------ # patch_hcorr_app.py applied again smoothing seapath1 data: ---Original Calibration--- - Water Track Calibration: median mean std amplitude 1.0000 0.9992 0.0074 phase 0.0450 0.0479 0.3976 ---New Calibration--- - Water Track Calibration: median mean std amplitude 1.0020 1.0012 0.0075 phase 0.0755 0.0612 0.4006 (4) small editing again: dataviewer.py -e quick_adcp.py --steps2rerun apply_edit:navsteps:calib --auto **watertrack** ------------ Number of edited points: 197 out of 210 amp = 1.0013 + 0.0001 (t - 101.3) phase = 0.06 + -0.0023 (t - 101.3) median mean std amplitude 1.0020 1.0013 0.0070 phase 0.0780 0.0649 0.4044 ------------ **transducer-gps offset** ------------ guessing ADCP (dx=starboard, dy=fwd) meters from GPS positions from aship.agt calculation done at 2019/01/02 23:11:35 xducer_dx = -0.295208 xducer_dy = 1.842284 signal = 6312.485595 ------------ (5) make web plots quick_web.py --interactive - view with a browser, look at webpy/index.html (6) Extract data: quick_adcp.py --steps2rerun matfiles --auto You can extract NetCDF files as well. adcp_nc.py adcpdb contour/os75nb_short nbp1802 os75nb adcp_nc.py --long adcpdb contour/os75nb_long nbp1802 os75nb To check netcdf file: ncdump -h contour/os75nb_short.nc # cat cruise_info.txt at the top of this file!