htrdrPy 1.1.6__tar.gz → 1.1.7__tar.gz

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/.gitlab-ci.yml

@@ -27,9 +27,9 @@ pages:
     - python3 -m venv venv
     - source venv/bin/activate
     - pip install -e .
-    - python3 -m pip install sphinx furo myst-parser sphinx_rtd_theme numpy matplotlib scipy tqdm netCDF4 vtk pyvista
+    - python3 -m pip install sphinx furo myst-parser sphinx_rtd_theme numpy matplotlib scipy tqdm netCDF4 vtk pyvista xarray
   script:
-    - d docs && python -c "import htrdPy; print('Module imported successfully')" && cd ..
+    - cd docs && python -c "import htrdrPy; print('Module imported successfully')" && cd ..
     - cd docs && make html
   after_script:
     - mv docs/build/html/ ./public/

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: htrdrPy
-Version: 1.1.6
+Version: 1.1.7
 Summary: htrdrPy is a library that aims to wrap the radiative transfer code htrdr-planets, aiming to make the use of htrdr-planets more user-friendly.
 Project-URL: Homepage, https://gitlab.obspm.fr/aarfaux/htrdr-wrapper.git
 Project-URL: Documentation, https://htrdr-wrapper-ebbb17.pages.obspm.fr

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/_modules/htrdrPy/data.html

@@ -1145,6 +1145,7 @@
 <span class="sd">        ----------</span>
 <span class="sd">        data : dict</span>
 <span class="sd">            Dictionnary with the following items: </span>
+
 <span class="sd">            - &quot;nLevel&quot; : int</span>
 <span class="sd">                Number of levels. Required only if sameGrid is ``False``.</span>
 <span class="sd">            - &quot;nCoeff&quot;: int</span>
@@ -1157,13 +1158,15 @@
 <span class="sd">                Number of angles in the phase functions (not required if</span>
 <span class="sd">                using the builtin Henyey-Greenstein phase function).</span>
 <span class="sd">            - &quot;altitude (m)&quot;: ``numpy.ndarray``</span>
-<span class="sd">                Array of altitudes (shape=(nLevel, nLat, nLon), [m]).</span>
+<span class="sd">                Array of altitudes (shape=(nLevel, [nLat, nLon]), [m]).</span>
 <span class="sd">            - &quot;latitude (°)&quot; : ``numpy.ndarray``, optional</span>
 <span class="sd">                List of latitudes (shape=(nLat), [°]). Required only if ``dim``</span>
 <span class="sd">                is 2 or 3.</span>
 <span class="sd">            - &quot;longitude (°)&quot; : ``numpy.ndarray``, optional</span>
 <span class="sd">                List of longitudes (shape=(nLon), [°]). Required only if ``dim``</span>
 <span class="sd">                is 3.</span>
+<span class="sd">            - &quot;temperature (K)&quot;: ``numpy.ndarray``</span>
+<span class="sd">                Array of temperatures (shape=(nLevel, [nLat, nLon]), [K]).</span>
 <span class="sd">            - &quot;nWavelength&quot;:  int</span>
 <span class="sd">                Number of wavelengths.</span>
 <span class="sd">            - &quot;wavelength&quot;: ``numpy.ndarray``</span>
@@ -1176,6 +1179,9 @@
 <span class="sd">            - &quot;band up&quot;:  ``numpy.ndarray``</span>
 <span class="sd">                Array of wavelength (shape = (nWavelength), [m]). The</span>
 <span class="sd">                values corresponds to the lower boundary of the band.</span>
+<span class="sd">            - &quot;weights&quot;: ``numpy.ndarray``</span>
+<span class="sd">                The weigths of the nCoeff quadrature points</span>
+<span class="sd">                (shape=(nCoeff)).</span>
 <span class="sd">            - &quot;scattering (m-1)&quot;: ``numpy.ndarray``</span>
 <span class="sd">                 Array of scattering coefficients</span>
 <span class="sd">                 (shape=(nWavelength, nLevel, [nLat, nLon]), [m-1]).</span>
@@ -1761,6 +1767,7 @@
 <span class="sd">            as a plan-parallel.</span>
 <span class="sd">        data : dict</span>
 <span class="sd">            Dictionnary with the following items: </span>
+
 <span class="sd">            - &quot;nAngle&quot;: int, optional</span>
 <span class="sd">                Number of angles in the phase functions (not required if</span>
 <span class="sd">                using the builtin Henyey-Greenstein phase function).</span>

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/_modules/htrdrPy/geometry.html

@@ -409,6 +409,13 @@
 <span class="sd">                    cells : ``numpy.ndarray``</span>
 <span class="sd">                        Array of the cells to be extracted (shape=(nCell,3))</span>
 <span class="sd">                        with the altitude, latitude and longitude indices.</span>
+
+<span class="sd">        Warning</span>
+<span class="sd">        -------</span>
+<span class="sd">        In the table of coordiantes that must be provided, the latitude should</span>
+<span class="sd">        be given by decreasing latitudes (North pole corresponding to the first</span>
+<span class="sd">        latitude index and South pole to the last), according to LMDZ PCM</span>
+<span class="sd">        convention.</span>
 <span class="sd">        &#39;&#39;&#39;</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">volrad</span><span class="p">[</span><span class="s2">&quot;sampling&quot;</span><span class="p">]</span> <span class="o">=</span> <span class="n">sampling</span>
         <span class="bp">self</span><span class="o">.</span><span class="n">GCMmesh</span> <span class="o">=</span> <span class="kc">False</span>
@@ -799,13 +806,19 @@
 <span class="w">        </span><span class="sd">&#39;&#39;&#39;</span>
 <span class="sd">        Make the mesh from the cell coordinates</span>
 <span class="sd">        # Input</span>
-<span class="sd">        - cellCoord (numpy 4-D array (shape=(nAltitudes,nLatitudes,nLongitudes,3), float [m])): coordinates of the cells centers in the atmosphere</span>
+<span class="sd">        - cellCoord (numpy 4-D array</span>
+<span class="sd">          (shape=(nAltitudes,nLatitudes,nLongitudes,3), float [m])): coordinates</span>
+<span class="sd">          of the cells centers or levels in the atmosphere</span>
 <span class="sd">        - radius (float [m]): radius of the planet</span>
 <span class="sd">        - poles (optional, default=False, bool): whether or not the first and</span>
 <span class="sd">          last latitudes correspond to the poles in the given coordinates</span>
 <span class="sd">        - onLevels (optional, default=False, bool): whether or not the altitudes</span>
 <span class="sd">          provided in the coordinates are given on the levels or in the center</span>
 <span class="sd">          of the cells</span>
+
+<span class="sd">        # Warning</span>
+<span class="sd">        The latitudes are provided from North pole to south pole according to</span>
+<span class="sd">        LMDZ PCM convention !</span>
 <span class="sd">        &#39;&#39;&#39;</span>
 
         <span class="bp">self</span><span class="o">.</span><span class="n">latitudes</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">unique</span><span class="p">(</span><span class="n">cellCoord</span><span class="p">[:,:,:,</span><span class="mi">1</span><span class="p">])</span>

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/_modules/htrdrPy/postprocess.html

