/*
* 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 examples;
import ucar.nc2.NetcdfFile;
import ucar.nc2.Variable;
import ucar.ma2.ArrayFloat;
import java.io.IOException;
/**
* User: yuanho
* Date: Nov 16, 2006
* Time: 11:59:38 AM
* To change this template use File | Settings | File Templates.
*/
public class Sfc_pres_temp_rd {
public static void main(String args[]) throws Exception
{
final int NLAT = 6;
final int NLON = 12;
// These are used to calculate the values we expect to find.
final float SAMPLE_PRESSURE = 900f;
final float SAMPLE_TEMP = 9.0f;
final float START_LAT = 25.0f;
final float START_LON = -125.0f;
// These will hold our pressure and temperature data.
float [][] presIn = new float[NLAT][NLON];
float [][] tempIn = new float[NLAT][NLON];
// These will hold our latitudes and longitudes.
float [] latsIn = new float[NLAT];
float [] lonsIn = new float[NLON];
// Open the file and check to make sure it's valid.
String filename = "sfc_pres_temp.nc";
NetcdfFile dataFile = null;
try {
dataFile = NetcdfFile.open(filename, null);
Variable latVar = dataFile.findVariable("latitude");
if (latVar == null) {
System.out.println("Cant find Variable latitude");
return;
}
Variable lonVar= dataFile.findVariable("longitude");
if (lonVar == null) {
System.out.println("Cant find Variable longitude");
return;
}
Variable presVar= dataFile.findVariable("pressure");
if (presVar == null) {
System.out.println("Cant find Variable pressure");
return;
}
Variable tempVar= dataFile.findVariable("temperature");
if (tempVar == null) {
System.out.println("Cant find Variable temperature");
return;
}
if(latVar.getDimensions().size() != 1) {
System.out.println(" fail to get the dimensions of variable latitude");
return;
}
if(presVar.getDimensions().size() != 2) {
System.out.println(" fail to get the dimensions of variable pressure");
return;
}
// Read the latitude and longitude coordinate variables into arrays
// latsIn and lonsIn.
ArrayFloat.D1 latArray;
ArrayFloat.D1 lonArray;
latArray = (ArrayFloat.D1)latVar.read();
lonArray = (ArrayFloat.D1)lonVar.read();
int[] shape = latArray.getShape();
for (int i=0; i<shape[0]; i++) {
latsIn[i] = latArray.get(i);
}
shape = lonArray.getShape();
for (int j=0; j<shape[0]; j++) {
lonsIn[j] = lonArray.get(j);
}
// Check the coordinate variable data.
for(int lat = 0; lat < NLAT; lat++)
if (latsIn[lat] != START_LAT + 5. * lat)
System.err.println("ERROR reading variable latitude");
// Check longitude values.
for (int lon = 0; lon < NLON; lon++)
if (lonsIn[lon] != START_LON + 5. * lon)
System.err.println("ERROR reading variable longitude");
// Read the data. Since we know the contents of the file we know
// that the data arrays in this program are the correct size to
// hold all the data.
ArrayFloat.D2 presArray, tempArray;
presArray = (ArrayFloat.D2)presVar.read();
tempArray = (ArrayFloat.D2)tempVar.read();
int [] shape1 = presArray.getShape();
for (int i=0; i<shape1[0]; i++) {
for (int j=0; j<shape1[1]; j++) {
presIn[i][j] = presArray.get(i,j);
tempIn[i][j] = tempArray.get(i,j);
}
}
// Check the data.
for (int lat = 0; lat < NLAT; lat++)
for (int lon = 0; lon < NLON; lon++)
if (presIn[lat][lon] != SAMPLE_PRESSURE + (lon * NLAT + lat)
|| tempIn[lat][lon] != SAMPLE_TEMP + .25 * (lon * NLAT + lat))
System.err.println("ERROR reading variable pressure or temperature");
// Each of the netCDF variables has a "units" attribute. Let's read
// them and check them.
if(!latVar.findAttributeIgnoreCase("units").getStringValue().equalsIgnoreCase("degrees_north"))
System.err.println("ERROR reading variable latitude units");
if(!lonVar.findAttributeIgnoreCase("units").getStringValue().equalsIgnoreCase("degrees_east"))
System.err.println("ERROR reading variable longitude units");
if(!presVar.findAttributeIgnoreCase("units").getStringValue().equalsIgnoreCase("hPa"))
System.err.println("ERROR reading variable pressure units");
if(!tempVar.findAttributeIgnoreCase("units").getStringValue().equalsIgnoreCase("celsius"))
System.err.println("ERROR reading variable temperature units");
} catch (java.io.IOException e) {
System.out.println(" fail = "+e);
e.printStackTrace();
} finally {
if (dataFile != null)
try {
dataFile.close();
} catch (IOException ioe) {
ioe.printStackTrace();
}
}
System.out.println("*** SUCCESS reading example file sfc_pres_temp.nc!");
}
}