#!/bin/bash

set -x

yvals="0,4000,...,4804000"
xvals="0,4000,...,6404000"

#xvals="2012000,2014000,...,4790000" 
#yvals="688000,690000,...,2486000"

infile=era5_2008.nc
outfile=era5_2008_intrp.nc
interpolation=coord_nearestneighbor
echo $xvals
echo $yvals
echo $outfile
echo $interpolation


fimex --interpolate.projString="+proj=stere +R=6371000.0 +lat_0=90.0 +lat_ts=60.0 +x_0=4180000.0 +y_0=2570000.0 +lon_0=58.0" --input.file=${infile}  --interpolate.xAxisUnit=m --interpolate.yAxisUnit=m --interpolate.xAxisValues=${xvals} --interpolate.yAxisValues=${yvals} --output.file=$outfile --interpolate.method=$interpolation

set +x
