{"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":70,"created_at":"2018-06-21T18:16:07Z","updated_at":"2018-06-22T09:12: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/399377412","html_url":"https://github.com/StrikeNP/trac_test/issues/211#issuecomment-399377412","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/211","id":399377412,"node_id":"MDEyOklzc3VlQ29tbWVudDM5OTM3NzQxMg==","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:12:52Z","updated_at":"2018-06-22T09:12:52Z","author_association":"OWNER","body":"**Trac update on 2009/12/16 at 22:40**\n**vlarson@uwm.edu commented:**\n\nThe results in comment:63 look promising.  \n\n1.  Leah suggested adding |dLscale/dz| to both the const1 term and the const2 term.  I think this is a good thing to try.\n\n2.  Of the results described in comment:63, my favorite was output_test5, which set const1 = 4.0.  We could try to boost const1 to 10.0; however, this might increase up2 in Wangara too much.\n\n3.  However, in output_test5, up2 was too large near the ground in BOMEX.  What part of term_pr2 is causing this?  Is it C4=1.0?  The const1 term?  The const2 term?\n\n4.  CLUBB doesn't pick up the maximum in up2 or especially vp2 at 125 m in gabls3_night.  Leah suggested that perhaps these maxima are related to dum/dz and dvm/dz respectively.  However, CLUBB doesn't simulate um and vm accurately, and so we might need to do an inputfields run in order to understand this better.\n\n5.  We might want to run the better configurations using inputfields, so that we can see how the parameterizations behave when the correct wind fields are fed in.\n\n6.  We might re-try adding the term\n\n```text\nconst3 * ( | d(upwp)/dz | + | d(vpwp)/dz | ) * ( |<d(um)/dz>_200m| + |<d(vm)/dz>_200m| )\n\n```\n\n"}}
{"actor":{"display_login":"StrikeNP"}}