====== IDL Installing 6.3 ======

  - Untar the **idl_6.3_x86_64_linux.tar.gz** program into **/usr/local** <code>
cd /usr/local
sudo tar -zxvf ~/down/idl/idl_6.3_x86_64_linux.tar.gz
</code>
  - Install the package
    - Run the install script <code>
cd /usr/local/rsi
sudo ./install
</code>
    - Answer the following questions:
      - Execute the install command now? (y/n): **y**
      - Do you accept all of the terms of the preceding license agreement? (y/n): **y**
      - Do you wish to continue? (y/n): **y**
      - Create the symbolic links described in option 1 above? (y/n): **y**
      - Do you want to start the DICOM Network Services at boot time?  Root permissions are required.? (y/n): **n** (no need for this)
  - Make symbolic links to this library to fool it into working <code>
sudo apt-get install libxpm4
cd /usr/lib/x86_64-linux-gnu
ls -l libXp*
sudo ln -s libXpm.so.4.11.0 libXp.so.6
</code>

===== Environment =====

<code>

#idl license and startup
export LM_LICENSE_FILE="$HOME/idl/license.dat"
export IDL_STARTUP="$HOME/idl/idl_startup.pro"

#idl path
export IDL_PATH="<IDL_DEFAULT>"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvis_current"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvis_utils"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvisf"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvis_classic_2019"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvis21"
export IDL_PATH="$IDL_PATH:$HOME/idl/lvis"
export IDL_PATH="$IDL_PATH:$HOME/idl/hdl"
export IDL_PATH="$IDL_PATH:$HOME/idl/local"
</code>

====== IDL Processing ======

===== 2017.04.07 =====

==== Center Image ====

<code>
cd /home/lvis/work/lvisF/idl
idl read_lvisf_merged.pro

# /mnt/nfs/internal
# sweep*.dat

v=where(c.system_status EQ 129)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3
</code>

===== Stationary RX =====

Roll 25000
<code>
@read_lvisf_merged.pro
; select file #1
v=where(c.system_status GT 100 and c.system_status LT 500)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=25000 RX=12104 File=131830'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
</code>

Roll 35000
<code>
@read_lvisf_merged.pro
; select file #2
v=where(c.system_status GT 100 and c.system_status LT 500)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=35000 RX=24712 File=131930'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
</code>

Roll 45000
<code>
@read_lvisf_merged.pro
; select file #3
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod LT 48080)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=45000 RX=35836 File=132030'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
</code>

Roll 55000
<code>
@read_lvisf_merged.pro
; select file #3
c1=c
@read_lvisf_merged.pro
; select file #4
c2=c
c=[c1,c2] & c1=0 & c2=0
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod GT 48080 and c.time_gps_sod LT 48120)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=55000 RX=47340 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
</code>

60000
<code>
@read_lvisf_merged.pro
; select file #3
c1=c
@read_lvisf_merged.pro
; select file #4
c2=c
c=[c1,c2] & c1=0 & c2=0
i=where(c.time_gps_sod GT 1)
plot,c(i).time_gps_sod,c(i).system_status
v=where(c.system_status GT 100 and c.system_status LT 500 and c.time_gps_sod GT 48120)
pa = peak(c.txawaveform)
pb = peak(c.txbwaveform)
pc = peak(c.txcwaveform)
mtitle='2017.03.31 CHA Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pa(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHB Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pb(v)),psym=2,symsize=0.3

mtitle='2017.03.31 CHC Roll=60000 RX=54054 File=132130'
plot,c(v).roll_beam(2),c(v).pitch_beam(2),/ynoz,/iso,/nodata,title=mtitle,xtitle='Roll (cnt)',ytitle='Pitch (cnt)'
plots,c(v).roll_beam(2),c(v).pitch_beam(2),color=bytscl(pc(v)),psym=2,symsize=0.3
</code>