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 :
- transformation de l’humidité spécifique de DRIAS ou SAFRAN donnée en g/kg en g/g $$QS=\frac{QS}{1000}$$
- calcul du rapport de mélange RM $$RM=\frac{QS}{1-QS}$$
- conversion des températures moyennes de °C en Kelvin, avec :
- 273,15°C = référence à la température absolue $$TM2K=TM2+273.15$$
- 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}$$
- 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$$
- 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})$$
- expression de la pression de vapeur saturante (en hPa) $$Es=Es/100$$
- calcul de la pression de vapeur (en hPa) $$Ea=(HR/100)*Es$$
- calcul du déficit de saturation de l’air (en hPa) $$DSAT=Es-Ea$$