#!/bin/sh # # qc_script # # Number of seconds to sleep between each check for new updates SLEEPTIME=1 TRACE=1 #The name of the LinAcq datafile is stored in 'which.file' #The LOGFILE variable holds the name of this datafile; it's #reloaded for every loop (see below). read LOGFILE < "which.file" #The FILE variable holds the name of the last shot file. See example usage below. PRESSUREFILE="/home/linacq/OBS2003/pressure.temp" LASTPRESSUREFILE="/home/linacq/OBS2003/pressure_last100.temp" PRESSPLOTCMD="/home/linacq/OBS2003/pressplot.cmd" PRESSPLOTPNGCMD="/home/linacq/OBS2003/pressplotpng.cmd" LASTPRESREAD="/var/www/html/agp/lastreading.txt" SHOWMODE="web" #echo "Processing $LOGDIRECTORY$LOGFILE" #rm nav.temp #rm $PRESSUREFILE #========================================================================================= # L O O P S T A R T #========================================================================================= while [ 0 ] do read LOGFILE < "which.file" FILE=data.su`./wqctrig ${SLEEPTIME}` #echo $FILE # if [ "$FILE" = "data.su0" ] # then # FILE=data.su1 # else # FILE=data.su0 # fi #echo "New file name: $FILE" # zap xwigb # zap ximage zap xgraph # killall gnuplot_x11 #--------------------------------------------------------------------------------------------------------------------- # Plot ship's trackline based on Easting, Northing header values (SEG-Y trace header words sx,sy) #--------------------------------------------------------------------------------------------------------------------- # sugethw < $FILE key=sx,sy | sed -n '/^./ p' | cut -c8-13,15,22-28 >> nav.temp # cat nav.temp # gnuplot -persist navplot.cmd #-------------------------------------------------------------------------------------------------------------- # Process whole data file #-------------------------------------------------------------------------------------------------------------- ###supef < $LOGFILE | sugain agc=1 tpow=0.7 | \ ### suximage xbox=10 ybox=10 wbox=830 cmap=hsv5 perc=95 \ ### title="Mini-streamer data" & ###supef < $LOGFILE | \ ###sufilter f=2000,2500,4000,4500 | \ # ####supef < $LOGFILE | \ # # sufilter < $LOGFILE f=7.5,15,40,60 | \ # sugain jon=1 | \ # suximage perc=95 xbox=10 ybox=10 wbox=830 cmap=hsv0 \ # label1="Two-way traveltime [s]" label2="Shot no." windowtitle="File: $LOGFILE" & # ##convert plot1.ps /var/www/html/mosby/plot1.png # # #-------------------------------------------------------------------------------------------------------------- # Process whole data file : Generate frequency plot vs. event no #-------------------------------------------------------------------------------------------------------------- # X version # suspecfx < $LOGFILE | \ # suximage xbox=850 ybox=10 wbox=410 hbox=460 cmap=hsv6 x1end=200 \ # label1="Frequency [Hz]" label2="Shot no." windowtitle="File: $LOGFILE" & #======================================================= ############## N A X Y S T E S T ################## #======================================================= # # Process whole file, display wiggle trace of one channel (selected by min/max value) # #suwind key=tracf min=1 max=1 < $LOGFILE | \ # #sugain agc=1 | \ # suxwigb xbox=10 ybox=10 wbox=730 hbox=500 perc=95 \ # label1="Time [s]" label2="Trace no." & # # Process whole file, display spectrum of all traces: # #suspecfx < $LOGFILE | \ # suximage xbox=10 ybox=600 wbox=730 hbox=450 cmap=hsv4 x1end=250 \ # label1="Frequency [Hz]" label2="Trace no." windowtitle="File: $LOGFILE" & # # # Geophone X,Y,Z time signal # #suwind key=tracf min=1 max=3 < $FILE | \ # suxgraph style=normal -geometry 730x490+10+10 \ # grid1=dot grid2=dot \ # title="Geophones time signal. X=red, Y=green, Z=blue" \ # label1="Time [s]" \ # label2="Amplitude" & # # # Hydrophone (Bottom sensor) time signal # #suwind key=tracf j=4 s=4 < $FILE | \ # #suxwigb xbox=1030 ybox=10 wbox=500 hbox=200 \ # suxgraph style=normal -geometry 500x200+1030+10 \ # grid1=dot grid2=dot \ # title="Hydrophone (bottom) time signal" \ # label1="Time [s]" \ # label2="Amplitude" & # # # Hydrophone (Bottom sensor) frequency domain # #suwind key=tracf j=4 s=4 < $FILE | \ # #suspecfx | \ # sufft | suamp mode=amp | \ # suop op=db | \ # suxgraph style=normal -geometry 500x200+1030+250 \ # grid1=dot grid2=dot \ # title="Hydrophone (bottom) spectrum" \ # label1="Frequency [Hz]" \ # label2="Amplitude [dB]" & # # # Channel #1 time display # suwind key=tracf min=1 max=1 < $FILE | \ suxgraph style=normal -geometry 500x200+1030+10 \ grid1=dot grid2=dot \ title="Ch #1: X" \ label1="Time [s]" \ label2="Amplitude" & # # Channel #2 time display # #suwind key=tracf min=2 max=2 < $FILE | \ # #suxwigb xbox=1030 ybox=10 wbox=500 hbox=200 \ # suxgraph style=normal -geometry 500x200+1030+250 \ # grid1=dot grid2=dot \ # x1beg=0 x1end=4.6 \ # x2beg=-2.5 x2end=2.5 \ # title="Ch #2: Y" \ # label1="Time [s]" \ # label2="Amplitude" & # # Channel #3 time display # #suwind key=tracf min=3 max=3 < $FILE | \ # suxgraph style=normal -geometry 500x200+1030+490 \ # grid1=dot grid2=dot \ # x1beg=0 x1end=4.6 \ # x2beg=-2.5 x2end=2.5 \ # title="Ch #3: Z" \ # label1="Time [s]" \ # label2="Amplitude" & # # Channel #4 time display # #suwind key=tracf min=4 max=4 < $FILE | \ # suxgraph style=normal -geometry 500x200+1030+730 \ # grid1=dot grid2=dot \ # x1beg=0 x1end=4.6 \ # x2beg=-2.5 x2end=2.5 \ # title="Ch #4: Streamer" \ # label1="Time [s]" \ # label2="Amplitude [Volt]" & # # # Hydrophone (NAXYS) frequency domain # #suwind key=tracf min=2 max=2 < $FILE | \ # #suspecfx | \ # sufft | suamp mode=amp | \ # suop op=db | \ # suxgraph style=normal -geometry 500x200+1030+750 \ # grid1=dot grid2=dot \ # title="Streamerspectrum" \ # label1="Frequency [Hz]" \ # label2="Amplitude [dB]" & # # # # GUN Trigger signal # #suwind key=tracf j=6 s=6 < $FILE | \ # suxgraph style=normal -geometry 500x150+10+820 \ # grid1=dot grid2=dot \ # title="Digital trigger signal - Ch.6" \ # label1="Digital trigger" & # # # # #-------------------------------------------------------------------------------------------------------------- # Process last shot : Plot amplitude vs. time #-------------------------------------------------------------------------------------------------------------- # # # # suxgraph < $FILE style=normal -geometry 415x250+5+740 \ # grid1=dot grid2=dot \ # windowtitle="Last shot time signal" \ # label1="Time [s]" \ # label2="Amplitude [V]" & # #-------------------------------------------------------------------------------------------------------------- # Process last shot : Plot frequency spectrum #-------------------------------------------------------------------------------------------------------------- # # sufft < $FILE | suamp mode=amp | \ # suop op=db | \ # suxgraph style=normal -geometry 415x250+430+740 \ # grid1=dot grid2=dot \ # windowtitle="Last shot frequency spectrum" \ # label1="Frequency [Hz]" \ # label2="Amplitude [dB]" & # # #-------------------------------------------------------------------------------------------------------------- # Process last shot : Plot instantanious frequency (Gabor spectrogram) #-------------------------------------------------------------------------------------------------------------- # # sugabor < $FILE | suximage xbox=850 ybox=740 wbox=410 hbox=495 cmap=hsv6 x2end=150 \ # label1="Two-way traveltime [s]" label2="Frequency [Hz]" windowtitle="Last shot instantaneous frequency (Gabor spec.)" & # # #-------------------------------------------------------------------------------------------------------------- # Process last shot : Plot dB amplitude vs. time #-------------------------------------------------------------------------------------------------------------- # ### sugain < $FILE | \ ### suattributes mode=amp | suop op=db | \ ### suxgraph -geometry 190x530+1790+5 \ ### label1="Traveltime [s]" label2="Amplitude [dB]" title="Amplitude" \ ### grid1=dot grid2=dot nTic2=2 verbose=0 & # #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Experimental stuff # #suxwigb < $FILE xbox=670 ybox=1 wbox=350 title="Wiggle plot last shot" & # #suspecfx < $FILE | \ # suximage xbox=1030 ybox=1 wbox=350 cmap=hsv4 x1end=500 \ # title="Spectrum last shot" & # #suwind key=tracf j=4 s=$TRACE < test.su | \ # suximage xbox=1030 ybox=10 wbox=500 cmap=hsv4 title="Plot of Ch. no. $TRACE" & # # #---------------------------------------------------------------------------------------------------------------------- # Plot Air Gun Pressure vs. Event number (TS, 220503) #---------------------------------------------------------------------------------------------------------------------- # #####sugethw < $FILE key=ep,sut | sed -n '/^./ p' | sed -n 's/[a-zA-Z\t\v\r\n\f ]*=//; s/[a-zA-Z\t\v\r\n\f ]*=//; p' \ ##### | awk '{ print $1" "$2 }' >> $PRESSUREFILE #####tail -n 100 $PRESSUREFILE > $LASTPRESSUREFILE #####tail -n 1 $PRESSUREFILE | awk '{ print $2 }' > $LASTPRESREAD ## sleep 3 #gnuplot -geometry 600x400+680+0 -persist -background white $PRESSPLOTCMD #####gnuplot $PRESSPLOTPNGCMD > /var/www/html/agp/agp.png # # #------------------------------------------------------------------------- # GUNCO test(!) #------------------------------------------------------------------------- # # Plot Gun status ####rm -f /var/www/html/agp/gun.status ####tail -n 1 /home/syncix/gunco.dat | awk -F "," -f /home/syncix/trond.awk > /var/www/html/gunco/gun.status # # TODO: Plot Error / Delay over time ... # #------------------------------------------------------------------------------- # H y d r o p h o n e p l o t s # Used on Yngve Kristoffersen's "Reflector" test April 2003 #------------------------------------------------------------------------------- # Time domain # # suwind key=tracf j=4 s=4 < $FILE | \ # #suxwigb xbox=1030 ybox=10 wbox=500 hbox=300 \ # suxgraph style=normal -geometry 500x300+1030+10 \ # grid1=dot grid2=dot \ # title="Hydrophone time signal" \ # label1="Time [s]" \ # label2="Amplitude" & # # Frequency domain # # suwind key=tracf j=4 s=4 < $FILE | \ # #suspecfx | \ # sufft | suamp mode=amp | \ # suop op=db | \ # suxgraph style=normal -geometry 500x300+1030+350 \ # grid1=dot grid2=dot \ # title="Hydrophone spectrum" \ # label1="Frequency [Hz]" \ # label2="Amplitude [dB]" & # # #=========================================================================== # L O O P E N D #=========================================================================== # done