david.rabine

User Tools

Site Tools

Trace: • idl
code:work:lvisf:processing:idl:start

Table of Contents

IDL Installing 6.3

  1. Untar the idl_6.3_x86_64_linux.tar.gz program into /usr/local 
    cd /usr/local
sudo tar -zxvf ~/down/idl/idl_6.3_x86_64_linux.tar.gz
  2. Install the package
    1. Run the install script 
      cd /usr/local/rsi
sudo ./install
    2. Answer the following questions:
      1. Execute the install command now? (y/n): y
      2. Do you accept all of the terms of the preceding license agreement? (y/n): y
      3. Do you wish to continue? (y/n): y
      4. Create the symbolic links described in option 1 above? (y/n): y
      5. Do you want to start the DICOM Network Services at boot time? Root permissions are required.? (y/n): n (no need for this)
  3. Make symbolic links to this library to fool it into working 
    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

Environment

#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"

IDL Processing

2017.04.07

Center Image

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

Stationary RX

Roll 25000 

@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

Roll 35000 

@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

Roll 45000 

@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

Roll 55000 

@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

60000 

@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/work/lvisf/processing/idl/start.txt · Last modified: 2022/03/09 15:29 by david