{"action":"created","issue":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/31","repository_url":"https://api.github.com/repos/StrikeNP/trac_test","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/31/labels{/name}","comments_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/31/comments","events_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/31/events","html_url":"https://github.com/StrikeNP/trac_test/issues/31","id":334599157,"node_id":"MDU6SXNzdWUzMzQ1OTkxNTc=","number":31,"title":"sounding.in: Allow input of pressure as the independent coordinate as well as height (Trac #31)","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":971577070,"node_id":"MDU6TGFiZWw5NzE1NzcwNzA=","url":"https://api.github.com/repos/StrikeNP/trac_test/labels/fasching@uwm.edu","name":"fasching@uwm.edu","color":"FFFFFF","default":false}],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":{"url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/4","html_url":"https://github.com/StrikeNP/trac_test/milestone/4","labels_url":"https://api.github.com/repos/StrikeNP/trac_test/milestones/4/labels","id":3442833,"node_id":"MDk6TWlsZXN0b25lMzQ0MjgzMw==","number":4,"title":"3. Refactor CLUBB","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":141,"closed_issues":3,"state":"open","created_at":"2018-06-21T18:01:51Z","updated_at":"2018-06-22T00:20:39Z","due_on":null,"closed_at":null},"comments":13,"created_at":"2018-06-21T18:03:20Z","updated_at":"2018-06-22T00:21:12Z","closed_at":null,"author_association":"OWNER","body":"Right now, a CLUBB user can input temperature in the initial sounding as either theta or theta_l.  However, all sounding data must be specified as a function of height (altitude), not pressure.  The drawback is that many GCSS cases specify data in terms of pressure, and we have to manually convert these data beforehand using MATLAB.\n\nTherefore, allow the height column to be entered as pressure values.  Then convert pressure to height within CLUBB.  This will require creating a new subroutine, inverse_hydrostatic, which takes thvm(pressure) as input and integrates in order to find z(pressure).  It would be the inverse of hydrostatic and would use the algorithm:\n\n1.  Compute exner and thvm on the zt_snd (given) sounding grid.\n\n2.  Interpolate exner to zm_snd grid.  Approximate using assumption of equal spacing:\n\nexner_zm(k) = 0.5 * ( exner( k ) + exner( k+1 ) )\n\n3.  Integrate hydrostatic equation to find altitudes of the zm_snd levels:\n\nzm_snd(k) =  zm_snd( k-1 ) - ( Cp / grav ) * thvm(k) * d_exner_zm(k) \n\nwhere\n\nd_exner_zm(k) = exner_zm(k) - exner_zm(k-1).\n\nWe need to iterate in the old hydrostatic code, because the hydrostatic code doesn't know pressure and therefore can't compute thvm.  Here we know pressure, and so we can compute thvm without approximation or interpolation.    \n\nThen we could set up M-PACE B using this initialize procedure, using a case-specific function to input profiles rather than sounding.in.  \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/31\n```json\n{\n    \"status\": \"closed\", \n    \"changetime\": \"2009-06-05T14:05:33\", \n    \"description\": \"Right now, a CLUBB user can input temperature in the initial sounding as either theta or theta_l.  However, all sounding data must be specified as a function of height (altitude), not pressure.  The drawback is that many GCSS cases specify data in terms of pressure, and we have to manually convert these data beforehand using MATLAB.\\n\\nTherefore, allow the height column to be entered as pressure values.  Then convert pressure to height within CLUBB.  This will require creating a new subroutine, inverse_hydrostatic, which takes thvm(pressure) as input and integrates in order to find z(pressure).  It would be the inverse of hydrostatic and would use the algorithm:\\n\\n1.  Compute exner and thvm on the zt_snd (given) sounding grid.\\n\\n2.  Interpolate exner to zm_snd grid.  Approximate using assumption of equal spacing:\\n\\nexner_zm(k) = 0.5 * ( exner( k ) + exner( k+1 ) )\\n\\n3.  Integrate hydrostatic equation to find altitudes of the zm_snd levels:\\n\\nzm_snd(k) =  zm_snd( k-1 ) - ( Cp / grav ) * thvm(k) * d_exner_zm(k) \\n\\nwhere\\n\\nd_exner_zm(k) = exner_zm(k) - exner_zm(k-1).\\n\\nWe need to iterate in the old hydrostatic code, because the hydrostatic code doesn't know pressure and therefore can't compute thvm.  Here we know pressure, and so we can compute thvm without approximation or interpolation.    \\n\\nThen we could set up M-PACE B using this initialize procedure, using a case-specific function to input profiles rather than sounding.in.  \", \n    \"reporter\": \"vlarson@uwm.edu\", \n    \"cc\": \"\", \n    \"resolution\": \"Verified by V. Larson\", \n    \"_ts\": \"1244210733000000\", \n    \"component\": \"clubb_src\", \n    \"summary\": \"sounding.in: Allow input of pressure as the independent coordinate as well as height\", \n    \"priority\": \"major\", \n    \"keywords\": \"pressure, height, altitude, sounding.in\", \n    \"time\": \"2009-05-14T15:00:06\", \n    \"milestone\": \"3. Refactor CLUBB\", \n    \"owner\": \"fasching@uwm.edu\", \n    \"type\": \"enhancement\"\n}\n```\n"},"comment":{"url":"https://api.github.com/repos/StrikeNP/trac_test/issues/comments/399282748","html_url":"https://github.com/StrikeNP/trac_test/issues/31#issuecomment-399282748","issue_url":"https://api.github.com/repos/StrikeNP/trac_test/issues/31","id":399282748,"node_id":"MDEyOklzc3VlQ29tbWVudDM5OTI4Mjc0OA==","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-22T00:21:12Z","updated_at":"2018-06-22T00:21:12Z","author_association":"OWNER","body":"**Trac update on 2009/06/03 at 21:23**\n\n* **vlarson@uwm.edu commented:**\n\nJoshua checked the code by running arm_0003 using height as the independent coordinate and then pressure.  The results almost overlapped in GrADS.\n\n\n* vlarson@uwm.edu changed resolution from \"fixed\" to \"Verified by V. Larson\""}}
{"actor":{"display_login":"StrikeNP"}}