27 | | - Problème de SPLIT precisol_ns corrigé. |
28 | | - Quelques bugs à corriger quand on active le routage: |
29 | | * mettre getin au lieu de getin_p dans CALL getin_p('RIVER_DESC', river_file) et CALL getin_p('RIVER_DESC_FILE', river_file_name) |
| 27 | * Problème de SPLIT precisol_ns corrigé. |
| 28 | * Quelques bugs à corriger quand on active le routage: |
| 29 | * mettre getin au lieu de getin_p dans CALL getin_p('RIVER_DESC', river_file) et CALL getin_p('RIVER_DESC_FILE', river_file_name) |
| 30 | * 31/01/2012: initaliser à FALSE le logique err_basin_number |
| 31 | |
| 32 | --------------------------- |
| 33 | 01/02/2012 |
| 34 | Essai sur 15ans d'une simulation SECHIBA+routage (sans floodplains ni irrigation) |
| 35 | * Le modèle tourne bien, avec tous les flags de vérification activés (waterbal, flags pour vérifier la construction des bassins versants ...) |
| 36 | * Problème de sous-estimation de l'évapotranspiration totale. En ordre de grandeur, voici les valeurs en moyennes sur le bassin amazonien de l'ET et de ses composantes (en mm/an): |
| 37 | - evap: 500 |
| 38 | - interception loss: 350 |
| 39 | - transpir: 10 |
| 40 | - evapnu: 140 |
| 41 | * Pour comparaison, voici les valeurs des mêmes variables pour une simulation que j'avais faite avec "ma" version d'ORCHIDEE lors de mon postdoc au LOCEAN (c.a.d. version avant le mergehydro). Rq: même forçage (NCC) pour les 2 simuls: |
| 42 | - evap: 1030 |
| 43 | - interception loss: 440 |
| 44 | - transpir: 460 |
| 45 | - evapnu: 130 |
| 46 | * L'évap est donc sous-estimée dans la version actuelle du mergehydro (evap "observée" de l'ordre de 1200 mm/an) à cause d'une transpiration quasi-inexistante. |
| 47 | * Le problème se situe au niveau du calcul des vbeta2, vbeta23 et vbeta3. |
| 48 | |
| 49 | Pour aller plus loin, j'ai fait 2 simulations sur 5 jours au pas de temps de 30min. Elles sont toutes deux issues de la même version (mergehydro). Seul le calcul des vbetas diffèrent entre les 2 simuls. |
| 50 | |
| 51 | Résultats donnés à Manaus. Voir figure: |
| 52 | [[Image(https://forge.ipsl.jussieu.fr/orchidee/attachment/wiki/Branches/MergeHydro/diag_transpir.png)]] |
| 53 | |
| 54 | Légende: |
| 55 | * "Merghy svn": avec le calcul original des vbetas dans la version mergehydro |
| 56 | * "Merghy Tristan": avec le calcul des vbetas d'après Tristan d'Orgeval dans la version mergehydro |
| 57 | |
| 58 | A priori, les 2 calculs sont sensés représenter un même concept. Mais ils divergent à mon sens: |
| 59 | * "Merghy svn": à chaque fois que le flux intégré de l'interception est > à qsintveg, la transpiration est déclenchée. Mais ce cas arrive très rarement voire jamais dans cet exemple. Voir dernier graphe de cette figure: |
| 60 | [[Image(https://forge.ipsl.jussieu.fr/orchidee/attachment/wiki/Branches/MergeHydro/diag_transpir2.png)]] |
| 61 | Le beta de l'interception est calculé sur la partie mouillée de la feuille. |
| 62 | |
| 63 | * "Merghy Tristan": à chaque fois que le flux maximal au niveau de la feuille n'est pas satisfait par le flux d'interception, la transpiration est déclenchée (c'est bien ce que l'on voit quand on compare les courbes bleues des 2 premiers graphes). Le beta de l'interception est calculé sur toute la feuille. |