Task 2.4: Correct the Areal Interpolation Plot's "save" that exports the tile plot and values instead of the polygonal values and Shapefile

[joellenb]

FastTilePlot and FastAreaTilePlot called the same code for menu items including exporting a Shapefile. Split the code such that the the FastTilePlot function will be called for tile plots and the FastAreaTilePlot function will be called for the areal interpolation plots. The new function to save a Shapefile for a FastAreaTilePlot is not yet written.

[systemsgo]

This is working in shapefile_mpas

[lizadams]

using VERDI_2.0_beta_linux64_20190623.tar.gz verified that this works correctly for the CCTM file using the following script.

#!/bin/csh -f

#script for testing command line options

echo 'running verdi_script_areal.csh'
setenv verdi_dir /proj/ie/apps/longleaf/VERDI_2.0_jun_23/VERDI_2.0_beta
cd $verdi_dir

foreach species ( O3 )
./verdi.sh \
#         -f $verdi_dir/data/model/history.2013-07-10.nc \
          -f $verdi_dir/data/model/CCTM_N1a_drydep_O3.20060701 \
         -s "${species}[1]" \
#          -quit
end

File > Save as Shapefile, then loaded it into QGIS to verify that the data matches the Areal Interpolation Plot.

However, If I try to do the same thing with an MPAS file, the data that is exported as a shapefile is the mpas cell data, not the areal interpolation data.

This shows the precipw variable in QGIS, and it doesn't match the areal interpolation plot. It matches the Tile Plot.

[lizadams]

The shapefile export for the MPAS areal interpolation plot doesn't work. The Shapefile that is saved is of the Tile Plot, but it does seem to use the extent of the Areal Interpolation plot (the huc03.shp covers only the SE region.) The data doesn't appear to match the areal interpolation plot because there are values over the atlantic ocean, and the areal interpolation plot only has values over the continental US. Areal interpolation plot created using the following script:

!/bin/csh -f

script for testing command line options

echo 'running verdi_script_areal.csh' setenv verdi_dir /applications/VERDI_2.1 cd $verdi_dir

foreach species ( qv ) ./verdi.sh \ -f $verdi_dir/data/model/history_qv_only.2016-07-04_00UTC.nc \ -s "${species}[1]" \

-quit

end Then you need to load the area file /applications/VERDI_2.1/data/hucRegion/huc03.shp

Then create an areal interpolation plot. Then save the areal interpolation plot as a shapefile.

Load saved shapefile into QGIS and then select Properties Change Simple Fill to Categorized Then Select Classify A plot is generated as follows:

I would expect that if the plot was of the areal interpolation values, then there wouldn't be any cells with color over the ocean.

@systemsgo - Please take a look and see if you agree that this still needs to be fixed.