Integration of theta-l dycore into EAM Hydrostatic (H) and - PowerPoint PPT Presentation

Integration of theta-l dycore into EAM Hydrostatic (H) and non-hydrostatic (NH) formulations for theta-l: u t + ( r s ⇥ u + 2 Ω ) ⇥ u + 1 s u s + 1 2 r s u 2 + ˙ hor. momentum (H/NH) κθ v r s Π + µ r s φ = 0 , vertical momentum (NH) w t + u · r s w + ˙ sw s + g (1 � µ ) = 0 , geopotential eqn (NH) φ t + u · r s φ + ˙ s φ s � gw = 0 , theta eqn (H/NH) Θ t + r s · ( u Θ ) + ( ˙ s Θ ) s = 0 , continuity eqn (H/NH) ( π s ) t + r s · ( u π s ) + ( ˙ s π s ) s = 0 , π � hydrostatic pressure , p � nonhydrostatic pressure , EOS : φ s = � ΘΠ p ; Θ = π s θ v , µ = p s , π s eqn. of state (NH)

<latexit sha1_base64="nSl0dlTWSIhicdtzJwD2hqgF8Tg=">ACJ3icbVDLSgMxFM3UV62vqks3wSKISJkRQRGUohuXtfQFnTpk0kwbmnmY3CmUoX/jxl9xI6iILv0T04dQWw8EDuecy809biS4AtP8MlILi0vLK+nVzNr6xuZWdnunqsJYUlahoQhl3SWKCR6wCnAQrB5JRnxXsJrbvRn6tR6TiodBGfoRa/qkHXCPUwJacrJXNnQYEKeHLy6x7UlCk9L9UXmQlOwiHxj+yEmLWyXp0IRd9TvlJPNmXlzBDxPrAnJoQmKTvbVboU09lkAVBClGpYZQTMhEjgVbJCxY8UiQrukzRqaBsRnqpmM7hzgA620sBdK/QLAI3V6IiG+Un3f1UmfQEfNekPxP68Rg3feTHgQxcACOl7kxQJDiIel4RaXjILoa0Ko5PqvmHaILgt0tRldgjV78jypnuQtze9Oc4XrSR1ptIf20SGy0BkqoFtURBVE0SN6Rm/o3XgyXowP43McTRmTmV30B8b3D6KZpTE=</latexit> <latexit sha1_base64="nSl0dlTWSIhicdtzJwD2hqgF8Tg=">ACJ3icbVDLSgMxFM3UV62vqks3wSKISJkRQRGUohuXtfQFnTpk0kwbmnmY3CmUoX/jxl9xI6iILv0T04dQWw8EDuecy809biS4AtP8MlILi0vLK+nVzNr6xuZWdnunqsJYUlahoQhl3SWKCR6wCnAQrB5JRnxXsJrbvRn6tR6TiodBGfoRa/qkHXCPUwJacrJXNnQYEKeHLy6x7UlCk9L9UXmQlOwiHxj+yEmLWyXp0IRd9TvlJPNmXlzBDxPrAnJoQmKTvbVboU09lkAVBClGpYZQTMhEjgVbJCxY8UiQrukzRqaBsRnqpmM7hzgA620sBdK/QLAI3V6IiG+Un3f1UmfQEfNekPxP68Rg3feTHgQxcACOl7kxQJDiIel4RaXjILoa0Ko5PqvmHaILgt0tRldgjV78jypnuQtze9Oc4XrSR1ptIf20SGy0BkqoFtURBVE0SN6Rm/o3XgyXowP43McTRmTmV30B8b3D6KZpTE=</latexit> <latexit sha1_base64="nSl0dlTWSIhicdtzJwD2hqgF8Tg=">ACJ3icbVDLSgMxFM3UV62vqks3wSKISJkRQRGUohuXtfQFnTpk0kwbmnmY3CmUoX/jxl9xI6iILv0T04dQWw8EDuecy809biS4AtP8MlILi0vLK+nVzNr6xuZWdnunqsJYUlahoQhl3SWKCR6wCnAQrB5JRnxXsJrbvRn6tR6TiodBGfoRa/qkHXCPUwJacrJXNnQYEKeHLy6x7UlCk9L9UXmQlOwiHxj+yEmLWyXp0IRd9TvlJPNmXlzBDxPrAnJoQmKTvbVboU09lkAVBClGpYZQTMhEjgVbJCxY8UiQrukzRqaBsRnqpmM7hzgA620sBdK/QLAI3V6IiG+Un3f1UmfQEfNekPxP68Rg3feTHgQxcACOl7kxQJDiIel4RaXjILoa0Ko5PqvmHaILgt0tRldgjV78jypnuQtze9Oc4XrSR1ptIf20SGy0BkqoFtURBVE0SN6Rm/o3XgyXowP43McTRmTmV30B8b3D6KZpTE=</latexit> <latexit sha1_base64="nSl0dlTWSIhicdtzJwD2hqgF8Tg=">ACJ3icbVDLSgMxFM3UV62vqks3wSKISJkRQRGUohuXtfQFnTpk0kwbmnmY3CmUoX/jxl9xI6iILv0T04dQWw8EDuecy809biS4AtP8MlILi0vLK+nVzNr6xuZWdnunqsJYUlahoQhl3SWKCR6wCnAQrB5JRnxXsJrbvRn6tR6TiodBGfoRa/qkHXCPUwJacrJXNnQYEKeHLy6x7UlCk9L9UXmQlOwiHxj+yEmLWyXp0IRd9TvlJPNmXlzBDxPrAnJoQmKTvbVboU09lkAVBClGpYZQTMhEjgVbJCxY8UiQrukzRqaBsRnqpmM7hzgA620sBdK/QLAI3V6IiG+Un3f1UmfQEfNekPxP68Rg3feTHgQxcACOl7kxQJDiIel4RaXjILoa0Ko5PqvmHaILgt0tRldgjV78jypnuQtze9Oc4XrSR1ptIf20SGy0BkqoFtURBVE0SN6Rm/o3XgyXowP43McTRmTmV30B8b3D6KZpTE=</latexit> <latexit sha1_base64="SmdrLGrfvBQiIblCZjc/F4Vukbc=">ACVnicbZFbS8MwFMfTetmct6qPvgSH4J1WBH0RH3wcYpTYZ0jTdMtmKY1ORVG6ZfUF/0ovojZVkWnB0J+Z9zcvknSAX4Lpvlj0xOTVdqc7UZufmFxadpeUbnWSKsiZNRKLuAqKZ4JI1gYNgd6liJA4Euw0ezgb52yemNE/kNfRT1o5JV/KIUwJG6jixHylC83MfegxI56nIz6E4/lrd53LbK/Au/hbkjv+YkXCs4Bj74XCfjav7rc2ReD2aivzKb/CSO07d3XOHgf+CV0IdldHoOM9+mNAsZhKoIFq3PDeFdk4UcCpYUfMzVJCH0iXtQxKEjPdzoe2FHjdKCGOEmWGBDxUf3bkJNa6HwemMibQ0+O5gfhfrpVBdNTOuUwzYJKODoygSHBA49xyBWjIPoGCFXc3BXTHjH+gPmJmjHBG3/yX7jZ3/MXx7UT05LO6poFa2hDeShQ3SCLlADNRFL+jdsq0J69X6sKfsyqjUtsqeFfQrbOcTEN6z6g=</latexit> <latexit sha1_base64="SmdrLGrfvBQiIblCZjc/F4Vukbc=">ACVnicbZFbS8MwFMfTetmct6qPvgSH4J1WBH0RH3wcYpTYZ0jTdMtmKY1ORVG6ZfUF/0ovojZVkWnB0J+Z9zcvknSAX4Lpvlj0xOTVdqc7UZufmFxadpeUbnWSKsiZNRKLuAqKZ4JI1gYNgd6liJA4Euw0ezgb52yemNE/kNfRT1o5JV/KIUwJG6jixHylC83MfegxI56nIz6E4/lrd53LbK/Au/hbkjv+YkXCs4Bj74XCfjav7rc2ReD2aivzKb/CSO07d3XOHgf+CV0IdldHoOM9+mNAsZhKoIFq3PDeFdk4UcCpYUfMzVJCH0iXtQxKEjPdzoe2FHjdKCGOEmWGBDxUf3bkJNa6HwemMibQ0+O5gfhfrpVBdNTOuUwzYJKODoygSHBA49xyBWjIPoGCFXc3BXTHjH+gPmJmjHBG3/yX7jZ3/MXx7UT05LO6poFa2hDeShQ3SCLlADNRFL+jdsq0J69X6sKfsyqjUtsqeFfQrbOcTEN6z6g=</latexit> <latexit sha1_base64="SmdrLGrfvBQiIblCZjc/F4Vukbc=">ACVnicbZFbS8MwFMfTetmct6qPvgSH4J1WBH0RH3wcYpTYZ0jTdMtmKY1ORVG6ZfUF/0ovojZVkWnB0J+Z9zcvknSAX4Lpvlj0xOTVdqc7UZufmFxadpeUbnWSKsiZNRKLuAqKZ4JI1gYNgd6liJA4Euw0ezgb52yemNE/kNfRT1o5JV/KIUwJG6jixHylC83MfegxI56nIz6E4/lrd53LbK/Au/hbkjv+YkXCs4Bj74XCfjav7rc2ReD2aivzKb/CSO07d3XOHgf+CV0IdldHoOM9+mNAsZhKoIFq3PDeFdk4UcCpYUfMzVJCH0iXtQxKEjPdzoe2FHjdKCGOEmWGBDxUf3bkJNa6HwemMibQ0+O5gfhfrpVBdNTOuUwzYJKODoygSHBA49xyBWjIPoGCFXc3BXTHjH+gPmJmjHBG3/yX7jZ3/MXx7UT05LO6poFa2hDeShQ3SCLlADNRFL+jdsq0J69X6sKfsyqjUtsqeFfQrbOcTEN6z6g=</latexit> <latexit sha1_base64="SmdrLGrfvBQiIblCZjc/F4Vukbc=">ACVnicbZFbS8MwFMfTetmct6qPvgSH4J1WBH0RH3wcYpTYZ0jTdMtmKY1ORVG6ZfUF/0ovojZVkWnB0J+Z9zcvknSAX4Lpvlj0xOTVdqc7UZufmFxadpeUbnWSKsiZNRKLuAqKZ4JI1gYNgd6liJA4Euw0ezgb52yemNE/kNfRT1o5JV/KIUwJG6jixHylC83MfegxI56nIz6E4/lrd53LbK/Au/hbkjv+YkXCs4Bj74XCfjav7rc2ReD2aivzKb/CSO07d3XOHgf+CV0IdldHoOM9+mNAsZhKoIFq3PDeFdk4UcCpYUfMzVJCH0iXtQxKEjPdzoe2FHjdKCGOEmWGBDxUf3bkJNa6HwemMibQ0+O5gfhfrpVBdNTOuUwzYJKODoygSHBA49xyBWjIPoGCFXc3BXTHjH+gPmJmjHBG3/yX7jZ3/MXx7UT05LO6poFa2hDeShQ3SCLlADNRFL+jdsq0J69X6sKfsyqjUtsqeFfQrbOcTEN6z6g=</latexit> <latexit sha1_base64="jxlMIRocSYQIH1xAFcAwb5QS8AM=">ACI3icbVBNS8NAFNzUr1q/qh69BItQPZREBEUQivbgsYpVoUnLZrtpl242YfelUEL+ixf/ihcPinjx4H9xm+agrQMLw8w83r7xIs4UWNaXUVhYXFpeKa6W1tY3NrfK2zv3KowloS0S8lA+elhRzgRtAQNOHyNJceBx+uANryb+w4hKxUJxB+OIugHuC+YzgkFL3fK540tMkoYDAwq4O0qTBqQXU7F62zk67Ig0uXWarJOINM3Dd1mqW65YNSuDOU/snFRQjma3/OH0QhIHVADhWKm2bUXgJlgCI5ymJSdWNMJkiPu0ranAVukt2Ymgda6Zl+KPUTYGbq74kEB0qNA08nAwDNetNxP+8dgz+mZswEcVABZku8mNuQmhOCjN7TFICfKwJpLpv5pkgHUPoGst6RLs2ZPnyf1xzdb85qRSv8zrKI9tI+qyEanqI6uURO1EFP6AW9oXfj2Xg1PozPabRg5DO76A+M7x/PEqWZ</latexit> <latexit sha1_base64="jxlMIRocSYQIH1xAFcAwb5QS8AM=">ACI3icbVBNS8NAFNzUr1q/qh69BItQPZREBEUQivbgsYpVoUnLZrtpl242YfelUEL+ixf/ihcPinjx4H9xm+agrQMLw8w83r7xIs4UWNaXUVhYXFpeKa6W1tY3NrfK2zv3KowloS0S8lA+elhRzgRtAQNOHyNJceBx+uANryb+w4hKxUJxB+OIugHuC+YzgkFL3fK540tMkoYDAwq4O0qTBqQXU7F62zk67Ig0uXWarJOINM3Dd1mqW65YNSuDOU/snFRQjma3/OH0QhIHVADhWKm2bUXgJlgCI5ymJSdWNMJkiPu0ranAVukt2Ymgda6Zl+KPUTYGbq74kEB0qNA08nAwDNetNxP+8dgz+mZswEcVABZku8mNuQmhOCjN7TFICfKwJpLpv5pkgHUPoGst6RLs2ZPnyf1xzdb85qRSv8zrKI9tI+qyEanqI6uURO1EFP6AW9oXfj2Xg1PozPabRg5DO76A+M7x/PEqWZ</latexit> <latexit sha1_base64="jxlMIRocSYQIH1xAFcAwb5QS8AM=">ACI3icbVBNS8NAFNzUr1q/qh69BItQPZREBEUQivbgsYpVoUnLZrtpl242YfelUEL+ixf/ihcPinjx4H9xm+agrQMLw8w83r7xIs4UWNaXUVhYXFpeKa6W1tY3NrfK2zv3KowloS0S8lA+elhRzgRtAQNOHyNJceBx+uANryb+w4hKxUJxB+OIugHuC+YzgkFL3fK540tMkoYDAwq4O0qTBqQXU7F62zk67Ig0uXWarJOINM3Dd1mqW65YNSuDOU/snFRQjma3/OH0QhIHVADhWKm2bUXgJlgCI5ymJSdWNMJkiPu0ranAVukt2Ymgda6Zl+KPUTYGbq74kEB0qNA08nAwDNetNxP+8dgz+mZswEcVABZku8mNuQmhOCjN7TFICfKwJpLpv5pkgHUPoGst6RLs2ZPnyf1xzdb85qRSv8zrKI9tI+qyEanqI6uURO1EFP6AW9oXfj2Xg1PozPabRg5DO76A+M7x/PEqWZ</latexit> <latexit sha1_base64="jxlMIRocSYQIH1xAFcAwb5QS8AM=">ACI3icbVBNS8NAFNzUr1q/qh69BItQPZREBEUQivbgsYpVoUnLZrtpl242YfelUEL+ixf/ihcPinjx4H9xm+agrQMLw8w83r7xIs4UWNaXUVhYXFpeKa6W1tY3NrfK2zv3KowloS0S8lA+elhRzgRtAQNOHyNJceBx+uANryb+w4hKxUJxB+OIugHuC+YzgkFL3fK540tMkoYDAwq4O0qTBqQXU7F62zk67Ig0uXWarJOINM3Dd1mqW65YNSuDOU/snFRQjma3/OH0QhIHVADhWKm2bUXgJlgCI5ymJSdWNMJkiPu0ranAVukt2Ymgda6Zl+KPUTYGbq74kEB0qNA08nAwDNetNxP+8dgz+mZswEcVABZku8mNuQmhOCjN7TFICfKwJpLpv5pkgHUPoGst6RLs2ZPnyf1xzdb85qRSv8zrKI9tI+qyEanqI6uURO1EFP6AW9oXfj2Xg1PozPabRg5DO76A+M7x/PEqWZ</latexit> Integration of hydrostatic model For hydrostatic model most of effort was spent on infrastructure and conversion of temperature tendencies. Theta model uses this particular definition of potential temperature: Potential temperature θ v := R ∗ T theta variable in code definition in homme, R Π , Θ v := π s θ v for conservation allows decoupling form of vapor in dycore There are a few options to convert temperature tendencies into theta tendencies: Dt = ( R ∗ ) n D θ v DT 1. R Π n Dt Implemented because it keeps theta definition consistent = ( R ∗ ) n +1 T n +1 D θ v Dt = θ n +1 θ n +1 2. − θ n v , v v R Π n +1

Integration of nonhydrostatic model Integration is done, some conceptual questions are work in progress.

Theta hydro, F case, 5 years climo (very preliminary) Theta model vs default model comparison:

Tests and consistency As we introduce more features into the dycore, including coupling options (ftypes) and dycore models, we get more tests into E3SM, too. Thanks to Jim Foucar now there is an option to reuse the same executable in tests to save build time (only namelists change). Potentially, there could be inconsistencies in definitions of potential temperature, Exner pressure, geopotential, energy, etc. To monitor for and to resolve these inconsistencies are on our list.