#!/bin/bash
#
# convert pmean sh model to grids removing mean
#
m=../models/pmean.31.m.ab
z=2850
rm depths.dat 2> /dev/null
p=3
i=1
while [ $z -ge 50 ];do
    echo $z >> depths.dat
    extract_layer $m $z | abconvert 0 8 | shsyn 0.5 0 dvp.$i
    if [ $p -eq 3 ];then
	grd2map dvp.$i.grd xxx "@~d@~v@-P@- [%]" 1 -JN170/7i "" "z = $z km" 0.5 -I    
	/usr/bin/convert -density 70 dvp.$i.ps dvp.$i.png
	p=0
    fi
    ((p=p+1))
    ((z=z-100))
    ((i=i+1))
done
