disort.h

Go to the documentation of this file.

 
#ifndef disort_h
#define disort_h
 
#include "agenda_class.h"
#include "matpack_data.h"
#include "mystring.h"
#include "optproperties.h"
#include "sun.h"
 
 
void add_normed_phase_functions(Tensor3View pftc1,
                                const MatrixView& sca1,
                                const MatrixView& pftc2,
                                const MatrixView& sca2);
 
void check_disort_input(  // Input
    const Index& cloudbox_on,
    const Index& atmfields_checked,
    const Index& atmgeom_checked,
    const Index& cloudbox_checked,
    const Index& scat_data_checked,
    const Index& atmosphere_dim,
    const Index& stokes_dim,
    const ArrayOfIndex& cloudbox_limits,
    const ArrayOfArrayOfSingleScatteringData& scat_data,
    ConstVectorView za_grid,
    const Index& nstreams);
 
void check_disort_irradiance_input(  // Input
    const Index& atmfields_checked,
    const Index& atmgeom_checked,
    const Index& scat_data_checked,
    const Index& atmosphere_dim,
    const Index& stokes_dim,
    const ArrayOfArrayOfSingleScatteringData& scat_data,
    const Index& nstreams);
 
void init_ifield(  // Output
    Tensor7& cloudbox_field,
    // Input
    const Vector& f_grid,
    const ArrayOfIndex& cloudbox_limits,
    const Index& n_za,
    const Index& n_aa,
    const Index& stokes_dim);
 
void get_disortsurf_props(  // Output
    Vector& albedo,
    Numeric& btemp,
    // Input
    ConstVectorView f_grid,
    const Numeric& surface_skin_t,
    ConstVectorView surface_scalar_reflectivity);
 
void run_cdisort(Workspace& ws,
                 Tensor7& cloudbox_field,
                 ArrayOfMatrix& disort_aux,
                 ConstVectorView f_grid,
                 ConstVectorView p_grid,
                 ConstVectorView z_profile,
                 const Numeric& z_surface,
                 ConstVectorView t_profile,
                 ConstMatrixView vmr_profiles,
                 ConstMatrixView pnd_profiles,
                 const ArrayOfArrayOfSingleScatteringData& scat_data,
                 const ArrayOfSun& suns,
                 const Agenda& propmat_clearsky_agenda,
                 const Agenda& gas_scattering_agenda,
                 const ArrayOfIndex& cloudbox_limits,
                 const Numeric& surface_skin_t,
                 const Vector& surface_scalar_reflectivity,
                 ConstVectorView za_grid,
                 ConstVectorView aa_grid,
                 ConstVectorView sun_rte_los,
                 const Index& gas_scattering_do,
                 const Index& suns_do,
                 const ArrayOfString& disort_aux_vars,
                 const Numeric& scale_factor,
                 const Index& nstreams,
                 const Index& Npfct,
                 const Index& only_tro,
                 const Index& quiet,
                 const Index& emission,
                 const Index& intensity_correction,
                 const Verbosity& verbosity);
 
void run_cdisort_flux(Workspace& ws,
                      Tensor5& spectral_irradiance_field,
                      ArrayOfMatrix& disort_aux,
                      ConstVectorView f_grid,
                      ConstVectorView p_grid,
                      ConstVectorView z_profile,
                      const Numeric& z_surface,
                      ConstVectorView t_profile,
                      ConstMatrixView vmr_profiles,
                      ConstMatrixView pnd_profiles,
                      const ArrayOfArrayOfSingleScatteringData& scat_data,
                      const ArrayOfSun& suns,
                      const Agenda& propmat_clearsky_agenda,
                      const Agenda& gas_scattering_agenda,
                      const ArrayOfIndex& cloudbox_limits,
                      const Numeric& surface_skin_t,
                      const Vector& surface_scalar_reflectivity,
                      ConstVectorView sun_rte_los,
                      const Index& gas_scattering_do,
                      const Index& suns_do,
                      const ArrayOfString& disort_aux_vars,
                      const Numeric& scale_factor,
                      const Index& nstreams,
                      const Index& Npfct,
                      const Index& only_tro,
                      const Index& quiet,
                      const Index& emission,
                      const Index& intensity_correction,
                      const Verbosity& verbosity);
 
