{"action":"created","issue":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/498","repository_url":"https://api.github.com/repos/StrikeNP/trac_test","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/498/labels{/name}","comments_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/498/comments","events_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/498/events","html_url":"https://github.com/StrikeNP/trac_test/issues/498","id":332177878,"node_id":"MDU6SXNzdWUzMzIxNzc4Nzg=","number":498,"title":"Latin hypercube: Test autoconversion and accretion formulas in Morrison microphysics (Trac #498)","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":964087254,"node_id":"MDU6TGFiZWw5NjQwODcyNTQ=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/janhft@uwm.edu","name":"janhft@uwm.edu","color":"FFFFFF","default":false},{"id":963997380,"node_id":"MDU6TGFiZWw5NjM5OTczODA=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/task","name":"task","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":48,"closed_issues":33,"state":"open","created_at":"2018-06-13T20:03:15Z","updated_at":"2018-06-15T06:11:59Z","due_on":null,"closed_at":null},"comments":17,"created_at":"2018-06-13T21:51:59Z","updated_at":"2018-06-15T06:12:38Z","closed_at":null,"author_association":"OWNER","body":"ticket:495 would test ice formulas in the KK module, but we'd also like to test the software interface between LH sampling and the Morrison microphysics.  This is partly feasible because Morrison microphysics uses the same [http://carson.math.uwm.edu/trac/sam_clubb/browser/trunk/SRC/MICRO_M2005/module_mp_graupel.F90?rev=822#L1833 autoconversion] and [http://carson.math.uwm.edu/trac/sam_clubb/browser/trunk/SRC/MICRO_M2005/module_mp_graupel.F90?rev=822#L1954 accretion] formulas as the  local KK scheme's [http://carson.math.uwm.edu/trac/clubb/browser/trunk/src/KK_micro/KK_local_means.F90?rev=5128#L89 autoconversion] and [http://carson.math.uwm.edu/trac/clubb/browser/trunk/src/KK_micro/KK_local_means.F90?rev=5128#L132 accretion].  \n\n * Run the dycoms2_rf02_ds case using the analytic upscaled KK microphysics.  In input/case_setups/dycoms2_rf02_ds_model.in, set micro_scheme = \"khairoutdinov_kogan\".  \n * Second simulation: In input/case_setups/dycoms2_rf02_ds_model.in, set micro_scheme = \"morrison\" and set l_input_fields = .true. in order to feed all the available KK profiles into a simulation of dycoms2_rf02_ds using Morrison microphysics.  For an example of the use of inputfields, see the [http://carson.math.uwm.edu/trac/clubb/browser/trunk/input/case_setups/bomex_model.in BOMEX case].  In this simulation, use LH sampling.  For an example of the settings, see [http://carson.math.uwm.edu/trac/clubb/log/trunk/input/case_setups/rico_LH_model.in  the RICO LH case].  Use LH_microphys_calls = 64.\n * In the output from the two simulations, compare 'rrainm_auto' and 'rrainm_accr'.  They should be similar.  Overplot them using GrADS, which is mentioned in the TWiki.  \n * In particular, the results should match better than if you do an inputfields simulation of dycoms2_rf02_ds with Morrison but without LH sampling.\n\n\nAttachments:\n\n[grads_test_rrainm_accr.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/grads_test_rrainm_accr.png)\n[grads_test_rrainm_auto.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/grads_test_rrainm_auto.png)\n[grads_test_wp3.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/grads_test_wp3.png)\n[grads_test_rrainm_mc.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/grads_test_rrainm_mc.png)\n[grads_test_rrainm.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/grads_test_rrainm.png)\n[kk_vs_mor_rrainm_accr.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/kk_vs_mor_rrainm_accr.png)\n[kk_vs_mor_rrainm_auto.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/kk_vs_mor_rrainm_auto.png)\n[kk_vs_mor_rrainm_accr_match.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/kk_vs_mor_rrainm_accr_match.png)\n[kk_vs_mor_rrainm_auto_match.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/kk_vs_mor_rrainm_auto_match.png)\n[12-04-17_kk_vs_mor_rrainm_accr_loc.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-17_kk_vs_mor_rrainm_accr_loc.png)\n[12-04-17_kk_vs_mor_rrainm_auto_loc.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-17_kk_vs_mor_rrainm_auto_loc.png)\n[12-04-17_kk_vs_mor_rrainm_accr_nonloc.2.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-17_kk_vs_mor_rrainm_accr_nonloc.2.png)\n[12-04-17_kk_vs_mor_rrainm_auto_nonloc.2.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-17_kk_vs_mor_rrainm_auto_nonloc.2.png)\n[12-04-24_kk_vs_mor_rrainm.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-24_kk_vs_mor_rrainm.png)\n[12-04-24_kk_vs_mor_rrainm_accr.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-24_kk_vs_mor_rrainm_accr.png)\n[12-04-24_kk_vs_mor_rrainm_auto.png](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/12-04-24_kk_vs_mor_rrainm_auto.png)\n[KK_dycoms2_rf02_ds_model.in](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/KK_dycoms2_rf02_ds_model.in)\n[MOR_dycoms2_rf02_ds_model.in](https://github.com/larson-group/trac_attachment_archive/blob/master/clubb/MOR_dycoms2_rf02_ds_model.in)\n\nMigrated from http://carson.math.uwm.edu/trac/clubb/ticket/498\n```json\n{\n    \"status\": \"closed\", \n    \"changetime\": \"2012-04-25T20:03:13\", \n    \"description\": \"ticket:495 would test ice formulas in the KK module, but we'd also like to test the software interface between LH sampling and the Morrison microphysics.  This is partly feasible because Morrison microphysics uses the same [http://carson.math.uwm.edu/trac/sam_clubb/browser/trunk/SRC/MICRO_M2005/module_mp_graupel.F90?rev=822#L1833 autoconversion] and [http://carson.math.uwm.edu/trac/sam_clubb/browser/trunk/SRC/MICRO_M2005/module_mp_graupel.F90?rev=822#L1954 accretion] formulas as the  local KK scheme's [http://carson.math.uwm.edu/trac/clubb/browser/trunk/src/KK_micro/KK_local_means.F90?rev=5128#L89 autoconversion] and [http://carson.math.uwm.edu/trac/clubb/browser/trunk/src/KK_micro/KK_local_means.F90?rev=5128#L132 accretion].  \\n\\n * Run the dycoms2_rf02_ds case using the analytic upscaled KK microphysics.  In input/case_setups/dycoms2_rf02_ds_model.in, set micro_scheme = \\\"khairoutdinov_kogan\\\".  \\n * Second simulation: In input/case_setups/dycoms2_rf02_ds_model.in, set micro_scheme = \\\"morrison\\\" and set l_input_fields = .true. in order to feed all the available KK profiles into a simulation of dycoms2_rf02_ds using Morrison microphysics.  For an example of the use of inputfields, see the [http://carson.math.uwm.edu/trac/clubb/browser/trunk/input/case_setups/bomex_model.in BOMEX case].  In this simulation, use LH sampling.  For an example of the settings, see [http://carson.math.uwm.edu/trac/clubb/log/trunk/input/case_setups/rico_LH_model.in  the RICO LH case].  Use LH_microphys_calls = 64.\\n * In the output from the two simulations, compare 'rrainm_auto' and 'rrainm_accr'.  They should be similar.  Overplot them using GrADS, which is mentioned in the TWiki.  \\n * In particular, the results should match better than if you do an inputfields simulation of dycoms2_rf02_ds with Morrison but without LH sampling.\\n\", \n    \"reporter\": \"vlarson@uwm.edu\", \n    \"cc\": \"vlarson@uwm.edu, charlass@uwm.edu, dschanen@uwm.edu\", \n    \"resolution\": \"Verified by V. Larson\", \n    \"_ts\": \"1335384193316760\", \n    \"component\": \"clubb_src\", \n    \"summary\": \"Latin hypercube: Test autoconversion and accretion formulas in Morrison microphysics\", \n    \"priority\": \"major\", \n    \"keywords\": \"\", \n    \"time\": \"2012-02-13T22:30:16\", \n    \"milestone\": \"Improve SILHS\", \n    \"owner\": \"janhft@uwm.edu\", \n    \"type\": \"task\"\n}\n```\n"},"comment":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/comments/397523780","html_url":"https://github.com/StrikeNP/trac_test/issues/498#issuecomment-397523780","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/498","id":397523780,"node_id":"MDEyOklzc3VlQ29tbWVudDM5NzUyMzc4MA==","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-15T06:12:38Z","updated_at":"2018-06-15T06:12:38Z","author_association":"OWNER","body":"**Trac update on 2012/04/10 at 22:05**\n**dschanen@uwm.edu commented:**\n\nReplying to [comment:7 dschanen@…]:\n> Replying to [comment:6 janhft@…]:\n> > Dave,\n>\n> > \n> > in r5653 you changed the ```module_mp_graupel.F90``` to output ```rrainm_accr``` and ```rrainm_auto``` for the Morrison microphysics, but it just gives me zeros. Does it work for you?\n>\n> > \n> > Thank you.\n>\n> The individual budgets (including rrainm_auto and rrainm_accr for both K&K and Morrison) will only work if latin hypercube is disabled using the current code.  This is because the accumulation of the budget stats in K&K was originally written with the assumption that the KK_microphys subroutine would be called only once per timestep.  Since the current code lacks the generality to deal with multiple samples per stats_tsamp interval, enabling it by letting l_stats_tsamp to be true in subroutine KK_microphys when latin hypercube sampling is being used would cause an error for LH_microphys_calls > 1 (e.g. if LH_microphys_calls was 2, and dt_main = stats_tsamp = 60 seconds, then we would have an extra sample of rrain_auto per minute of stats), so the code would have to be modified to make a greater number of samples average properly.  I'm not sure how much work this would be; in general I've found the stats_accumulate code to be unnecessarily complex and hard to understand due to its age.\n\nAfter examining the stats code a little closer I think I might have a fairly elegant solution to this problem.  Currently we assume all files should contain ```stats_tout / stats_tsamp``` samples and output our latin hypercube statistics to the zt file.  We could simply add a new file with the name <fname_prefix>_LH_zt that contains ``` LH_micropphys_calls * stats_tout / stats_tsamp```.  The added benefit of this generalization is the we could disable latin hypercube stats when LH is disable entirely, and make the set of ```standard_stats``` smaller for non-LH simulations.\n\n"}}
{"actor":{"display_login":"StrikeNP"}}