{"action":"created","issue":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/755","repository_url":"https://api.github.com/repos/StrikeNP/trac_test","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/755/labels{/name}","comments_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/755/comments","events_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/755/events","html_url":"https://github.com/StrikeNP/trac_test/issues/755","id":330468988,"node_id":"MDU6SXNzdWUzMzA0Njg5ODg=","number":755,"title":"New Methods to Calculate Hydrometeor Component Means and Component Precipitation Fraction. (Trac #755)","user":{"login":"StrikeNP","id":2577059,"node_id":"MDQ6VXNlcjI1NzcwNTk=","avatar_url":"https://avatars0.githubusercontent.com/u/2577059?v=4","gravatar_id":"","url":"https://api.github.com/users/StrikeNP","html_url":"https://github.com/StrikeNP","followers_url":"https://api.github.com/users/StrikeNP/followers","following_url":"https://api.github.com/users/StrikeNP/following{/other_user}","gists_url":"https://api.github.com/users/StrikeNP/gists{/gist_id}","starred_url":"https://api.github.com/users/StrikeNP/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/StrikeNP/subscriptions","organizations_url":"https://api.github.com/users/StrikeNP/orgs","repos_url":"https://api.github.com/users/StrikeNP/repos","events_url":"https://api.github.com/users/StrikeNP/events{/privacy}","received_events_url":"https://api.github.com/users/StrikeNP/received_events","type":"User","site_admin":false},"labels":[{"id":958019852,"node_id":"MDU6TGFiZWw5NTgwMTk4NTI=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/Migrated%20from%20Trac","name":"Migrated from Trac","color":"FFFFFF","default":false},{"id":958043400,"node_id":"MDU6TGFiZWw5NTgwNDM0MDA=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/bmg2@uwm.edu","name":"bmg2@uwm.edu","color":"FFFFFF","default":false},{"id":958042830,"node_id":"MDU6TGFiZWw5NTgwNDI4MzA=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/clubb_src","name":"clubb_src","color":"FFFFFF","default":false},{"id":958019578,"node_id":"MDU6TGFiZWw5NTgwMTk1Nzg=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/enhancement","name":"enhancement","color":"a2eeef","default":true}],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":86,"created_at":"2018-06-07T23:22:25Z","updated_at":"2018-06-10T09:52:38Z","closed_at":null,"author_association":"OWNER","body":"Our current method of calculating ```hm_1``` and ```hm_2``` is based on vertically-integrated cloud water mixing ratio in each component.  In highly skewed cases, 1st component cloud water mixing ratio greatly dwarfs 2nd component cloud water mixing ratio.  The ratio of ```hm_1``` to ```hm_2``` becomes huge.  This produces a hydrometeor PDF that looks disjointed between the two components.  We could provide a better match to SAM.  Additionally, having almost all the rain in one component allows SILHS to draw extreme sample points, resulting in the need for hole-filling.\n\nHowever, RICO (and perhaps other highly skewed cases) requires that ```mu_hm_1``` remains greater (by many multiples) than ```mu_hm_2```.  When ```mu_hm_1``` is set equal to ```mu_hm_2```, the case over-evaporates and under-accretes, losing almost all the rain.  So, for a highly skewed case, we need a method that will set ```mu_hm_1``` > ```mu_hm_2``` without such an extreme overkill as found in the vertically-integrated cloud water method.\n\n\nAdditionally, the vertically-integrated cloud water method produces a PDF realizability issue.  In RICO, near or below cloud base, the turbulent fluxes of hydrometeor fields (```wprrp``` and ```wpNrp```) are negative.  Yet, the vertically-integrated component rc method produces ```hm_1``` >> ```hm_2```.  Since ```w_1``` is always greater than ```w_2```, we are naturally set up for a positive overall correlation (and covariance) between w and the hydrometeor.  However, when w'hm' is negative (fairly substantially at some levels), the within-component in-precip correlation of w and the hydrometeor (```corr_w_hm_i```) needs to be unrealizably negative (less than -1) in order to match ```w_1```, ```hm_1```, ```w_2```, and ```hm_2``` with w'hm'.  The PDF becomes unrealizable.  This is why we can't turn on the ```l_calc_w_corr``` code with the current setup.\n\n\nTo top it all off, \"the ratio method\" of setting precipitation fraction ultimately ends up setting ```mu_hm_1``` = ```mu_hm_2```, regardless of ```hm_1```/```hm_2```.  The only reason this isn't observed in the current standalone model is because the ```hm_1```/```hm_2``` ratio is so great that ```precip_frac_1```/```precip_frac_2``` = ```hm_1```/```hm_2``` results in ```precip_frac_1``` > 1.  Since this must be truncated, ```precip_frac_1``` = 1, ```precip_frac_2``` is increased a little bit to compensate, and the result is ```mu_hm_1``` is substantially larger than ```mu_hm_2```, saving us.  Please see the svn log for r7379 for more details.\n\n\nWe need to rework ```hm_1``` and ```hm_2```, as well as ```precip_frac_1``` and ```precip_frac_2```.\n\nI have made progress on new ideas in r7361, r7364, r7365, r7366, r7369, r7378, and r7379.\n\nAttachments:\n\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/RICO_rr_Nr_z1500_t4260_inptflds_sing_log.png\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/RICO_rr_Nr_z1500_t4260_inptflds_sing_dlog.png\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/RICO_accr_ov_evap_ratio.png\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/RICO_evap_accr_incr_voms.png\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/RICO_evap_accr_new.png\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/agu2014_sgs_hydrometeor.jpg\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/mu_hm_i_method_description.pdf\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/plot_norm_corr.maff\nhttp://carson.math.uwm.edu/trac/clubb/attachment/ticket/755/plot_standard_new_mus.maff\n\nMigrated from http://carson.math.uwm.edu/trac/clubb/ticket/755\n```json\n{\n    \"status\": \"closed\", \n    \"changetime\": \"2016-09-06T19:29:35\", \n    \"description\": \"Our current method of calculating {{{hm_1}}} and {{{hm_2}}} is based on vertically-integrated cloud water mixing ratio in each component.  In highly skewed cases, 1st component cloud water mixing ratio greatly dwarfs 2nd component cloud water mixing ratio.  The ratio of {{{hm_1}}} to {{{hm_2}}} becomes huge.  This produces a hydrometeor PDF that looks disjointed between the two components.  We could provide a better match to SAM.  Additionally, having almost all the rain in one component allows SILHS to draw extreme sample points, resulting in the need for hole-filling.\\n\\nHowever, RICO (and perhaps other highly skewed cases) requires that {{{mu_hm_1}}} remains greater (by many multiples) than {{{mu_hm_2}}}.  When {{{mu_hm_1}}} is set equal to {{{mu_hm_2}}}, the case over-evaporates and under-accretes, losing almost all the rain.  So, for a highly skewed case, we need a method that will set {{{mu_hm_1}}} > {{{mu_hm_2}}} without such an extreme overkill as found in the vertically-integrated cloud water method.\\n\\nAdditionally, the vertically-integrated cloud water method produces a PDF realizability issue.  In RICO, near or below cloud base, the turbulent fluxes of hydrometeor fields ({{{wprrp}}} and {{{wpNrp}}}) are negative.  Yet, the vertically-integrated component rc method produces {{{hm_1}}} >> {{{hm_2}}}.  Since {{{w_1}}} is always greater than {{{w_2}}}, we are naturally set up for a positive overall correlation (and covariance) between w and the hydrometeor.  However, when w'hm' is negative (fairly substantially at some levels), the within-component in-precip correlation of w and the hydrometeor ({{{corr_w_hm_i}}}) needs to be unrealizably negative (less than -1) in order to match {{{w_1}}}, {{{hm_1}}}, {{{w_2}}}, and {{{hm_2}}} with w'hm'.  The PDF becomes unrealizable.  This is why we can't turn on the {{{l_calc_w_corr}}} code with the current setup.\\n\\nTo top it all off, \\\"the ratio method\\\" of setting precipitation fraction ultimately ends up setting {{{mu_hm_1}}} = {{{mu_hm_2}}}, regardless of {{{hm_1}}}/{{{hm_2}}}.  The only reason this isn't observed in the current standalone model is because the {{{hm_1}}}/{{{hm_2}}} ratio is so great that {{{precip_frac_1}}}/{{{precip_frac_2}}} = {{{hm_1}}}/{{{hm_2}}} results in {{{precip_frac_1}}} > 1.  Since this must be truncated, {{{precip_frac_1}}} = 1, {{{precip_frac_2}}} is increased a little bit to compensate, and the result is {{{mu_hm_1}}} is substantially larger than {{{mu_hm_2}}}, saving us.  Please see the svn log for r7379 for more details.\\n\\nWe need to rework {{{hm_1}}} and {{{hm_2}}}, as well as {{{precip_frac_1}}} and {{{precip_frac_2}}}.\\n\\nI have made progress on new ideas in r7361, r7364, r7365, r7366, r7369, r7378, and r7379.\", \n    \"reporter\": \"bmg2@uwm.edu\", \n    \"cc\": \"vlarson@uwm.edu, weberjk@uwm.edu, raut@uwm.edu\", \n    \"resolution\": \"fixed\", \n    \"_ts\": \"1473190175562854\", \n    \"component\": \"clubb_src\", \n    \"summary\": \"New Methods to Calculate Hydrometeor Component Means and Component Precipitation Fraction.\", \n    \"priority\": \"major\", \n    \"keywords\": \"\", \n    \"time\": \"2014-11-12T04:51:30\", \n    \"milestone\": \"\", \n    \"owner\": \"bmg2@uwm.edu\", \n    \"type\": \"enhancement\"\n}\n```\n"},"comment":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/comments/396035941","html_url":"https://github.com/StrikeNP/trac_test/issues/755#issuecomment-396035941","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/755","id":396035941,"node_id":"MDEyOklzc3VlQ29tbWVudDM5NjAzNTk0MQ==","user":{"login":"StrikeNP","id":2577059,"node_id":"MDQ6VXNlcjI1NzcwNTk=","avatar_url":"https://avatars0.githubusercontent.com/u/2577059?v=4","gravatar_id":"","url":"https://api.github.com/users/StrikeNP","html_url":"https://github.com/StrikeNP","followers_url":"https://api.github.com/users/StrikeNP/followers","following_url":"https://api.github.com/users/StrikeNP/following{/other_user}","gists_url":"https://api.github.com/users/StrikeNP/gists{/gist_id}","starred_url":"https://api.github.com/users/StrikeNP/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/StrikeNP/subscriptions","organizations_url":"https://api.github.com/users/StrikeNP/orgs","repos_url":"https://api.github.com/users/StrikeNP/repos","events_url":"https://api.github.com/users/StrikeNP/events{/privacy}","received_events_url":"https://api.github.com/users/StrikeNP/received_events","type":"User","site_admin":false},"created_at":"2018-06-10T09:52:38Z","updated_at":"2018-06-10T09:52:38Z","author_association":"OWNER","body":"Trac update at 2015/03/11 at 03:33\n\n\n* bmg2@uwm.edu changed attachment from \"\" to \"plot_norm_corr.maff\"\n* bmg2@uwm.edu commented:\n\nBefore-and-after plot of all cases for CLUBB r7531.\n\n"}}
{"actor":{"display_login":"StrikeNP"}}