htrdr

commit 08d94d80ee0ac36195c553e622dea997a00c8a45
parent b07d94977d941f26e024ac87881df0d84ad548bc
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu, 26 Jul 2018 10:02:50 +0200

Add the HG g parameter getter on the sky.

Diffstat:
M
src/htrdr_sky.c
 | 
15
+++++++++++++++
M
src/htrdr_sky.h
 | 
5
+++++

2 files changed, 20 insertions(+), 0 deletions(-)
diff --git a/src/htrdr_sky.c b/src/htrdr_sky.c
@@ -96,6 +96,9 @@ struct htrdr_sky {
   double Ca_avg_sw;
   double Cs_avg_sw;
 
+  /* Average asymmetry parameter in Short Wave, in [380, 780] nanometers */
+  double g_avg_sw;
+
   ref_T ref;
   struct htrdr* htrdr;
 };
@@ -501,6 +504,9 @@ htrdr_sky_create
     (sky->htmie, band_sw, HTMIE_FILTER_LINEAR);
   sky->Cs_avg_sw = htmie_compute_xsection_scattering_average
     (sky->htmie, band_sw, HTMIE_FILTER_LINEAR);
+  sky->g_avg_sw = htmie_compute_asymmetry_parameter_average
+    (sky->htmie, band_sw, HTMIE_FILTER_LINEAR);
+  ASSERT(sky->Ca_avg_sw > 0 && sky->Cs_avg_sw > 0 && sky->g_avg_sw > 0);
 
   res = setup_clouds(sky);
   if(res != RES_OK) goto error;
@@ -531,6 +537,15 @@ htrdr_sky_ref_put(struct htrdr_sky* sky)
 }
 
 double
+htrdr_sky_fetch_particle_phase_function_asymmetry_parameter
+  (const struct htrdr_sky* sky, const double wavelength)
+{
+  ASSERT(sky && wavelength > 0);
+  (void)wavelength;
+  return sky->g_avg_sw;
+}
+
+double
 htrdr_sky_fetch_raw_property
   (const struct htrdr_sky* sky,
    const enum htrdr_sky_property prop,
diff --git a/src/htrdr_sky.h b/src/htrdr_sky.h
@@ -63,6 +63,11 @@ htrdr_sky_ref_put
   (struct htrdr_sky* sky);
 
 extern LOCAL_SYM double
+htrdr_sky_fetch_particle_phase_function_asymmetry_parameter
+  (const struct htrdr_sky* sky,
+   const double wavelength);
+
+extern LOCAL_SYM double
 htrdr_sky_fetch_raw_property
   (const struct htrdr_sky* sky,
    const enum htrdr_sky_property prop,