/*
* Copyright 1998-2009 University Corporation for Atmospheric Research/Unidata
*
* Portions of this software were developed by the Unidata Program at the
* University Corporation for Atmospheric Research.
*
* Access and use of this software shall impose the following obligations
* and understandings on the user. The user is granted the right, without
* any fee or cost, to use, copy, modify, alter, enhance and distribute
* this software, and any derivative works thereof, and its supporting
* documentation for any purpose whatsoever, provided that this entire
* notice appears in all copies of the software, derivative works and
* supporting documentation. Further, UCAR requests that the user credit
* UCAR/Unidata in any publications that result from the use of this
* software or in any product that includes this software. The names UCAR
* and/or Unidata, however, may not be used in any advertising or publicity
* to endorse or promote any products or commercial entity unless specific
* written permission is obtained from UCAR/Unidata. The user also
* understands that UCAR/Unidata is not obligated to provide the user with
* any support, consulting, training or assistance of any kind with regard
* to the use, operation and performance of this software nor to provide
* the user with any updates, revisions, new versions or "bug fixes."
*
* THIS SOFTWARE IS PROVIDED BY UCAR/UNIDATA "AS IS" AND ANY EXPRESS OR
* IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
* WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
* DISCLAIMED. IN NO EVENT SHALL UCAR/UNIDATA BE LIABLE FOR ANY SPECIAL,
* INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING
* FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
* NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
* WITH THE ACCESS, USE OR PERFORMANCE OF THIS SOFTWARE.
*/
package ucar.nc2.dt.grid;
import junit.framework.*;
import ucar.nc2.iosp.grib.GribGridServiceProvider;
import ucar.nc2.dataset.CoordinateAxis;
import ucar.nc2.dataset.NetcdfDataset;
import ucar.nc2.constants.AxisType;
import ucar.nc2.TestAll;
import ucar.nc2.ncml.NcMLReader;
import ucar.nc2.ncml.TestNcML;
import java.util.List;
import java.util.Iterator;
import java.io.ByteArrayInputStream;
import java.io.IOException;
/** Count geogrid objects - sanity check when anything changes. */
public class TestReadandCount extends TestCase {
private static boolean show = false, showCount = true;
public TestReadandCount( String name) {
super(name);
}
private String griddir = TestAll.cdmUnitTestDir +"conventions/";
private String grib1dir = TestAll.cdmUnitTestDir +"formats/grib1/";
private String grib2dir = TestAll.cdmUnitTestDir +"formats/grib2/";
public void testRead1() throws Exception {
doOne(griddir+"avhrr/","amsr-avhrr-v2.20040729.nc", 0, 1, 4, 0);
}
public void testRead() throws Exception {
/* doOne(griddir+"atd-radar/","rgg.20020411.000000.lel.ll.nc", 5, 1, 4, 1);
doOne(griddir+"atd-radar/","SPOL_3Volumes.nc", 3, 1, 4, 1);
doOne(griddir+"gief/","coamps.wind_uv.nc", 2, 1, 4, 1);
//// coards derived
doOne(griddir+"coards/", "cldc.mean.nc", 1, 1, 3, 0);
doOne(griddir+"coards/","inittest24.QRIDV07200.nc", -1, -1, -1, -1); // no "positive" att on level
*/
doOne(griddir+"coards/","inittest24.QRIDV07200.ncml", 1, 1, 3, 1);
// doOne(griddir+"avhrr/","amsr-avhrr-v2.20040729.nc", 4, 1, 4, 1);
doOne(griddir+"csm/","o3monthly.nc", 4, 1, 7, 2);
doOne(griddir+"csm/","ha0001.nc", 35, 3, 5, 2); //
doOne(griddir+"cf/","cf1.nc", 1,1,5,2);
doOne(griddir+"cf/","ccsm2.nc", 107, 3, 5, 2); //
doOne(griddir+"cf/","tomw.nc", 19, 3, 4, 1);
doOne(griddir+"cf/","cf1_rap.nc", 11, 2, 4, 0); // not getting x, y
// doOne("C:/data/conventions/cf/signell/","signell_july2_03.nc", -1, -1, -1, -1); // 2D lat, lon; no x,y
//** doOne(griddir+"cf/","feb2003_short.nc", 14, 4, 4, 1);
doOne(griddir+"cf/","feb2003_short2.nc", 22, 9, 8, 1);
doOne(griddir+"cf/","temperature.nc", 4, 3, 5, 1);
doOne(griddir+"gdv/","testGDV.nc", 30, 1, 4, 1);
doOne(griddir+"gdv/","OceanDJF.nc", 15, 1, 4, 1);
// uses GDV as default
doOne(griddir+"mars/","temp_air_01082000.nc", 1, 1, 4, 1); // uses GDV
//doOne("C:/data/conventions/mm5/","n040.nc", -1, -1, -1, -1); // no Conventions code
doOne(griddir+"m3io/","agg.cctmJ3fx.b312.nh3c1.dep_wa.annual.2001base.nc", 13, 1, 5, 1); // m3io
doOne(griddir+"m3io/","19L.nc", 23, 1, 4, 1); // M3IOVGGrid
//// the uglies
doOne(griddir+"nuwg/", "avn-x.nc", 31, 4, 8, 4);
doOne(griddir+"nuwg/", "2003021212_avn-x.nc", 30, 5, 7, 4);
doOne(griddir+"nuwg/", "avn-q.nc", 22, 7, 9, 6);
doOne(griddir+"nuwg/", "eta.nc", 28,9,11,8);
doOne(griddir+"nuwg/", "ocean.nc", 5, 1, 3, 0);
doOne(griddir+"nuwg/", "ruc.nc", 31,5, 6, 3);
doOne(griddir+"nuwg/", "CMC-HGT.nc", 1, 1, 3, 0); // */
doOne(griddir+"wrf/","wrfout_v2_Lambert.nc", 57, 11, 8, 3);
doOne(griddir+"wrf/","wrf2-2005-02-01_12.nc", 60, 11, 8, 3);
doOne(griddir+"wrf/","wrfout_d01_2006-03-08_21-00-00", 70, 11, 8, 3);
doOne(griddir+"wrf/","wrfrst_d01_2002-07-02_12_00_00.nc", 162, 11, 8, 3);
doOne(griddir+"awips/","19981109_1200.nc", 36, 13, 14, 11);
doOne(griddir+"awips/","awips.nc", 38, 12, 13, 10); //
doOne(griddir+"ifps/","HUNGrids.netcdf", 26, 26, 29, 0); // *
// our grib reader */
doOne( grib1dir ,"AVN.wmo", 22, -1, -1, -1);
doOne( grib1dir,"RUC_W.wmo", 44,-1, -1, -1);
doOne( grib1dir,"NOGAPS-Temp-Regional.grib", 1, -1, -1, -1); // */
doOne( grib1dir,"eta.Y.Q.wmo", 25, -1, -1, -1);
doOne( grib2dir,"ndfd.wmo", 1, -1, -1, -1);
// radar mosaic
doOne( grib1dir,"radar_national.grib", 1, 1, 3, 0);
doOne( grib1dir,"radar_regional.grib", 1, 1, 3, 0);
// redo grib files, forcing new index
GribGridServiceProvider.forceNewIndex( true );
doOne( grib1dir,"AVN.wmo", 22, -1, -1, -1);
doOne( grib1dir,"RUC_W.wmo", 44, -1, -1, -1);
doOne( grib1dir,"eta.Y.Q.wmo", 25, -1, -1, -1);
doOne( grib2dir,"ndfd.wmo", 1, -1, -1, -1);
GribGridServiceProvider.forceNewIndex( false );
}
static void doOne(String dir, String filename, int ngrids, int ncoordSys, int ncoordAxes, int nVertCooordAxes) throws Exception {
System.out.println("test read GridDataset = " + dir + filename);
ucar.nc2.dt.grid.GridDataset gridDs = GridDataset.open(dir + filename);
int countGrids = gridDs.getGrids().size();
int countCoordAxes = gridDs.getNetcdfDataset().getCoordinateAxes().size();
int countCoordSys = gridDs.getNetcdfDataset().getCoordinateSystems().size();
// count vertical axes
int countVertCooordAxes = 0;
List axes = gridDs.getNetcdfDataset().getCoordinateAxes();
for (int i = 0; i < axes.size(); i++) {
CoordinateAxis axis = (CoordinateAxis) axes.get(i);
AxisType t = axis.getAxisType();
if ((t == AxisType.GeoZ) || (t == AxisType.Height) || (t == AxisType.Pressure) )
countVertCooordAxes++;
}
Iterator iter = gridDs.getGridsets().iterator();
while (iter.hasNext()) {
GridDataset.Gridset gridset = (GridDataset.Gridset) iter.next();
GridCoordSys gcs = gridset.getGeoCoordSys();
//if (gcs.hasTimeAxis())
// System.out.println(" "+gcs.isDate()+" "+gcs.getName());
}
if (showCount) {
System.out.println(" grids=" + countGrids + ((ngrids < 0) ? " *" : ""));
System.out.println(" coordSys=" + countCoordSys + ((ncoordSys < 0) ? " *" : ""));
System.out.println(" coordAxes=" + countCoordAxes + ((ncoordAxes < 0) ? " *" : ""));
System.out.println(" vertAxes=" + countVertCooordAxes + ((nVertCooordAxes < 0) ? " *" : ""));
}
if (ngrids >= 0)
assert ngrids == countGrids : "Grids " + ngrids + " != " + countGrids;
if (ncoordSys >= 0)
assert ncoordSys == countCoordSys : "CoordSys " + ncoordSys + " != " + countCoordSys;
if (ncoordAxes >= 0)
assert ncoordAxes == countCoordAxes : "CoordAxes " + ncoordAxes + " != " + countCoordAxes;
if (nVertCooordAxes >= 0)
assert nVertCooordAxes == countVertCooordAxes : "VertAxes" + nVertCooordAxes + " != " + countVertCooordAxes;
gridDs.close();
}
public void utestProblem() throws Exception {
doOne( grib1dir,"eta.Y.Q.wmo", 25, -1, -1, -1);
doOne( grib2dir,"ndfd.wmo", 1, -1, -1, -1);
}
public void utestReadNcMLInputStream() throws Exception {
String ncmlLoc = "file:C:\\data\\work\\margolis\\test.ncml";
GridDataset fullDataset = GridDataset.open( ncmlLoc);
System.out.printf("full size= %d%n", fullDataset.getGrids().size());
// real ncml through a InputStream.
String ncml =
"<?xml version='1.0' encoding='UTF-8'?>\n" +
"<netcdf xmlns='http://www.unidata.ucar.edu/namespaces/netcdf/ncml-2.2'>\n" +
" <variable name='time'>\n" +
" <attribute name='ncmlAdded' value='timeAtt'/>\n" +
" </variable>\n" +
" <aggregation dimName='time' type='joinExisting'>\n" +
" <netcdf location='file:C:\\data\\work\\margolis\\f_0000000.nc'/>\n" +
" <netcdf location='file:C:\\data\\work\\margolis\\f_0032400.nc'/>\n" +
" </aggregation>\n" +
"</netcdf>";
NetcdfDataset aggregatedDataset = NcMLReader.readNcML( new ByteArrayInputStream(ncml.getBytes()), null );
GridDataset emptyDataset = new GridDataset( aggregatedDataset );
System.out.printf("empty= %s%n", emptyDataset.getGrids().size());
assert emptyDataset.getGrids().size() == fullDataset.getGrids().size();
}
public static void main( String arg[]) throws Exception {
// new TestReadandCount("dummy").doOne("C:/data/conventions/wrf/","wrf.nc", 33, 5, 7, 7); // missing TSLB
//new TestReadandCount("dummy").testRead(); // missing TSLB
//doOne(TestAll.testdataDir + "grid/grib/grib1/data/","AVN.wmo", 22, -1, -1, -1);
}
}