private void BindDensityLayer(ComputeShader compute, DensityProfileLayer layer)
{
compute.SetFloat(layer.Name + "_width", (float)(layer.Width / LengthUnitInMeters));
compute.SetFloat(layer.Name + "_exp_term", (float)layer.ExpTerm);
compute.SetFloat(layer.Name + "_exp_scale", (float)(layer.ExpScale * LengthUnitInMeters));
compute.SetFloat(layer.Name + "_linear_term", (float)(layer.LinearTerm * LengthUnitInMeters));
compute.SetFloat(layer.Name + "_constant_term", (float)layer.ConstantTerm);
}
public void Start()
{
MonoBehaviour.print("TURefs Start: " + GetHashCode());
TUFXScatteringResources.PrecomputeShader = precompute;
TUFXScatteringResources.ScatteringShader = scateringShader;
TUFXScatteringResources.WaterShader = waterShader;
sun = GameObject.Find("Sun");
//create model here
foreach (Model m in TUFXScatteringResources.Models)
{
m.Release();
}
TUFXScatteringResources.Models.Clear();
float kSunAngularRadius = 0.00935f / 2.0f;
float kBottomRadius = 6360000.0f;
float kLengthUnitInMeters = 100000.0f;
bool UseConstantSolarSpectrum = false;
bool UseOzone = true;
bool UseCombinedTextures = true;
bool UseHalfPrecision = false;
LUMINANCE UseLuminance = LUMINANCE.NONE;
// Values from "Reference Solar Spectral Irradiance: ASTM G-173", ETR column
// (see http://rredc.nrel.gov/solar/spectra/am1.5/ASTMG173/ASTMG173.html),
// summed and averaged in each bin (e.g. the value for 360nm is the average
// of the ASTM G-173 values for all wavelengths between 360 and 370nm).
// Values in W.m^-2.
int kLambdaMin = 360;
int kLambdaMax = 830;
double[] kSolarIrradiance = new double[]
{
1.11776, 1.14259, 1.01249, 1.14716, 1.72765, 1.73054, 1.6887, 1.61253,
1.91198, 2.03474, 2.02042, 2.02212, 1.93377, 1.95809, 1.91686, 1.8298,
1.8685, 1.8931, 1.85149, 1.8504, 1.8341, 1.8345, 1.8147, 1.78158, 1.7533,
1.6965, 1.68194, 1.64654, 1.6048, 1.52143, 1.55622, 1.5113, 1.474, 1.4482,
1.41018, 1.36775, 1.34188, 1.31429, 1.28303, 1.26758, 1.2367, 1.2082,
1.18737, 1.14683, 1.12362, 1.1058, 1.07124, 1.04992
};
// Values from http://www.iup.uni-bremen.de/gruppen/molspec/databases/
// referencespectra/o3spectra2011/index.html for 233K, summed and averaged in
// each bin (e.g. the value for 360nm is the average of the original values
// for all wavelengths between 360 and 370nm). Values in m^2.
double[] kOzoneCrossSection = new double[]
{
1.18e-27, 2.182e-28, 2.818e-28, 6.636e-28, 1.527e-27, 2.763e-27, 5.52e-27,
8.451e-27, 1.582e-26, 2.316e-26, 3.669e-26, 4.924e-26, 7.752e-26, 9.016e-26,
1.48e-25, 1.602e-25, 2.139e-25, 2.755e-25, 3.091e-25, 3.5e-25, 4.266e-25,
4.672e-25, 4.398e-25, 4.701e-25, 5.019e-25, 4.305e-25, 3.74e-25, 3.215e-25,
2.662e-25, 2.238e-25, 1.852e-25, 1.473e-25, 1.209e-25, 9.423e-26, 7.455e-26,
6.566e-26, 5.105e-26, 4.15e-26, 4.228e-26, 3.237e-26, 2.451e-26, 2.801e-26,
2.534e-26, 1.624e-26, 1.465e-26, 2.078e-26, 1.383e-26, 7.105e-27
};
// From https://en.wikipedia.org/wiki/Dobson_unit, in molecules.m^-2.
double kDobsonUnit = 2.687e20;
// Maximum number density of ozone molecules, in m^-3 (computed so at to get
// 300 Dobson units of ozone - for this we divide 300 DU by the integral of
// the ozone density profile defined below, which is equal to 15km).
double kMaxOzoneNumberDensity = 300.0 * kDobsonUnit / 15000.0;
// Wavelength independent solar irradiance "spectrum" (not physically
// realistic, but was used in the original implementation).
double kConstantSolarIrradiance = 1.5;
//double kTopRadius = 6420000.0;
double kTopRadius = 6560000.0;
double kRayleigh = 1.24062e-6;
double kRayleighScaleHeight = 8000.0;
double kMieScaleHeight = 1200.0;
double kMieAngstromAlpha = 0.0;
double kMieAngstromBeta = 5.328e-3;
double kMieSingleScatteringAlbedo = 0.9;
double kMiePhaseFunctionG = 0.8;
double kGroundAlbedo = 0.1;
double max_sun_zenith_angle = (UseHalfPrecision ? 102.0 : 120.0) / 180.0 * Mathf.PI;
DensityProfileLayer rayleigh_layer = new DensityProfileLayer("rayleigh", 0.0, 1.0, -1.0 / kRayleighScaleHeight, 0.0, 0.0);
DensityProfileLayer mie_layer = new DensityProfileLayer("mie", 0.0, 1.0, -1.0 / kMieScaleHeight, 0.0, 0.0);
// Density profile increasing linearly from 0 to 1 between 10 and 25km, and
// decreasing linearly from 1 to 0 between 25 and 40km. This is an approximate
// profile from http://www.kln.ac.lk/science/Chemistry/Teaching_Resources/Documents/Introduction%20to%20atmospheric%20chemistry.pdf (page 10).
List <DensityProfileLayer> ozone_density = new List <DensityProfileLayer>();
ozone_density.Add(new DensityProfileLayer("absorption0", 25000.0, 0.0, 0.0, 1.0 / 15000.0, -2.0 / 3.0));
ozone_density.Add(new DensityProfileLayer("absorption1", 0.0, 0.0, 0.0, -1.0 / 15000.0, 8.0 / 3.0));
List <double> wavelengths = new List <double>();
List <double> solar_irradiance = new List <double>();
List <double> rayleigh_scattering = new List <double>();
List <double> mie_scattering = new List <double>();
List <double> mie_extinction = new List <double>();
List <double> absorption_extinction = new List <double>();
List <double> ground_albedo = new List <double>();
for (int l = kLambdaMin; l <= kLambdaMax; l += 10)
{
double lambda = l * 1e-3; // micro-meters
double mie = kMieAngstromBeta / kMieScaleHeight * Math.Pow(lambda, -kMieAngstromAlpha);
wavelengths.Add(l);
if (UseConstantSolarSpectrum)
{
solar_irradiance.Add(kConstantSolarIrradiance);
}
else
{
solar_irradiance.Add(kSolarIrradiance[(l - kLambdaMin) / 10]);
}
rayleigh_scattering.Add(kRayleigh * Math.Pow(lambda, -4));
mie_scattering.Add(mie * kMieSingleScatteringAlbedo);
mie_extinction.Add(mie);
absorption_extinction.Add(UseOzone ? kMaxOzoneNumberDensity * kOzoneCrossSection[(l - kLambdaMin) / 10] : 0.0);
ground_albedo.Add(kGroundAlbedo);
}
Model model = new Model();
model.HalfPrecision = UseHalfPrecision;
model.CombineScatteringTextures = UseCombinedTextures;
model.UseLuminance = UseLuminance;
model.Wavelengths = wavelengths;
model.SolarIrradiance = solar_irradiance;
model.SunAngularRadius = kSunAngularRadius;
model.BottomRadius = kBottomRadius;
model.TopRadius = kTopRadius;
model.RayleighDensity = rayleigh_layer;
model.RayleighScattering = rayleigh_scattering;
model.MieDensity = mie_layer;
model.MieScattering = mie_scattering;
model.MieExtinction = mie_extinction;
model.MiePhaseFunctionG = kMiePhaseFunctionG;
model.AbsorptionDensity = ozone_density;
model.AbsorptionExtinction = absorption_extinction;
model.GroundAlbedo = ground_albedo;
model.MaxSunZenithAngle = max_sun_zenith_angle;
model.LengthUnitInMeters = kLengthUnitInMeters;
model.SunDirection = -sun.transform.forward;
int numScatteringOrders = 6;
//model.Init(TUFXScatteringResources.PrecomputeShader, numScatteringOrders);
//TUFXScatteringResources.Models.Add(model);
kRayleighScaleHeight *= 0.8;
kMieScaleHeight *= 0.8;
rayleigh_layer = new DensityProfileLayer("rayleigh", 0.0, 1.0, -1.0 / kRayleighScaleHeight, 0.0, 0.0);
mie_layer = new DensityProfileLayer("mie", 0.0, 1.0, -1.0 / kMieScaleHeight, 0.0, 0.0);
model = new Model();
model.HalfPrecision = UseHalfPrecision;
model.CombineScatteringTextures = UseCombinedTextures;
model.UseLuminance = UseLuminance;
model.Wavelengths = wavelengths;
model.SolarIrradiance = solar_irradiance;
model.SunAngularRadius = kSunAngularRadius;
model.BottomRadius = 600000;
model.TopRadius = 670000;
model.RayleighDensity = rayleigh_layer;
model.RayleighScattering = rayleigh_scattering;
model.MieDensity = mie_layer;
model.MieScattering = mie_scattering;
model.MieExtinction = mie_extinction;
model.MiePhaseFunctionG = kMiePhaseFunctionG;
model.AbsorptionDensity = ozone_density;
model.AbsorptionExtinction = absorption_extinction;
model.GroundAlbedo = ground_albedo;
model.MaxSunZenithAngle = max_sun_zenith_angle;
model.LengthUnitInMeters = 100;
model.SunDirection = -Vector3.forward;
model.PlanetCenter = Vector3.forward * 7005;
numScatteringOrders = 6;
model.Init(TUFXScatteringResources.PrecomputeShader, numScatteringOrders);
TUFXScatteringResources.Models.Add(model);
}
