{"action":"created","issue":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211","repository_url":"https://api.github.com/repos/StrikeNP/trac_test","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211/labels{/name}","comments_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211/comments","events_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211/events","html_url":"https://github.com/StrikeNP/trac_test/issues/211","id":334603195,"node_id":"MDU6SXNzdWUzMzQ2MDMxOTU=","number":211,"title":"Modify damping for up2 and vp2 equations (Trac #211)","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":971577014,"node_id":"MDU6TGFiZWw5NzE1NzcwMTQ=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/Migrated%20from%20Trac","name":"Migrated from Trac","color":"FFFFFF","default":false},{"id":971577274,"node_id":"MDU6TGFiZWw5NzE1NzcyNzQ=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/clubb_src","name":"clubb_src","color":"FFFFFF","default":false},{"id":971576851,"node_id":"MDU6TGFiZWw5NzE1NzY4NTE=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/enhancement","name":"enhancement","color":"a2eeef","default":true},{"id":971588502,"node_id":"MDU6TGFiZWw5NzE1ODg1MDI=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/schemena@uwm.edu","name":"schemena@uwm.edu","color":"FFFFFF","default":false}],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":{"url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/1","html_url":"https://github.com/StrikeNP/trac_test/milestone/1","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/1/labels","id":3442829,"node_id":"MDk6TWlsZXN0b25lMzQ0MjgyOQ==","number":1,"title":"Improve the accuracy of CLUBB's horizontal winds","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":6,"closed_issues":9,"state":"open","created_at":"2018-06-21T18:01:23Z","updated_at":"2018-06-22T09:10:26Z","due_on":null,"closed_at":null},"comments":40,"created_at":"2018-06-21T18:16:07Z","updated_at":"2018-06-22T09:11:52Z","closed_at":null,"author_association":"OWNER","body":"up2 and vp2 are damped too much near the inversion and the ground.  We conjecture that this happens because they are damped using the same expression as for wp2.  To fix this, let's try to modify the damping terms on up2 and vp2 (up2_pr1 and vp2_pr1).  \n\nThat is, instead of feeding in C14 to all the functions where it's fed in, we'll feed in\n\n```text\n  C14_fnc(k) = C14 * sqrt( wp2 ) / \n                   sqrt( min( up2_max, 0.5*(up2+vp2) ) )  \n```\n\nwhere up2_max = 1e4 m^2^s^-2^ and C14_fnc is defined on the zm grid.  We'll also need to replace C14 by C14_fnc(k) in C4_C14_1d.\n\nThis change will decrease the damping on up2 and vp2; therefore, it might cause numerical instability.  \n\nWe also might want to consider making a similar change on C4, but let's start with C14 for now.\n\nTests:\n\n1)  Make sure that when C14_fnc(k) = C14, we recover the standard results for all cases.\n2)  Make sure that C14 is no longer fed into any terms, but rather is replaced by C14_fnc.\n\nAttachments:\n\n[plot_explicit_ta_configs.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_explicit_ta_configs.maff)\n[plot_new_pdf_config_1_plot_2.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_new_pdf_config_1_plot_2.maff)\n[plot_combo_pdf_run_3.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_combo_pdf_run_3.maff)\n[plot_input_fields_rtp3_thlp3_1.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_input_fields_rtp3_thlp3_1.maff)\n[plot_new_pdf_20180522_test_1.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_new_pdf_20180522_test_1.maff)\n[plot_attempts_8_10.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_attempts_8_10.maff)\n[plot_attempt_8_only.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_attempt_8_only.maff)\n[plot_beta_1p3.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_beta_1p3.maff)\n[plot_beta_1p3_all.maff](https://github.com/larson-group/trac_attachment_archive/blob/master/trac_test/822/plot_beta_1p3_all.maff)\n\nMigrated from http://carson.math.uwm.edu/trac/clubb/ticket/211\n```json\n{\n    \"status\": \"assigned\", \n    \"changetime\": \"2014-12-18T22:52:36\", \n    \"description\": \"up2 and vp2 are damped too much near the inversion and the ground.  We conjecture that this happens because they are damped using the same expression as for wp2.  To fix this, let's try to modify the damping terms on up2 and vp2 (up2_pr1 and vp2_pr1).  \\n\\nThat is, instead of feeding in C14 to all the functions where it's fed in, we'll feed in\\n{{{\\n  C14_fnc(k) = C14 * sqrt( wp2 ) / \\n                   sqrt( min( up2_max, 0.5*(up2+vp2) ) )  \\n}}}\\n\\nwhere up2_max = 1e4 m^2^s^-2^ and C14_fnc is defined on the zm grid.  We'll also need to replace C14 by C14_fnc(k) in C4_C14_1d.\\n\\nThis change will decrease the damping on up2 and vp2; therefore, it might cause numerical instability.  \\n\\nWe also might want to consider making a similar change on C4, but let's start with C14 for now.\\n\\nTests:\\n\\n1)  Make sure that when C14_fnc(k) = C14, we recover the standard results for all cases.\\n2)  Make sure that C14 is no longer fed into any terms, but rather is replaced by C14_fnc.\", \n    \"reporter\": \"vlarson@uwm.edu\", \n    \"cc\": \"\", \n    \"resolution\": \"\", \n    \"_ts\": \"1418943156030795\", \n    \"component\": \"clubb_src\", \n    \"summary\": \"Modify damping for up2 and vp2 equations\", \n    \"priority\": \"minor\", \n    \"keywords\": \"\", \n    \"time\": \"2009-08-31T22:01:34\", \n    \"milestone\": \"Improve the accuracy of CLUBB's horizontal winds\", \n    \"owner\": \"schemena@uwm.edu\", \n    \"type\": \"enhancement\"\n}\n```\n"},"comment":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/comments/399377089","html_url":"https://github.com/StrikeNP/trac_test/issues/211#issuecomment-399377089","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211","id":399377089,"node_id":"MDEyOklzc3VlQ29tbWVudDM5OTM3NzA4OQ==","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-22T09:11:52Z","updated_at":"2018-06-22T09:11:52Z","author_association":"OWNER","body":"**Trac update on 2009/11/11 at 20:19**\n**ldgrant@uwm.edu commented:**\n\nReplying to [comment:36 vlarson@…]:\n> Replying to [comment:35 ldgrant@…]:\n> > Replying to [comment:31 vlarson@…]:\n>\n> > > I am interested in whether a term like \n>\n>>\n> > > \n>\n\n```text\n> > >             1.0 * < dum/dz >_400m * d(upwp)/dz\n>\n>>\n> > > ```\n>\n> > > could help us generate up2 just beneath the inversion.  Here <>_400m means a running average over 400 m (i.e. +- 200 m).\n>\n>>\n> > I posted some plots [http://www.larson-group.com/ldgrant/protected/ here].  These are results (both without and with inputfields) of using a 400-m average of d(um)/dz in the pr2 shear term.  Compared to the set of plots without inputfields posted below these plots on my protected page (the plots that say, \"RF01 best results from 10/30/09 experiments - using d(upwp)/dz rather than upwp in the pr2 shear term\"), it looks like there's some improvement, and the inputfields run looks good.  In the non-inputfields run, the inversion is still not high enough.\n>\n> This looks like a promising start.  Perhaps it is worth fiddling with the simulations in order to see if it is possible to better match the LES benchmark results.  Plotting the budgets for both inputfields and non-inputfields runs might reveal problems.  It might also be worth adjusting parameters to increase up2 near the inversion and ground.  For instance, you could increase the magnitude of the parameter in front of the dupwp/dz term.  \n\nI attached a printscreen of plotgen results (above) after I increased the dupwp/dz term - I increased \"length\", multiplied by the pr2 dupwp/dz shear terms, from 10,000 to 100,000.  It seems like the shape just isn't quite right yet.  I also tried a tuning run, but that didn't seem to help either.  \n\nI recall something that you mentioned a while ago when we were fiddling with the shear term - the up2 and vp2 plots seem to follow the shape of the pr2 term in the budget plots.  It appears that pr2 is just too large in the ~100-400 m range and peaks too low, below 800 m, when compared to the inputfields run budget plots.\n\n"}}
{"actor":{"display_login":"StrikeNP"}}