# -----------------------------------
# Makefile - Created by Brad Artman
# Wed Mar 14 18:23:55 2001.
# -----------------------------------
# Usage:
# Parameters:
# Description:manipulate WSM data found at
# http://www-wsm.physik.uni-karlsruhe.de/
# -----------------------------------
include ${SEPINC}/SEP.top
HTMLDIR = /net/sepwww/sepwww/pub/sep/brad/Stress
RESDIR = ./Fig
SRCDIR = ./src
PARDIR = ./par
BINDIR = ./bin/${MTYPE}
OBJDIR = ./obj/${MTYPE}
INCDIR = ./inc/${MTYPE}
DATDIR= ./Data
R = ${RESDIR}
B = ${BINDIR}
P = ${PARDIR}
D = ${DATDIR}
#UF90FLAGS = -Mbounds
UF90INCLUDES = -I${INCDIR}
ifeq (${ARCH},LINUX) # LINUX
F90C = pgf90 -v
F90LD = pgf90 -v
UF90LIBS = -lsepgeef90 -lsep2df90
else # SGI64
F90C = f90 -64
F90LD = f90
UF90LIBS = -lgeef90 -lsep2df90 ${GEELIB}
endif
RESULTSCR =
RESULTSNR = wsm_nolegend hist_rand regime
RESULTSER = world plates coast var_0 hist_0 hist_af plot_0 plot_af plot_nContNFSp1 plot_rand_af plot_cont plot_Namerica plot_nNamerica plot_ocean plot_nocean plot_rand_0 plot_nNFSpacific plot_TFSpacific
.PRECIOUS = ${R}/regime.ps ${D}/coast.dat ${D}/plates.dat ${D}/wsm2000-1.csv ${D}/wsm.ref.codes.txt
GP = Graph min1=-180 max1=180 min2=-5 > /dev/null
GA = Graph min1=0 max1=180 > /dev/null
GW = Graph min1=-180 max1=180 min2=-90 max2=90 symbol="." title=" " > /dev/null
G = Grey min1=90 max1=-90 min2=-180 max2=180 color=v title=" " > /dev/null
#compile%: ${B}/variogram.x ${B}/mapmaker.x
# echo "stress"
# echo $*
#___________________________________________________________________________________________
#------ Make seplib format from World Stress Map database format
#------ original data file awk'ed with script: ./src/get.awk
#------ original data file in: ./Data/wsm2000-1.csv
data_0.HH:${B}/wsm2seplib.x ${D}/wsm2000-1.csv
gawk -f src/get.awk Data/wsm2000-1.csv > ws.dat
src/changereg.sed > ${D}/wsm.dat
${B}/wsm2seplib.x fname=${D}/wsm.dat > $@
ws.dat
#_____________________________________________________________________________________________
#---- mapping density of intersecting data azimuths, where (biggest value)=(fewest intersections)
map%.HH: # data_0.HH ${B}/mapmaker.x ${P}/model%.P
< data_af.HH ${B}/mapmaker.x par=${P}/model$*.P datafile=actualdata$*.H > $@
touch par/plot$*.P
#______________________________________________________________________________________________
#---- Plots colorscale intersection map w/ world geography overlays
${R}/plot%.v ${R}/plot%.v3: #${P}/plot%.P ${R}/world.v map%.HH
<map$*.HH Reverse which=1 | ${G} title=$* par=${P}/plot$*.P out=${R}/plotsc$*.v wantscalebar=y
<map$*.HH Reverse which=1 | ${G} title=$* par=${P}/plot$*.P out=${R}/plot2$*.v
<map$*.HH Reverse which=1 | ${G} title=$* par=${P}/plot$*.P out=${R}/plot3$*.v wantscalebar=y color=a
< actualdata$*.H Window n2=1 f2=0 > lon.H
< actualdata$*.H Window n2=1 f2=1 > lat.H
Cmplx lon.H lat.H | ${GW} plotcol=0 plotfat=7 title=" " out=${R}/loc$*.v
vp_Overlay ${R}/plot2$*.v ${R}/loc$*.v > ${R}/plotlocs$*.v
vp_Overlay ${R}/plot2$*.v ${R}/world.v ${R}/loc$*.v> ${R}/plotworld$*.v
vp_SideBySideAniso ${R}/plotsc$*.v ${R}/plotworld$*.v > ${R}/plot$*.v
vp_Movie ${R}/plotsc$*.v ${R}/plot3$*.v ${R}/plotlocs$*.v ${R}/plotworld$*.v > ${R}/plot$*.v3
Rm lat.H lon.H
#___________________________________________________________________________________________
#------ Make variogram of azimuth dotproducts (almost) vs. angular distance
var_0.H: #data_0.HH ${B}/variogram.x ${P}/model%.P
<data_0.HH ${B}/variogram.x par=${P}/model$*.P > $@
${R}/var%.v:# var.H
< var%.H Graph title="Variogram of WSM data"$* min1=0 min2=-10 out=${R}/$@ >/dev/null
#____________________________________________________________________________________________
#------ Make histogram of az's and location to show data bunching
${R}/hist%.v: data_0.H
<data_0.HH Window n2=1 f2=2 | Histogram | ${GA} title="Azimuth Histogram" out=az$*.v
<data_0.HH Window n2=1 f2=0 | Histogram | ${GP} title="Latitude Histogram" out=lat$*.v
<data_0.HH Window n2=1 f2=1 | Histogram | ${GP} title="Longitude Histogram" out=lon$*.v
vp_OverUnderAniso az$*.v lat$*.v lon$*.v > temp$*.v
vp_annotate <temp$*.v text=${D}/georef.txt > $@
Rm temp$*.v lat$*.v lon$*.v az$*.v
#${R}/hist_af.v: data_af.HH
# <data_af.HH Window n2=1 f2=2 | Histogram | ${GA} title="Azimuth Histogram" out=az.v
# <data_af.HH Window n2=1 f2=0 | Histogram | ${GP} title="Latitude Histogram" out=lat.v
# <data_af.HH Window n2=1 f2=1 | Histogram | ${GP} title="Longitude Histogram" out=lon.v
# vp_OverUnderAniso az.v lat.v lon.v > temp.v
# vp_annotate <temp.v text=${D}/afghan.txt > ${R}/$@
# Rm az.v lat.v lon.v temp.v
#___________________________________________________________________________________________
#---- insert random azimuths at data locations to test
mmap_rand%.HH:data%.HH ${B}/mapmaker.x ${P}/model%.P
< data$*.HH ${B}/mapmaker.x par=${P}/model$*.P randaz=randaz.H datafile=data$*.H > $@
Mv $@ map_rand$*.H
Rm $@
#_____________________________________________________________________________________
#----make data subset(s)
#---- first 100 data points happen to be around afghanistan (40, 20)(lat,lon)
#data_af.HH: data_0.HH
# < data_0.HH Window n1=100 f1=0 > $@
#______________________________________________________________________________________
#------plate boundaries and coastlines
plates.H: ${D}/plates.dat ${B}/plabo2seplib.x
${B}/plabo2seplib.x fname=${D}/plates.dat npick=5535> $@
${R}/world.v: #coast.HH plates.H
< plates.H Window n2=1 f2=0 > lon.H
< plates.H Window n2=1 f2=1 > lat.H
Cmplx lon.H lat.H | ${GW} plotcol=0 plotfat=4 title=" " out=${R}/plates.v
< coast.HH Window n2=1 f2=0 > lat.H
< coast.HH Window n2=1 f2=1 > lon.H
Cmplx lon.H lat.H | ${GW} plotcol=0 plotfat=4 title=" " out=${R}/coast.v
vp_Overlay ${R}/coast.v ${R}/plates.v title=" " > $@
Rm lat.H lon.H
done: paper.tex
gmake paper.ps
rm -fR ${HTMLDIR}
mkdir ${HTMLDIR}
gmake html
# cp paper.pdf ${HTMLDIR}/.
# to make an adobe presentation type:
#make LATOPTS=option=web aread
clean: jclean
xclean:
cd bin/${MTYPE}; *.x
cd obj/${MTYPE}; *.o
include ${SEPINC}/SEP.bottom