<?xml version="1.0" encoding="utf-8"?><Code_Saturne_GUI case="case1" solver_version="6.0" study="test_atmospheric_dispersion" version="2.0">
  <additional_scalars>
    <users/>
    <variable label="concentration1" name="concentration1" type="user">
      <blending_factor>1</blending_factor>
      <max_value>1e+12</max_value>
      <min_value>-1e+12</min_value>
      <property choice="constant" label="Dscal1" name="concentration1_diffusivity">
        <initial_value>1.83e-05</initial_value>
      </property>
      <rhs_reconstruction>1</rhs_reconstruction>
      <turbulent_flux_model>SGDH</turbulent_flux_model>
    </variable>
  </additional_scalars>
  <analysis_control>
    <output>
      <listing_printing_frequency>1</listing_printing_frequency>
      <mesh id="-1" label="Fluid domain" type="cells">
        <all_variables status="on"/>
        <location>all[]</location>
        <writer id="-1"/>
      </mesh>
      <mesh id="-2" label="Boundary" type="boundary_faces">
        <all_variables status="on"/>
        <location>all[]</location>
        <writer id="-1"/>
      </mesh>
      <probe name="1" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>34</probe_z>
      </probe>
      <probe name="2" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>84</probe_z>
      </probe>
      <probe name="3" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>134</probe_z>
      </probe>
      <probe name="4" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>184</probe_z>
      </probe>
      <probe name="5" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>234</probe_z>
      </probe>
      <probe name="6" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>284</probe_z>
      </probe>
      <probe name="7" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>334</probe_z>
      </probe>
      <probe name="8" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>384</probe_z>
      </probe>
      <probe name="9" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>434</probe_z>
      </probe>
      <probe name="10" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>484</probe_z>
      </probe>
      <probe name="11" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>534</probe_z>
      </probe>
      <probe name="12" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>584</probe_z>
      </probe>
      <probe name="13" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>634</probe_z>
      </probe>
      <probe name="14" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>684</probe_z>
      </probe>
      <probe name="15" status="on">
        <probe_x>193.341</probe_x>
        <probe_y>617.53</probe_y>
        <probe_z>734</probe_z>
      </probe>
      <probe name="16" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>1</probe_z>
      </probe>
      <probe name="17" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>51</probe_z>
      </probe>
      <probe name="18" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>101</probe_z>
      </probe>
      <probe name="19" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>151</probe_z>
      </probe>
      <probe name="20" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>201</probe_z>
      </probe>
      <probe name="21" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>251</probe_z>
      </probe>
      <probe name="22" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>301</probe_z>
      </probe>
      <probe name="23" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>351</probe_z>
      </probe>
      <probe name="24" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>401</probe_z>
      </probe>
      <probe name="25" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>451</probe_z>
      </probe>
      <probe name="26" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>501</probe_z>
      </probe>
      <probe name="27" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>551</probe_z>
      </probe>
      <probe name="28" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>601</probe_z>
      </probe>
      <probe name="29" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>651</probe_z>
      </probe>
      <probe name="30" status="on">
        <probe_x>420</probe_x>
        <probe_y>650</probe_y>
        <probe_z>701</probe_z>
      </probe>
      <probe name="31" status="on">
        <probe_x>255</probe_x>
        <probe_y>649</probe_y>
        <probe_z>54.35</probe_z>
      </probe>
      <probe name="32" status="on">
        <probe_x>1550</probe_x>
        <probe_y>641.5</probe_y>
        <probe_z>122</probe_z>
      </probe>
      <probe_format choice="CSV"/>
      <probe_recording_frequency>1</probe_recording_frequency>
      <writer id="-1" label="results">
        <directory name="postprocessing"/>
        <format name="ensight" options="separate_meshes"/>
        <frequency period="time_step">100</frequency>
        <output_at_end status="on"/>
        <time_dependency choice="fixed_mesh"/>
      </writer>
    </output>
    <profiles/>
    <scalar_balances/>
    <time_averages/>
    <time_parameters>
      <iterations>500</iterations>
      <max_courant_num>1.2</max_courant_num>
      <max_fourier_num>5</max_fourier_num>
      <property label="CourantNb" name="courant_number">
        <postprocessing_recording status="off"/>
      </property>
      <property label="FourierNb" name="fourier_number">
        <postprocessing_recording status="off"/>
      </property>
      <property label="LocalTime" name="local_time_step">
        <postprocessing_recording status="off"/>
      </property>
      <relaxation_coefficient>0.3</relaxation_coefficient>
      <thermal_time_step status="off"/>
      <time_passing>2</time_passing>
      <time_step_max_factor>1000</time_step_max_factor>
      <time_step_min_factor>0.1</time_step_min_factor>
      <time_step_ref>0.05</time_step_ref>
      <time_step_var>0.005</time_step_var>
    </time_parameters>
  </analysis_control>
  <boundary_conditions>
    <boundary label="BC_1" name="1" nature="inlet">inlet</boundary>
    <boundary label="BC_2" name="2" nature="inlet">outlet</boundary>
    <boundary label="BC_3" name="3" nature="inlet">source_polluants</boundary>
    <boundary label="BC_4" name="4" nature="inlet">symmetry</boundary>
    <boundary label="BC_5" name="5" nature="wall">w_bottom</boundary>
    <inlet field_id="none" label="BC_4">
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="concentration1" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <meteo_automatic status="off"/>
        <meteo_data status="on"/>
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="BC_2">
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="concentration1" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <meteo_automatic status="on"/>
        <meteo_data status="on"/>
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="BC_1">
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>0</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="concentration1" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <meteo_automatic status="on"/>
        <meteo_data status="on"/>
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="None">
      <scalar choice="dirichlet" label="concentration1" name="concentration1" type="user">
        <dirichlet>0</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <inlet field_id="none" label="BC_3">
      <convective_inlet status="on"/>
      <scalar choice="dirichlet" name="concentration1" type="user">
        <dirichlet>0.001</dirichlet>
      </scalar>
      <scalar choice="dirichlet" name="temperature" type="thermal">
        <dirichlet>293.15</dirichlet>
      </scalar>
      <turbulence choice="hydraulic_diameter">
        <hydraulic_diameter>1</hydraulic_diameter>
      </turbulence>
      <velocity_pressure choice="norm" direction="normal">
        <meteo_automatic status="off"/>
        <meteo_data status="off"/>
        <norm>1</norm>
      </velocity_pressure>
    </inlet>
    <wall field_id="none" label="BC_5">
      <scalar choice="neumann_formula" name="temperature" type="thermal">
        <neumann_formula>flux = 0.0;</neumann_formula>
      </scalar>
      <scalar choice="neumann" name="concentration1" type="user">
        <neumann>0</neumann>
      </scalar>
      <velocity_pressure choice="off">
        <dirichlet component="0" name="velocity">0</dirichlet>
        <dirichlet component="1" name="velocity">0</dirichlet>
        <dirichlet component="2" name="velocity">0</dirichlet>
        <roughness>0.05</roughness>
      </velocity_pressure>
    </wall>
  </boundary_conditions>
  <calculation_management>
    <block_io/>
    <partitioning/>
    <start_restart>
      <frozen_field status="off"/>
      <restart_rescue>0</restart_rescue>
      <restart_with_auxiliary status="on"/>
    </start_restart>
  </calculation_management>
  <lagrangian model="off"/>
  <numerical_parameters>
    <gradient_reconstruction choice="2"/>
    <gradient_transposed status="on"/>
    <hydrostatic_pressure status="on"/>
    <pressure_relaxation>1</pressure_relaxation>
    <velocity_pressure_algo choice="simplec">
      <piso_sweep_number>1</piso_sweep_number>
    </velocity_pressure_algo>
    <velocity_pressure_coupling status="off"/>
    <wall_pressure_extrapolation>1</wall_pressure_extrapolation>
  </numerical_parameters>
  <physical_properties>
    <fluid_properties>
      <material choice="user_material"/>
      <method choice="user_properties"/>
      <property choice="predefined_law" label="Density" name="density">
        <initial_value>1.17862</initial_value>
      </property>
      <property choice="constant" label="DiffDyn" name="dynamic_diffusion">
        <initial_value>0.01</initial_value>
        <listing_printing status="off"/>
        <postprocessing_recording status="off"/>
      </property>
      <property choice="user_law" label="LamVisc" name="molecular_viscosity">
        <formula># Utilisation de la loi de Sutherland
