#!/bin/sh

n=$(gawk '$1=='$1'{print $2}' MAXTIMES)
cdo seltimestep,$n -divc,1e7 tmp/phi500_dom.ano.ERA5.nc toto.nc
d=$(cdo showdate toto.nc)
tit="Phi 500 Anomaly $d"

f=toto.nc

# pars="-R-11/25/62/69r -Jd20/55/22/66/0.34"; tloc="-48 63"
pars="-R-30/24/40/59r -Jd-3/55/22/66/0.34"; tloc="-52 60"
makepalette11 palette.warm11 -2.5 -1.8 -1.4 -1.0 -0.6 -0.2 0.2 0.6 1.0 1.4 1.8 2.5

cdo outputtab,lon,lat,value $f | gawk 'FNR>1&&$3<1e10{print}' > DATA

gmt pscontour DATA ${pars} -I -CPAL.cpt -P -K -X0.5 -Y10 > .ps
gmt pscoast ${pars} -Bf20g20 -A8000 -W1 -Dl -N1 -O -K >> .ps
gmt psscale -D15.3/7/14/0.5 -CPAL.cpt -L -O -K >> .ps
gmt psxy ANALOG ${pars} -W3 -O -K >> .ps
echo ${tloc} 20 0 1 5 $tit | gmt pstext ${pars} -W -N -O >> .ps
mv .ps figscore.ps
gs figscore.ps