#include "cs_defs.h"
#include <assert.h>
#include <errno.h>
#include <ctype.h>
#include <stdio.h>
#include <stdarg.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <float.h>
#include <ple_locator.h>
#include "bft_mem.h"
#include "bft_error.h"
#include "bft_printf.h"
#include "fvm_nodal_extract.h"
#include "cs_atmo_profile_std.h"
#include "cs_base.h"
#include "cs_boundary_conditions.h"
#include "cs_boundary_zone.h"
#include "cs_domain.h"
#include "cs_field.h"
#include "cs_field_default.h"
#include "cs_field_pointer.h"
#include "cs_halo.h"
#include "cs_halo_perio.h"
#include "cs_log.h"
#include "cs_math.h"
#include "cs_mesh.h"
#include "cs_mesh_location.h"
#include "cs_mesh_quantities.h"
#include "cs_parall.h"
#include "cs_equation_iterative_solve.h"
#include "cs_physical_constants.h"
#include "cs_physical_model.h"
#include "cs_prototypes.h"
#include "cs_post.h"
#include "cs_restart.h"
#include "cs_selector.h"
#include "cs_thermal_model.h"
#include "cs_turbulence_model.h"
#include "cs_volume_zone.h"
#include "cs_balance.h"
#include "cs_blas.h"
#include "cs_convection_diffusion.h"
#include "cs_gradient.h"
#include "cs_parameters.h"
#include "cs_porous_model.h"
#include "cs_timer.h"
#include "cs_matrix_building.h"
#include "cs_matrix_default.h"
#include "cs_sles.h"
#include "cs_sles_default.h"
#include "cs_face_viscosity.h"
#include "cs_divergence.h"
#include "cs_velocity_pressure.h"
#include "cs_atmo.h"
#include "cs_atmo_aerosol.h"

Functions
void	cs_hydrostatic_pressure_atmo_compute (cs_real_3_t f_ext, cs_real_3_t dfext, int f_id, cs_real_t i_massflux, cs_real_t b_massflux, cs_real_t i_viscm[], cs_real_t b_viscm[], cs_real_t dam, cs_real_t xam, cs_real_t dpvar, cs_real_t rhs)
	This auxiliary function solve a Poisson equation for hydrostatic pressure : $\divs ( \grad P ) = \divs ( f )$ . More...

void	cs_atmo_init_meteo_profiles (void)
	Initialize meteo profiles if no meteo file is given. More...

void	cs_atmo_compute_meteo_profiles (void)
	Compute meteo profiles if no meteo file is given. More...

void	cs_atmo_z_ground_compute (void)
	Compute the gruond elevation. More...

void	cs_atmo_hydrostatic_profiles_compute (void)
	Compute hydrostatic profiles of density and pressure. More...

void	cs_atmo_set_meteo_file_name (const char *file_name)
	This function set the file name of the meteo file. More...

void	cs_atmo_chemistry_set_spack_file_name (const char *file_name)
	This function set the file name of the SPACK file. More...

void	cs_atmo_chemistry_set_aerosol_file_name (const char *file_name)
	This function sets the file name to initialize the aerosol library. More...

void	cs_atmo_declare_chem_from_spack (void)
	This function declares additional transported variables for atmospheric module for the chemistry defined from SPACK. More...

void	cs_atmo_compute_solar_angles (cs_real_t latitude, cs_real_t longitude, cs_real_t squant, cs_real_t utc, int sea_id, cs_real_t albedo, cs_real_t muzero, cs_real_t omega, cs_real_t fo)
	1D Radiative scheme - Solar data + zenithal angle) More...

void	cs_atmo_log_setup (void)
	Print the atmospheric module options to setup.log. More...

void	cs_atmo_chemistry_log_setup (void)
	Print the atmospheric chemistry options to setup.log. More...

void	cs_atmo_aerosol_log_setup (void)
	Print the atmospheric aerosols options to setup.log. More...

void	cs_atmo_finalize (void)
	Deallocate arrays for atmo module. More...

Function Documentation

◆ cs_atmo_aerosol_log_setup()

void cs_atmo_aerosol_log_setup ( void )

Print the atmospheric aerosols options to setup.log.

◆ cs_atmo_chemistry_log_setup()

void cs_atmo_chemistry_log_setup ( void )

Print the atmospheric chemistry options to setup.log.

◆ cs_atmo_chemistry_set_aerosol_file_name()

void cs_atmo_chemistry_set_aerosol_file_name ( const char * file_name )

This function sets the file name to initialize the aerosol library.

Parameters

[in] file_name name of the file.

◆ cs_atmo_chemistry_set_spack_file_name()

void cs_atmo_chemistry_set_spack_file_name ( const char * file_name )

This function set the file name of the SPACK file.

Parameters

[in] file_name name of the file.

◆ cs_atmo_compute_meteo_profiles()

void cs_atmo_compute_meteo_profiles ( void )

Compute meteo profiles if no meteo file is given.

◆ cs_atmo_compute_solar_angles()

