+ SG_LOG( SG_GENERAL, SG_DEBUG, "Wind : " << dir_from << "@" << speed);
+ }
+
+ // Now handle the visibility. We convert both visibility values
+ // to X-values, then interpolate from there, then back to real values.
+ // The length_scale is fixed to 1000m, so the visibility changes by
+ // by MaxVisChangePercentSec or 1000m X MaxVisChangePercentSec,
+ // whichever is more.
+ double currentvis =
+ fgGetDouble("/environment/config/boundary/entry/visibility-m");
+ double metarvis = fgGetDouble("/environment/metar/min-visibility-m");
+ double currentxval = log(1000.0 + currentvis);
+ double metarxval = log(1000.0 + metarvis);
+
+ currentxval = interpolate_val(currentxval, metarxval, MaxVisChangePercentSec);
+
+ // Now convert back from an X-value to a straightforward visibility.
+ vis = exp(currentxval) - 1000.0;
+
+ pressure = interpolate_prop(
+ "/environment/config/boundary/entry/pressure-sea-level-inhg",
+ "/environment/metar/pressure-inhg",
+ MaxPressureChangeInHgSec);
+
+ dir_to = fgGetDouble("/environment/metar/base-wind-range-to");
+ gust = fgGetDouble("/environment/metar/gust-wind-speed-kt");
+ temp = fgGetDouble("/environment/metar/temperature-degc");
+ dewpoint = fgGetDouble("/environment/metar/dewpoint-degc");
+
+ // Set the cloud layers by interpolating over the METAR versions.
+ SGPropertyNode * clouds = fgGetNode("/environment/metar/clouds");
+
+ vector<SGPropertyNode_ptr> layers = clouds->getChildren("layer");
+ vector<SGPropertyNode_ptr>::const_iterator layer;
+ vector<SGPropertyNode_ptr>::const_iterator layers_end = layers.end();
+
+ const char *cl = "/environment/clouds/layer[%i]";
+ double aircraft_alt = fgGetDouble("/position/altitude-ft");
+ char s[128];
+ int i;
+
+ for (i = 0, layer = layers.begin(); layer != layers_end; ++layer, i++) {
+ double currentval;
+ double requiredval;
+
+ // In the case of clouds, we want to avoid writing if nothing has
+ // changed, as these properties are tied to the renderer and will
+ // cause the clouds to be updated, reseting the texture locations.
+
+ // We don't interpolate the coverage values as no-matter how we
+ // do it, it will be quite a sudden change of texture. Better to
+ // have a single change than four or five.
+ snprintf(s, 128, cl, i);
+ strncat(s, "/coverage", 128);
+ const char* coverage = (*layer)->getStringValue("coverage", "clear");
+ if (strncmp(fgGetString(s), coverage, 128) != 0)
+ fgSetString(s, coverage);
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/elevation-ft", 128);
+ double current_alt = fgGetDouble(s);
+ double required_alt = (*layer)->getDoubleValue("elevation-ft");
+
+ if (current_alt < -9000 || required_alt < -9000 ||
+ fabs(aircraft_alt - required_alt) > MaxCloudInterpolationHeightFt ||
+ fabs(current_alt - required_alt) > MaxCloudInterpolationDeltaFt) {
+ // We don't interpolate any layers that are
+ // - too far above us to be visible
+ // - too far below us to be visible
+ // - with too large a difference to make interpolation sensible
+ // - to or from -9999 (used as a placeholder)
+ // - any values that are too high above us,
+ snprintf(s, 128, cl, i);
+ strncat(s, "/elevation-ft", 128);
+ if (current_alt != required_alt)
+ fgSetDouble(s, required_alt);
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/thickness-ft", 128);
+ if (fgGetDouble(s) != (*layer)->getDoubleValue("thickness-ft"))
+ fgSetDouble(s, (*layer)->getDoubleValue("thickness-ft"));
+
+ } else {
+ // Interpolate the other values in the usual way
+ if (current_alt != required_alt) {
+ current_alt = interpolate_val(current_alt,
+ required_alt,
+ MaxCloudAltitudeChangeFtSec);
+ fgSetDouble(s, current_alt);
+ }
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/thickness-ft", 128);
+ currentval = fgGetDouble(s);
+ requiredval = (*layer)->getDoubleValue("thickness-ft");
+
+ if (currentval != requiredval) {
+ currentval = interpolate_val(currentval,
+ requiredval,
+ MaxCloudThicknessChangeFtSec);
+ fgSetDouble(s, currentval);
+ }
+ }
+ }
+
+ } else {
+ // We haven't already loaded a METAR, so apply it immediately.
+ dir_from = fgGetDouble("/environment/metar/base-wind-range-from");
+ dir_to = fgGetDouble("/environment/metar/base-wind-range-to");
+ speed = fgGetDouble("/environment/metar/base-wind-speed-kt");
+ gust = fgGetDouble("/environment/metar/gust-wind-speed-kt");
+ vis = fgGetDouble("/environment/metar/min-visibility-m");
+ pressure = fgGetDouble("/environment/metar/pressure-inhg");
+ temp = fgGetDouble("/environment/metar/temperature-degc");
+ dewpoint = fgGetDouble("/environment/metar/dewpoint-degc");
+
+ // Set the cloud layers by copying over the METAR versions.
+ SGPropertyNode * clouds = fgGetNode("/environment/metar/clouds", true);
+
+ vector<SGPropertyNode_ptr> layers = clouds->getChildren("layer");
+ vector<SGPropertyNode_ptr>::const_iterator layer;
+ vector<SGPropertyNode_ptr>::const_iterator layers_end = layers.end();
+
+ const char *cl = "/environment/clouds/layer[%i]";
+ char s[128];
+ int i;
+
+ for (i = 0, layer = layers.begin(); layer != layers_end; ++layer, i++) {
+ snprintf(s, 128, cl, i);
+ strncat(s, "/coverage", 128);
+ fgSetString(s, (*layer)->getStringValue("coverage", "clear"));
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/elevation-ft", 128);
+ fgSetDouble(s, (*layer)->getDoubleValue("elevation-ft"));
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/thickness-ft", 128);
+ fgSetDouble(s, (*layer)->getDoubleValue("thickness-ft"));
+
+ snprintf(s, 128, cl, i);
+ strncat(s, "/span-m", 128);
+ fgSetDouble(s, 40000.0);
+ }
+
+ // Force an update of the 3D clouds
+ fgSetDouble("/environment/rebuild-layers", 1.0);
+ }
+
+ fgSetupWind(dir_from, dir_to, speed, gust);
+ fgDefaultWeatherValue("visibility-m", vis);
+ set_temp_at_altitude(temp, station_elevation_ft);
+ set_dewpoint_at_altitude(dewpoint, station_elevation_ft);
+ fgDefaultWeatherValue("pressure-sea-level-inhg", pressure);
+
+ // We've now successfully loaded a METAR into the configuration
+ metar_loaded = true;
+}
+
+double FGMetarEnvironmentCtrl::interpolate_prop(const char * currentname,
+ const char * requiredname,
+ double dt)
+{
+ double currentval = fgGetDouble(currentname);
+ double requiredval = fgGetDouble(requiredname);
+ return interpolate_val(currentval, requiredval, dt);
+}
+
+double FGMetarEnvironmentCtrl::interpolate_val(double currentval,
+ double requiredval,
+ double dt)
+{
+ double dval = EnvironmentUpdatePeriodSec * dt;
+
+ if (fabs(currentval - requiredval) < dval) return requiredval;
+ if (currentval < requiredval) return (currentval + dval);
+ if (currentval > requiredval) return (currentval - dval);
+ return requiredval;