private void showDates1D(CoordinateAxis1D axis1D, CalendarDateUnit cdu) {
if (!axis1D.isInterval()) {
for (double val : axis1D.getCoordValues()) {
if (Double.isNaN(val)) infoTA.appendLine(" N/A");
else infoTA.appendLine(" " + cdu.makeCalendarDate(val));
}
} else { // is interval
Formatter f = new Formatter();
double[] b1 = axis1D.getBound1();
double[] b2 = axis1D.getBound2();
for (int i = 0; i < b1.length; i++)
f.format(
" (%f, %f) == (%s, %s)%n",
b1[i], b2[i], cdu.makeCalendarDate((b1[i])), cdu.makeCalendarDate((b2[i])));
infoTA.appendLine(f.toString());
}
}
private void showValueDiffs(CoordinateAxis axis) {
if (!axis.isNumeric()) return;
try {
if (axis instanceof CoordinateAxis1D) {
CoordinateAxis1D axis1D = (CoordinateAxis1D) axis;
double[] mids = axis1D.getCoordValues();
double[] diffs = new double[mids.length];
for (int i = 0; i < mids.length - 1; i++) diffs[i] = mids[i + 1] - mids[i];
printArrays("midpoint differences", mids, diffs);
} else if (axis instanceof CoordinateAxis2D) {
CoordinateAxis2D axis2D = (CoordinateAxis2D) axis;
ArrayDouble.D2 mids = axis2D.getCoordValuesArray();
int[] shape = mids.getShape();
ArrayDouble.D2 diffx =
(ArrayDouble.D2) Array.factory(DataType.DOUBLE, new int[] {shape[0], shape[1] - 1});
for (int j = 0; j < shape[0]; j++) {
for (int i = 0; i < shape[1] - 1; i++) {
double diff = mids.get(j, i + 1) - mids.get(j, i);
diffx.set(j, i, diff);
}
}
infoTA.appendLine(NCdumpW.toString(diffx, "diff in x", null));
ArrayDouble.D2 diffy =
(ArrayDouble.D2) Array.factory(DataType.DOUBLE, new int[] {shape[0] - 1, shape[1]});
for (int j = 0; j < shape[0] - 1; j++) {
for (int i = 0; i < shape[1]; i++) {
double diff = mids.get(j + 1, i) - mids.get(j, i);
diffy.set(j, i, diff);
}
}
infoTA.appendLine("\n\n\n");
infoTA.appendLine(NCdumpW.toString(diffy, "diff in y", null));
}
} catch (Exception e1) {
e1.printStackTrace();
infoTA.appendLine(e1.getMessage());
}
}
public AxisBean(CoordinateAxis v) {
this.axis = v;
setName(v.getFullName());
setCoordVar(v.isCoordinateVariable());
setDescription(v.getDescription());
setUnits(v.getUnitsString());
// collect dimensions
StringBuilder lens = new StringBuilder();
StringBuilder names = new StringBuilder();
java.util.List dims = v.getDimensions();
for (int j = 0; j < dims.size(); j++) {
ucar.nc2.Dimension dim = (ucar.nc2.Dimension) dims.get(j);
if (j > 0) {
lens.append(",");
names.append(",");
}
String name = dim.isShared() ? dim.getShortName() : "anon";
names.append(name);
lens.append(dim.getLength());
}
setDims(names.toString());
setShape(lens.toString());
AxisType at = v.getAxisType();
if (at != null) setAxisType(at.toString());
String p = v.getPositive();
if (p != null) setPositive(p);
if (v instanceof CoordinateAxis1D) {
CoordinateAxis1D v1 = (CoordinateAxis1D) v;
if (v1.isRegular()) setRegular(Double.toString(v1.getIncrement()));
}
isLayer = (null != axis.findAttribute(_Coordinate.ZisLayer));
isInterval = axis.isInterval();
}
private void showValues(CoordinateAxis axis) {
try {
if (axis instanceof CoordinateAxis1D && axis.isNumeric()) {
CoordinateAxis1D axis1D = (CoordinateAxis1D) axis;
printArray("midpoints=", axis1D.getCoordValues());
if (!axis1D.isInterval()) {
printArray("edges=", axis1D.getCoordEdges());
} else {
printArray("bound1=", axis1D.getBound1());
printArray("bound2=", axis1D.getBound2());
Formatter f = new Formatter();
double[] mid = axis1D.getCoordValues();
double[] b1 = axis1D.getBound1();
double[] b2 = axis1D.getBound2();
for (int i = 0; i < b1.length; i++) {
f.format("%f (%f,%f) = %f%n", mid[i], b1[i], b2[i], b2[i] - b1[i]);
}
infoTA.appendLine(f.toString());
}
} else if (axis instanceof CoordinateAxis2D && axis.isNumeric()) {
infoTA.appendLine(NCdumpW.printVariableData(axis, null));
showValues2D((CoordinateAxis2D) axis);
} else {
infoTA.appendLine(NCdumpW.printVariableData(axis, null));
}
} catch (IOException e1) {
e1.printStackTrace();
infoTA.appendLine(e1.getMessage());
}
}