| 24 | | Je propose un plan en 4 étapes ou nous laissons quelques semaines entre chaque étape afin que tout le monde puisse faire des runs et des tests avec cette nouvelle version. Ceux ceux qui ont piloté les changement pourront écrire une documentation technique pendant ce temps. J'ai essayé aussi de séparer les changements structuraux des améliorations des paramétrisations actuelles. Donc voici mes 4 étapes : |
| 25 | | |
| 26 | | 1) implémentation de la discrétisation verticale commune entre T et W. L'important est de garder la profondeur en T et W comme des paramètres indépendants afin qu'on soit libre de choisir plus tard ce qui convient. |
| 27 | | |
| 28 | | 2) Enlever la neige de Thermosoil et la mettre en sandwich entre enerbil et thermosoil. Ceci est donc la partie thermodynamique de la neige que Thao vient d’implémenter dans MICT mais avec la l'interface neige/sol finalisée. |
| 29 | | |
| 30 | | Après ces deux changements on aura une version qui a la nouvelle structure et les nouvelles approches numériques. Les paramètre auront peu changés et donc à priori le comportement du modèle sera similaire au trunk. Un environnement idéal pour voir si tout marche bien techniquement. |
| 31 | | |
| 32 | | 3) On rajoute les nouvelles capacité et diffusion thermiques de Fuxing avec le transport de la chaleur par l'eau. |
| 33 | | |
| | 24 | Je propose un plan en 4 étapes ou nous laissons quelques semaines entre chaque étape afin que tout le monde puisse faire des runs et des tests avec cette nouvelle version. Ceux ceux qui ont piloté les changement pourront écrire une documentation technique pendant ce temps. J'ai essayé aussi de séparer les changements structuraux des améliorations des paramétrisations actuelles. Donc voici mes 4 étapes :[[BR]] |
| | 25 | 1) implémentation de la discrétisation verticale commune entre T et W. L'important est de garder la profondeur en T et W comme des paramètres indépendants afin qu'on soit libre de choisir plus tard ce qui convient.[[BR]] |
| | 26 | 2) Enlever la neige de Thermosoil et la mettre en sandwich entre enerbil et thermosoil. Ceci est donc la partie thermodynamique de la neige que Thao vient d’implémenter dans MICT mais avec la l'interface neige/sol finalisée.[[BR]] |
| | 27 | Après ces deux changements on aura une version qui a la nouvelle structure et les nouvelles approches numériques. Les paramètre auront peu changés et donc à priori le comportement du modèle sera similaire au trunk. Un environnement idéal pour voir si tout marche bien techniquement.[[BR]] |
| | 28 | 3) On rajoute les nouvelles capacité et diffusion thermiques de Fuxing avec le transport de la chaleur par l'eau.[[BR]] |
| | 91 | '''Jan Polcher, 26/01/15''' to Agnès 23/01/15 |
| | 92 | |
| | 93 | thank you for this analysis ... which I totally share. I will give you my feeling on a few of the points you raise : |
| | 94 | |
| | 95 | The challenges : |
| | 96 | - I have put the snow second because I find the interaction with Tao more difficult and thus I wanted to give us more time for it. |
| | 97 | |
| | 98 | The depth : |
| | 99 | This should be configurable by the user as different applications will have different needs. We should code things so that each one can test the configuration which suits him/her best. The limiting conditions are : |
| | 100 | - The top layers need to have the fine vertical resolution of the soil moisture. |
| | 101 | - The total depth for soil moisture should be equal or smaller than for temperature (The total depth of soil moisture and temperature need to be parameters in the run.def). |
| | 102 | - After that we can debate if some upper and lower bounds for the depth should be introduced. |
| | 103 | |
| | 104 | As for the technicalities : |
| | 105 | - soil moisture below its total depth will be take as constant and equal to that of the deepest layer. |
| | 106 | - The smoothing of the soil temperature profile at depth of less than 5cm needs to be done in thermosoil once the new profile is calculated. That should be enough to avoid numerical instabilities. |
| | 107 | |
| | 108 | For the snow, I am still waiting to see the equations to be used to coupled the bottom temperature of snow pack with the one of the soil. Once we have all the equations of the temperature diffusion (from the skin temperature to the bottom of the soil and through the snow pack) in implicit form, things should flow smoothly. |
| | 109 | |
| | 110 | '''Fuxing Wang, 26/01/15''' |
| | 111 | |
| | 112 | (1) The depth. This should be configurable by the user as different applications will have different needs. We should code things so that each one can test the configuration which suits him/her best. |
| | 113 | |
| | 114 | In the new discretization, |
| | 115 | the HYDROL_SOIL_DEPTH (=dpu_max), THERMOSOIL_NBLEV (=ngrnd) can be defined in run.def. |
| | 116 | I also added a flag (VERDIS_FLAG, can be defined in run.def ) for different discretizations: verdis_flag = 1, the standard discretization; verdis_flag = 2, the revised discretization. |
| | 117 | But 'nbdl' (Number of diagnostic levels) and 'nslm' (number of CWRR levels) in 'constantes_soil_var' (=11 in standard discretization) was not defined in run.def. |
| | 118 | Is this kind of configuration OK ? |
| | 119 | Do the variables 'nbdl' and 'nslm' also need to be defined in run.def ? Because they also change when using the revised discretization. |
| | 120 | |
| | 121 | (2) The technicalities : The smoothing of the soil temperature profile at depth of less than 5cm needs to be done in thermosoil once the new profile is calculated. |
| | 122 | |
| | 123 | In this case, we smooth the top 5 soil layers for soil temperature ? |
| | 124 | The bottom positions from top 4st to 6th layers are : 2.151cm, 4.497cm, 9.189cm. |
| | 125 | |
| | 126 | '''Frédérique Cheruy, 27/01/15''' |
| | 127 | |
| | 128 | ''About the optimal vertical discretization for W and T''[[BR]] |
| | 129 | 1- We still need to verify that the 8m deep soil, is a problem for the runoff and the drainage.[[BR]] |
| | 130 | 2-Afterwards, as already discussed last week, and reported in Agnes mail, we might want to |
| | 131 | have different depths for moisture and temperature but we want to have the dependance of the soil thermal properties with |
| | 132 | the soil texture (3 or 12?), and the same vertical discretization for moisture and temperature to avoid interpolations.[[BR]] |
| | 133 | I am not sure we want to externalise nslm or ngrnd, and let the user play with it. If it is an advanced user he can go into the code |
| | 134 | and change the value of the parameter, if not may be it's better he does not play with it at all.[[BR]] |
| | 135 | But we need to test upper and lower bound that we want to adopt for CMIP6 at least.[[BR]] |
| | 136 | Once we get the answer to point 1, |
| | 137 | Is the 8m T, 2m W with same vertical discretization for the upper 2 m and moisture constant and equal to the one of the deepest layer |
| | 138 | a good point to start with? Can we agree on that?[[BR]] |
| | 139 | 3- Is someone able to explain to me (sorry ), why a smoothing is a good way to handle with numerical instabilities, to me it looks like an artificial way of masking a potential problem.[[BR]] |
| | 140 | 4- Do we want to consider the problem of the soil freezing now? |
| | 141 | |
| | 142 | ''About the albedo''[[BR]] |
| | 143 | Whoever will work on that subject, it would be good to know if the bare soil albedo maps prepared by Sebastian and team can be directly pluged into Orchidee-trunk or if we need to see with ISBA people, or to redo the work ourselves. |
| | 144 | |
| | 145 | '''Jan Polcher, 27/01/15''' to F. Chéruy, 27/01/15 |
| | 146 | |
| | 147 | Je suis d'accord avec toi qu'il faudra définir la "configurabilité" des niveaux verticaux en fonction des différentes contraintes. Le nombre de nivaux sera sans doute un résultat plutôt qu'un paramètre. |
| | 148 | |
| | 149 | Sinon pour le filtre. Tu as une onde thermique qui a une vitesse de propagation dans le sol. Il faut avec la discrétisation satisfaire le critère de Courrant-Levy : ne pas sauter plus d'une maille par pas de temps. Cela n'est pas possible avec les niveau de l'hydrologie. Au lieu de filtrer les ondes rapides (On ne veut pas de cela !) il faudra filtrer le résultat et aboutir à une discrétisation effective (au plan numérique) plus lâche que celle de l'hydrologie. |
| | 150 | |
| | 151 | C'est la même chose qui est fait dans LMDZ au pôle ! |
| | 152 | |
