Nearest: Add test for ellipsoid lambert

This commit is contained in:
Shahram Najm 2020-06-20 14:20:25 +01:00
parent c438c13c84
commit 93911a8830
1 changed files with 4 additions and 4 deletions

View File

@ -52,7 +52,6 @@ ${tools_dir}/grib_ls -l 50,0 $tempGrib
cat > $tempFilter <<EOF
set gridType="lambert";
set numberOfDataPoints=294000;
set shapeOfTheEarth=6;
set Nx=588;
set Ny=500;
set latitudeOfFirstGridPoint=40442000;
@ -63,7 +62,7 @@ cat > $tempFilter <<EOF
set Dy=2499000;
set Latin1=46401000;
set Latin2=46401000;
set shapeOfTheEarth=2;
set shapeOfTheEarth=2; # oblate earth
set numberOfValues=294000;
write;
EOF
@ -76,10 +75,11 @@ if [ ! -f "$tempGrib" ]; then
fi
grib_check_key_equals $tempGrib 'earthIsOblate,earthMinorAxisInMetres,earthMajorAxisInMetres' '1 6356775 6378160'
# Invoke Geoiterator on the newly created GRIB file
# Invoke Geoiterator on the oblate lambert GRIB
${tools_dir}/grib_get_data $tempGrib > $tempOut
# Nearest neighbour on the oblate lambert GRIB
${tools_dir}/grib_ls -l 40.44,353.56 $tempGrib
# Clean up