7.7.31. quick_enrpy_help.txtΒΆ

This 'help' refers to the adcp_py3demos  ps0918 cruise.
(manual steps to process ENR data, via uhdas-style directory)

=======================================================
(1) find out information about the data

Make a **project** directory (ps0918_vmdas)
From here:

        vmdas_info.py --logfile lta_manual_info.txt  os ../../vmdas_data/ps0918/*LTA
        vmdas_info.py --logfile enr_manual_info.txt  os ../../vmdas_data/ps0918/*ENR



The output says:

      - 3 files, all have data
      - averaging length: 300sec
      - instrument frequency 75
      - pingtypes: bb pings
      - bottom track was on for part of the cruise
      - beam angle 30 (normal for OS75)
      - transducer angle (EA) 1.18
      - messages available and what was used:

      N1R      HEHDT           <-- VmDAS primary heading

      N2R      GPGGA
      N2R      GPHDT
      N2R      PASHR,AT2       <-- VmDAS backup heading
      N2R      PASHR,ATT
      N2R      PASHR,POS

      N3R      GPGGA           <-- VmDAS primary position
      N3R      GPGLL


Look at the info files and note:
    - what kind of pings
    - averaging length
    - are filenames sorted in chronological order?
    - what was the heading source?


(2)  Start a text file "os75bb_ENR_manualproc.txt" adjacent to the processing directory,
     "os75bb_lta" and put this information in it:

       --> these are 'bb' pings
       --> averaged over 300sec
       --> frequency is 75kHz, so "sonar" os "os75bb"
       --> primary heading was HEHDT
       --> there is an Ashtech, used as secondary heading


(3) convert vmdas data to uhdas-style data

    mkdir config
    reform_vmdas.py --project_dir_path ./ --vmdas_dir_path ../../vmdas_data/ps0918 --uhdas_style_dir ../../uhdas_style_data  --cruisename ps0918_manual --start ..
    cd  config
    python3 vmdas2uhdas.py

(4) create the processing configuration file

   proc_starter.py reform_defs.py

   cd ..

(5a) set up processing directory:

     adcptree.py os75bb_manual --datatype uhdas --cruisename ps0918_manual
     cd os75bb_manual
     # make q_py.cnt

(5b) create a quick_adcp.py control file

 create a quick_adcp.py control file

    windows:
        - you can create q_py.cnt with an editor to have the
          parts between (or including) the comments
        - it may be easier to call it "q_py.txt" to avoid
          fighting with your editor, which may try to set the suffix

    unix: this example  is written as a bash "heredoc" for
          cut-and-past ease, but you can create q_py.cnt with
          an editor to have the parts between (or including)
          the comments


cat << EOF > q_py.cnt

    ##### q_py.cnt is

    --yearbase 2009             ## required, for decimal day conversion
                                ##     (year of first data)
    --cruisename ps0918_manual  ## *must* match prefix in config dir
    --dbname aship              ## database name; in adcpdb.  eg. a0918
                                ##
    --datatype uhdas            ## datafile type
    --sonar os75bb              ## specify instrument letters, frequency,
                                ##     (and ping type for ocean surveyors)
    --ens_len  300              ## averages of 300sec duration
                                ##
    --update_gbin               ## required for this kind of processing
    --configtype  python        ## file used in config/ dirctory is python
                                ##
    --ping_headcorr             ## ps0918_manual_proc.py says use HDT first,
                                ##    correct to ashtech
    --max_search_depth 1500     ## try to identify the bottom and eliminate
                                ##    data below the bottom IF topo says
                                ##    the bottom is shallower than 1000m

    #### end of q_py.cnt

EOF


(6) run quick_adcp.py:

   quick_adcp.py --cntfile q_py.cnt --auto



(7) Check calibrations

   **watertrack**
------------
Number of edited points:   2 out of   2
   amp   = 1.0030  + 0.1228 (t - 264.8)
   phase =  -0.01  + 3.2807 (t - 264.8)
            median     mean      std
amplitude   1.0030   1.0030   0.0099
phase      -0.0150  -0.0150   0.2645
------------


   **bottomtrack**
------------
unedited: 21 points
edited:   18 points, 2.0 min speed, 2.5 max dev
            median     mean      std
amplitude   0.9996   1.0004   0.0044
phase       0.3090   0.3392   0.1713

------------


(8)  apply rotation correction (no scale factor necessary)
     (all in one line)


     quick_adcp.py --steps2rerun rotate:navsteps:calib --rotate_angle 0.32  --auto



   **bottomtrack**
------------
unedited: 20 points
edited:   16 points, 2.0 min speed, 2.5 max dev
            median     mean      std
amplitude   1.0002   1.0000   0.0024
phase      -0.0380  -0.0183   0.1113

------------


(9)  edit bad data

   dataviewer.py -e


 # apply editing:


   quick_adcp.py --steps2rerun apply_edit:navsteps:calib  --auto


(10) make web plots

       quick_web.py --interactive

 - view with a browser, look at webpy/index.html

(11) extract data

   quick_adcp.py --steps2rerun matfiles --auto


    adcp_nc.py adcpdb  contour/os75bb  ps0918_demo os75bb --ship_name Point Sur



done with processing.
Add  cruise_info.txt to the top of this file.
Edit with correct info