CST = 111;
T0 = 273;
mu_ref = 17.6e-6;

if ( temperature &gt; 0 &amp;&amp; temperature &lt; 555) {
molecular_viscosity=mu_ref * ((temperature/T0)^(3./2.))*(T0+CST)/(temperature+CST);
} else {
molecular_viscosity = -999.0;
}

#molecular_viscosity = 1.83e-05;</formula>
        <initial_value>1.83e-05</initial_value>
      </property>
      <property choice="user_law" label="SpecHeat" name="specific_heat">
        <formula># Loi polynomiale en temperature
specific_heat = 1057.5-0.4489*temperature+0.0011407*temperature^2. -0.0000007999*temperature^3.+0.00000000019327*temperature^4.;

#specific_heat = 1017.24;</formula>
        <initial_value>1017.24</initial_value>
      </property>
      <property choice="user_law" label="ThermalCond" name="thermal_conductivity">
        <formula># Utilisation de la loi de Sutherland
CST = 114;
T0 = 273;
lambda_ref = 0.0241;

if ( temperature &gt; 0 &amp;&amp; temperature &lt; 555) {
thermal_conductivity=lambda_ref * ((temperature/T0)^(3./2.))*(T0+CST)/(temperature+CST);
} else {
thermal_conductivity = -999.0;
}