void cs_atmo_compute_solar_angles	(	cs_real_t	latitude,
		cs_real_t	longitude,
		cs_real_t	squant,
		cs_real_t	utc,
		int	sea_id,
		cs_real_t *	albedo,
		cs_real_t *	muzero,
		cs_real_t *	omega,
		cs_real_t *	fo
	)

1D Radiative scheme - Solar data + zenithal angle)

Compute:

zenithal angle
solar contant (with correction for earth - solar length)

albedo if above the sea (Use analytical formulae of Paltrige and Platt dev.in atm. science no 5)

Parameters

[in]	latitude	latitude
[in]	longitude	longitude
[in]	squant	start day in the year
[in]	utc	Universal time (hour)
[in]	sea_id	sea index
[out]	albedo	albedo
[out]	muzero	cosin of zenithal angle
[out]	omega	solar azimut angle
[out]	fo	solar constant

◆ cs_atmo_declare_chem_from_spack()

void cs_atmo_declare_chem_from_spack ( void )

This function declares additional transported variables for atmospheric module for the chemistry defined from SPACK.

This function declare additional transported variables for atmospheric module for the chemistry defined from SPACK.

◆ cs_atmo_finalize()

void cs_atmo_finalize ( void )

Deallocate arrays for atmo module.

◆ cs_atmo_hydrostatic_profiles_compute()

void cs_atmo_hydrostatic_profiles_compute ( void )

Compute hydrostatic profiles of density and pressure.

This function solves the following transport equation on $\varia$ :

$\divs \left( \grad \varia \right) = \divs \left( \dfrac{\vect{g}}{c_p \theta} \right)$

where $\vect{g}$ is the gravity field and $\theta$ is the potential temperature.

The boundary conditions on $\varia$ read:

$\varia = \left(\dfrac{P_{sea}}{p_s}\right)^{R/C_p} \textrm{on the ground}$

and Neumann elsewhere.

◆ cs_atmo_init_meteo_profiles()

void cs_atmo_init_meteo_profiles ( void )

Initialize meteo profiles if no meteo file is given.

(DOXYGEN_SHOULD_SKIP_THIS)

◆ cs_atmo_log_setup()

void cs_atmo_log_setup ( void )

Print the atmospheric module options to setup.log.

◆ cs_atmo_set_meteo_file_name()

void cs_atmo_set_meteo_file_name ( const char * file_name )

This function set the file name of the meteo file.

Parameters

[in] file_name name of the file.

◆ cs_atmo_z_ground_compute()

void cs_atmo_z_ground_compute ( void )

Compute the gruond elevation.

This function computes the ground elevation.

This function solves the following transport equation on $\varia$ :

$\dfrac{\partial \varia}{\partial t} + \divs \left( \varia \vect{g} \right) - \divs \left( \vect{V} \right) \varia = 0$

where $\vect{g}$ is the gravity field

The boundary conditions on $\varia$ read:

$\varia = z \textrm{ on walls}$

$\dfrac{\partial \varia}{\partial n} = 0 \textrm{ elsewhere}$

Remarks:

a steady state is looked for.

◆ cs_hydrostatic_pressure_atmo_compute()

void cs_hydrostatic_pressure_atmo_compute	(	cs_real_3_t *	f_ext,
		cs_real_3_t *	dfext,
		int	f_id,
		cs_real_t *	i_massflux,
		cs_real_t *	b_massflux,
		cs_real_t	i_viscm[],
		cs_real_t	b_viscm[],
		cs_real_t *	dam,
		cs_real_t *	xam,
		cs_real_t *	dpvar,
		cs_real_t *	rhs
	)

This auxiliary function solve a Poisson equation for hydrostatic pressure : $\divs ( \grad P ) = \divs ( f )$ .

Parameters

[in]	f_ext	Total (NOT incremental) external forcing
[in]	dfext	external forcing increment
[in,out]	i_massflux	Internal mass flux
[in,out]	b_massflux	Boundary mass flux
[in]	f_id	Field id from the field to update
[in,out]	i_viscm	Internal face viscosity
[in,out]	i_viscm	Boundary face viscosity
[in,out]	dam	Working array
[in,out]	xam	Working array
[in,out]	dpvar	Pressure increment at each sweep
[in,out]	rhs	Working array

Functions

Function Documentation

◆ cs_atmo_aerosol_log_setup()

◆ cs_atmo_chemistry_log_setup()

◆ cs_atmo_chemistry_set_aerosol_file_name()

◆ cs_atmo_chemistry_set_spack_file_name()

◆ cs_atmo_compute_meteo_profiles()

◆ cs_atmo_compute_solar_angles()

◆ cs_atmo_declare_chem_from_spack()

◆ cs_atmo_finalize()

◆ cs_atmo_hydrostatic_profiles_compute()

◆ cs_atmo_init_meteo_profiles()

◆ cs_atmo_log_setup()

◆ cs_atmo_set_meteo_file_name()

◆ cs_atmo_z_ground_compute()

◆ cs_hydrostatic_pressure_atmo_compute()