Read data variables and attributes from the file
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;
final float SAMPLE_PRESSURE = 900f;
final float SAMPLE_TEMP = 9.0f;
final float START_LAT = 25.0f;
final float START_LON = -125.0f;
float[][] presIn = new float[NLAT][NLON];
float[][] tempIn = new float[NLAT][NLON];
float[] latsIn = new float[NLAT];
float[] lonsIn = new float[NLON];
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;
}
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);
}
for (int lat = 0; lat < NLAT; lat++)
if (latsIn[lat] != START_LAT + 5. * lat)
System.err.println("ERROR reading variable latitude");
for (int lon = 0; lon < NLON; lon++)
if (lonsIn[lon] != START_LON + 5. * lon)
System.err.println("ERROR reading variable longitude");
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);
}
}
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");
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!");
}
}