#thermal_conductivity = 0.02495;</formula>
        <initial_value>0.02495</initial_value>
      </property>
      <reference_pressure>101325</reference_pressure>
      <reference_temperature>293.15</reference_temperature>
    </fluid_properties>
    <gravity>
      <gravity_x>0</gravity_x>
      <gravity_y>0</gravity_y>
      <gravity_z>-9.81</gravity_z>
    </gravity>
    <notebook/>
    <omega>
      <omega_x>0</omega_x>
      <omega_y>0</omega_y>
      <omega_z>0</omega_z>
    </omega>
  </physical_properties>
  <solution_domain>
    <extrusion/>
    <faces_cutting status="off"/>
    <joining/>
    <mesh_smoothing status="off"/>
    <meshes_list>
      <mesh name="Mesh_1.med"/>
    </meshes_list>
    <periodicity/>
    <thin_walls/>
    <volumic_conditions>
      <zone groundwater_law="off" head_losses="off" id="1" initialization="on" label="all_cells" mass_source_term="off" momentum_source_term="off" porosity="off" scalar_source_term="off" thermal_source_term="off">all[]</zone>
    </volumic_conditions>
  </solution_domain>
  <thermophysical_models>
    <ale_method/>
    <atmospheric_flows model="dry">
      <meteo_data>meteo</meteo_data>
      <property label="RealTemp" name="real_temperature"/>
      <read_meteo_data status="on"/>
    </atmospheric_flows>
    <compressible_model model="off"/>
    <conjugate_heat_transfer>
      <external_coupling/>
    </conjugate_heat_transfer>
    <gas_combustion model="off"/>
    <groundwater_model model="off"/>
    <hgn_model model="off"/>
    <joule_effect model="off"/>
    <porosities/>
    <radiative_transfer model="off">
      <absorption_coefficient type="constant">0</absorption_coefficient>
      <restart status="on"/>
    </radiative_transfer>
    <reference_values>
      <length/>
    </reference_values>
    <solid_fuels model="off"/>
    <source_terms/>
    <thermal_scalar model="potential_temperature">
      <property label="Dimensionless Thermal flux" name="boundary_layer_nusselt" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Boundary temperature" name="boundary_temperature" support="boundary"/>
      <property label="Thermal flux" name="thermal_flux" support="boundary"/>
      <property label="Tplus" name="tplus" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <variable label="PotTemp" name="temperature" type="thermal">
        <blending_factor>1</blending_factor>
        <formula zone_id="1">temperature = 300.;</formula>
        <max_value>1e+12</max_value>
        <min_value>0</min_value>
        <rhs_reconstruction>1</rhs_reconstruction>
        <turbulent_flux_model>SGDH</turbulent_flux_model>
      </variable>
    </thermal_scalar>
    <turbomachinery model="off">
      <joining/>
    </turbomachinery>
    <turbulence model="k-epsilon-PL">
      <gravity_terms status="on"/>
      <initialization choice="reference_value" zone_id="1"/>
      <property label="TurbVisc" name="turbulent_viscosity"/>
      <reference_velocity>5</reference_velocity>
      <variable label="epsilon" name="epsilon">
        <blending_factor>1</blending_factor>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <variable label="k" name="k">
        <blending_factor>1</blending_factor>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
      <wall_function>3</wall_function>
    </turbulence>
    <velocity_pressure>
      <initialization>
        <formula zone_id="1">velocity[0] = 0.;
velocity[1] = 0.;
velocity[2] = 0.;</formula>
      </initialization>
      <property label="Stress" name="stress" support="boundary"/>
      <property label="Stress, normal" name="stress_normal" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="Stress, tangential" name="stress_tangential" support="boundary">
        <postprocessing_recording status="off"/>
      </property>
      <property label="total_pressure" name="total_pressure"/>
      <property label="Yplus" name="yplus" support="boundary"/>
      <variable label="Pressure" name="pressure">
        <rhs_reconstruction>2</rhs_reconstruction>
      </variable>
      <variable dimension="3" label="Velocity" name="velocity">
        <blending_factor>1</blending_factor>
        <rhs_reconstruction>1</rhs_reconstruction>
      </variable>
    </velocity_pressure>
  </thermophysical_models>
</Code_Saturne_GUI>