{"action":"created","issue":{"url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/165","repository_url":"https://api.github.com/repos/NCAR/ccpp-physics","labels_url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/165/labels{/name}","comments_url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/165/comments","events_url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/165/events","html_url":"https://github.com/NCAR/ccpp-physics/pull/165","id":366054833,"node_id":"MDExOlB1bGxSZXF1ZXN0MjE5ODI1Mzcz","number":165,"title":"Aerosol-aware Thompson microphysics","user":{"login":"climbfuji","id":8006981,"node_id":"MDQ6VXNlcjgwMDY5ODE=","avatar_url":"https://avatars3.githubusercontent.com/u/8006981?v=4","gravatar_id":"","url":"https://api.github.com/users/climbfuji","html_url":"https://github.com/climbfuji","followers_url":"https://api.github.com/users/climbfuji/followers","following_url":"https://api.github.com/users/climbfuji/following{/other_user}","gists_url":"https://api.github.com/users/climbfuji/gists{/gist_id}","starred_url":"https://api.github.com/users/climbfuji/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/climbfuji/subscriptions","organizations_url":"https://api.github.com/users/climbfuji/orgs","repos_url":"https://api.github.com/users/climbfuji/repos","events_url":"https://api.github.com/users/climbfuji/events{/privacy}","received_events_url":"https://api.github.com/users/climbfuji/received_events","type":"User","site_admin":false},"labels":[{"id":942730116,"node_id":"MDU6TGFiZWw5NDI3MzAxMTY=","url":"https://api.github.com/repos/NCAR/ccpp-physics/labels/do%20not%20merge","name":"do not merge","color":"d93f0b","default":false},{"id":942729356,"node_id":"MDU6TGFiZWw5NDI3MjkzNTY=","url":"https://api.github.com/repos/NCAR/ccpp-physics/labels/work%20in%20progress","name":"work in progress","color":"fbca04","default":false}],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":9,"created_at":"2018-10-02T19:54:50Z","updated_at":"2018-10-05T20:18:10Z","closed_at":null,"author_association":"COLLABORATOR","pull_request":{"url":"https://api.github.com/repos/NCAR/ccpp-physics/pulls/165","html_url":"https://github.com/NCAR/ccpp-physics/pull/165","diff_url":"https://github.com/NCAR/ccpp-physics/pull/165.diff","patch_url":"https://github.com/NCAR/ccpp-physics/pull/165.patch"},"body":"This PR and associated PRs listed below add the Thompson aerosol-aware microphysics scheme to CCPP and enable its use in NEMSfv3gfs.\r\n\r\nThe scheme is based on the version in WRFV4.0 with bugfixes, updates and recommendations from @gthompsnWRF. Additional infrastructure changes/improvements are made to speed up the writing/reading of the lookup tables. This requires building the host model with SIONlib support (for future scalability from single to thousands of processors, see http://www.fz-juelich.de/ias/jsc/EN/Expertise/Support/Software/SIONlib/_node.html).\r\n\r\n**Notes:**\r\n\r\n- The Thompson MP scheme can be run without aerosols by setting the corresponding namelist option to .false.\r\n- If aerosols are switched on, the user can supply 3-dim aerosol concentrations for CCN and IN as well as 2-dimensional CCN and IN surface emission rates as part of the initial conditions. This is host-model dependent. In the case of NEMSfv3gfs, the aerosol climatology that is supplied with WRF (as intermediate files) can be interpolated into the horizontal and vertical grid of the NGGPS initial conditions and injected into these files using a set of Python scripts developed by @climbfuji. - The user can supply all four 3-dim/2-dim aerosol fields as initial conditions, only some of them or none.\r\n- If aerosols are switched on and no initial 3-dim aerosol concentrations data are supplied, a default vertical profile is calculated for every grid column following @gthompsnWRF's original approach. The main difference to the original code in WRF/HRRR is that the calculation of these profiles has been moved out of the scheme and placed into the _pre interstitial scheme. In WRFV4.0, this code exists in the actual scheme and is also copied into the \"real\" part of WRF, which is redundant.\r\n- If aerosols are switched on and no initial 2-dim aerosol emission rates are supplied, these are calculated from the lowest level of the 3-dim aerosol concentrations following @gthompsnWRF's original approach (see previous bullet point for details).\r\n- The RAP/HRRR version of the Thompson microphysics scheme used limits on temperature tendencies produced by the scheme. This code has been placed in the mp_thompson_hrrr_post scheme, but is currently active only for grid spacings <50km (for >50km, the limit is set to \"huge\").\r\n\r\n**Limitations/Todo:**\r\n\r\n- The scheme was tested to run and produce comparable results to the reference (Zhao-Carr microphysics) and the old (partially wrong non-aerosol-aware) implementation of the Thompson scheme in FV3. These tests were conducted on a coarse C96 mesh for 6h-long forecasts (cold starts) on a Macbook and on a C768 mesh for 6h-long forecasts (cold starts) on Theia.\r\n- Todo: Investigate the need/best value for temperature limits for the different grid spacings.\r\n- Todo: Test/develop restarts (warm starts) with aerosol concentrations / aerosol surface emissions.\r\n- Todo: Speedup calculation of lookup tables using OpenMP (or MPI - but this is against the CCPP standards).\r\n- Todo: Test/improve SIONlib writing and reading of lookup tables on Theia and Cheyenne.\r\n"},"comment":{"url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/comments/427486855","html_url":"https://github.com/NCAR/ccpp-physics/pull/165#issuecomment-427486855","issue_url":"https://api.github.com/repos/NCAR/ccpp-physics/issues/165","id":427486855,"node_id":"MDEyOklzc3VlQ29tbWVudDQyNzQ4Njg1NQ==","user":{"login":"climbfuji","id":8006981,"node_id":"MDQ6VXNlcjgwMDY5ODE=","avatar_url":"https://avatars3.githubusercontent.com/u/8006981?v=4","gravatar_id":"","url":"https://api.github.com/users/climbfuji","html_url":"https://github.com/climbfuji","followers_url":"https://api.github.com/users/climbfuji/followers","following_url":"https://api.github.com/users/climbfuji/following{/other_user}","gists_url":"https://api.github.com/users/climbfuji/gists{/gist_id}","starred_url":"https://api.github.com/users/climbfuji/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/climbfuji/subscriptions","organizations_url":"https://api.github.com/users/climbfuji/orgs","repos_url":"https://api.github.com/users/climbfuji/repos","events_url":"https://api.github.com/users/climbfuji/events{/privacy}","received_events_url":"https://api.github.com/users/climbfuji/received_events","type":"User","site_admin":false},"created_at":"2018-10-05T20:18:10Z","updated_at":"2018-10-05T20:18:10Z","author_association":"COLLABORATOR","body":"Update: test results for total precipitation and 2m temperature after 6h of integration on a C768 mesh for a reference run (Zhao-Carr MP), the old Thompson GFS version without aerosols, and the new Thompson CCPP version without aerosols, default aerosol profiles and aerosol climatology.\r\n\r\n<img width=\"772\" alt=\"reference_zhao_carr_t2m_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558149-65568f80-c8a9-11e8-995c-c54e563ed13c.png\">\r\n<img width=\"771\" alt=\"reference_zhao_carr_tprcp_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558150-65568f80-c8a9-11e8-915c-3baf2745976c.png\">\r\n<img width=\"772\" alt=\"thompson_gfs_no_aerosols_t2m_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558151-65568f80-c8a9-11e8-9e8d-6a256cf57aa9.png\">\r\n<img width=\"768\" alt=\"thompson_gfs_no_aerosols_tprcp_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558152-65ef2600-c8a9-11e8-9735-de595d9ea5ab.png\">\r\n<img width=\"777\" alt=\"thompson_hrrr_aerosol_climatology_t2m_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558153-65ef2600-c8a9-11e8-80f4-89617b7440d5.png\">\r\n<img width=\"773\" alt=\"thompson_hrrr_aerosol_climatology_tprcp_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558154-65ef2600-c8a9-11e8-969a-da7e1264b29b.png\">\r\n<img width=\"771\" alt=\"thompson_hrrr_default_aerosols_t2m_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558155-65ef2600-c8a9-11e8-880a-5b00e11a9e13.png\">\r\n<img width=\"775\" alt=\"thompson_hrrr_default_aerosols_tprcp_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558156-65ef2600-c8a9-11e8-88b7-84e808014644.png\">\r\n<img width=\"775\" alt=\"thompson_hrrr_no_aerosols_t2m_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558157-65ef2600-c8a9-11e8-96ca-80f2590baa6b.png\">\r\n<img width=\"778\" alt=\"thompson_hrrr_no_aerosols_tprcp_c768_6h\" src=\"https://user-images.githubusercontent.com/8006981/46558159-65ef2600-c8a9-11e8-9745-1cab589e0690.png\">\r\n"}}
{"actor":{"display_login":"climbfuji"}}