| 2 | |
| 3 | = Update (March 22nd 2016) = |
| 4 | An implicit snow scheme was introduced by Fuxing Wang and Tao Wang in revision [3059]. However there is a problem with the snowmelt, which is much too low in this revision.[[BR]] |
| 5 | The problem was analysed by Catherine and Fabienne, working with Sylvie Charbit and Christophe Dumas from LSCE, who are applying the explicit snow scheme on Greenland.[[BR]] |
| 6 | A first commit was done in revision [3269] to fix the snowmelt problem.[[BR]] |
| 7 | |
| 8 | Some explanations (soon in English):[[BR]] |
| 9 | '''Révision du modèle de neige dans Orchidee'''[[BR]] |
| 10 | Les modifications effectuées par Tao et Fuxing pour rendre le schéma de neige implicite, ont consisté en : |
| 11 | - déplacement des routines de calcul de la température des 3 couches de neige (explicit-snow-profile) et des coefficients thermiques (explicit-snow-coef) dans thermosoil |
| 12 | - non prise en compte de temp_sol_add (excès de température disponible pour la fonte de la neige à l'interface avec l'atmosphère, calculée à partir du bilan d'énergie de surface) dans le calcul de snow-temp |
| 13 | - remise à la température de fusion (0°C) de temp_sol_new (température de surface issue du bilan d'énergie recalculée par enerbil_fusion), en présence de neige après le calcul de la fonte et avant le calcul du profil de température dans la neige et le sol |
| 14 | Ces modifications entrainaient que la neige ne pouvait jamais fondre puisque temp_sol_add n’était pas considérée dans le calcul du snowmelt. Par ailleurs, temp_sol_new étant toujours maintenue à 0 après la fonte de la neige et avant thermosoil, la température de la neige et du sol (en présence de neige) était calculée sans jamais prendre en compte l'énergie disponible à la surface.[[BR]] |
| 15 | Pour corriger le schéma, on propose : |
| 16 | - de remettre le calcul du profil de température dans la neige dans explicit-snow en prenant en compte temp_sol_add pour pouvoir calculer correctement le snowmelt mais de laisser le calcul des coefficients thermiques (pour le prochain pas de temps) dans thermosoil pour qu'ils soient bien calculés avec les températures du sol et de la neige calculées au pas de temps t |
| 17 | - de transférer les lignes de code initialement dans enerbil_fusion qui ramènent la température de surface de la neige (temp_sol_new) à la température de fusion (0°C), après le calcul des températures de la neige et du sol, c’est-à-dire à la fin de thermosoil. |
| 18 | Ainsi la résolution des températures dans la neige et le sol restent bien implicites. |
| 19 | |
| 20 | |
| 21 | |