Format Biljou

Un fichier de données météo pour le modèle de bilan hydrique forestier Biljou est un fichier texte comportant les données d’une maille. Les données sont séparées par un point-virgule, un blanc ou une tabulation. Le séparateur utilisé par SICLIMA est la tabulation.

Conversions des variables climatiques

Le modèle DRIAS/IPSL ne fournissant pas de rayonnement global, le format Biljou n’est pas disponible.

Colonne du fichier Description Format Formule pour ALADIN Formule pour SAFRAN
1 annee Année entier, entre 1000 et 9999
2 jour Jour de l’année entier, entre 1 et 366
3 vent Vent à 2 m [m/s] réel sfcwind * 0.7 $$ff_q * \frac{4.87}{log(10 * 67.8 - 5.42)}$$
4 pluie Pluie [mm] réel, entre 0 et 300 rstr + prsnls preliq_q + prenei_q
5 tsec Température moyenne de l’air [°C] réel, entre 0 et 300 tmean t_q
6 dsat Déficit de saturation de l’air [hPa] réel, entre 0 et 50 voir ci-dessous voir ci-dessous
7 rgl Rayonnement global [J/cm²] réel, entre 0 et 3500 $$rsds * \frac{60 * 60 * 24}{10000}$$ ssi_q

Calcul du déficit de saturation de l’air

Le calcul de dsat (Es - Ea) procède ainsi :

  1. transformation de l’humidité spécifique de DRIAS ou SAFRAN donnée en g/kg en g/g $$QS=\frac{QS}{1000}$$
  2. calcul du rapport de mélange RM $$RM=\frac{QS}{1-QS}$$
  3. conversion des températures moyennes de °C en Kelvin, avec :
    • 273,15°C = référence à la température absolue $$TM2K=TM2+273.15$$
  4. calcul de la pression atmosphérique (en Pa) à l’altitude de la maille SAFRAN selon la formule de Nadeau & Puiggali / Jean-Christophe Calvet / Mathieu Regimbeau avec :
    • 101325 = pression standard au niveau de la mer (Pa)
    • 288.15 K = 15°C = température de référence au niveau de la mer, Jean-Christophe Calvet utilise la température de référence de 293 : pas d’incidence sur HR
    • ALT, l’altitude de la maille SAFRAN $$PA=101325*(\frac{288.15-(0.0065*ALT)}{288.15})^{5.257}$$
  5. calcul de la pression de vapeur saturante selon Goff & Gratch 1957 (en Pa) préconisation WMO-2000 dans sa formulation appliquée par Météo-France, avec :
    • 273,16°C = référence au point triple où la vapeur, l’eau liquide et la glace peuvent coexister. $$Es=10^{10.79574*(1-\frac{273.16}{TM2K}) -5.028 * {\log10 \frac{TM2K}{273.16}} + 0.000150475 * (1-10^{-8.2969*\frac{TM2K}{273.16-1}}) + 0.00042873 * (10^{4.76955 * (1-273.16/TM2K)}-1) + 0.78614}*100$$
  6. calcul de l’humidité relative (en %) avec :
    • Masse molaire de l’air sec Ma=28,9644 g/mol
    • Masse molaire de la vapeur d’eau Mw=18,01534 g/mol
    • Epsilon = Mw/Ma = 0,621982157 = 0,622 $$HR=MIN(100,\frac{100 * RM * PA}{Es * 0.622+ Es * RM})$$
  7. expression de la pression de vapeur saturante (en hPa) $$Es=Es/100$$
  8. calcul de la pression de vapeur (en hPa) $$Ea=(HR/100)*Es$$
  9. calcul du déficit de saturation de l’air (en hPa) $$DSAT=Es-Ea$$