cdm/src/main/java/ucar/unidata/geoloc/projection/proj4/LambertConformalConicEllipse.java
private double computeTheta(double lon) {
    double dlon = LatLonPointImpl.lonNormal(lon - lon0deg);
    return n * Math.toRadians(dlon);
  }
private double computeTheta(double lon) {
    double dlon = LatLonPointImpl.lonNormal(lon - lon0deg);
    return n * Math.toRadians(dlon);
  }
/*
public ProjectionPoint latLonToProj(LatLonPoint latLon,
ProjectionPointImpl result) {
double toX, toY;
double fromLat = latLon.getLatitude();
double fromLon = latLon.getLongitude();

fromLat = Math.toRadians(fromLat);
double dlon = LatLonPointImpl.lonNormal(fromLon - lon0Degrees);
double theta = n * Math.toRadians(dlon);
double tn = Math.pow(Math.tan(PI_OVER_4 + fromLat / 2), n);
double r = earthRadiusTimesF / tn;
toX = r * Math.sin(theta);
toY = rho - r * Math.cos(theta);

result.setLocation(toX + falseEasting, toY + falseNorthing);
return result;
}

/*
public Point2D.Double project(double x, double y, Point2D.Double out) {
double rho;
if (Math.abs(Math.abs(y) - MapMath.HALFPI) < 1e-10)
rho = 0.0;
else
rho = c * (spherical ?
Math.pow(Math.tan(MapMath.QUARTERPI + .5 * y), -n) :
Math.pow(MapMath.tsfn(y, Math.sin(y), e), n));
out.x = scaleFactor * (rho * Math.sin(x *= n));
out.y = scaleFactor * (rho0 - rho * Math.cos(x));
return out;
}