#!/bin/bash
#
# convert smean sh model to grids removing mean
#
m=../models/smean.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 dvs.$i
    if [ $p -eq 3 ];then
	grd2map dvs.$i.grd xxx "@~d@~v@-S@- [%]" 1 -JN170/7i "" "z = $z km" 0.5 -I    
	/usr/bin/convert -density 70 dvs.$i.ps dvs.$i.png
	p=0
    fi
    ((p=p+1))
    ((z=z-100))
    ((i=i+1))
done

tar --gzip --verbose --create --file ~/tmp/smean_grd.tgz dvs.*grd convert README depths.dat dvs*png
