{"action":"created","issue":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/682","repository_url":"https://api.github.com/repos/StrikeNP/trac_test","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/682/labels{/name}","comments_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/682/comments","events_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/682/events","html_url":"https://github.com/StrikeNP/trac_test/issues/682","id":332199120,"node_id":"MDU6SXNzdWUzMzIxOTkxMjA=","number":682,"title":"Create mirrored sampling option in SILHS (Trac #682)","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":963997231,"node_id":"MDU6TGFiZWw5NjM5OTcyMzE=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/Migrated%20from%20Trac","name":"Migrated from Trac","color":"FFFFFF","default":false},{"id":963997516,"node_id":"MDU6TGFiZWw5NjM5OTc1MTY=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/clubb_src","name":"clubb_src","color":"FFFFFF","default":false},{"id":963997043,"node_id":"MDU6TGFiZWw5NjM5OTcwNDM=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/enhancement","name":"enhancement","color":"a2eeef","default":true},{"id":964076735,"node_id":"MDU6TGFiZWw5NjQwNzY3MzU=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/raut@uwm.edu","name":"raut@uwm.edu","color":"FFFFFF","default":false}],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":{"url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/8","html_url":"https://github.com/StrikeNP/trac_test/milestone/8","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/8/labels","id":3423184,"node_id":"MDk6TWlsZXN0b25lMzQyMzE4NA==","number":8,"title":"Improve SILHS","description":null,"creator":{"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},"open_issues":33,"closed_issues":48,"state":"open","created_at":"2018-06-13T20:03:15Z","updated_at":"2018-06-15T18:00:06Z","due_on":null,"closed_at":null},"comments":2,"created_at":"2018-06-13T23:21:12Z","updated_at":"2018-06-15T22:18:35Z","closed_at":null,"author_association":"OWNER","body":"\n**Introduction**\n\nCurrently, WRF-CLUBB-SILHS is crashing (https://github.com/larson-group/wrf/issues/79) and so is CAM-CLUBB-SILHS (https://github.com/larson-group/cpt/issues/189).   Perhaps the reason is that the samples are too noisy.  SILHS tries to reduce noise by sampling preferentially within cloud at the k_lh_start grid level (```l_lh_cloud_weighted_sampling```).  However, that doesn't guarantee preferential sampling within cloud at other altitudes.\n\nPerhaps instead we need an option that tends to produce samples whose sample mean is close to CLUBB's grid mean.  One approximate way to do this is, for every subcolumn of uniformly distributed points, choose the mirror image subcolumn.  For a symmetric function like a single Gaussian, this would guarantee that the sample mean equals the grid mean.  For a mixture of normal or lognormals, this is not guaranteed, but I still think we might achieve good variance reduction at all grid levels, not just k_lh_start.\n\n**Technical spec**\n\nLet's call the option ```l_lh_mirrored_sampling```.  When this is true, we'll generate a full vertical subcolumn of ''uniformly'' distributed points (\"```x```\"), and then we'll generate a mirror subcolumn simply by computing ```1-x``` for all variates, including the d+1 and d+2 components.  Then we'll proceed with choosing the mixture component and transforming to normal/lognormal space as before.  \n\nWe'll need to idiot-proof the code by making sure that a user cannot set both \n```l_lh_mirrored_sampling``` and ```l_lh_cloud_weighted_sampling``` to true.  \n\n**Tests**\n\nShould we test convergence to the analytic solution?\n\nAttachments:\n\n[plout.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/plout.maff)\n[plout.diff](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/plout.diff)\n[rrainm.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/rrainm.png)\n[precip_frac.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/precip_frac.png)\n[rrain_2.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/rrain_2.png)\n[mixt_frac_imp.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/mixt_frac_imp.maff)\n[mixt_frac_imp.diff](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/mixt_frac_imp.diff)\n\nMigrated from http://carson.math.uwm.edu/trac/clubb/ticket/682\n```json\n{\n    \"status\": \"closed\", \n    \"changetime\": \"2014-09-01T20:47:30\", \n    \"description\": \"\\n'''Introduction'''\\n\\nCurrently, WRF-CLUBB-SILHS is crashing (wrf:ticket:79) and so is CAM-CLUBB-SILHS (cpt:ticket:189).   Perhaps the reason is that the samples are too noisy.  SILHS tries to reduce noise by sampling preferentially within cloud at the k_lh_start grid level ({{{l_lh_cloud_weighted_sampling}}}).  However, that doesn't guarantee preferential sampling within cloud at other altitudes.\\n\\nPerhaps instead we need an option that tends to produce samples whose sample mean is close to CLUBB's grid mean.  One approximate way to do this is, for every subcolumn of uniformly distributed points, choose the mirror image subcolumn.  For a symmetric function like a single Gaussian, this would guarantee that the sample mean equals the grid mean.  For a mixture of normal or lognormals, this is not guaranteed, but I still think we might achieve good variance reduction at all grid levels, not just k_lh_start.\\n\\n'''Technical spec'''\\n\\nLet's call the option {{{l_lh_mirrored_sampling}}}.  When this is true, we'll generate a full vertical subcolumn of ''uniformly'' distributed points (\\\"{{{x}}}\\\"), and then we'll generate a mirror subcolumn simply by computing {{{1-x}}} for all variates, including the d+1 and d+2 components.  Then we'll proceed with choosing the mixture component and transforming to normal/lognormal space as before.  \\n\\nWe'll need to idiot-proof the code by making sure that a user cannot set both \\n{{{l_lh_mirrored_sampling}}} and {{{l_lh_cloud_weighted_sampling}}} to true.  \\n\\n'''Tests'''\\n\\nShould we test convergence to the analytic solution?\", \n    \"reporter\": \"vlarson@uwm.edu\", \n    \"cc\": \"vlarson@uwm.edu\", \n    \"resolution\": \"wontfix\", \n    \"_ts\": \"1409604450401269\", \n    \"component\": \"clubb_src\", \n    \"summary\": \"Create mirrored sampling option in SILHS\", \n    \"priority\": \"major\", \n    \"keywords\": \"\", \n    \"time\": \"2014-05-08T21:44:49\", \n    \"milestone\": \"Improve SILHS\", \n    \"owner\": \"raut@uwm.edu\", \n    \"type\": \"enhancement\"\n}\n```\n"},"comment":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/comments/397756326","html_url":"https://github.com/StrikeNP/trac_test/issues/682#issuecomment-397756326","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/682","id":397756326,"node_id":"MDEyOklzc3VlQ29tbWVudDM5Nzc1NjMyNg==","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-15T22:18:35Z","updated_at":"2018-06-15T22:18:35Z","author_association":"OWNER","body":"**Trac update on 2014/05/22 at 20:11**\n\n* raut@uwm.edu changed attachment from \"\" to \"plout.diff\"\n* **raut@uwm.edu commented:**\n\nDifferences from CLUBB r6890 used to create raw-attachment:plout.maff.\n\n"}}
{"actor":{"display_login":"StrikeNP"}}