Valeurs négative dans cvltr_noscav

[Laurent Fairhead]

Envoyé par O. Boucher:

Didier et Anne m'ont fait remonter un cas un peu curieux où les sels marins augmentent progressivement d'un facteur 10⁶ sur un point de grille puis redescendent, le tout en quelques jours. Après enquête, j'observe que le coeur du problème semble venir de la routine cvltr_noscav qu'on utilise quand on est en couplé avec INCA car le lessivage reste celui qui est fait dans INCA. On a un profil de sels marins très concentrés vers la surface (normal pour de gros aérosols), et cvltr_noscav conserve bien la masse sur la colonne mais créée une valeur assez négative sur le niveau n°4. Puis on remet cette valeur négative à 0 dans phytrac (ce qui crée pas mal de masse), et ainsi de suite.

Ce n'est pas très physique de découpler transport convectif et lessivage, surtout pour des aérosols de sels marins. Et une solution pourrait être de désactiver cvltr_noscav pour les modes grossiers des aérosols. Néanmoins ça vaut peut-être le coup de comprendre pourquoi cvltr_noscav n'est pas suffisamment robuste sur ce cas de figure. J'ai vérifié qu'on utilise bien la dernière version.

[Laurent Fairhead]

Jean-Yves: Le problème arrive avec les aérosols grossiers qui présentent des gradients de concentration élevés en bas de la couche limite. Le schéma explicite utilisé entraîne alors des concentrations négatives. Développement à prévoir: schéma implicite pour lermes d'échanges locaux.