cs_atmo.h

Go to the documentation of this file.

#ifndef __CS_ATMO_H__
#define __CS_ATMO_H__
 
/*============================================================================
 * Main for atmospheric related functions
 *============================================================================*/
 
/*
  This file is part of code_saturne, a general-purpose CFD tool.
 
  Copyright (C) 1998-2025 EDF S.A.
 
  This program is free software; you can redistribute it and/or modify it under
  the terms of the GNU General Public License as published by the Free Software
  Foundation; either version 2 of the License, or (at your option) any later
  version.
 
  This program is distributed in the hope that it will be useful, but WITHOUT
  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  details.
 
  You should have received a copy of the GNU General Public License along with
  this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
  Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
 
/*----------------------------------------------------------------------------*/
 
/*----------------------------------------------------------------------------
 *  Local headers
 *----------------------------------------------------------------------------*/
 
#include "base/cs_base.h"
 
/*----------------------------------------------------------------------------*/
 
BEGIN_C_DECLS
 
/*============================================================================
 * Local Macro definitions
 *============================================================================*/
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Atmospheric models
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_OFF = -1,
  CS_ATMO_CONSTANT_DENSITY = 0,
  CS_ATMO_DRY = 1,
  CS_ATMO_HUMID = 2,
 
} cs_atmo_model_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric nucleation models
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_NUC_OFF = 0,
  CS_ATMO_NUC_PRUPPACHER_KLETT = 1,
  CS_ATMO_NUC_COHARD = 2,
  CS_ATMO_NUC_ABDUL_RAZZAK = 3
 
} cs_atmo_nucleation_type_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric universal functions
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_UNIV_FN_CHENG = 0,
  CS_ATMO_UNIV_FN_HOGSTROM = 1,
  CS_ATMO_UNIV_FN_BUSINGER = 2,
  CS_ATMO_UNIV_FN_HARTOGENSIS = 3,
  CS_ATMO_UNIV_FN_CARL = 4,
 
} cs_atmo_universal_functions_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric ground model
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_GROUND_5_CAT = 0,
  CS_ATMO_GROUND_7_CAT = 1,
  CS_ATMO_GROUND_23_CAT = 2,
 
} cs_atmo_ground_cat_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric ground micro-scale options
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_GROUND_GENUINE = 0,
  CS_ATMO_GROUND_PHOTOVOLTAICS = 1,
  CS_ATMO_GROUND_VEGETATION = 2,
 
} cs_atmo_ground_meb_model_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric ground micro-scale options
 *----------------------------------------------------------------------------*/
 
typedef enum {
 
  CS_ATMO_PROJ_UNDEF = -1,
  CS_ATMO_PROJ_WGS84 = 0,
  /* Lambert93 projection type */
  CS_ATMO_PROJ_LAMBERT_93 = 1,
  /* Universal transverse mercator, zone  needs to be indicated*/
  CS_ATMO_PROJ_UTM = 2,
  /* Tangantial */
  CS_ATMO_PROJ_TAN = 3,
 
} cs_atmo_projection_t;
 
/*============================================================================
 * Type definitions
 *============================================================================*/
 
/*----------------------------------------------------------------------------
 * Atmospheric model options descriptor
 *----------------------------------------------------------------------------*/
 
typedef struct {
  /* Space and time reference of the run */
  int syear;
  int squant;
  int shour;
  int smin;
  cs_real_t ssec;
  cs_real_t longitude;
  cs_real_t latitude;
  cs_real_t altitude;
  cs_real_t x_l93; // conserved for retrocompatibility
  cs_real_t y_l93; //
  cs_real_t x_origin; // similar to l93 but more general
  cs_real_t y_origin; //
 
  cs_atmo_projection_t projection_type;
  int utm_zone; // If using utm, a zone needs to be specified
 
  union {
    int met_1d_nlevels_d;
    int nbmetd; /* deprecated */
  };
  union {
    int met_1d_nlevels_t;
    int nbmett; /* deprecated */
  };
  union {
    int met_1d_ntimes;
    int nbmetm;
  };
 
  union {
    int met_1d_nlevels_max_t;
    int nbmaxt;
  };
 
  cs_real_t domain_orientation;
 
  bool compute_z_ground;
 
  int open_bcs_treatment;
  int theo_interp;
 
  /* Model options */
  int sedimentation_model;
  int deposition_model;
  int nucleation_model;
  int subgrid_model;
  int distribution_model;
  int meteo_profile;
 
  char *meteo_file_name;
 
  cs_real_t meteo_dlmo;
  cs_real_t meteo_z0;
  cs_real_t meteo_zref;
  cs_real_t meteo_zi;
  cs_real_t meteo_zu1;
  cs_real_t meteo_zu2;
  cs_real_t meteo_zt1;
  cs_real_t meteo_zt2;
  cs_real_t meteo_uref;
  cs_real_t meteo_u1;
  cs_real_t meteo_u2;
  cs_real_t meteo_ustar0;
  cs_real_t meteo_wstar0;
  cs_real_t meteo_angle;
  cs_real_t meteo_t0;
  cs_real_t meteo_t1;
  cs_real_t meteo_t2;
  cs_real_t meteo_tstar;
  cs_real_t meteo_psea;
 
  cs_real_t meteo_qw0;
  cs_real_t meteo_qwstar;
  cs_real_t meteo_qw1;
  cs_real_t meteo_qw2;
  cs_real_t meteo_ql0;
  cs_real_t meteo_evapor;
  cs_real_t meteo_sensi;
  int meteo_phim_s;
  int meteo_phih_s;
  int meteo_phim_u;
  int meteo_phih_u;
 
  /* 1D meteo profiles */
  cs_real_t *xyp_met;
  cs_real_t *u_met;
  cs_real_t *v_met;
  cs_real_t *w_met;
  cs_real_t *ek_met;
  cs_real_t *ep_met;
  cs_real_t *temp_met;
  cs_real_t *rho_met;
  cs_real_t *qw_met;
  cs_real_t *ndrop_met;
  cs_real_t *z_dyn_met;
  cs_real_t *z_temp_met;
  cs_real_t *time_met;
  cs_real_t *hyd_p_met;
  cs_real_t *pot_t_met;
  cs_real_t *dpdt_met;
  cs_real_3_t *mom_met;
  cs_real_3_t *mom_cs;
  int hydrostatic_pressure_model;
  int qv_profile;
 
  int ground_model;
  cs_atmo_ground_cat_t ground_cat;
  int ground_zone_id;
  cs_atmo_ground_meb_model_t ground_meb_model;
 
  bool rain;
  int cloud_type;
  bool accretion;
  bool autoconversion;
  bool autocollection_cloud;
  bool autocollection_rain;
  bool precipitation;
  bool evaporation;
  bool rupture;
 
  cs_real_t ground_surf_temp;
  cs_real_t ground_temperature;
  cs_real_t ground_humidity;
  cs_real_t ground_w1_ini;
  cs_real_t ground_w2_ini;
  cs_real_t *ground_cat_thermal_inertia;
  cs_real_t *ground_cat_roughness;
  cs_real_t *ground_cat_thermal_roughness;
  cs_real_t *ground_cat_albedo;
  cs_real_t *ground_cat_emissi;
  cs_real_t *ground_cat_vegeta;
  cs_real_t *ground_cat_w1;
  cs_real_t *ground_cat_w2;
  cs_real_t *ground_cat_r1;
  cs_real_t *ground_cat_r2;
  cs_real_t sigc;
  int infrared_1D_profile;
  int solar_1D_profile;
  int profiles_grid_id;
 
  cs_real_t aod_o3_tot;
 
  cs_real_t aod_h2o_tot;
 
} cs_atmo_option_t;
 
