private void showValues2D(CoordinateAxis2D axis2D) {
Formatter f = new Formatter();
if (axis2D.isInterval()) {
ArrayDouble.D2 coords = axis2D.getCoordValuesArray();
ArrayDouble.D3 bounds = axis2D.getCoordBoundsArray();
if (bounds == null) {
infoTA.appendLine("No bounds for interval " + axis2D.getFullName());
return;
}
IndexIterator coordIter = coords.getIndexIterator();
IndexIterator boundsIter = bounds.getIndexIterator();
while (coordIter.hasNext()) {
double coordValue = coordIter.getDoubleNext();
if (!boundsIter.hasNext()) break;
double bounds1 = boundsIter.getDoubleNext();
if (!boundsIter.hasNext()) break;
double bounds2 = boundsIter.getDoubleNext();
f.format("%f (%f,%f) = %f%n", coordValue, bounds1, bounds2, bounds2 - bounds1);
}
} else {
ArrayDouble.D2 coords = axis2D.getCoordValuesArray();
IndexIterator coordIter = coords.getIndexIterator();
while (coordIter.hasNext()) {
double coordValue = coordIter.getDoubleNext();
f.format("%f%n", coordValue);
}
}
infoTA.appendLine(f.toString());
}
private void printArrays(String title, double vals[], double vals2[]) {
Formatter sbuff = new Formatter();
sbuff.format(" %s%n", title);
for (int i = 0; i < vals.length; i++) {
sbuff.format(" %3d: %10.2f %10.2f%n", i, vals[i], vals2[i]);
}
sbuff.format("%n");
infoTA.appendLine(sbuff.toString());
}
private void printArray(String title, double vals[]) {
Formatter sbuff = new Formatter();
sbuff.format(" %s=", title);
for (double val : vals) {
sbuff.format(" %f", val);
}
sbuff.format("%n");
infoTA.appendLine(sbuff.toString());
}
private String showMissing(Variable v) {
if (!(v instanceof VariableDS)) return "";
VariableDS ve = (VariableDS) v;
Formatter buff = new Formatter();
buff.format("%s:", v.getFullName());
EnumSet<NetcdfDataset.Enhance> enhanceMode = ve.getEnhanceMode();
buff.format("enhanceMode= %s%n", enhanceMode);
ve.showScaleMissingProxy(buff);
return buff.toString();
}
public void summaryInfo(Formatter f) {
if (ds == null) return;
f.format("%s%n", ds.getLocation());
int ngrids = 0;
for (Object varo : varTable.getBeans()) {
VariableBean varBean = (VariableBean) varo;
if (varBean.getDataType().trim().equalsIgnoreCase("grid")) ngrids++;
}
int ncoordSys = csTable.getBeans().size();
int ncoords = axisTable.getBeans().size();
f.format(" ngrids=%d, ncoords=%d, ncoordSys=%d%n", ngrids, ncoords, ncoordSys);
}
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 showValuesAsDates(CoordinateAxis axis) {
String units = axis.getUnitsString();
String cal = getCalendarAttribute(axis);
CalendarDateUnit cdu = CalendarDateUnit.of(cal, units);
try {
infoTA.appendLine(units);
infoTA.appendLine(NCdumpW.printVariableData(axis, null));
if (axis.getDataType().isNumeric()) {
if (axis instanceof CoordinateAxis2D) {
showDates2D((CoordinateAxis2D) axis, cdu);
} else if (axis instanceof CoordinateAxis1D) { // 1D
showDates1D((CoordinateAxis1D) axis, cdu);
} else { // > 2D
Array data = axis.read();
IndexIterator ii = data.getIndexIterator();
while (ii.hasNext()) {
double val = ii.getDoubleNext();
infoTA.appendLine(makeCalendarDateStringOrMissing(cdu, val));
}
}
} else { // must be iso dates
Array data = axis.read();
Formatter f = new Formatter();
if (data instanceof ArrayChar) {
ArrayChar dataS = (ArrayChar) data;
ArrayChar.StringIterator iter = dataS.getStringIterator();
while (iter.hasNext()) f.format(" %s%n", iter.next());
infoTA.appendLine(f.toString());
} else if (data instanceof ArrayObject) {
IndexIterator iter = data.getIndexIterator();
while (iter.hasNext()) f.format(" %s%n", iter.next());
infoTA.appendLine(f.toString());
}
}
} catch (Exception ex) {
ex.printStackTrace();
infoTA.appendLine(ex.getMessage());
}
}
private void writeHeader(ProfileFeature profile, PointFeature obs) throws IOException {
Formatter coordNames = new Formatter().format("%s %s %s", profileTimeName, latName, lonName);
List<VariableSimpleIF> coords = new ArrayList<>();
if (useAlt) {
coords.add(
VariableSimpleImpl.makeScalar(
altitudeCoordinateName, "obs altitude", altUnits, DataType.DOUBLE)
.add(new Attribute(CF.STANDARD_NAME, "altitude"))
.add(
new Attribute(
CF.POSITIVE,
CF1Convention.getZisPositive(altitudeCoordinateName, altUnits))));
coordNames.format(" %s", altitudeCoordinateName);
}
super.writeHeader(
coords, profile.getFeatureData(), obs.getFeatureData(), coordNames.toString());
}
@Override
public void getDetailInfo(Formatter f) {
super.getDetailInfo(f);
if (runMatcher != null) f.format(" runMatcher=%s%n", runMatcher);
/* if (forecastMatcher != null)
f.format(" forecastMatcher=%s%n", forecastMatcher);
if (offsetMatcher != null)
f.format(" offsetMatcher=%s%n", offsetMatcher); */
}
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());
}
}
private void showDates2D(CoordinateAxis2D axis2D, CalendarDateUnit cdu) {
Formatter f = new Formatter();
if (axis2D.isInterval()) {
ArrayDouble.D2 coords = axis2D.getCoordValuesArray();
ArrayDouble.D3 bounds = axis2D.getCoordBoundsArray();
if (bounds == null) {
infoTA.appendLine("No bounds for interval " + axis2D.getFullName());
return;
}
IndexIterator coordIter = coords.getIndexIterator();
IndexIterator boundsIter = bounds.getIndexIterator();
while (coordIter.hasNext()) {
double coordValue = coordIter.getDoubleNext();
if (!boundsIter.hasNext()) break;
double bounds1 = boundsIter.getDoubleNext();
if (!boundsIter.hasNext()) break;
double bounds2 = boundsIter.getDoubleNext();
f.format(
"%s (%s,%s)%n",
makeCalendarDateStringOrMissing(cdu, coordValue),
makeCalendarDateStringOrMissing(cdu, bounds1),
makeCalendarDateStringOrMissing(cdu, bounds2));
}
} else {
ArrayDouble.D2 coords = axis2D.getCoordValuesArray();
IndexIterator coordIter = coords.getIndexIterator();
while (coordIter.hasNext()) {
double coordValue = coordIter.getDoubleNext();
f.format("%s%n", makeCalendarDateStringOrMissing(cdu, coordValue));
}
}
infoTA.appendLine(f.toString());
}
private String tryGrid(VariableEnhanced v) {
Formatter buff = new Formatter();
buff.format("%s:", v.getFullName());
List<CoordinateSystem> csList = v.getCoordinateSystems();
if (csList.size() == 0) buff.format(" No Coord System found");
else {
for (CoordinateSystem cs : csList) {
buff.format("%s:", cs.getName());
if (GridCoordSys.isGridCoordSys(buff, cs, v)) {
buff.format("GRID OK%n");
} else {
buff.format(" NOT GRID");
}
}
}
return buff.toString();
}