#!/bin/sh #------------------------------------------------------------- MDL="${HOME}/hysplit" OUT="${MDL}/working" MET="${HOME}/Tutorial/volcano" cd $OUT if [ ! -f ASCDATA.CFG ]; then echo "-90.0 -180.0" >ASCDATA.CFG echo "1.0 1.0" >>ASCDATA.CFG echo "180 360" >>ASCDATA.CFG echo "2" >>ASCDATA.CFG echo "0.2" >>ASCDATA.CFG echo "'$MDL/bdyfiles/'" >>ASCDATA.CFG fi echo "### $0 ###" #-------------------------------------------------------------- syr=10 smo=04 sda=14 shr=00 run=36 ztop=32000.0 data="apr1420.bin" #---------------------------------------------------------- echo "$syr $smo $sda $shr " >CONTROL echo "2 ">>CONTROL echo "63.63 -19.62 100.0 ">>CONTROL echo "63.63 -19.62 6000.0 ">>CONTROL echo "$run " >CONTROL1 echo "0 ">>CONTROL1 echo "$ztop ">>CONTROL1 echo "1 ">>CONTROL1 echo "$MET/ ">>CONTROL1 echo "$data ">>CONTROL1 echo "4 ">>CONTROL1 echo "p006 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "p020 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "p060 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "p200 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "1 ">>CONTROL1 echo "45.0 0.0 ">>CONTROL1 echo "0.25 0.25 ">>CONTROL1 echo "90.0 360.0 ">>CONTROL1 echo "$OUT/ ">>CONTROL1 echo "volcano.bin ">>CONTROL1 echo "1 ">>CONTROL1 echo "10000 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "00 00 00 00 00 ">>CONTROL1 echo "00 06 00 ">>CONTROL1 echo "4 ">>CONTROL1 echo "0.6 2.5 1.0 ">>CONTROL1 echo "0.0 0.0 0.0 0.0 0.0 ">>CONTROL1 echo "0.0 8.0E-05 8.0E-05 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "2.0 2.5 1.0 ">>CONTROL1 echo "0.0 0.0 0.0 0.0 0.0 ">>CONTROL1 echo "0.0 8.0E-05 8.0E-05 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "6.0 2.5 1.0 ">>CONTROL1 echo "0.0 0.0 0.0 0.0 0.0 ">>CONTROL1 echo "0.0 8.0E-05 8.0E-05 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "20.0 2.5 1.0 ">>CONTROL1 echo "0.0 0.0 0.0 0.0 0.0 ">>CONTROL1 echo "0.0 8.0E-05 8.0E-05 ">>CONTROL1 echo "0.0 ">>CONTROL1 echo "0.0 ">>CONTROL1 cat CONTROL1 >>CONTROL echo "&SETUP " >SETUP.CFG echo "kmsl=0, " >>SETUP.CFG echo "ninit=1, " >>SETUP.CFG echo "ndump=12, " >>SETUP.CFG echo "ncycl=12, " >>SETUP.CFG echo "numpar=5000, " >>SETUP.CFG echo "maxpar=50000, " >>SETUP.CFG echo "pinpf = 'PARINIT', " >>SETUP.CFG echo "poutf = 'PARDUMP', " >>SETUP.CFG echo "efile = 'EMITIMES', " >>SETUP.CFG echo "/ " >>SETUP.CFG echo "YYYY MM DD HH DURATION(hhhh) #RECORDS " >EMITIMES echo "YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w) " >>EMITIMES echo "2010 04 14 00 0012 8 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 100.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 100.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 100.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 100.0 0.67 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 6000.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 6000.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 6000.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 14 00 00 1200 63.63 -19.62 6000.0 0.67 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 0012 8 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 100.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 100.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 100.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 100.0 0.67 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 12000.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 12000.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 12000.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 14 12 00 1200 63.63 -19.62 12000.0 0.67 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 0012 8 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 100.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 100.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 100.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 100.0 0.67 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 6000.0 0.01 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 6000.0 0.07 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 6000.0 0.25 0.0 0.0 " >>EMITIMES echo "2010 04 15 00 00 1200 63.63 -19.62 6000.0 0.67 0.0 0.0 " >>EMITIMES #---------------------------------------------------------- rm -f volcano.bin ${MDL}/exec/hycs_std echo "'TITLE&','### $0: Eruption using EMITIMES ### &'" >LABELS.CFG rm -f parxplot.ps ${MDL}/exec/parxplot -iPARDUMP -k1 -z80 -j${MDL}/graphics/arlmap if [[ "$OSTYPE" == "darwin"* ]]; then open parxplot.ps else gs parxplot.ps fi echo "Enter to continue ..."; read x rm -f concplot.ps ${MDL}/exec/concplot -ivolcano.bin -j${MDL}/graphics/arlmap -s0 -z80 \ -c4 -uug -v100+50+20+10 -x1.0E+16 if [[ "$OSTYPE" == "darwin"* ]]; then open concplot.ps else gs concplot.ps fi