/*----------------------------------------------------------------------------
 * Atmospheric model constants descriptor
 *----------------------------------------------------------------------------*/
 
typedef struct {
  /* Space and time reference of the run */
  cs_real_t ps;
 
} cs_atmo_constants_t;
 
/*============================================================================
 * Static global variables
 *============================================================================*/
 
/* Pointer to atmo options structure */
extern cs_atmo_option_t *cs_glob_atmo_option;
 
/* Pointer to atmo constants structure */
extern cs_atmo_constants_t *cs_glob_atmo_constants;
 
/*============================================================================
 * Public function definitions
 *============================================================================*/
 
void
cs_atmo_phyvar_update(void);
 
/*----------------------------------------------------------------------------
 * initialize fields, stage 0
 *----------------------------------------------------------------------------*/
 
void
cs_atmo_fields_init0(void);
 
/*----------------------------------------------------------------------------
 * Automatic boundary condition for cooling towers
 *----------------------------------------------------------------------------*/
 
void
cs_atmo_bcond(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_init_meteo_profiles(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_compute_meteo_profiles(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_z_ground_compute(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_hydrostatic_profiles_compute(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_read_meteo_profile(int mode);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_mo_phim(cs_real_t              z,
           cs_real_t              dlmo);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_mo_phih(cs_real_t              z,
           cs_real_t              dlmo,
           cs_real_t              prt);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_mo_psim(cs_real_t              z,
           cs_real_t              z0,
           cs_real_t              dlmo);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
cs_real_t
cs_mo_psih(cs_real_t              z,
           cs_real_t              z0,
           cs_real_t              dlmo,
           cs_real_t              prt);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_mo_compute_from_thermal_diff(cs_real_t   z,
                                cs_real_t   z0,
                                cs_real_t   du,
                                cs_real_t   buoyant_param,
                                cs_real_t   dt,
                                cs_real_t   *dlmo,
                                cs_real_t   *ustar);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_mo_compute_from_thermal_flux(cs_real_t   z,
                                cs_real_t   z0,
                                cs_real_t   du,
                                cs_real_t   buoyant_param,
                                cs_real_t   flux,
                                cs_real_t   *dlmo,
                                cs_real_t   *ustar);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_set_meteo_file_name(const char *file_name);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_compute_solar_angles(cs_real_t xlat,
                             cs_real_t xlong,
                             cs_real_t jour,
                             cs_real_t heurtu,
                             int       imer,
                             cs_real_t *albe,
                             cs_real_t *za,
                             cs_real_t *muzero,
                             cs_real_t *omega,
                             cs_real_t *fo);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_log_setup(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_finalize(void);
 
/*----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------*/
 
void
cs_user_ground_model(void);
 
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_get_ground_zone(cs_lnum_t         *n_elts,
                        int               *n_ground_cat,
                        const cs_lnum_t  **elt_ids);
 
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_ground_init_arrays(int        *n_ground_cat,
                           cs_real_t  **csol,
                           cs_real_t  **rugdyn,
                           cs_real_t  **rugthe,
                           cs_real_t  **albedo,
                           cs_real_t  **emissi,
                           cs_real_t  **vegeta,
                           cs_real_t  **c1w,
                           cs_real_t  **c2w,
                           cs_real_t  **r1,
                           cs_real_t  **r2);
 
/*----------------------------------------------------------------------------*/
 
void
cs_atmo_projection
(
  cs_atmo_projection_t origin_projection,
  cs_atmo_projection_t target_projection,
  cs_real_t            x_origin,
  cs_real_t            y_origin,
  cs_real_t*           x_target,
  cs_real_t*           y_target,
  int                  utm_fixed_zone
);
 
/*----------------------------------------------------------------------------*/
 
END_C_DECLS
 
#endif /* __CS_ATMO_H__ */