The script is given below. It is will commented. The only difference between this version and that for GMT Versions 5 and up is that DOSOLUTION is replaced by DOSOLUTION5.
#!/bin/bash set -x ##### # set up directory structure for an event that we get from the IRIS ##### ##### # test usage ##### if [ $# -ne 16 ] then echo DOSETUP YEAR MO DY HR MN SC MSC LAT LON DEP MAG REG NEIC FELTID STATE MYPWD echo DOSETUP 2006 02 05 08 17 04 000 66.25 -143.02 18.60 5.1 CUS iuae Xiuae_06 Alaska /backup/rbh/MOMENT_TENSOR/MECH.NA echo $* exit 0 fi ##### # define internal variables: the quotes are required since some fields can have # blank spacing within" ##### YEAR="${1}" MO="${2}" DY="${3}" HR="${4}" MN="${5}" SEC="${6}" MSEC="${7}" LAT="${8}" LON="${9}" DEP="${10}" MAG="${11}" REG="${12}" NEIC="${13}" FELTID="${14}" STATE="${15}" MYPWD="${16}" cd ${MYPWD} if [ "${FELTID}" != "NONE" ] then ISFELT="y" else ISFELT="NO" fi ##### # now run three scripts that set up everything # DOSETUP - sets up the directory structure for the inversion ##### ##### # define the DefaultVelocity Model ##### case ${REG} in WUS) GMODEL="WUS.REG";VMODEL="WUS";FIDREG="imw";REGIONTEXT=" West " ;; CUS) GMODEL="CUS.REG";VMODEL="CUS";FIDREG="cus";REGIONTEXT="Central and Southeastern US" ;; GSKAN) GMODEL="GSKAN.REG";VMODEL="GSKAN";FIDREG="ok" ;REGIONTEXT="S. KansasS" ;; *) echo Improper REGION specification echo Select on of CA PNW IMW AK NE CUS ECAN echo DOSETUP YEAR MO DY HR MN SC MSC LAT LON DEP MAG REG NEIC FELTID STATE/COUNTRY echo DOSETUP 2006 02 05 08 17 04 000 66.25 -143.02 18.60 5.1 CUS iuae Xiuae_06 "Yukon" ;; esac echo ${REG} ${GMODEL} ${VMODEL} ${FIDREG} ${REGIONTEXT} echo $* XXXDIR="${MYPWD}/0XXXREG" echo $MYPWD $PROTODIR ${XXXDIR} ##### # now define the directory name ##### DIR=${YEAR}${MO}${DY}${HR}${MN}${SEC} echo DIR $DIR echo PROTODIR $PROTODIR echo 0XXXDIR $XXXDIR ##### # Create and populate the directory structure ##### if [ ! -d ${DIR} ] then mkdir ${DIR} fi ##### # indicate the Velocity model in a file ##### echo ${VMODEL} > ${MYPWD}/${DIR}/VMODEL.used cd ${XXXDIR} cp -pr * ${MYPWD}/${DIR} cd ${MYPWD} ##### # NOW THIS IS VERY VERY SPECIFIC - MODIFY THE SCRIPTS FOR THIS REGION ##### # get correct model into the scripts ##### cd ${MYPWD}/${DIR} for i in */DOSTA */DOGRID do echo Editing $i ed $i << EOH g/GMODEL/s//${GMODEL}/gp w q EOH done for i in */DOHTML */DOPLTRAD do echo Editing $i ed $i << EOH g/VMODEL.mod/s//${VMODEL}.mod/gp g/GMODEL/s//${GMODEL}/gp g/VMODEL/s//${VMODEL}/gp w q EOH done ##### # go to the HTML directory to define the processing parameters ##### cd ${MYPWD}/${DIR}/HTML.REG cat > html.tmp << EOF REGION="${FIDREG}" FELTEVID="${FELTID}" TITLE="" EVID="${NEIC}" DATE="${YEAR}/${MO}/${DY}" TIME="${HR}:${MN}:${SEC}" LAT="${LAT}" LON="${LON}" DEPTH="${DEP}" MAG="${MAG}" FELT="${ISFELT}" DIR="${DIR}" STATE="${STATE}" YEAR="${YEAR}" MO="${MO}" DY="${DY}" HR="${HR}" MN="${MN}" SEC="${SEC}" MSEC="${MSEC}" EOF set -x ed DOHTML << EOF /IDHERE/ .r html.tmp w q EOF echo Processing ${DIR}
Last changed May 21, 2022