@@ -70,6 +70,7 @@
              
   <h1>Source code for htrdrPy.postprocess</h1><div class="highlight"><pre>
 <span></span><span class="kn">import</span><span class="w"> </span><span class="nn">numpy</span><span class="w"> </span><span class="k">as</span><span class="w"> </span><span class="nn">np</span>
+<span class="kn">import</span><span class="w"> </span><span class="nn">xarray</span><span class="w"> </span><span class="k">as</span><span class="w"> </span><span class="nn">xr</span>
 
 <span class="kn">from</span><span class="w"> </span><span class="nn">htrdrPy.helperFunctions</span><span class="w"> </span><span class="kn">import</span> <span class="o">*</span>
 <span class="kn">from</span><span class="w"> </span><span class="nn">htrdrPy.data</span><span class="w"> </span><span class="kn">import</span> <span class="o">*</span>
@@ -94,7 +95,7 @@
 
 <span class="sd">    &gt;&gt;&gt; pp = htrdrPy.Postprocess(script)</span>
 
-<span class="sd">    This will generate the output files (&#39;.json&#39;, images, ...) in a</span>
+<span class="sd">    This will generate the output files (&#39;.nc&#39;, &#39;.json&#39;, images, ...) in a</span>
 <span class="sd">    &#39;result_{name}/&#39; folder, where the &#39;name&#39; is the name of the</span>
 <span class="sd">    :py:class:`.Data` instance passed to the :py:class:`.Script`.</span>
 <span class="sd">    &quot;&quot;&quot;</span>
@@ -162,7 +163,7 @@
 
 <span class="sd">        Notes</span>
 <span class="sd">        -----</span>
-<span class="sd">        Im most cases, the module recognize the kind of script and automatically</span>
+<span class="sd">        In most cases, the module recognize the kind of script and automatically</span>
 <span class="sd">        apply the required post-process function. This is not the case for</span>
 <span class="sd">        script based on &quot;startMultipleObsGeometry&quot;, as the post-processing</span>
 <span class="sd">        routine to be used is not trivial. A bunch of additional methods are</span>
@@ -623,6 +624,56 @@
         <span class="n">deviationCellDiffuse</span><span class="p">[:,</span><span class="o">-</span><span class="mi">1</span><span class="p">,:]</span> <span class="o">=</span> <span class="n">deviationCellDiffuse</span><span class="p">[:,</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span><span class="mi">0</span><span class="p">][:,</span><span class="kc">None</span><span class="p">]</span>
         <span class="n">varianceCellDiffuse</span><span class="p">[:,</span><span class="o">-</span><span class="mi">1</span><span class="p">,:]</span> <span class="o">=</span> <span class="n">varianceCellDiffuse</span><span class="p">[:,</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span><span class="mi">0</span><span class="p">][:,</span><span class="kc">None</span><span class="p">]</span>
 
+        <span class="n">dicNC</span> <span class="o">=</span> <span class="p">{</span>
+                <span class="s1">&#39;coords&#39;</span><span class="p">:</span> <span class="p">{</span>
+                    <span class="s1">&#39;longitude&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:</span> <span class="s2">&quot;longitude&quot;</span><span class="p">,</span> <span class="s1">&#39;data&#39;</span><span class="p">:</span><span class="n">longitudes</span><span class="p">,</span>
+                                  <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Longitude&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;°&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;latitude&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:</span> <span class="s2">&quot;latitude&quot;</span><span class="p">,</span> <span class="s1">&#39;data&#39;</span><span class="p">:</span><span class="n">latitudes</span><span class="p">,</span>
+                                 <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Latitude&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;°&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;altitude&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:</span> <span class="p">(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                                 <span class="s1">&#39;data&#39;</span><span class="p">:</span><span class="n">altitudes</span><span class="p">,</span>
+                                 <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Altitude&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;m&quot;</span><span class="p">}},</span>
+                    <span class="p">},</span>
+                <span class="s1">&#39;data_vars&#39;</span><span class="p">:</span> <span class="p">{</span>
+                    <span class="s1">&#39;mean&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">meanCellTotal</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;variance&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">varianceCellTotal</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Variance of the radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;std_dev&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellTotal</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Standard deviation of the radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;mean_direct&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">meanCellDirect</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Direct component of Radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;variance_direct&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">varianceCellDirect</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Variance of the Direct component of radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;std_dev_direct&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellDirect</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Standard deviation of the Direct component of radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;mean_diffuse&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">meanCellDiffuse</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Diffuse component of Radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;variance_diffuse&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">varianceCellDiffuse</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Variance of the Diffuse component of radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="s1">&#39;std_dev_diffuse&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                             <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellDiffuse</span><span class="p">,</span>
+                             <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Standard deviation of the Diffuse component of radiative budget&quot;</span><span class="p">,</span>
+                                      <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;W/m3&quot;</span><span class="p">}},</span>
+                    <span class="p">}</span>
+                <span class="p">}</span>
+
         <span class="n">dic</span> <span class="o">=</span> <span class="p">{</span>
                 <span class="s1">&#39;latitudes&#39;</span><span class="p">:</span> <span class="n">latitudes</span><span class="o">.</span><span class="n">tolist</span><span class="p">(),</span>
                 <span class="s1">&#39;longitudes&#39;</span><span class="p">:</span> <span class="n">longitudes</span><span class="o">.</span><span class="n">tolist</span><span class="p">(),</span>
@@ -650,11 +701,31 @@
             <span class="n">dic</span><span class="p">[</span><span class="s1">&#39;heating rate diffuse&#39;</span><span class="p">]</span> <span class="o">=</span> <span class="n">dT_dt_diffuse</span><span class="o">.</span><span class="n">tolist</span><span class="p">()</span>
             <span class="n">dic</span><span class="p">[</span><span class="s1">&#39;heating rate diffuse deviation&#39;</span><span class="p">]</span> <span class="o">=</span> <span class="p">(</span><span class="n">deviationCellDiffuse</span> <span class="o">/</span> <span class="p">(</span><span class="n">rho</span> <span class="o">*</span> <span class="n">heatCapacity</span><span class="p">))</span><span class="o">.</span><span class="n">tolist</span><span class="p">()</span>
 
-            <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;</span><span class="si">{</span><span class="bp">self</span><span class="o">.</span><span class="n">resultsPath</span><span class="si">}</span><span class="s2">radiativeBudget.json&quot;</span><span class="p">,</span> <span class="s1">&#39;w&#39;</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
-                <span class="n">f</span><span class="o">.</span><span class="n">write</span><span class="p">(</span><span class="n">json</span><span class="o">.</span><span class="n">dumps</span><span class="p">(</span><span class="n">dic</span><span class="p">,</span> <span class="n">indent</span><span class="o">=</span><span class="mi">4</span><span class="p">))</span>
-        <span class="k">else</span><span class="p">:</span>
-            <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;</span><span class="si">{</span><span class="bp">self</span><span class="o">.</span><span class="n">resultsPath</span><span class="si">}</span><span class="s2">radiativeBudget.json&quot;</span><span class="p">,</span> <span class="s1">&#39;w&#39;</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
-                <span class="n">f</span><span class="o">.</span><span class="n">write</span><span class="p">(</span><span class="n">json</span><span class="o">.</span><span class="n">dumps</span><span class="p">(</span><span class="n">dic</span><span class="p">,</span> <span class="n">indent</span><span class="o">=</span><span class="mi">4</span><span class="p">))</span>
+            <span class="n">dicNC</span><span class="p">[</span><span class="s1">&#39;data_vars&#39;</span><span class="p">]</span><span class="o">.</span><span class="n">update</span><span class="p">({</span>
+                <span class="s1">&#39;dtsw&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">dT_dt_total</span><span class="p">,</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Heating Rate&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="s1">&#39;dtsw_dev&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellTotal</span> <span class="o">/</span> <span class="p">(</span><span class="n">rho</span> <span class="o">*</span> <span class="n">heatCapacity</span><span class="p">),</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Heating Rate deviation&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="s1">&#39;dtsw_direct&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">dT_dt_direct</span><span class="p">,</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Direct component of Heating Rate&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="s1">&#39;dtsw_direct_dev&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellDirect</span> <span class="o">/</span> <span class="p">(</span><span class="n">rho</span> <span class="o">*</span> <span class="n">heatCapacity</span><span class="p">),</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Direct component of Heating Rate deviation&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="s1">&#39;dtsw_diffuse&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">dT_dt_diffuse</span><span class="p">,</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Diffuse component of Heating Rate&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="s1">&#39;dtsw_diffuse_dev&#39;</span><span class="p">:</span> <span class="p">{</span><span class="s1">&#39;dims&#39;</span><span class="p">:(</span><span class="s2">&quot;altitude&quot;</span><span class="p">,</span><span class="s2">&quot;latitude&quot;</span><span class="p">,</span><span class="s2">&quot;longitude&quot;</span><span class="p">),</span>
+                         <span class="s1">&#39;data&#39;</span><span class="p">:</span> <span class="n">deviationCellDiffuse</span> <span class="o">/</span> <span class="p">(</span><span class="n">rho</span> <span class="o">*</span> <span class="n">heatCapacity</span><span class="p">),</span>
+                         <span class="s1">&#39;attrs&#39;</span><span class="p">:{</span><span class="s1">&#39;long_name&#39;</span><span class="p">:</span><span class="s2">&quot;Diffuse component of Heating Rate deviation&quot;</span><span class="p">,</span> <span class="s1">&#39;units&#39;</span><span class="p">:</span><span class="s2">&quot;K/s&quot;</span><span class="p">}},</span>
+                <span class="p">})</span>
+
+        <span class="n">xr</span><span class="o">.</span><span class="n">Dataset</span><span class="o">.</span><span class="n">from_dict</span><span class="p">(</span><span class="n">dicNC</span><span class="p">)</span><span class="o">.</span><span class="n">to_netcdf</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;</span><span class="si">{</span><span class="bp">self</span><span class="o">.</span><span class="n">resultsPath</span><span class="si">}</span><span class="s2">radiativeBudget.nc&quot;</span><span class="p">)</span>
+
+        <span class="k">with</span> <span class="nb">open</span><span class="p">(</span><span class="sa">f</span><span class="s2">&quot;</span><span class="si">{</span><span class="bp">self</span><span class="o">.</span><span class="n">resultsPath</span><span class="si">}</span><span class="s2">radiativeBudget.json&quot;</span><span class="p">,</span> <span class="s1">&#39;w&#39;</span><span class="p">)</span> <span class="k">as</span> <span class="n">f</span><span class="p">:</span>
+            <span class="n">f</span><span class="o">.</span><span class="n">write</span><span class="p">(</span><span class="n">json</span><span class="o">.</span><span class="n">dumps</span><span class="p">(</span><span class="n">dic</span><span class="p">,</span> <span class="n">indent</span><span class="o">=</span><span class="mi">4</span><span class="p">))</span>
         <span class="k">return</span>
 
     <span class="k">def</span><span class="w"> </span><span class="nf">__twoStreamV1</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">heatCapacity</span><span class="p">):</span>

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/htrdrPy.data.html

@@ -137,14 +137,15 @@ to name the input folders.</p></li>
 <dd class="field-odd"><ul class="simple">
 <li><p><strong>dim</strong> (<em>int</em>) – Dimension of the original data. If set to 0, the case is considered
 as a plan-parallel.</p></li>
-<li><p><strong>data</strong> (<em>dict</em>) – <p>Dictionnary with the following items:
-- “nAngle”: int, optional</p>
-<blockquote>
-<div><p>Number of angles in the phase functions (not required if
-using the builtin Henyey-Greenstein phase function).</p>
-</div></blockquote>
+<li><p><strong>data</strong> (<em>dict</em>) – <p>Dictionnary with the following items:</p>
 <ul>
 <li><dl class="simple">
+<dt>”nAngle”: int, optional</dt><dd><p>Number of angles in the phase functions (not required if
+using the builtin Henyey-Greenstein phase function).</p>
+</dd>
+</dl>
+</li>
+<li><dl class="simple">
 <dt>”scattering (m-1)”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>Array of scattering coefficients
 (shape=(nWavelength, nLevel, [nLat, nLon]), [m-1]).</p>
 </dd>
@@ -1087,13 +1088,14 @@ instance are not used and instead are defined from the length of the
 <dl class="field-list simple">
 <dt class="field-odd">Parameters<span class="colon">:</span></dt>
 <dd class="field-odd"><ul class="simple">
-<li><p><strong>data</strong> (<em>dict</em>) – <p>Dictionnary with the following items:
-- “nLevel” : int</p>
-<blockquote>
-<div><p>Number of levels. Required only if sameGrid is <code class="docutils literal notranslate"><span class="pre">False</span></code>.</p>
-</div></blockquote>
+<li><p><strong>data</strong> (<em>dict</em>) – <p>Dictionnary with the following items:</p>
 <ul>
 <li><dl class="simple">
+<dt>”nLevel”<span class="classifier">int</span></dt><dd><p>Number of levels. Required only if sameGrid is <code class="docutils literal notranslate"><span class="pre">False</span></code>.</p>
+</dd>
+</dl>
+</li>
+<li><dl class="simple">
 <dt>”nCoeff”: int</dt><dd><p>Number of quadrature points for the k-coeff.</p>
 </dd>
 </dl>
@@ -1115,7 +1117,7 @@ using the builtin Henyey-Greenstein phase function).</p>
 </dl>
 </li>
 <li><dl class="simple">
-<dt>”altitude (m)”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>Array of altitudes (shape=(nLevel, nLat, nLon), [m]).</p>
+<dt>”altitude (m)”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>Array of altitudes (shape=(nLevel, [nLat, nLon]), [m]).</p>
 </dd>
 </dl>
 </li>
@@ -1132,6 +1134,11 @@ is 3.</p>
 </dl>
 </li>
 <li><dl class="simple">
+<dt>”temperature (K)”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>Array of temperatures (shape=(nLevel, [nLat, nLon]), [K]).</p>
+</dd>
+</dl>
+</li>
+<li><dl class="simple">
 <dt>”nWavelength”:  int</dt><dd><p>Number of wavelengths.</p>
 </dd>
 </dl>
@@ -1156,6 +1163,12 @@ values corresponds to the lower boundary of the band.</p>
 </dl>
 </li>
 <li><dl class="simple">
+<dt>”weights”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>The weigths of the nCoeff quadrature points
+(shape=(nCoeff)).</p>
+</dd>
+</dl>
+</li>
+<li><dl class="simple">
 <dt>”scattering (m-1)”: <code class="docutils literal notranslate"><span class="pre">numpy.ndarray</span></code></dt><dd><p>Array of scattering coefficients
 (shape=(nWavelength, nLevel, [nLat, nLon]), [m-1]).</p>
 </dd>

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/htrdrPy.geometry.html

@@ -505,6 +505,13 @@ with the altitude, latitude and longitude indices.</p>
 </ul>
 </dd>
 </dl>
+<div class="admonition warning">
+<p class="admonition-title">Warning</p>
+<p>In the table of coordiantes that must be provided, the latitude should
+be given by decreasing latitudes (North pole corresponding to the first
+latitude index and South pole to the last), according to LMDZ PCM
+convention.</p>
+</div>
 </dd></dl>
 
 </dd></dl>

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/docs/build/html/htrdrPy.postprocess.html

@@ -107,7 +107,7 @@ provide the <a class="reference internal" href="htrdrPy.script.html#htrdrPy.scri
 <div class="doctest highlight-default notranslate"><div class="highlight"><pre><span></span><span class="gp">&gt;&gt;&gt; </span><span class="n">pp</span> <span class="o">=</span> <span class="n">htrdrPy</span><span class="o">.</span><span class="n">Postprocess</span><span class="p">(</span><span class="n">script</span><span class="p">)</span>
 </pre></div>
 </div>
-<p>This will generate the output files (‘.json’, images, …) in a
+<p>This will generate the output files (‘.nc’, ‘.json’, images, …) in a
 ‘result_{name}/’ folder, where the ‘name’ is the name of the
 <a class="reference internal" href="htrdrPy.data.html#htrdrPy.data.Data" title="htrdrPy.data.Data"><code class="xref py py-class docutils literal notranslate"><span class="pre">Data</span></code></a> instance passed to the <a class="reference internal" href="htrdrPy.script.html#htrdrPy.script.Script" title="htrdrPy.script.Script"><code class="xref py py-class docutils literal notranslate"><span class="pre">Script</span></code></a>.</p>
 <dl class="field-list simple">
@@ -207,7 +207,7 @@ divergence.</p>
 </dd>
 </dl>
 <p class="rubric">Notes</p>
-<p>Im most cases, the module recognize the kind of script and automatically
+<p>In most cases, the module recognize the kind of script and automatically
 apply the required post-process function. This is not the case for
 script based on “startMultipleObsGeometry”, as the post-processing
 routine to be used is not trivial. A bunch of additional methods are

htrdrpy-1.1.7/docs/build/html/searchindex.js

@@ -0,0 +1 @@
+Search.setIndex({"alltitles":{"Arguments":[[4,"id2"]],"Building method":[[4,"id1"]],"Data":[[1,"data"]],"Examples":[[0,null]],"Generating input files":[[1,"generating-input-files"]],"Imports":[[1,"imports"]],"Installation":[[8,"installation"]],"Introduction":[[8,"introduction"]],"Observation geometry":[[1,"observation-geometry"]],"Post-process":[[1,"post-process"]],"Recover a Script":[[9,null]],"Scripts":[[1,"scripts"]],"Spectra: comparison between plan-parallel and spherical":[[1,null]],"Summary":[[8,"summary"]],"htrdrPy":[[8,null]],"htrdrPy documentation":[[2,null]],"htrdrPy.data module":[[3,null]],"htrdrPy.geometry module":[[4,null]],"htrdrPy.helperFunctions module":[[5,null]],"htrdrPy.postprocess module":[[6,null]],"htrdrPy.script module":[[7,null]],"kwargs":[[6,"id1"]]},"docnames":["Examples","example_comparison_PP_sphere","htrdrPy","htrdrPy.data","htrdrPy.geometry","htrdrPy.helperFunctions","htrdrPy.postprocess","htrdrPy.script","index","loadScript"],"envversion":{"sphinx":65,"sphinx.domains.c":3,"sphinx.domains.changeset":1,"sphinx.domains.citation":1,"sphinx.domains.cpp":9,"sphinx.domains.index":1,"sphinx.domains.javascript":3,"sphinx.domains.math":2,"sphinx.domains.python":4,"sphinx.domains.rst":2,"sphinx.domains.std":2,"sphinx.ext.viewcode":1},"filenames":["Examples.rst","example_comparison_PP_sphere.rst","htrdrPy.rst","htrdrPy.data.rst","htrdrPy.geometry.rst","htrdrPy.helperFunctions.rst","htrdrPy.postprocess.rst","htrdrPy.script.rst","index.rst","loadScript.rst"],"indexentries":{"addparticle() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.addParticle",false]],"bandintegratedimage() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.bandIntegratedImage",false]],"cleaninputs() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.cleanInputs",false]],"compositerbg() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.compositeRBG",false]],"data (class in htrdrpy.data)":[[3,"htrdrPy.data.Data",false]],"exportgeometry() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.exportGeometry",false]],"extractmeanradiancefromoutput() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.extractMeanRadianceFromOutput",false]],"extractmeanradiances() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.extractMeanRadiances",false]],"geometry (class in htrdrpy.geometry)":[[4,"htrdrPy.geometry.Geometry",false]],"getimage() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.getImage",false]],"getimages() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.getImages",false]],"htrdrpy.data":[[3,"module-htrdrPy.data",false]],"htrdrpy.geometry":[[4,"module-htrdrPy.geometry",false]],"htrdrpy.postprocess":[[6,"module-htrdrPy.postprocess",false]],"htrdrpy.script":[[7,"module-htrdrPy.script",false]],"imageratio() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.imageRatio",false]],"loadscript() (in module htrdrpy.script)":[[7,"htrdrPy.script.loadScript",false]],"makeatmospherefrom1d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeAtmosphereFrom1D",false]],"makeatmospherefrom1d_pp() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeAtmosphereFrom1D_PP",false]],"makeatmospherefrom2d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeAtmosphereFrom2D",false]],"makeatmospherefrom3d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeAtmosphereFrom3D",false]],"makegeomfromapie() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.makeGeomFromAPIE",false]],"makegroundfrom1d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeGroundFrom1D",false]],"makegroundfrom1d_pp() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeGroundFrom1D_PP",false]],"makegroundfrom2d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeGroundFrom2D",false]],"makegroundfrom3d() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeGroundFrom3D",false]],"makemixture() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.makeMixture",false]],"module":[[3,"module-htrdrPy.data",false],[4,"module-htrdrPy.geometry",false],[6,"module-htrdrPy.postprocess",false],[7,"module-htrdrPy.script",false]],"monochromaticimage() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.monochromaticImage",false]],"plotgeometry() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.plotGeometry",false]],"postprocess (class in htrdrpy.postprocess)":[[6,"htrdrPy.postprocess.Postprocess",false]],"processimages() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.processImages",false]],"processsinglearrayobssw() (htrdrpy.postprocess.postprocess method)":[[6,"htrdrPy.postprocess.Postprocess.processSingleArrayObsSW",false]],"reflectancespectrum() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.reflectanceSpectrum",false]],"script (class in htrdrpy.script)":[[7,"htrdrPy.script.Script",false]],"setcamera() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.setCamera",false]],"setimage() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.setImage",false]],"setsource() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.setSource",false]],"setspectralcumuldist() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.setSpectralCumulDist",false]],"setvolrad() (htrdrpy.geometry.geometry method)":[[4,"htrdrPy.geometry.Geometry.setVolrad",false]],"spectrum() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.spectrum",false]],"startmultipleobsgeometry() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.startMultipleObsGeometry",false]],"startradbudgetgcm() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.startRadBudgetGCM",false]],"visibleimage() (htrdrpy.script.script method)":[[7,"htrdrPy.script.Script.visibleImage",false]],"writeinputatmosphere() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeInputAtmosphere",false]],"writeinputground() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeInputGround",false]],"writeinputs() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeInputs",false]],"writevtkfiles() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeVTKfiles",false]],"writevtkfilesatmosphere() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeVTKfilesAtmosphere",false]],"writevtkfilesground() (htrdrpy.data.data method)":[[3,"htrdrPy.data.Data.writeVTKfilesGround",false]]},"objects":{"htrdrPy":[[3,0,0,"-","data"],[4,0,0,"-","geometry"],[5,0,0,"-","helperFunctions"],[6,0,0,"-","postprocess"],[7,0,0,"-","script"]],"htrdrPy.data":[[3,1,1,"","Data"]],"htrdrPy.data.Data":[[3,2,1,"","addParticle"],[3,2,1,"","cleanInputs"],[3,2,1,"","makeAtmosphereFrom1D"],[3,2,1,"","makeAtmosphereFrom1D_PP"],[3,2,1,"","makeAtmosphereFrom2D"],[3,2,1,"","makeAtmosphereFrom3D"],[3,2,1,"","makeGroundFrom1D"],[3,2,1,"","makeGroundFrom1D_PP"],[3,2,1,"","makeGroundFrom2D"],[3,2,1,"","makeGroundFrom3D"],[3,2,1,"","makeMixture"],[3,2,1,"","writeInputAtmosphere"],[3,2,1,"","writeInputGround"],[3,2,1,"","writeInputs"],[3,2,1,"","writeVTKfiles"],[3,2,1,"","writeVTKfilesAtmosphere"],[3,2,1,"","writeVTKfilesGround"]],"htrdrPy.geometry":[[4,1,1,"","Geometry"]],"htrdrPy.geometry.Geometry":[[4,2,1,"","exportGeometry"],[4,2,1,"","makeGeomFromAPIE"],[4,2,1,"","plotGeometry"],[4,2,1,"","setCamera"],[4,2,1,"","setImage"],[4,2,1,"","setSource"],[4,2,1,"","setSpectralCumulDist"],[4,2,1,"","setVolrad"]],"htrdrPy.helperFunctions":[[5,3,1,"","cart2sphere"],[5,3,1,"","combineEstimates"],[5,3,1,"","dplanck_dT"],[5,3,1,"","planck"],[5,3,1,"","sphere2cart"]],"htrdrPy.postprocess":[[6,1,1,"","Postprocess"]],"htrdrPy.postprocess.Postprocess":[[6,2,1,"","extractMeanRadianceFromOutput"],[6,2,1,"","extractMeanRadiances"],[6,2,1,"","getImage"],[6,2,1,"","getImages"],[6,2,1,"","processImages"],[6,2,1,"","processSingleArrayObsSW"]],"htrdrPy.script":[[7,1,1,"","Script"],[7,3,1,"","loadScript"]],"htrdrPy.script.Script":[[7,2,1,"","bandIntegratedImage"],[7,2,1,"","compositeRBG"],[7,2,1,"","imageRatio"],[7,2,1,"","monochromaticImage"],[7,2,1,"","reflectanceSpectrum"],[7,2,1,"","spectrum"],[7,2,1,"","startMultipleObsGeometry"],[7,2,1,"","startRadBudgetGCM"],[7,2,1,"","visibleImage"]]},"objnames":{"0":["py","module","Python module"],"1":["py","class","Python class"],"2":["py","method","Python method"],"3":["py","function","Python function"]},"objtypes":{"0":"py:module","1":"py:class","2":"py:method","3":"py:function"},"terms":{"":[4,6],"0":[1,3,4,6],"00":[],"04":[],"1":[1,3,5,6],"10":1,"100":[],"1000":1,"1000000":3,"1024":3,"11":[],"1323":3,"15":1,"1500000":3,"180":3,"181":[],"1d":1,"1d_spetral_data":1,"1e":[3,4],"1e6":3,"1e9":1,"2":[1,3,4,5,6,7],"20":[3,4],"2744":[],"2d":3,"2e7":3,"3":[1,3,4,7],"30":3,"300":3,"338541":[],"360":3,"3d":[1,3,4],"4":[1,3],"4802":3,"5":3,"50":3,"500":3,"512":3,"547428":[],"5586":[],"575e6":1,"5800":1,"5e5":3,"64":1,"64827":3,"69it":[],"7e8":1,"8":3,"80":1,"9":1,"90":3,"91":1,"95it":[],"98":3,"9e7":3,"A":[1,6,7],"At":1,"For":[1,7],"If":[3,4,5,6,7],"In":[1,3,4,6,9],"It":[4,9],"NOT":4,"Not":[3,4],"The":[1,3,4,5,6,7,9],"Then":1,"These":1,"To":1,"_":9,"_________________________________________________________________________________________________________________________________________________________":[],"absor":3,"absorpt":[1,3],"accord":4,"accordingli":1,"activ":7,"actual":1,"add":[1,3],"addit":[1,6],"additon":1,"addparticl":[2,3,8],"aerosol":1,"affect":4,"after":[7,9],"aim":[1,3,4,6,7],"albedo":[1,3],"all":[1,3,4,6,7,9],"along":[1,3],"alreadi":[1,3,4,7],"also":[1,3],"altern":3,"altitud":[1,3,4,5],"an":[0,1,3,4,5,7,9],"angl":[1,3,4],"ani":1,"append":1,"appli":6,"ar":[1,3,4,5,6,7,9],"arg":4,"argument":6,"around":[4,7],"arrai":[1,3,4,5,6,7],"array_lik":4,"aspect":4,"assess":3,"assign":3,"associ":6,"assum":4,"assym":[],"assymetri":[],"asymetri":[],"asymm":3,"asymmetri":3,"ath":7,"atmospher":[1,3,4,6,7],"au":1,"auto":6,"automat":[1,6],"averag":6,"ax":[1,4],"axi":4,"azimut":[1,4],"band":[3,7],"bandintegratedimag":[2,6,7,8],"bandslow":3,"bandsup":3,"bandwidth":7,"base":[3,4,6,7],"becaus":1,"bee":[],"been":[4,6,7,9],"befor":[1,7,9],"being":[1,7],"below":6,"between":[0,3,4,8,9],"bin":[3,9],"binari":[3,7,9],"black":4,"blank":3,"blue":4,"bool":[3,4,6,7],"both":[1,5],"bound":7,"boundari":[3,4,7],"brdf":[1,3],"budget":[4,7],"built":[],"builtin":3,"bunch":[1,6],"c":[1,4,6,9],"calcul":[1,3,4,5,6,7,9],"call":[1,6,7,9],"callabl":7,"calul":[1,4],"camdist":4,"camera":[1,4],"cameradist":1,"camerafov":[],"camfov":4,"can":[1,3,4,6,7,8,9],"capac":6,"carlo":5,"cart2spher":[2,5,8],"cartesian":[4,5],"case":[1,3,4,6,7],"casenam":7,"cell":[3,4,7],"cellcoord":4,"center":[3,4],"central":[],"chanel":7,"chang":1,"chosen":4,"cie_xyz":7,"class":[3,4,6,7],"cleaninput":[2,3,8],"cloud":3,"cmap":6,"co":[],"code":9,"coeff":3,"coeffic":[],"coefficeint":[],"coeffici":[3,4],"color":6,"column":[1,3,4,6],"com":8,"combin":5,"combineestim":[2,5,8],"command":[1,7,8],"compar":1,"comparison":[0,8,9],"complet":[3,4],"compos":[],"composit":7,"compositerbg":[2,6,7,8],"compris":1,"comput":[3,6,7],"concern":1,"conduct":4,"consid":[3,6],"constant":1,"constraint":1,"contain":[1,4,5,6,7],"context":1,"contina":4,"convent":4,"convert":[1,5],"coordiant":4,"coordin":[4,5],"core":7,"correct":1,"correctli":7,"correl":4,"correspond":[1,3,4,6,7],"could":7,"count":4,"counter":3,"crash":9,"creat":[1,3,4,6,7,9],"creatin":6,"creation":[1,9],"cst":1,"cumul":4,"current":[1,4],"d":[3,4,5,6],"data":[0,2,4,6,7,8,9],"datapp":1,"dataset":3,"datasp":1,"decreas":4,"default":[3,4,6,7],"defin":[1,3,7],"definit":[1,4,6],"defint":3,"denomin":7,"densiti":[4,6],"depend":[1,3,5,6,7],"deriv":5,"desctruct":4,"detail":1,"determin":1,"deviat":[1,5,6],"dict":[1,3,4,6,7],"dictionari":6,"dictionnari":[1,3,4,6,7],"differ":[1,7],"differenti":7,"dim":[1,3],"dim_cel":[],"dim_nod":[],"dimens":[1,3,4,5],"direct":4,"directli":[1,4],"directori":6,"discret":3,"disctionnari":1,"disk":3,"distanc":[1,4,5],"distinct":4,"distribut":[1,4],"diverg":[4,6],"dlongitud":4,"do":[3,9],"document":[1,6,7,8],"don":9,"done":1,"down":1,"download":1,"dplanck_dt":[2,5,8],"drax":4,"dtype":1,"e":[3,4],"each":[1,4,5,6,7],"earlier":1,"either":[4,5],"emerg":[1,4],"emiss":[5,7],"emit":5,"empti":7,"enough":[1,3],"errorbar":1,"estim":5,"etc":1,"evalu":4,"everi":[],"exampl":[1,3,6,7,8,9],"example_comparison_pp_spher":[],"example_spectra_data":1,"execut":[6,7,9],"exist":1,"expans":3,"export":4,"exportgeometri":[1,2,4,8],"exposur":6,"express":6,"extens":1,"extern":7,"extinct":[],"extract":[3,4,6],"extractfromdata":4,"extractmeanradi":[2,6,8],"extractmeanradiancefromoutput":[2,6,8],"f":[1,4,6,9],"fals":[3,4,5,6],"field":4,"fieldofview":4,"fig":1,"file":[0,3,4,6,7,8,9],"filenam":[3,7],"final":[1,4,5,7],"find":[1,6],"first":[1,4,6,7],"flaot":5,"float":[1,3,4,5,6,7],"flux":[4,6],"folder":[3,6,7],"follow":[3,4,6,7,8,9],"fontsiz":1,"form":7,"format":4,"found":1,"frac":1,"from":[1,3,4,5,6,7,8],"full":1,"fulli":1,"func":[],"function":[1,3,4,6,7],"further":1,"futur":[3,7],"g":[],"ga":3,"gase":1,"gaussian":[],"gcm":7,"gener":[0,3,4,6,7,8],"geometri":[0,2,6,7,8],"geometry_":4,"geometryfromapi":[],"geompp":1,"geomsp":1,"geopotenti":3,"geopotentialsurf":[],"getimag":[2,6,8],"given":[1,4,5,6,7,9],"gometri":4,"graviti":3,"green":[4,7],"greenstein":3,"grid":[3,6],"ground":[3,4],"gv":[],"ha":[1,4,5,7,9],"handl":[1,3,4],"happen":9,"happi":1,"have":[1,4,6,9],"haze":3,"heat":6,"heatcapac":6,"helperfunct":[2,8],"henyei":3,"here":1,"heteorgen":1,"heterogen":[],"hexahedron":3,"hg":[],"homogen":[],"horizont":[1,3,4],"how":1,"html":8,"htrdpy":[],"htrdr":[1,3,4,6,7,8,9],"htrdrpy":[0,1,9],"http":8,"hwidth":4,"i":[0,1,3,4,5,6,7,8,9],"identifi":7,"if_mcc":1,"if_ob":1,"if_shdom":1,"im":[],"imag":[1,4,6,7],"imageratio":[2,6,7,8],"import":[0,8],"incid":[1,4],"includ":8,"incom":1,"inconsist":1,"index":[4,6,8],"indic":[1,4,6],"inferno":6,"inform":[1,3,4,6,7,9],"initialis":3,"input":[0,3,7,8],"input_":[],"inputs_":[1,3],"insatnc":[1,7],"instanc":[1,3,4,6,7,9],"instead":3,"int":[3,4],"integr":7,"interfac":[3,4],"issu":[1,4,9],"istanc":3,"item":[1,3,4,6,7],"its":[1,9],"j":6,"json":[1,4,6],"jupyt":9,"k":[1,3,4,5,6],"keep":3,"kei":[1,3,4],"keyslmdztitan_vi":[],"keyword":3,"kg":[3,6],"kilo":1,"kind":[1,3,6,7],"kv":[],"l":1,"label":1,"lambda":[],"lambertian":[1,3],"larg":[1,3],"larger":1,"last":[1,4],"lat":[],"latitud":[1,3,4,5],"launch":3,"layer":3,"least":4,"left":7,"legend":1,"len":1,"length":[3,5],"level":[3,4],"librari":1,"like":[4,6,7],"limit":[3,6],"line":[1,4],"linspac":3,"list":[1,3,4,5,6,7],"lmdz":[3,4],"lmdzinput":[],"lmdzouput":[],"loaction":4,"load":[1,3,7,9],"loadscript":[2,7,8,9],"localpath":6,"locat":[3,4,6,9],"log":[],"lon":[],"longitud":[1,3,4,5],"longitudin":3,"low":[3,7],"lower":[3,7],"lw":[1,7],"m":[1,3,4,5,7],"m2":[4,5],"m3":6,"mai":[4,7],"make":[1,3,9],"makeatmospherefrom1d":[2,3,8],"makeatmospherefrom1d_pp":[2,3,8],"makeatmospherefrom2d":[2,3,8],"makeatmospherefrom3d":[2,3,8],"makecolumnpp":4,"makefromcellcoord":4,"makefromlmdz":[],"makegeomfromapi":[1,2,4,8],"makeground":1,"makegroundfrom1d":[1,2,3,8],"makegroundfrom1d_pp":[1,2,3,8],"makegroundfrom2d":[2,3,8],"makegroundfrom3d":[2,3,8],"makemixtur":[1,2,3,8],"makeslicealtlat":4,"manag":[4,8],"manatori":3,"mandatori":[3,4],"mani":[3,7],"map":[1,6],"marker":1,"mass":[3,6],"master":3,"matplotlib":[1,4],"max":6,"maxim":3,"maximal":6,"maximum":7,"mccsphe":1,"md":[],"mean":[1,5,6],"meant":[3,8],"memori":3,"merg":3,"mesh":[3,4],"meso":8,"method":[1,6,7],"micro":1,"mixtur":[1,3],"mode":[3,4],"modul":[1,2,8],"monochromat":7,"monochromaticimag":[2,6,7,8],"mont":5,"more":4,"most":6,"mostli":7,"moveaxi":3,"mpi":7,"mpicmd":[1,7],"mpirun":[1,7],"multipl":[1,6,7],"must":[1,3,4,8],"mutlipl":7,"n":5,"n_hexahedron":3,"n_octahedron":3,"n_rectangl":3,"n_triangl":3,"naltitud":4,"name":[1,3,6,7,9],"nangl":[1,3],"nc":6,"ncell":4,"ncoeff":[1,3,4],"ndarrai":[3,4,5,6,7],"necessari":1,"need":[1,6],"netcdf":[],"new":9,"next":1,"nlat":3,"nlatitud":4,"nlayer":4,"nlevel":[1,3,4],"nlon":3,"nlongitud":4,"nm":4,"nnode":[],"node":3,"noeds":3,"non":4,"none":[3,4,6,7],"normal":4,"north":4,"note":[1,3,4,6],"notebook":9,"now":[1,9],"np":[1,3],"nphi":[1,3],"npix":6,"npixel":6,"npixelx":6,"npz":1,"nspectralel":4,"ntheta":[1,3],"nthoctre":3,"num":7,"number":[1,3,4,5,7],"numer":7,"numpi":[1,3,4,5,6,7],"nwavelength":[1,3,4,5,7],"nweight":[],"nz":3,"o":1,"ob":1,"obj":3,"object":[1,3,4,6,7],"observ":[0,4,7,8],"obslist":7,"obtain":3,"obvious":1,"octahedron":3,"octre":3,"octree_def":3,"octreefil":3,"omit":[],"onc":[1,9],"one":[1,4,6],"ones":[1,3],"ones_lik":3,"onlevel":4,"onli":[1,3,4,6],"open":1,"oper":6,"opthick":3,"optic":[1,3],"option":[3,4,5,6,7],"order":[1,6],"orient":4,"origin":[3,4],"other":[3,6],"otherwis":[6,9],"output":[1,3,6,7,9],"outputs_":[1,3,9],"outputs_planet":3,"outputs_spher":9,"over":[6,7],"p":[],"packag":[0,8],"panic":9,"parallel":[0,3,4,8,9],"paramet":[1,3,4,5,6,7],"paraview":3,"particl":3,"pass":[1,3,6,7],"path":[4,6,7],"pausibl":1,"pcm":[3,4],"pdf":4,"pe":1,"per":[1,6],"perpendicular":4,"phase":[1,3,4],"phasefunc":3,"physic":3,"pip":8,"pixel":[1,4,6],"plan":[0,3,4,8,9],"planck":[2,4,5,8],"plane":4,"planet":[1,3,4,6,7,8,9],"plot":[1,4],"plotgeometri":[2,4,8],"plt":1,"plu":5,"png":1,"point":[1,3,4],"pole":4,"posit":[1,4],"possibl":[1,6,7],"post":[0,6,8,9],"postprocess":[1,2,8,9],"pp":[1,6],"pphi":[],"precalcul":3,"predefin":[1,6,7],"present":1,"pressur":[],"previou":9,"previous":[1,6,7],"print":1,"prioriti":4,"probabl":4,"process":[0,3,6,8,9],"processimag":[2,6,8],"processsinglearrayobssw":[2,6,8],"prococtre":3,"produc":[1,7],"project":[4,8],"properti":[1,3,4],"provid":[1,3,4,5,6,7],"put":3,"py":[],"pyplot":[1,4],"python":[1,8,9],"quadratur":3,"quick":3,"r":[1,3],"r_d":5,"r_max":3,"r_min":3,"rad":[],"radi":[4,7],"radianc":[4,5,6],"radiancefil":4,"radiat":1,"radiu":[1,3,4,5],"rai":4,"rang":[3,7],"rate":6,"ratio":[4,7],"raw":1,"read":1,"readabl":4,"readlin":1,"realis":6,"realiz":[3,4,5],"receiv":5,"recogn":6,"recov":[0,1,6,8],"rectangl":3,"red":[4,7],"refer":[],"reflect":[1,7],"reflectance_spectrum_spectrum":1,"reflectancespectrum":[1,2,6,7,8],"reflectancespectrum_sphere_spectrum":9,"reflexion":3,"regard":5,"rel":6,"relat":4,"remov":3,"repositori":[1,4,9],"repres":4,"reproduc":1,"request":[3,6],"requir":[1,3,4,6],"requirt":1,"reshap":3,"resolut":3,"respect":[3,5],"respp":1,"ressp":1,"restart":9,"result":[1,5,6,8],"result_":6,"results_pp":1,"results_spher":1,"return":[5,6],"review":1,"rgb":7,"rho":6,"right":[1,4,9],"role":[],"roll":4,"root":6,"rotat":1,"routin":[1,6],"run":[3,4,6,7],"runner":7,"s2":3,"same":[1,3,6,7],"samegrid":3,"sampl":[1,4],"save":[3,7],"savefig":1,"scale":4,"scater":[],"scatt":3,"scatter":3,"scenario":[1,9],"scene":1,"scipi":1,"scirpt":1,"script":[0,2,6,8],"scriptpp":1,"scriptsp":[1,9],"second":[4,6,7],"see":[6,7],"sen":4,"send":1,"separ":[1,8],"separt":[],"serv":1,"set":[1,3,4,5,7],"set_xlabel":1,"set_ylabel":1,"set_yscal":[],"setcamera":[2,4,8],"setimag":[1,2,4,8],"setsourc":[1,2,4,8],"setspectralcumuldist":[2,4,8],"setup":4,"setvolrad":[2,4,8],"shape":[1,3,4,5,6,7],"share":3,"shdompp":1,"shift":1,"should":[4,7],"sight":4,"similar":6,"simplfi":8,"simpli":6,"simplifi":1,"sinc":3,"singl":[1,3,6],"size":1,"skip":6,"slice":[1,4],"slurm":7,"socket":1,"softwar":3,"solar":4,"soldist":1,"solrad":1,"soltemp":1,"some":[1,4],"sourc":[1,3,4,5,6,7],"sourcedist":[],"sourcerad":[],"sourcetemp":[],"south":4,"space":3,"specif":[],"specifi":3,"specrtum":1,"spectra":[0,8,9],"spectra_exampl":[],"spectral":[1,7],"spectrum":[1,2,4,6,7,8],"specular":3,"sphere":[1,4],"sphere2cart":[2,5,8],"spheric":[0,3,5,8,9],"split":1,"squar":[4,5],"sr":[4,5],"srcdist":4,"srcpo":4,"srcrad":4,"srctemp":[1,4],"srun":7,"ssa":[],"standard":[5,6,9],"star":8,"start":[1,3,7],"startmultiplegeometri":6,"startmultipleobsgeometri":[2,6,7,8],"startradbudgetgcm":[2,6,7,8],"std":1,"steal":4,"stellar":1,"step":[1,7,9],"storag":3,"store":[1,3,4,7,9],"str":[3,4,6,7],"string":7,"strore":3,"structur":9,"sub":4,"subplot":1,"subset":6,"suit":7,"sum":5,"sumx":5,"sumxsquar":5,"supercalcul":7,"sure":[1,3],"surfac":[1,3,4,5],"surfacetemperatur":3,"surfalb":3,"surftemp":3,"sw":[1,7],"system":[3,7],"t":[1,3,5,7,9],"tabl":[1,3,4,6],"take":7,"taken":4,"target":[1,4],"targetposit":4,"temp":[],"temperatur":[1,3,4,5],"termin":9,"test":[4,7],"tetrahedron":[3,4],"textmu":[],"tgtpo":4,"than":[4,6],"thei":[3,4],"therefor":[1,3,6],"theta":[],"thi":[0,1,3,4,6,7,9],"thick":3,"those":[1,4],"thread":[3,7],"threadflag":7,"threshold":[3,6],"through":1,"ti":[],"tile":3,"time":[1,3,6],"titan":[],"too":1,"tool":6,"treat":[1,7],"treatement":8,"triangl":3,"trivial":6,"true":[3,4,6,7],"tsurf":[],"tupl":4,"turn":4,"two":1,"txt":1,"type":[1,5,6,7,9],"uniform":3,"unit":6,"up":[3,7],"updat":1,"upper":7,"upward":4,"us":[0,1,3,4,6,7,8],"useful":7,"useless":3,"user":[3,6],"v_23":[],"valu":[1,3,4,6],"vari":3,"varianc":5,"vec":5,"vector":4,"verbos":7,"veri":3,"version":7,"vertic":4,"via":[1,4],"view":[3,4],"virtual":7,"visibl":7,"visibleimag":[2,6,7,8],"visual":1,"visualis":3,"volrad":4,"volum":[4,6],"vtk":[1,3],"w":[4,5],"wa":[6,9],"want":4,"wavelength":[1,3,4,5,6,7],"wavelength_vi":[],"wavelengthblu":7,"wavelengthden":7,"wavelengthgreen":7,"wavelengthlow":7,"wavelengthnum":7,"wavelengthr":7,"wavelengthup":7,"wavlength":1,"we":1,"weight":[3,5],"weigth":3,"well":[1,4],"when":[1,3,4,9],"where":[1,3,6,7],"wherea":4,"wheteh":3,"whether":[4,6,7],"which":[1,3,4,6],"width":4,"within":[1,4,9],"won":1,"work":[1,7],"would":1,"wrapper":8,"write":3,"writeinput":[1,2,3,8],"writeinputatmospher":[2,3,8],"writeinputground":[2,3,8],"writevtkfil":[1,2,3,8],"writevtkfilesatmospher":[2,3,8],"writevtkfilesground":[2,3,8],"written":9,"wv":[],"wvl":[1,5],"www":8,"x":[1,4,5],"y":[1,5],"yerr":1,"yet":7,"you":[1,6,7,9],"your":7,"z":[4,5],"zero":[3,6],"\u00b5":6,"\u00b5m":1},"titles":["Examples","Spectra: comparison between plan-parallel and spherical","htrdrPy documentation","htrdrPy.data module","htrdrPy.geometry module","htrdrPy.helperFunctions module","htrdrPy.postprocess module","htrdrPy.script module","htrdrPy","Recover a Script"],"titleterms":{"argument":4,"between":1,"build":4,"calcul":[],"comparison":1,"creat":[],"data":[1,3],"document":2,"exampl":0,"file":1,"gener":1,"geometri":[1,4],"helperfunct":5,"htrdrpy":[2,3,4,5,6,7,8],"import":1,"input":1,"instal":8,"introduct":8,"kwarg":6,"method":4,"modul":[3,4,5,6,7],"object":[],"observ":1,"parallel":1,"plan":1,"plot":[],"post":1,"postprocess":6,"process":1,"read":[],"recov":9,"result":[],"script":[1,7,9],"spectra":1,"spheric":1,"start":[],"summari":8,"vtk":[],"write":[]}})

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "htrdrPy"
-version = "1.1.6"
+version = "1.1.7"
 authors = [
   { name="Anthony Arfaux", email="anthony.arfaux@obspm.fr" },
 ]

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/src/htrdrPy/data.py

