#!/bin/sh ##### # get the A0 and ellipticity from the ROCK/SREGN.ASC and SOIL/SREGN.ASC ##### # SREGN.ASC # RMODE NFREQ PERIOD(S) FREQUENCY(Hz) C(KM/S) U(KM/S) AR GAMMA(1/KM) ELLIPTICITY # 0 1 0.10000 10.000 3.3869 3.3814 0.81972E-01 0.20280E-01 0.67126 # 0 2 0.10476 9.5456 3.3868 3.3814 0.78247E-01 0.19358E-01 0.67126 # 0 3 0.10975 9.1116 3.3867 3.3814 0.74689E-01 0.18478E-01 0.67126 ##### CNT=0 CNTMAX=90 rm -f t.txt while [ $CNT -lt $CNTMAX ] do CNT=`echo $CNT | awk '{print $1 + 1}' ` PER=`grep " ${CNT} " SOIL/SREGN.ASC | awk '{print $3}' ` C=`grep " ${CNT} " SOIL/SREGN.ASC | awk '{print $5}' ` AR=`grep " ${CNT} " SOIL/SREGN.ASC | awk '{print $7}' ` E=`grep " ${CNT} " SOIL/SREGN.ASC | awk '{print $9}' ` PERH=`grep " ${CNT} " ROCK/SREGN.ASC | awk '{print $3}' ` CH=`grep " ${CNT} " ROCK/SREGN.ASC | awk '{print $5}' ` ARH=`grep " ${CNT} " ROCK/SREGN.ASC | awk '{print $7}' ` EH=`grep " ${CNT} " ROCK/SREGN.ASC | awk '{print $9}' ` echo $CNT $PER $AR $E $PERH $ARH $EH $C $CH>> t.txt done cat > t.f << EOF 1000 continue read(5,*,end=2000)i,per,ar,e,perh,arh,eh,c,ch ratior = (e/eh)*sqrt(ar/arh) ratioz = sqrt(ar/arh) ratiotr = (e/eh)*sqrt((c*ar)/(ch*arh)) ratiotz = sqrt((c*ar)/(ch*arh)) write(6,*)i,per,ar,e,perh,arh,eh,ratior,ratioz,ratiotz,ratiotr go to 1000 2000 continue end EOF gfortran t.f a.out < t.txt > t.out rm a.out cat t.out | awk '{print 1./$2, $10}' > tz.xy cat t.out | awk '{print 1./$2, $11}' > tr.xy cat > lcmdfil << EOF 'tz.xy' 1095 0.05 'NO' 0.05 'Rz' 'tr.xy' 1005 0.05 'NO' 0.05 'Rr' EOF genplt -XMIN 0.1 -XMAX 20 -XLIN -YMIN -40 -YMAX 40 -L lcmdfil -TX 'Frequency (Hz)' -TY ' ' -LPOS 'BL' mv GENPLT.PLT AMP.PLT ##### # plot absolute values ##### cat t.out | awk '{if($10 < 0) print 1./$2, -1*$4; else print 1./$2, $10}' > tz.xy cat t.out | awk '{if($11 < 0) print 1./$2, -1*$4; else print 1./$2, $11}' > tr.xy cat > lcmdfil << EOF 'tz.xy' 1095 0.05 'NO' 0.05 "Rz" 'tr.xy' 1005 0.05 'NO' 0.05 "Rr" EOF genplt -XMIN 0.1 -XMAX 20 -XLIN -YMIN -0 -YMAX 30 -L lcmdfil -TX 'Frequency (Hz)' -TY ' ' -LPOS 'BL' mv GENPLT.PLT ABSAMP.PLT DOPLTPNG AMP.PLT ABSAMP.PLT