void get_gasoptprop(Workspace& ws,
                    MatrixView ext_bulk_gas,
                    const Agenda& propmat_clearsky_agenda,
                    ConstVectorView t_profile,
                    ConstMatrixView vmr_profiles,
                    ConstVectorView p_grid,
                    ConstVectorView f_grid);
 
void get_gas_scattering_properties(Workspace& ws,
                                   MatrixView sca_coeff_gas,
                                   MatrixView sca_coeff_gas_level,
                                   MatrixView pfct_gas,
                                   const ConstVectorView& f_grid,
                                   const VectorView& p,
                                   const VectorView& t,
                                   const MatrixView& vmr,
                                   const Agenda& gas_scattering_agenda);
 
void get_paroptprop(MatrixView ext_bulk_par,
                    MatrixView abs_bulk_par,
                    const ArrayOfArrayOfSingleScatteringData& scat_data,
                    ConstMatrixView pnd_profiles,
                    ConstVectorView t_profile,
                    ConstVectorView DEBUG_ONLY(p_grid),
                    const ArrayOfIndex& cloudbox_limits,
                    const Index f_index);
 
void get_dtauc_ssalb(MatrixView dtauc,
                     MatrixView ssalb,
                     ConstMatrixView ext_bulk_gas,
                     ConstMatrixView ext_bulk_par,
                     ConstMatrixView abs_bulk_par,
                     ConstVectorView z_profile);
 
void get_angs(Vector& pfct_angs,
              const ArrayOfArrayOfSingleScatteringData& scat_data,
              const Index& Npfct);
 
void get_parZ(Tensor3& pha_bulk_par,
              const ArrayOfArrayOfSingleScatteringData& scat_data,
              ConstMatrixView pnd_profiles,
              ConstVectorView t_profile,
              ConstVectorView pfct_angs,
              const ArrayOfIndex& cloudbox_limits,
              const Index f_index);
 
void get_pfct(Tensor3& pfct_bulk_par,
              ConstTensor3View pha_bulk_par,
              ConstMatrixView ext_bulk_par,
              ConstMatrixView abs_bulk_par,
              const ArrayOfIndex& cloudbox_limits);
 
void get_pmom(Tensor3View pmom,
              ConstTensor3View pfct_bulk_par,
              ConstVectorView pfct_angs,
              const Index& Nlegendre);
 
void get_scat_bulk_layer(MatrixView sca_bulk_layer,
                         const MatrixView& ext_bulk,
                         const MatrixView& abs_bulk);
 
void reduced_1datm(Vector& p,
                   Vector& z,
                   Vector& t,
                   Matrix& vmr,
                   Matrix& pnd,
                   ArrayOfIndex& cboxlims,
                   Index& ncboxremoved,
                   ConstVectorView p_grid,
                   ConstVectorView z_profile,
                   const Numeric& z_surface,
                   ConstVectorView t_profile,
                   ConstMatrixView vmr_profiles,
                   ConstMatrixView pnd_profiles,
                   const ArrayOfIndex& cloudbox_limits);
 
void surf_albedoCalc(Workspace& ws,
                     //Output
                     VectorView albedo,
                     Numeric& btemp,
                     //Input
                     const Agenda& surface_rtprop_agenda,
                     ConstVectorView f_grid,
                     ConstVectorView scat_za_grid,
                     const Numeric& surf_alt,
                     const Verbosity& verbosity);
 
void surf_albedoCalcSingleAngle(Workspace& ws,
                                //Output
                                VectorView albedo,
                                Numeric& btemp,
                                //Input
                                const Agenda& surface_rtprop_agenda,
                                ConstVectorView f_grid,
                                const Numeric& surf_alt,
                                const Numeric& inc_angle);
 
#endif /* disort_h */