@@ -872,7 +872,7 @@ class Data:
         absorption = absorption.reshape((self.nWavelength, self.nNodes, self.nCoeff))
         temperatureNodes = self.temperatureCells.reshape(self.nNodes)
 
-        cell_ids = Delaunay(node_coord)
+        cell_ids = Delaunay(node_coord, qhull_options='QJ QbB Qc Qz')
 
         ##############################################################
         # generating the atmosphere dictionnary used by self.prepareInput()
@@ -1022,7 +1022,7 @@ class Data:
         absorption = absorption.reshape((self.nWavelength, self.nNodes, self.nCoeff))
         temperatureNodes = self.temperatureCells.reshape(self.nNodes)
 
-        cell_ids = Delaunay(node_coord)
+        cell_ids = Delaunay(node_coord, qhull_options='QJ QbB Qc Qz')
 
         ##############################################################
         # generating the atmosphere dictionnary used by self.prepareInput()
@@ -1046,6 +1046,7 @@ class Data:
         ----------
         data : dict
             Dictionnary with the following items: 
+
             - "nLevel" : int
                 Number of levels. Required only if sameGrid is ``False``.
             - "nCoeff": int
@@ -1664,6 +1665,7 @@ class Data:
             as a plan-parallel.
         data : dict
             Dictionnary with the following items: 
+
             - "nAngle": int, optional
                 Number of angles in the phase functions (not required if
                 using the builtin Henyey-Greenstein phase function).

{htrdrpy-1.1.6 → htrdrpy-1.1.7}/src/htrdrPy/geometry.py

@@ -467,11 +467,15 @@ class Geometry:
             pass
         elif phase is not None:
             ap = phase * cst.degree
-            aa = np.acos(
-                         (np.cos(ai) * np.cos(ae) - np.cos(ap))
-                         / (np.sin(ai) * np.sin(ae))
-                         )
+            if ai == 0 or ae == 0:
+                aa = 0
+            else:
+                aa = np.acos(
+                             (np.cos(ai) * np.cos(ae) - np.cos(ap))
+                             / (np.sin(ai) * np.sin(ae))
+                             )
             azimut = aa / cst.degree
+            print(azimut)
         else:
             raise ValueError("At least one of the azimut or the phase angle must be provided")