diff --git a/.DS_Store b/.DS_Store index 59899e9..5008ddf 100644 Binary files a/.DS_Store and b/.DS_Store differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Data.pic b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Data.pic index bd1b2b0..9cc743b 100755 --- a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Data.pic +++ b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Data.pic @@ -46,7 +46,7 @@ Rp18 tp19 g11 I00 -S'\xfc\x14\xd84*S\x8f?\xa4\x0bi\x88\x1cE\x8f?i,3*\x93t\x8f?)\xfe-\x07\xfd\x8f\x8f?\x80V(\x9b\xdc2\x8f?\x19y_\x80\x05\x7f\x8f?|\xd6\x1d[2D\x8f?/&:b\x9c\\\x8f?.i\x0e\xc86.\x8f?"\xa2\x86}\xeeS\x8f?\x08E>\xf0\x8aq\x8f?Q\x0ct\xfc\xfc}\x8f?\x8d\xc8i\xe5dX\x8f?pAF\xae.W\x8f?qLg\xd2\xf4-\x8f?\xbf\x02\xf5~>T\x8f?\xdbL\xda\xd7\xd3c\x8f?\x11\xe3~F\xd5|\x8f?\xa9\x93]?:v\x8f?\x7f\xe8\x88\xbb2\x8e\x8f?\x98l\xb2\x15pG\x8f?\x0ekpe8\x8c\x8f?[\xe4\x9e\x83&\x92\x8f?\xeb\xb1\x1f\x13B\xbc\x8f?\x0bv\xce\xb0\xbc\x9a\x8f?\xffuo\xb1V\xa1\x8f?:<\x03\x06\xa4\xb3\x8f?m/\x18\'\x99\xc4\x8f?\xbd\xeb\x8d\xbe\x0b\xb1\x8f?:c\xcb\xac\x06\xdc\x8f?\xa85G\xa9\xfe\xe7\x8f?\x07\x02\x8d\x07\xce\xa7\x8f?\xc7\x80\xc8)Z\xf9\x8f?\xca\xc3R\x07\x87\xc5\x8f?/<\xf1Lo\xdc\x8f?\x1e\xe7dF\x08\xe9\x8f?w8\x9cg\x0c\xe2\x8f?\x9a\xfdQ\xa4j\xe1\x8f?\x95>\x17\x18\xc2\xe0\x8f?\xc1)\x1f\x08y\xe9\x8f?\x8d7\xe5\\m\xf9\x8f?$\xa8\x1bi\x1d\xea\x8f?S\x04\xfe.\x0e\xde\x8f?\x0c\x18\x15"f\xfa\x8f?o9\xfd\x91\x84\x00\x90?]\xb4\x83\xa6B\xe3\x8f?*k\xde\xce\xe6\x05\x90?D\xce\x96a\x04\r\x90?\xa0P\x13\xfc\xdc\x02\x90?\x95&\xb7\t\x17\xf8\x8f?p\x12gYk\x1e\x90?\xe1J\xee\x9c\n\x18\x90?\x11\xb4\x18\xce\x9a\xca\x8f?\x84\xf2\xe3>\x02\xc9\x8f?!\xd3\xc7\xe8(\xe7\x8f?\xb6*\xc7\xb5\x12\xd0\x8f?\x16\xfc7\xdcj\x91\x8f?\x16\xd5>\xceow\x8f?q\xc4\xc2\x0c\x08\xca\x8f?\x1f\xd54\xa0\x16\xc3\x8f?\xb3\xe7\xe7%s\xad\x8f?Q\xc1\x1fS[\xf4\x8f?B\x1b\xa2\r\x06\xeb\x8f?LO\xcfV~\xd2\x8f?\x17\xbfY\xee\xa9\xc7\x8f?\xb32\x19\x927\xea\x8f?\x8d\xa4C\x86\xa2\xc2\x8f?J\x9c\xeeL\xdf\xde\x8f?B\xf3\xbdx\xdb\xf3\x8f?\x11\x86\x87\x9c\x8b\xef\x8f?d\xea}+\x7f\xfa\x8f?\x10\xd6:\t\xfc\xbb\x8f?\x1d\x0f\xa7P`\xb6\x8f?Vqo\xb8=\xd8\x8f?\xafV\'\xd0"\xc2\x8f?syP\xfe\t\xeb\x8f?%&hA\xf1\xc5\x8f?\xf0\xa3b\x14\xad\xa4\x8f?\xb1\x17\x0b\xd0\xc6\xe9\x8f?2{\xea;D\xc9\x8f?\x11\xff{\xe3\xca\xe6\x8f?o\x90\x05\xcd;\xcc\x8f?X\xd6\xf6\x1aK\xf1\x8f?\xedV\xa5J\x8d\xe8\x8f?@\x8d\x85H\xca\xf1\x8f?\xef\xbf\xbf\x9a\xdd\xc5\x8f?\x82\x0b\xb7\x8a6\xf6\x8f?\xffO \x82\xd0\xd6\x8f?\xda%\x00\x8f\x87\xf9\x8f?\x1c\xb1f\x05G\x05\x90?\xaa\xa8\x89\xa0f\xff\x8f?\xc6i\x8c\xdd\\\xf4\x8f?\xb95\x99\xdb\xe0\x0b\x90?k\xa3\xe9MG\x02\x90?D\tp-`\xff\x8f?\x1d\xd1.\x87\xfe\xa5\x8f?p/\xb3\xa8\xc4\x93\x8f?\x91\xf1\xff\xad\xd3m\x8f?B6@[\xc3<\x8f?\xd7\xd2:\x90\xc5:\x8f?\x8e\x88c\xad4\'\x8f?\xb7\xb5\xdd\xdd\xa0\x19\x8f?\x8e\xca\xfacY2\x8f?v9S_O\x0c\x8f?\xd1\xc7]\xf6\xf3\xfc\x8e?4ds\x92X\xd8\x8e?\x95^\xe1\x80\\"\x8f?\x00\xf6\xe4\x0b\x940\x8f?68\xb9\xf9\xef\x19\x8f?\xdf",!\x8aO\x8f?q\xd7\xfe^\x1d\x85\x8f?9&\xe7\xe1\x08D\x8f?\xe9\xdd\x94\xc8\x06<\x8f?\x8a7\xdaY\xa7]\x8f?G\xef\xcb \x97c\x8f?VX\xc104r\x8f?\xcfo\xb0`\xc2`\x8f?\x94\xff(\xa2\x8d\x89\x8f?\xeb.#bZ\x93\x8f?\x94\x04\x10\xbb\xd5\xb5\x8f?\xbd\xa1\xa8\xca\xb5z\x8f?\x1c\r{\x08NY\x8f?]\xe8\xe4\x19\xcfu\x8f?\xb3\xf8\xedj:\xa1\x8f?\xfc\xa8\xa8\xd3\xc2\xa4\x8f?V\r\xf6\xe2\x8d\xe8\x8f?~^\xbe\xcei\xbe\x8f?\xac!f_\x0c\xc3\x8f?\x9f\x06\xde\xd2Q\xd4\x8f?D\xc1\xd5l\x8f\xf7\x8f?\xce\xdf\xc1t\x05\xf3\x8f?m]\xb4\x9f\x88\x03\x90?\xbf\x0e\xb5\x15\xfe\xfd\x8f?\xe6>\xad(\xe2\xfb\x8f?\xd0qUc\xa3\x0b\x90?e\xdf\x1c\xd3\x99\x1a\x90?i\xdf\xf7\xd5\xa1\x07\x90?\xc2Xv\x19\x81!\x90?-!\xaa\xc9]\x07\x90?\x01\x1f\x1e\xb5\x03\x10\x90?S\x1f\x19\xf3i\x1a\x90?s\xe5\xc6\xe9w\x1a\x90?p\xefVw\xfe\x15\x90?bz\xee\x91\xac\xd8\x8f?\xca6^\xcc\x91\x18\x90?\xc0]\xa0O\x0f\x07\x90?\xe8\xcb\xf8c\xf2\xfa\x8f?w\x1c#\xcd\xeb\x16\x90?+\x0c\x9a\x07q\x06\x90?\x98\x1eX\x8f$\xd8\x8f?\x8b\x07\xc9DX\xe3\x8f?\xcc\xad\xc6I\xce\xbe\x8f?\x85n\x8e\x94\xc5\x83\x8f?Q)\xbej\xf6\x89\x8f?\xbf=\xb6\xff\xe0_\x8f?\x8c\xad\xd1<\xba\x81\x8f?\x0e-n\x95\xf4b\x8f?\x13\xfdMZ\xbdc\x8f?\xben\xd4\x1f\xa63\x8f?\xc3a_#\xc68\x8f?e"\x92Zz\x13\x8f?\xdc\xf9:h\xa3\x03\x8f?\x07\x0c\x0c\xc8o\x10\x8f?\xc5U\x99=\xb76\x8f?\x80>z\r2:\x8f?>\x88Q+\xea\x06\x8f?`@4XWC\x8f?Y\x88\x8f\xc1\xe5V\x8f?)o\x88\xba\xa6P\x8f?\x05\xb8\xd0\x99XP\x8f?\xbb\xd1\xca\x1b\x86\x91\x8f?\x85k\x14\xb5\xb9\x97\x8f?\x14~\xf9_\x10\xe6\x8f?\xdf%p(\x11\xa1\x8f?\xdd\x1d\xb0\xdd\xae\x98\x8f?\xbd\xa1\x84\xae\x82\xb6\x8f?\x06\xb6\x02\xa9\xb0\xc1\x8f?\xa2\xe7\x15u]\xba\x8f?\x95w\xf0\xccI\xd4\x8f?\xbb\xf1%\x83\x14\xda\x8f?,\xb6z\xa2\xac\xac\x8f?\x84\xc3i\xf9j\xb4\x8f?\x93D&P\x8dz\x8f?\xa6\x00E\x87\xa8\xbc\x8f?B\xb9wazR\x8f?\xbd\x01R\x01\xf34\x8f?\x08\x8a4@\x9d\x12\x8f?K\x97W\xda\xa9\x0f\x8f?z\xad\xb6W\x18 \x8f?\xba\xa1=\x04V\x12\x8f?3\x15~\xd0\xe5\xf4\x8e?\x10\x0c\x1e\xe9P\xc4\x8e?b\xbcj\xd5`)\x8f?\xf2\xc6\x85/A.\x8f?F\xb9z\x9cW\x19\x8f?\xb7\x99\xed\xcet"\x8f?\xfa\x8c\xf2\t\xb8X\x8f?\x95\xc9\x1bW\xa7\x88\x8f?*q\xea\xbe\x9dy\x8f?\xf3\xc8[\x06J\x89\x8f?|\xb1\xf4\x9c\x9f\xa3\x8f?v\x19C\x1f\xe5\xc9\x8f?x+\xdb\xabg\x7f\x8f?g\xb1\xd6\x86\xf9\xa6\x8f?\x7f\x94\x1b?\xb8\xcc\x8f?\x0f\xd2R`r\xbe\x8f?R9\x8f\xaf\xd9\xa4\x8f?u\xcf\xbapf\xc5\x8f?\x06\x8c\x14\x8d\xfa\xb6\x8f?d\xe8\xb3Y|\xc3\x8f?\xd4\xbd\x10\xef\xbf\xe5\x8f?,7+,z\xbd\x8f?\x0fL\xb0N\x14W\x8f?\x9fz\xbdT\xd3\xe1\x8e?\xc2\xfd\xf4\x92-\xe7\x8e?\xd5\xfbWP_\x00\x8f?8\x86\xd9\x0c9\xdc\x8e?\xaeG}\xb9j\xb9\x8e?\xd7Zs\x00\x98\xb2\x8e?h\xe2]\x10T\x13\x8f?\x1d\xf0\xae~\xa8\xd6\x8e?\xe3\x85\xae\xa5\x01\x11\x8f?\xc4\xc4]\xde\x8f\xfe\x8e?:[\xf4\xa0\xc7\'\x8f?\x93Gw\xf8\xe6!\x8f?H\x1d\x81\x0f\x96[\x8f?\x96\xf6\xb5\xffaY\x8f?\xedpK\x8d\x93M\x8f?\xf8\x1e\xfcm\xa7\x94\x8f?\xb7*3\x0b"Q\x8f?oZ]\xd3\xb2m\x8f?I\xb0\x9c\xf9\xd5\xfe\x8e?I\xb5m)|\xc2\x8e?\x9b\x91\xb8y\x0f\xc6\x8e?\x9b\xbc\x0b\xcf2\xc0\x8e?\xaf4\rTR\xc1\x8e?\nI4!\x0f\xa6\x8e?\xa0\x04\xf7x\xcff\x8e?*2\xea\x80U\xd7\x8e?\xb0\x8f\x12\x83Y\xee\x8e?qOK\x99c\xdf\x8e?\x9a\xc7S\x1a\xc9\xf4\x8e?,\xe6\xf1Z\x07\n\x8f?\xed\xd6\xb8\xda\xde0\x8f?c\x18\xbb%E(\x8f?y{f9\xa9I\x8f?\xe1\x11\x08t\x8b\x15\x8f?\xf9,\xf8\n)&\x8f?\xf9P\xf3m7\x13\x8f?\xa3\xaa*"\xd9\xbc\x8e?i\xc5I\xb4\xb5\xfa\x8e?\xa9\xfc4\xc2\xe2\xe4\x8e?\xb8\xdb\xd5\x89\x16\xa5\x8e?8\xd4(\xaal\xc6\x8e?\x10\xb7\xe7\xa0\xc6\xaf\x8e?\xb5\xca\x17\xae\xdc\xdc\x8e?7@\x0e(8\x9d\x8e?3\xdf\x97\xccS\xb4\x8e?\xd9\x974\x86A\xdb\x8e?9\x9b78\xf4\xff\x8e?S#<\xc2\xe9\xf1\x8e?\x96\xe2d\xb4\xfb\xee\x8e?\x8d{5\xcaY\x1d\x8f?{\x00\x0f\xe5\x12\x0b\x8f?q\xfc\xe7\x87\xf6\xeb\x8e?5\xf7G\x9ef1\x8f?\xeaX\xcb\xf5}r\x8f??\x85\x81(\x87\xc5\x8f?\xdc\xcb\x8e\x85\xbdO\x8f?\x00\xcbw\x85\x8aQ\x8f?k.\x8a\xf1U\r\x8f?\x83V\xf7i\x19\x12\x8f?\xc8\x12\xc8AZ\xf5\x8e?\x0cm1\\6\xf9\x8e?>8\xe9\xfd\x9c\xdf\x8e?`\xa2.\xc0o\xfd\x8e?<\xe0\x03\xf3\x0f\xe9\x8e?\xe5\x06\xfbd6\xdd\x8e?\x8cW\x87\x974\xb5\x8e?\xae\x1d8\xa3\'\xcd\x8e?R\xc1\xab\x03{\xc2\x8e?\xf5\x12k\xc2\xe0\xea\x8e?\x11\x07\n\xed\x81\x97\x8e?>rO\xa6Q\xfa\x8e?Q\xf4S\xc5#\xf1\x8e?' +S'\xda|\x91\xcb\xac\x1b\xa3?o\xdc\xc9lc\x05\xa3?\xc3\x17}4\xbf\x15\xa3?\x03\x04\x82=\x19\x14\xa3?\x88v=\xfc0\x07\xa3?\x91\x04<\xaa\x98\x07\xa3?)g\xce\xd2\xd3\x0b\xa3?f\xa6\xda\xa9\xe3\xfe\xa2?d\xb0\xaf/&\x0b\xa3?~\xaa\xa1G\x0f\x02\xa3?\x11\x96&c\xee(\xa3?|&\x0e\xacB\x00\xa3?\xec\x9a\xbd}\x84\x15\xa3?\xdbUi\xc8\xc61\xa3?\x89dn|\xb6\x11\xa3?K\xb72P\xeb8\xa3?y`\xa8&\x922\xa3?\xe9;\x82^\xe4P\xa3?7\xbe\xc3\xe8\xb2N\xa3?\'\xca\x04\xeb\x18R\xa3?R>\xdb\x0fRx\xa3?\x01L$$\x01`\xa3?Q\x92ly\xf4\xa1\xa3?\x93^\xb31Q\xa3\xa3?\xb6^\xbfIc\x9b\xa3?\x9fv\x04\x01\xe3\xc1\xa3?|\'\xff\xb8\xe7\xaa\xa3?9\xa8\xcc!(\xd2\xa3?4\x15\x9a\x04\xaf\xc8\xa3?\x1d\x8dY\xb9\x0c\xd6\xa3?Q\xbc\x1d\xb1{\xd5\xa3?Dg\xe4D\xfb\xe0\xa3?[\xa4"\xa8\xd7\xd7\xa3?_\xee\xb4\xfc\x88\xfb\xa3?\xa6\xe5\xb2\x96.\xc9\xa3?\xe4\x8f\xf9\xd4}\xe8\xa3?\x95\x16;\xfa\x1e\xdd\xa3?Q4\x97\xc1\xde\xda\xa3?\xda\xc2Fr)\xdd\xa3?\xe7\x88\r\x16\xe7\xd8\xa3?`]\xa6\x10\xaa\xdb\xa3?~\xeaX\x98\xcd\xde\xa3?\x9b\x04\xfd\xd4m\xdf\xa3?\xc0\x1cK\xea\xf6\xfa\xa3?\x86h\'\xe5\xbb\xff\xa3?\xf1\xf3R\xf4=\xf3\xa3?b\xf8J{\xaf\x0f\xa4?\xa4\x98\xb7y\x8e4\xa4?U\x1bZ\xd9\x83/\xa4?%C\x9fyP\x13\xa4?\xf0\xc1j\x04\xeb_\xa4?\xa69\x88\x8e\x99J\xa4?\x00\x82Fud\xf6\xa3?R\x85\xaf\xd9\x05\xd0\xa3?\xaef\x11\xc7\x14w\xa3?\xd1;d\xfc\xc4N\xa3?P\xd8\x8a\x96JX\xa3?^\xcbi\xc5EN\xa3?\xc8\xbf\xe3ri_\xa3?\x84!\xc1\x19\x91p\xa3?\xaaf\xc4?%\x9f\xa3?\xbe\x990\x8es\xa4\xa3?J\xd1\xad\xef\xfe\xae\xa3?\xa9\xaa\x8c\x05-\xb1\xa3?\x80\xfdnD\x02\xcb\xa3?\x8d\x0e]\xb9\xc6\xc1\xa3?\xc2\xe7\x9c\xa4\xbc\xbe\xa3?\xb3Q\x8ez\xe4\xc4\xa3?\'\t\\\x90\xd2\xd7\xa3?i\xf0\xe3\xbf\x92\xca\xa3?\x13ME\xc9\x1d\xc7\xa3?\xaa<\x05\x04m\xd6\xa3?P\x98\xf7\x02\xb8\xc1\xa3?\x95\xd9=m\x96\xbe\xa3?)c\xbd\x00n\xb9\xa3?\x92e\xc5\xc2\xf0\xb3\xa3?\x1d\x9ccp\x13\xc0\xa3?c\xc6\x9c\xd7\xdf\xaf\xa3?\xca\xa1rf\xb9\xb5\xa3?\xae]h8y\x98\xa3?\xb7%u\xbe\xaf\x8f\xa3?\x805?\xb1\xc3\x93\xa3?~%Mk\xba\x7f\xa3?}\x88\x125\xc8\x96\xa3?\xb4m\xa8ca\x8b\xa3?W\xed\xe50\x16\x81\xa3?\xcbs\xd1Y\xf2\x85\xa3?`N\x1b\xbbE|\xa3?\x18v\x1e\xfc(\x8c\xa3?\xf9\xa0%\xf0E\x8f\xa3?\x84Wg\xa3\xfa\x99\xa3?\xa5x\x12\x94\xe1\xa3\xa3?\x9b\xa6\xf6\x10\xcb\x99\xa3?M\x14\x13\x84t\x86\xa3?b\xb8~]A\x80\xa3?\x82R\x1c=\xaeT\xa3?^\x8d\x9d|\x1d\x08\xa3?\x04\xc7\xff\xdc\x96\xd6\xa2?\xd6q\xa3\xc3\x9e\xbc\xa2?P~Z\xc6\xcc\xb1\xa2?\x9b\x11\xe0\x13g\x82\xa2?\xa4\xe04g\xb9g\xa2?\xcf\xac\xadH\x8e_\xa2?R\xf1\xe7\xec\x1fM\xa2?I\xe5\xc1\x03rD\xa2?\xcc0w\xcf\x05Z\xa2?`\xc4\x85\xe4\xefX\xa2?\xb6c\xbb\t*r\xa2?ra\x13J\xa6\x8e\xa2?\xfe\xbaW\x1a\xf7\x9f\xa2?\xde\xe8\xc9\x18>\xba\xa2?\x12P\'\\I\xce\xa2?\xca\xb5\x1a\x0b\x02\xd8\xa2?\x81\x04\x1aT\xa2\xe5\xa2?\x16\x85D\xd2\xe1\xf3\xa2?`\xcb\xea\x97\xcc\x14\xa3?\xe2~\xb5\xe7r\x07\xa3?\xbd\x0696\xff+\xa3?#\xff\t{\x16+\xa3?\xa9\x80\xa0\xe3\x8dA\xa3?\xd4d&AuW\xa3?\xc3\x171E\xc2i\xa3?\xf94\xa7-\x85p\xa3?3\xf4\xe0\xf0R\x86\xa3?\x03U\n\x8fR\xb9\xa3?BMQ\xfa\x02\xcf\xa3?t\xdb<:g\xdf\xa3?\xdd\x04c\x86\xb0\xf3\xa3?(K\x1f\xcc\xf7\xfe\xa3?\xc0\xe5\xc8\xd7\xc6$\xa4?\xae\xd9\x01\xa8_1\xa4?\x1d\xbd\x13\t\x80K\xa4?\xa7z\x9a\xf0\xce5\xa4?\x01\xc9\xc9\x88`j\xa4?*Y&!\xc91\xa4?b\xbcQ\xe0\x8dB\xa4?\xd0\x82\x98\x16cB\xa4?\xbb\x06\xccY&8\xa4?IZ\x0f\x88,S\xa4?\xfb\xa3\xca\xee\x1d9\xa4?Jj#\xe4@0\xa4?\x10\x8f\xfbl\x0c\x1a\xa4?[IN\x9d\x0e\xf5\xa3?b\xdb\t\xbd\xd3\x02\xa4?\x9d\xf7;\xb0\xab\x0c\xa4?\xc6@\x8bG\x00\x18\xa4?b\xe5\x10\x07; \xa4?\xe4\xa4\x0e\xc5\x1b\x0e\xa4?)>H\x87l\xe9\xa3?\x139p\xa4\xa9\xbe\xa3?\xee_\x98\x93\xc0\x85\xa3?c \x81\xd86^\xa3?\xef\x08\xad\xdc\x91\x16\xa3?;h\xc8\x98\xc9\xe4\xa2?\x9c\tHs\x0e\xc6\xa2?\tn\xe5\x12C\xb2\xa2?\x1c\x14\x04\xca\x19\x88\xa2?\xf8S\x83d\xede\xa2?]\xce\x9c\xaaPM\xa2?\xc0\xf5\xe73-N\xa2?\xf5\x06\xdf\x8a\xfbV\xa2?L\x9b\x8b\x8b\xfd[\xa2?\xd7\xf3\xd7o\xedb\xa2?T\xdb\xbfqHd\xa2?y\xe9H\xc6\xfc\x89\xa2?\xf7K\x80\x9e\x05\xb2\xa2?\x07\xa8/\x10\xc4\xc7\xa2?\x89\xe8=\x1e\xef\x02\xa3?\xe3\x8c\xe9\xe21&\xa3?\xbc\x0b\xd4\x9a\x15Q\xa3?O\xb5o\xc6\x80f\xa3?\x18\x18\x108\xe0s\xa3?\xc4JP\xb2\xd9o\xa3?|h.\xaa\xade\xa3?\x93\x12Qvx`\xa3?\x1a\xd7 \x86\xdcf\xa3?\x01\\=\x9e\xe8^\xa3?\x95,\x8cB2\x85\xa3?\xd4Mu\x7f\x05\xae\xa3?\xd5a\xef\x03E\x83\xa3?p\xe6\xcf\x1f*a\xa3?N\xb1[\xb6\x1aL\xa3?\xb3\x99\xf7\x03;7\xa3?E,\x82\x8e\x14\x05\xa3?>9\x98\xd1\xeb\xcc\xa2?OC\xf6\x81\xb3\x9d\xa2?\xbdM<<\x11\x83\xa2?\xaapdZ\x96=\xa2?\x10\x1f<^K\'\xa2?bCp\xd8\x1e#\xa2?t\xa4Rb\xb8/\xa2?p\x81%\x9b\xafM\xa2?\x01M\xa0\xe4#b\xa2?jZ\x11\x15\xb0\x85\xa2?\x90\x99\x8fg\xcb\x93\xa2?\xd9\xd712\xab\xb5\xa2?\xe1\\\xf0\x08\x95\xe7\xa2?\x0e\x8f\xf4\x05\x11\xf8\xa2?]>\xbbf\xbf\x06\xa3?\xf7lQ\xf2A+\xa3?\xe4,\xd8\x86\xb1P\xa3?\x8e\x1f\xdc|\xad]\xa3?\x14\x94\xa1-\xf6m\xa3?\x07\xceW\xc8Qh\xa3?!C/\xec|\x8f\xa3?\x18,\xd6\xd6\x7f\xa6\xa3?\xa2\xb8\x87\x1c\xe2\x91\xa3?\x154j\x92\xc7m\xa3?H\xedj\x82\x07`\xa3?\xa7\x1e\x16\x892u\xa3?TO\xc2\x85@s\xa3?\x1dV\xf5x\x949\xa3?\xeb\xb4+\xa9\xc8{\xa2?\xade\x136\xd2\x15\xa2?O\xe6v\xf32\xe8\xa1?Q\xc04c8\xec\xa1?1\xf0;\xc8>\xe2\xa1?\xdf\x10\x8f\xbf\xe9\xe7\xa1?\xa9v^J\xe7\xf0\xa1?r\xafi\xd1x\xf0\xa1?\tk\x1bX\xd8\x0b\xa2?\xffq\x9e6_A\xa2?"d\xbfT\x18{\xa2?\x8d\x9a\xb6Sy\xbb\xa2?\x93\xb2{\x85/\xc0\xa2?k\xf2\x7f\xca\xab\x8a\xa2?\xfb\tGu/r\xa2?\xaf\xbd;\x0f\x13u\xa2?c\xf4d1S\x84\xa2?\xbc\xd0.\xc3\xa9|\xa2?\xf69\xce\xd0\xbe{\xa2?u\xea\xb7/\xe7L\xa2?\x95\x06\xc7\x90\x94\x0e\xa2?\x0c\xce\xb5m\x02/\xa2?)\x99\xa2\xd7_\x96\xa2?(-\x1cZ\xb3\xf7\xa2?V\x10\xca\xb1lc\xa2?N\xfa\xd4\xd3\x9a\x02\xa2?\xd9\xb4\xbc\xb04\x1a\xa2?\xe7\x922\xb5)&\xa2?\x84\x86\x03Jd6\xa2?\xdc_5I\xfd_\xa2?\x11N%\xfb\xc7\x7f\xa2?\xa4\x1b\xadF\x02\xc3\xa2?\x88\x04K\xe6\x1d\xf1\xa2?)\x98\xaa?\xcc\xf6\xa2?T\x10\x9c\x02w\xb1\xa2?\x04\x12\xf2@\xaf\x8f\xa2?w\x99%@\xc6\xa6\xa2?\xc5\xb9}/H+\xa2?\xbf<\xc1u\xff)\xa2?\x8a\x12\xa7\xca\xa4n\xa2?B\xdfb4\xc9\xd3\xa2?v\x0e\xf0\x97\xaf\xf4\xa2?\x81\x1e\xab\xc7>\x16\xa3?\xa4\xc60\xddm0\xa3?\x88om?"\x0c\xa3?\xc8\xe2:\x83\x8aI\xa3?\x0e9\xf0A!\x81\xa2?\xe8\xe5?\xbbdx\xa2?Q\xa2\xd00[\x83\xa2?\xf2\xdei\xdcV\xa2\xa2?d\x0e\x99\xd6\x84&\xa3?\x15O\xd6(B\xb8\xa2?\xc0"\xc8\xc9X\x7f\xa2?K\xc1\x1d\xea\xc2\xca\xa2?\xb9Q?\x86q?\xa3?\xeeC\xc0\x18A\xde\xa3?\xd7\x81\xe7z*\xea\xa2?N\xe7\xb3\xc5\xe2\x04\xa2?\x16-\x04\x7f\xaf\xfe\xa1?LU\xa6\xbf=*\xa2?3\x0f\xdb\x01\xf5\xb8\xa2?\xa9/\xfa\xb4\x1d\xd7\xa3?Lz\x9a\xeb\x91\xfe\xa4?\x8d\xce5/\x15\xdc\xa4?\xd1\xbe\xe1\xf2y\xdc\xa2?\xf2\xae#\x17h\xbb\xa2?\xbd\x9cd\xa4\x99}\xa3?\x8a\xfd\xcbw\xe0L\xa4?\x90\x94\x16A\xd3\xf1\xa4?\xf7B\x93\xaa\xeb\xc4\xa4?w\x04#\x12&\xaf\xa3?S\xc9\x8a)\xe7;\xa3?\xfaT\xce\xe1\xf4U\xa3?' p20 tp21 bag1 diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Model.pic b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Model.pic index a7749e4..ca5093b 100755 --- a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Model.pic +++ b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/Model.pic @@ -46,7 +46,7 @@ Rp18 tp19 g11 I00 -S'\xb3r\xbd\xf0\xf2h\x8f?\x1f@\xd8\x85\xb2S\x8f?\x91\x9dK\xb05T\x8f?9\xd3\x9c\xe6%\\\x8f?\xcc1vY\x9b[\x8f?h]DE\x0co\x8f?\x9a\x94\x9clma\x8f?|\xf9\x07\xd0\x01F\x8f?8\xe4\x1d\x87\xcbP\x8f?N\xe0\xf7)2P\x8f?\xb2\xe9;\xac\xaem\x8f?\xae\x98\x9e#\x1dV\x8f?\x98\x00\xa7$CX\x8f?\xc2\xf2\xeb>Hj\x8f?2\x03\x0e\xc7\xcc\\\x8f?\x1c\x14\xb9k\xbcv\x8f?jk\xd8)Df\x8f?\xdb\xd1"#\xfbw\x8f?\x94\\\x9d\xf0\x89r\x8f?\xcf\x01\xe4\x91D\x7f\x8f?\xdf\x1bP\xbb\xe2\x8a\x8f?\x14:1/\xf7z\x8f?\x96Z\xa5\xa49\x9b\x8f?/\x96^l\xdf\x98\x8f?\\\xa3\xd50\xad\x92\x8f?I\xb8\xa7\x84X\xa9\x8f?\xa9\x85:\xc4\x19\xae\x8f?0c,\xb8\xfc\xbd\x8f?\xf9kEf2\xbe\x8f?8H\xf6I8\xcb\x8f?q\xc9\xd7\xfa\xb2\xc9\x8f?v\xf9\x8cS\xca\xd5\x8f?\xb4\x8b;\xc4\n\xce\x8f?\x8a$\xf7N+\xf4\x8f?v~\xa5aE\xcd\x8f?5\xce\x98\xaa\xcb\xe7\x8f?\x8f"3k\x1b\xd7\x8f?H\xf4\xd4\x02b\xe9\x8f?\xc3\xf3\xabJ[\xde\x8f?y~\xa1\x98\xd2\xde\x8f?\x91\xa1\xec\x04\x99\xe4\x8f?$Sk\xbc4\xde\x8f?\xb8L\xd3\xce)\xde\x8f?\xfe\xd7\x8f\xd4E\xec\x8f?\xbb1\x84 \xaf\xf4\x8f?n\xb1\xfb\xf8\xee\xe3\x8f?\xfdb\x9b\x9dJ\xf1\x8f?C`\xfb\x14\xbd\x06\x90?\xcc\xe6\x0e,\xda\x04\x90?& in^\xfc\x8f?c\xc0\x8d\xab\xd8\x14\x90?\x19\xd3\xd1\xea\xc4\x13\x90?\x8a\x07L\xc6D\xf8\x8f?\x95\xed\x96j\x8c\xe8\x8f?\xf8\xa0u\xec.\xbe\x8f?\x9c\x0f\xbaQ\x1b\x99\x8f?\x08\xb40\xc6\x1e\x9b\x8f?h1\x96\xc77\xa1\x8f?\r\x9c7\xbc\xd3\xa6\x8f?\x9a[#3b\xb9\x8f?u\x05\xcae\xc9\xcd\x8f?\xee\x02C\xad\xcd\xd3\x8f?\xa2J\xc4\xc8\xde\xd9\x8f?\xdfB\x8c\x0c\xcb\xe1\x8f?\x81Z\xca\x19\x8b\xdc\x8f?\xbb\xaf\x18\xc9\xa0\xe4\x8f?\x9dES\xad\x1d\xd5\x8f?\xdc\x045\xd9\xc0\xd7\x8f?\xf7\x86\xa8\xfe\x04\xe3\x8f?k\x92\xba\x17\xd7\xd7\x8f?\x94\xb2SP\xd0\xd1\x8f?\x07\x01\xd5h~\xe0\x8f?\xf40x\xe2|\xc9\x8f?_wKG?\xd1\x8f?\xaadE\x8ed\xc5\x8f?P\xa1\xff\x9f\xec\xc8\x8f?\x95]\xc0}}\xcc\x8f?\x89\xf8\xb3_\x8d\xc7\x8f?\x9b\x1aw\x88\xb6\xce\x8f?\x03\xa3@\xa1\xc4\xcb\x8f?J\xf4\xf6{y\xcc\x8f?\xd7\xd6\xbe\xf3\x86\xcf\x8f?\xaf\xb79"\x82\xcf\x8f?+\x83S}\x02\xd5\x8f?s\xb9\xb9\xc1\xde\xcf\x8f?\x8a\xb4\x7f#\x08\xcc\x8f?x\x94t%\x1e\xcf\x8f?\xdfR\x8e\x02\xdd\xd3\x8f?\xa2:*\x97D\xdc\x8f?N\xf1<\xdd\x07\xe8\x8f?T\xb6\xe7@\x81\xf1\x8f?\x88\xdf\x88i\xe6\xfc\x8f?\xe6Y\x1ask\xff\x8f?\xc48\xd0\x02z\xf7\x8f?\xefy#\xf7\x0c\xef\x8f?e\xd9\xe5\xcbS\xd0\x8f?\xab\xd0.\x84\x7f\x93\x8f?\xc8d\xc1V\x8d]\x8f?\xcb\x03\xbc\x1chD\x8f?\xceh\x8ejY;\x8f?\x8ba\xf9d\x8f%\x8f?1\xde\xdf\xdcD\x18\x8f?\x89\x9e\x03%\xa5\r\x8f?6h\xfadl\x07\x8f?\x9e\x16\x9a\xd3\xfb\x07\x8f?\xcc\xf6\x0c\x8aZ\x06\x8f?\xaf\x8f\x93\x01\x9c\x16\x8f?\xea\xe3\xb4&\x82$\x8f?`A\x0c\x0b\x04/\x8f?A\xdaqt\xa67\x8f?$o\xe2\x13\xb7G\x8f?\xda\x02\xf6r\x10N\x8f?\xb5K\x80\xfdMS\x8f?\x1fxM4p\\\x8f?\xe2\xda\xc5K\x8bb\x8f?\x07bqdUn\x8f?\x89[\x0f!\x93p\x8f?\xa8\x8b\xb4\xf3Ox\x8f?o\x98\xe7\xe4\x93}\x8f?\xc3\x80\xce\xd7\xe0\x8c\x8f?\r\'\xbe\xdb\x8f\x91\x8f?\x88\xfcf9\xe4\x89\x8f?^\xcd\xee\'\xd6\x8e\x8f?\x94\x94\\\x13C\xa0\x8f?BK\xa4\xc7\xb2\xb3\x8f?\x04T]\x17\xfe\xb8\x8f?8\xd1\xe5(\x8f\xc2\x8f?]\x01{\x04\xa4\xce\x8f?"\x94\xe0\xb5\x1e\xda\x8f?\\vp\x0e\xc4\xeb\x8f?X\xd2\xdeG/\xf6\x8f?GA\xef\xf2B\x03\x90?b\xbcKP\x8f\x04\x90?\xbd^\xcd\xd01\x11\x90?a0\x98N\x06\x08\x90?\xd6\xe6i\xa3\xab\r\x90?z\xe6\xb3n\xd7\x10\x90?\x89v\t\xdc\x1a\x0e\x90?~/9\xffB\x13\x90?\xf0\xe2W83\x0e\x90?L\xdb\x08\xdb9\x13\x90?\x0b\xdc-\xb3(\x14\x90?\xec\x1cJh\xc1\x08\x90?E\xaf\xfd\x9e\xc2\xff\x8f?qM\x14nM\x02\x90?j\xa4\xdf\x8c\t\x0b\x90?\x07S6[o\t\x90?\xc3\xbfs\xed\x81\x03\x90?\xe4I\x08\x01\xde\xf4\x8f?\x94U\xdde\x17\xd6\x8f?& a\x80\xfc\xb6\x8f?f\x07\xe6\t6\x9f\x8f?\x14L\xe7\xdc\xc3\x81\x8f?\x8bUL\xc5wq\x8f?\x17(\x15\xc5\x9dd\x8f?\xed\xa5\xbd\xe9\x89Z\x8f?\x8d\xc8q\x02\xc8O\x8f?\xdc\xf1\xa2c<@\x8f?\xc4KU\xfc9(\x8f?\xf0@!\xd9\xc5"\x8f?\xd1Z\xf0\xff\x0e"\x8f?\xea3\xc6\xde\xd9\x1e\x8f?\xe1\x91\x86\x8a\x8e\x1a\x8f?\x84Z\xeb\xf4\xb6\x18\x8f?\xeb\x8c\xd0\x81\x92\x1e\x8f?\x06\xbe\xf5\\4.\x8f?F\x05\xc8i,A\x8f?\xdc\xb7Z0iV\x8f?9xc\x91\x0ck\x8f?\xae\xd9\xd9a\xd6\x83\x8f?\x92\xd4+5*\x9d\x8f?\x92\x15gm\xa1\xb2\x8f?K]z\xde\xc2\xb8\x8f?\xf3\xb5z\x967\xbf\x8f?\xaa_bj\x96\xbc\x8f?r\xc6c\x86`\xab\x8f?\xcc\x8e\x1c\xfd\x9b\xad\x8f?E\\\x90\xddb\xbf\x8f?x\xf4\x99s\xd8\xd5\x8f?\xe8H\xad\x97\xc5\xbf\x8f?\xa4\x0cm\x08\xab\xb0\x8f?\x14\x0c\xeb\xbf\xd0\xa9\x8f?Q\xebZ\x0c\x01\x8e\x8f?\xac\xb0\xa3<\xdfz\x8f?\xb3YV\x90\xe4T\x8f?mk\xac.\xb7=\x8f?HS\xba\xfdL$\x8f?y\xd6\x16\x8e\x8d\xf9\x8e?6n>!=\xe6\x8e?\xc4yy\xed+\xe8\x8e?&\xa87\xb8\x85\xf1\x8e?v\x1fK\xb1\x1b\x01\x8f?\xba\x92\x05o\xed\x17\x8f?\xafh\x8a\x13\xda,\x8f? w\xd7\x98\xb59\x8f?\x0eq\x98\xb3\xdcH\x8f?\xf9\x81\xb7\xa6\xddc\x8f?\x16\xd9\t\xdc\xd4j\x8f?\x92\x0b\x1cT\xfdv\x8f?\x15\x15\xda\xf3\x88\x86\x8f?\xd9R!\x05p\x99\x8f?\xbc\xe6\x12*\xfa\xa8\x8f?\xecn\x81\xf3\x92\xb3\x8f?+vp\x824\xb9\x8f?\x82\x89\xd7\x00P\xd8\x8f?\xc8\xc12!\x0f\xe0\x8f?\x9a\x15\xc5\xad}\xd5\x8f?\xe0\x91\xdf\xa8\xb4\xc6\x8f?\xc0\x9f\x0c\x8e5\xc0\x8f?\x01w\xaa<)\xcc\x8f?\x04\x84\x06dx\xcc\x8f?\xda\x05\xfdY\xa5\xa9\x8f?\xd5N\x912\xe9?\x8f?\xc7/\xf2\xec\xea\t\x8f?3\x95,S\xdd\xf4\x8e?\xa2\xccfT\x88\xe4\x8e?\xb2\xd4\xfdZr\xd8\x8e?l;\xefo\x19\xce\x8e?\n$Z*\xa9\xc7\x8e?<\xe9\xaf_5\xd5\x8e?\xea\xdb\xfb\x8b}\xe5\x8e?z\xa1\x01`\xaa\xf5\x8e?\x04\x8a\xfd\xa2\x08\x06\x8f?\xcfg=!\xb5 \x8f?\xc2p\xa7\xf5g6\x8f?\x07\xf8"\x96XB\x8f?\x03\xc2\xa5\t\x11J\x8f?\x0c\xb2\x87\x95\xcdO\x8f?H\x1d\xb8\xb3\t`\x8f?\x0c\x0e\x03\x13%W\x8f?=!\xb5\xeb\x96J\x8f?\xdd\xd2o\xe7\xa1!\x8f?\x1a)h(\x12\xd9\x8e?\x7f\x87\xa9v\xa1\xa9\x8e?\x861,C\x8f\xae\x8e?\x01\xd2\x88p\xa3\xb0\x8e?\x0b\x1f\x97k"\xb2\x8e?\xe5\x9f\xcc\xe3\x92\xb4\x8e?q|\xfe\xd5\xcd\xc8\x8e?i\xbf\xa2V/\xdb\x8e?\xd25\x9c\xa0\xd5\xef\x8e?\xa8\xb5\xfa}\xab\x04\x8f?7c\x08\xa0\x14\x15\x8f?\xc3\x9el\xf2\x99\x1f\x8f?\x80\x94\xb8\xcf]"\x8f?\xde\xb8\x83\xd4\xae$\x8f?\xc2KY\xd2\x94$\x8f?L\x86\xef\xc0\x0b.\x8f?\xba\xdd\xe6\x0f \x1e\x8f?\xb4U\xc1\x983\xd7\x8e?\xf8~\'\xea\x9d\xc2\x8e?1\xb7\xb6\x07\xc8\xcf\x8e?\xa3\x02\xc9\xb0\xa2\xd6\x8e?\x08O\xc9k\r\xbc\x8e?\x02\x92\x10\xa3\xe8\xbd\x8e?\x17\xc4P\xbc{\xbc\x8e?\xe2\xa5\xe0\tp\xbf\x8e?7\xe0\'`\x83\xcc\x8e?p\x0f\xb0\x1e\xc9\xd6\x8e?\x9e\x12\x03*\xb3\xe1\x8e?\xd9!e\x8c\xd6\xeb\x8e?J\xfc*\x06B\xf9\x8e?\x8c\xf4\xdc\xe5M\xfd\x8e?\xc5:\xa4\x96\xae\x02\x8f?Q{\xb1\x1b\x87\x15\x8f?\x8f\x10\x19_\xb11\x8f?H\x95\xab\xb8\xf9^\x8f?\xe3\xe6\xb6l\xbc\xdf\x8f?\xed\x98\xbb\'|r\x8f?\xc7\xcd\'\xf0\xf9;\x8f?\x97\x10\x97K\xb5 \x8f?7\xa7\xa5\\\xef\x0f\x8f?\x10Z:\xb7\xc4\x02\x8f?\xdb\r\xdd\xa4\xc7\xf8\x8e?\xf7\xe1\xe9\x95\x7f\xf0\x8e?]\x1f\xbdS\x06\xdc\x8e?\x8e\xa1\x17P\x98\xd4\x8e?\xcc\x19\x17\x8a\t\xd4\x8e?\x9e\xdf\x80\xa6Q\xd0\x8e?iG\x1e\xd55\xce\x8e?"\x02\xf5\xf8\xcd\xd2\x8e?\xeaJ\xa0}\xef\xd8\x8e?\x94\xf5.\x1cA\xe0\x8e?\xbc\x9c\x97\x17\xb0\xe8\x8e?i\x16-\x8e\x86\xf2\x8e?' +S'\xc6\x96\x02E\x96\x1f\xa3?\xdf\x98qP\xab\x06\xa3?\x17\x1e\xf4\xe33\r\xa3?\xaba\xbe\x08\xcc\x12\xa3?\xfc\xee\xebj\xed\x07\xa3?\xd9\xf0"\xd1\xa4\x0e\xa3?\xf4\xa8kdi\t\xa3?2F\x1aQ\xb4\xf7\xa2?\xaa\x1d\x8b\x1cI\r\xa3?\x97\x8e\xfe\x05\x8a\xfb\xa2?p\xfa\xd0\xa0\xe5"\xa3?\xa1\xe5\xbbf"\x07\xa3?\xcb{y\xe8\xe6\x0f\xa3?\xa9\xd4\xf4\xfdy)\xa3?\x02sq\xc8\xb1\x14\xa3?\xd0\xfd\xcfh\xa4;\xa3?\xf8\xbf\\\xae\x1d4\xa3?\x80,\x9ef\x80U\xa3?\x0b\xe0\x16\xd8\xcaH\xa3?\x98\x02+\xd8\xdaU\xa3?\xc8/\x8a\n\xd2y\xa3?\x17\x98\n\\\xd6g\xa3?.\xfd2\x07%\x9b\xa3?\xf6r-\x81\xe5\xa0\xa3?rJ\x7fC\xd6\x97\xa3?\x9cU?\x0f\xe4\xb7\xa3?|(J\x03\xb3\xbd\xa3?TZ\x87\xc5\xe1\xda\xa3?\x8c"~O(\xcd\xa3?\x9f9V\xbec\xd8\xa3?\xc3Ki\xa6\xe1\xd0\xa3?\x91)CJ\x97\xd8\xa3?o\x19La\x08\xd0\xa3?\x82\x04\x14\xaef\xfa\xa3?\x07\xc7\xf4N:\xb8\xa3?\x00?o\xd5\xc2\xe0\xa3?i[\x14\xa9\xb8\xd9\xa3?\xc7}BV\x8a\xde\xa3?\x10\xf9Q\x13\xe9\xd7\xa3?\xea\xd5\xda\x82\x8a\xd7\xa3?\xa7N\xb6\x9dI\xe4\xa3?K:\xb7\x89^\xe1\xa3?+}&\xc0J\xe0\xa3?\x12\x0c\xbe""\xfb\xa3?G\x8d8\x11\xc1\x05\xa4?\xdc\xf9\x1d\xaeQ\xf2\xa3?\xf4\x9e\xe9\xb01\x11\xa4?\x0eZ\xe5\xf2q-\xa4?`x\x08\n\xa90\xa4?\x9e\x10\xa8I\xbe\x15\xa4?\x11\x8dYc\xcfe\xa4?\xc9\x11\x95\x9enA\xa4?\x03\xa0|J\xc1\xee\xa3?B\x90vZ\xd7\xc6\xa3?\\\xbcya\xbet\xa3?\xa1\xa8\x8d\x90\x9eO\xa3?y\xfa\xfb\xf1\x08W\xa3?\xe0\xf6\xf9\xe57R\xa3?\x18\xea\xf17g_\xa3?\xf0\xe5\xb4\x06\x8ew\xa3?\x9d\xbe\xb6\xf1\xe6\xa2\xa3?X\xf7H\xebr\xa7\xa3?\xfd\xac+\xcd\xf1\xad\xa3?\xea\xe9cM\xb6\xbc\xa3?\xdak/\x0b*\xbb\xa3?\x97\xb6\xe2\xbe%\xca\xa3?\xb7}\xd7\xb7\xc1\xbc\xa3?o\xa5L\x10\x88\xc3\xa3?\xbc7\xcf\xc8\xe2\xd7\xa3?0N8\x16\x91\xca\xa3?\x90\xb5\x8f\xb0\x17\xc4\xa3?\xeb\xcdC\x11\x1b\xd7\xa3?8\xc1\x98K\x0c\xbe\xa3?\x80\n\xa0x\xfc\xc4\xa3?\xf2\xac\xdf\xe0z\xb6\xa3?\x1a\xc8\xd0\xae\x98\xb8\xa3?\xee J\x86\xd3\xb3\xa3?H\x85\xab\xf5\xb9\xac\xa3?\xa1.\x91\x0c\\\xb2\xa3?\xe2t\xb0Y\x00\xa1\xa3?\x05\xf2\x8d\x92u\x8e\xa3?x\xfc\x16\x12"\x8e\xa3?\x946\xcbc\x06\x8a\xa3?\xfa\xe3\xc0\xa5f\x90\xa3?\n|\xe6\xe4}\x82\xa3?\r\x15\xa0Vu\x84\xa3?.1\x98o\xae\x8f\xa3?\x08u\x90E\xff\x81\xa3?\xa2\x84~/A\x8b\xa3?B\x1a6\xe7\xed\x95\xa3?\x9d\xac\xfb\x89l\x9a\xa3?\x19\x127[\xd7\xa0\xa3?\xd2\x89\xb1\xd9\x94\x9b\xa3?\xaf\xb6f\x99\xa7\x83\xa3?\xf0 G\xc1\xabt\xa3?\x95\x8f\xa2\x0f\x97W\xa3?\x8cN\xa1p\xf2\x0c\xa3?\xbb\xba\x04v\xdf\xd6\xa2?\x9b\x17\xbd\xff\xe1\xb8\xa2?\xc4dl\xad\x1a\xb2\xa2?G_\xea3\t\x89\xa2?\xb0\x7f\x991~l\xa2?\x06l\x0e\x8c\xac[\xa2?J\x9am\xc7wM\xa2?A(\xa9\x93\xa8K\xa2?\x93>>\xc2\x92T\xa2?\x9eW\x11R\xfcg\xa2?p\xf6\xa4\xd1\x0fw\xa2?E\xe5{+\xe4\x8d\xa2?\x1a\xb3\xa8\xbb_\x9f\xa2?\xa8b\xfc\xf1\x08\xbf\xa2?o\xef\x98V\xe8\xc9\xa2?b\xc0\xb65\x87\xd7\xa2?\xb0L\xf0\x16\x9f\xe5\xa2?\xda\x06\xcf(\x86\xf8\xa2?\xc4\x8as;.\x08\xa3?\xc6\x85\x173-\n\xa3?\x8e\x08\x186j!\xa3?f\x9b\xcd\xd6\xac1\xa3?1\x16"YrA\xa3?\x9b]\x0c+\xd5K\xa3?\xb1\x8eP8\x83e\xa3?\xa83D\xbbov\xa3?\xdc\x1c\xa7\xf0\xdf\x91\xa3?\x08\xb4\xeam\xf0\xba\xa3?.\xeb\xc1Q\xff\xc8\xa3?HA\xda\xb8]\xdf\xa3?>~e4?\xf1\xa3?\xa9\xc6z\xc8u\x06\xa4?\xc4\xa3\xd0my#\xa4?K\xf6[{\x08.\xa4?\xd8+2\x0f\xfbE\xa4?\x8el\x7f.\x85@\xa4?T\x80\xf0\'\xf7e\xa4?\xb2\xb8\xfb\x92\xdf6\xa4?V$\xb9\xe2\xd1H\xa4?\x97/(\x9d\x80T\xa4?_\xad&Im@\xa4?mK\xb8\x8a-J\xa4?\xc7b\x8eP?7\xa4?u\xa1\\\xf0^0\xa4?\xe7&\xc8\xf5\x89\x1c\xa4?\x11\xad\xb87\x1b\xf7\xa3?*J\xb1A\xdd\xfd\xa3?\x16\x86\x11EZ\x01\xa4?\x9c;a\xfaV\x1b\xa4?\xf3\x97\xe4|\x0b\x18\xa4?\xbejR\xed\xd5\x07\xa4?\x93\x1c.<)\xf1\xa3?\x0c\x08\x8b\x15\xde\xc1\xa3?\xaa<}J\xc5\x88\xa3?O*\x13x\xd3X\xa3?\x99\x07\x08\xd0\xc6\x19\xa3?\xce\xc6\xba}$\xe9\xa2?\xd4,\x84$\xde\xc4\xa2?\x10\x9c\xf7\xe4\xb0\xad\xa2?\x8f\x0c\xdd{n\x8f\xa2?\x88\xf2\x88\xbato\xa2?\xb2\xcb\xe2.\x87J\xa2?78\xda3lH\xa2?(q\x1e\x00UR\xa2?\xdb\xd6\xaa\x8d\x85X\xa2?LO\x88\x8f\x9de\xa2?%\xe4\xee\xe6\xd9o\xa2?\xf2\xc5\x0c\x81C\x8c\xa2?\xaa\x1c\xe9\xe7\xc4\xaf\xa2?-\xd6\x99\x89\xa8\xd7\xa2?%\xd7%$R\x01\xa3?C\x95\x18K\xfc%\xa3?\x80\xa3\x9b<\x9bG\xa3?\xf8\xd6p\xc0\xdde\xa3?\xac\xe2\x82\x8c\x88v\xa3?\x80\xe7NH\x0bo\xa3?\xc2\xb1\xdd\x88Ui\xa3?\xc8\x02\xd3\x9e\xc8h\xa3?\x14\x0c\x8aI\xe6^\xa3?H\xdb\xbf\xf9!d\xa3?\x84\xaa\x03\xf1@\x81\xa3?\x07\xb8\xf4\xadm\xa9\xa3?v\x97\xb5\x86\x02\x85\xa3?\x86$\xae\xc9\xfd`\xa3?\xfc*\xbc\x84\xcfV\xa3?\xf9\x82\x0c\xc6\xd6-\xa3?\x15\xbeh6\xdd\n\xa3?C"G\xca\xf9\xce\xa2?~Hb\xc9\xd9\xa4\xa2?\xd5\x972\xa2\x9c\x83\xa2?\xb3]k6c?\xa2?\xd6\xa3J\xc0 !\xa2?\x83&\xff\xb3%&\xa2?\xd9\xc1>\xc6\xa0:\xa2?\x1b${ \nP\xa2?\xf1Ty\xd0\x8bk\xa2?2\xf8\xc1\x0c\xf2\x85\xa2?F\xbc\x04=\x0e\x9d\xa2?\xf3_4?\xc1\xb0\xa2?\x85\xa6\xc4\xc2\xd7\xdc\xa2?]\xb4\xe7\xd34\xf5\xa2?\x1dR8\xd7j\x10\xa3?\x06\xea_\xc9$*\xa3?\x8b\x9b\n\xa6\xb2H\xa3?\xe4I\x88J\xfe_\xa3?"\x1a3?\xa8o\xa3?\xa7$\x01cgp\xa3?H\xa0J\xe7Y\x98\xa3?\x89<\xb0\xf3\x95\xa4\xa3?hhs\xaaZ\x89\xa3?\x18\xbaH\xfa\xcbn\xa3?U\xa9x\xde-[\xa3?\\\xe8\xd2\x8b\x1en\xa3?\xae\xa1\xe9\xdc\x9bi\xa3?\xcf>}\x88a=\xa3?\xc2]3JS}\xa2?p\xe9\xe4\xb8\\\x0f\xa2?o\xa9>q\xae\xf2\xa1?;)G\xa2g\xe8\xa1?WX\r\n4\xe3\xa1?\xc7\xa6\xcf2M\xe7\xa1?\x10\x04A\xf6\x87\xe7\xa1?\xd1\xea\x8d\'\x82\xf6\xa1?\x03;z\xb4\xb9\x0b\xa2?\xba#\n_\xe78\xa2?\xd8\x9f\x89\xe8\x10u\xa2?\xc4e\xa5\xc2\x1e\xb8\xa2?{\x1b\xa7-_\xc1\xa2?\xf2\xb3CS\x7f\x83\xa2?iq@\xc9\xd9w\xa2?\x00\x82\xc01\xb8{\xa2?\xca\xbe\xeeq\xf2\x86\xa2?\x16\x1f1+5y\xa2?;\x84\xffB\xcbx\xa2?DG\\\nYO\xa2?\xbf\xdc\xe5\x9b\xf6\x11\xa2?7\x9a\xb83b+\xa2?\x89\xeb\x85\xfd6\x9e\xa2?&\xac\xd5\x9f\x00\xef\xa2?O\xb8sF\x8d]\xa2?\xc5\x12\x1bl\xc3\x07\xa2?Z@5ft\x1d\xa2?2\x1e\xd7XV!\xa2?\xce\xc7\x85\x8bb:\xa2?2Ny\xd2\x1b_\xa2?\xb7X\xf6\xe7Py\xa2?\xcc\xfa\xd3\xfa\x0f\xbc\xa2?\x86\x9a6\x9d\xb0\xf9\xa2?T7\xcd\xd7\xa4\xef\xa2?\xbd\xbc\x8f2o\xab\xa2?3\xd6ZQ\x91\xa0\xa2?S\x8d\x1c5\xa8\xa0\xa2?\xe0q\x05\x1a,,\xa2?\xe1\xf7/\x97\xc7&\xa2?\x19;\xb67\x93l\xa2?\xcb`R\x8e\xb0\xd0\xa2?\x8c\xfdp\xa2\xe4\xf5\xa2?\xe9\xce\x07O\xab\x1b\xa3?\xbe#\xda\x94\x80\'\xa3?\x97\xe3&\x05\xd5\t\xa3?\xe1\x7f4F\x86A\xa3?\x80+?Z\xd9\x83\xa2?\xc3i\x9cx\x00|\xa2?B\xb5g\xee\xee\x8c\xa2?yZH\xb9`\x9f\xa2?\xc0*\xc5\xd5\xbf!\xa3?Z\xf4D\xa6\xd5\xc3\xa2?\x82?\xabg\x15|\xa2?\x91\x98\xa8\xa2\n\xce\xa2?[\xdaG:\xcf9\xa3?d\xef\xa1`u\xd9\xa3?\xf0\x90!\x83^\xf2\xa2?z\xf6X\xe9\xc9\x00\xa2?j\xe1\x12\xf3\xdd\xf8\xa1?+\x9f#R$(\xa2?\xd2\xf7L\xe9i\xbb\xa2?\xe2\x07\x85\x92\x87\xda\xa3?\xc5k\xb9C\x92\x05\xa5?K0\x0cg\xee\xdd\xa4?F\xf4\xb8\x7f\x98\xdc\xa2?C\xddx%\xea\xc3\xa2?\x948\xc6\x01\xebv\xa3?\r%B`\xf2H\xa4?\x80\xaa+d\x84\xf0\xa4?\xf9\x8a\x8a*\x00\xb8\xa4?f/\xb7>\xcb\xa5\xa3?\xdb,E7\'3\xa3?\xda\x07\xd0x\xd3Y\xa3?' p20 tp21 ba. \ No newline at end of file diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/_tran_module.so b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/_tran_module.so index fb2a2bc..a466138 100755 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/_tran_module.so and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/_tran_module.so differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau.o b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau.o old mode 100755 new mode 100644 index 91857bb..83896e1 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau.o and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau.o differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau_wrap.o b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau_wrap.o old mode 100755 new mode 100644 index 6d1536e..5284938 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau_wrap.o and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_tau_wrap.o differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_xsects.o b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_xsects.o old mode 100755 new mode 100644 index 03f9e5f..509435a Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_xsects.o and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/build/temp.macosx-10.6-x86_64-2.7/_xsects.o differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/cea2.x b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/cea2.x index e1c71a5..e25a997 100755 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/cea2.x and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/cea2.x differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/ctran.pyc b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/ctran.pyc old mode 100755 new mode 100644 index 13498d0..1e8ab96 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/ctran.pyc and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/ctran.pyc differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/fm.pyc b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/fm.pyc old mode 100755 new mode 100644 index f36b2f5..30dc87f Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/fm.pyc and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/fm.pyc differ diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/make_spec.py b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/make_spec.py index 4c0bed3..9d30c0b 100755 --- a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/make_spec.py +++ b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/make_spec.py @@ -4,143 +4,165 @@ rc('font',family='serif') -#the parameters -#planet/star system params--typically not free parameters in retrieval -Rp= 1.359# Planet radius in Jupiter Radii--this will be forced to be 10 bar radius--arbitrary (scaling to this is free par) -Rstar=1.155 #Stellar Radius in Solar Radii -M = 0.690 #Mass in Jupiter Masses -#TP profile params (3--Guillot 2010, Parmentier & Guillot 2013--see Line et al. 2013a for implementation) -Tirr=1200. #terminator **isothermal** temperature--if full redistribution this is equilibrium temp -logKir=-1.5 #TP profile IR opacity controlls the "vertical" location of the gradient -logg1=-1 #single channel Vis/IR opacity. Controls the delta T between deep T and TOA T -#Composition parameters---assumes "chemically consistnat model" described in Kreidberg et al. 2015 -logMet=0. #. #Metallicity relative to solar log--solar is 0, 10x=1, 0.1x = -1 used -1.01*log10(M)+0.6 -logCtoO=-0.26 #log C-to-O ratio: log solar is -0.26 -logPQCarbon=-5 #CH4, CO, H2O Qunech pressure--forces CH4, CO, and H2O to constant value at quench pressure value -logPQNitrogen=-5 #N2, NH3 Quench pressure--forces N2 and NH3 to "" --ad hoc for chemical kinetics--reasonable assumption -#generic cloud parameters--Rayleigh Haze(amplitude, slope) and Hard Cloud Top Pressure -RayAmp=0. #log Rayleigh Haze Amplitude (relative to H2) -RaySlp=4. #haze slope--4 is Rayeigh, 0 is "gray" or flat. -logPc=1.5 #log of the hard gray cloud top pressure - - - -#seting up input state vector. Must be in this order as indicies are hard wired in fx inside fm - # 0 1 2 3 4 5 6 7 8 9 10 11 12 -x=np.array([Tirr, logKir,logg1, logMet, logCtoO, logPQCarbon,logPQNitrogen, Rp, Rstar, M, RayAmp, RaySlp,logPc]) - - -#calling forward model -#thermochemical gas profile scaling factors -# 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 -#H2O CH4 CO CO2 NH3 N2 HCN H2S PH3 C2H2 C2H6 Na K TiO VO FeH H H2 He -gas_scale=np.array([1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.]) #can be made free params if desired (won't affect mmw) -y_mod,wnocrop,atm=fx(x,gas_scale) #returns model spectrum, wavenumber grid, and vertical abundance profiles from chemistry - - -#read in from noise model here -err=np.zeros(len(wnocrop)) -err[:]=50*1E-6 -y_meas=np.zeros(len(err)) - -#adding gaussian noise -for i in range(y_meas.shape[0]): y_meas[i]=y_mod[i]+np.random.randn(1)*err[i] - -#computing chi-square of random noise instance -print np.sum((y_meas-y_mod)**2/err**2)/len(y_meas) - - -#dumping pickles--model then noised up data data -output=[1E4/wnocrop, y_mod] -pickle.dump(output,open("Model.pic","wb")) #spectral model to be noised up by instrument noise model - -output=[1E4/wnocrop, y_meas,err] #noised up "synthetic" spectrum -pickle.dump(output,open("Data.pic","wb")) - - -#plotting------------------------------------------------------------------- -ymin=0.98*np.min(y_mod)*100 -ymax=1.02*np.max(y_mod)*100 -#ymin=2.9 -#ymax=3.05 - - -#plotting structure (TP, mixing ratios, cloud)----------------------- -# atm=np.array([P,T, H2Oarr, CH4arr, COarr,CO2arr, NH3arr, N2arr, H2Sarr, HCNarr, C2H2arr, C2H6arr, H2arr, Naarr, Karr]) -P=atm[0] -T=atm[1] -H2O=atm[2] -CH4=atm[3] -CO=atm[4] -CO2=atm[5] -NH3=atm[6] -Cld=atm[-1] - -#mixing ratios -fig2, ax1=subplots() -ax1.semilogx(H2O,P,'b',ls='--',lw=2,label='H2O') -ax1.semilogx(CH4,P,'black',ls='--',lw=2,label='CH4') -ax1.semilogx(CO,P,'g',ls='--',lw=2,label='CO') -ax1.semilogx(CO2,P,'orange',ls='--',lw=2,label='CO2') -ax1.semilogx(NH3,P,'darkblue',ls='--',lw=2,label='NH3') -ax1.set_xlabel('Mixing Ratio',fontsize=20) -ax1.axis([1E-9,1,100,1E-6]) - -#TP -ax2=ax1.twiny() -ax2.semilogy(T,P,'r-',lw='2',label='TP') -ax2.set_xlabel('Temperature [K]',color='r',fontsize=20) -ax2.axis([0,3000,100,1E-6]) -for tl in ax2.get_xticklabels(): tl.set_color('r') -ax1.set_ylabel('Pressure [bar]',fontsize=20) -ax1.legend(loc=4,frameon=False) - - -#plotting spectrum------------------------------------- -fig1, ax=subplots() -ax.plot(1E4/wnocrop, y_mod*100.) -xlabel('$\lambda$ ($\mu$m)',fontsize=18) -ylabel('(R$_{p}$/R$_{*}$)$^{2} \%$',fontsize=18) -minorticks_on() -ax.errorbar(1E4/wnocrop, y_meas*100, yerr=err*100,xerr=None,fmt='D',color='black',alpha=0.5) -ax.set_xscale('log') -ax.set_xticks([0.6,0.8,1,3,5,7,9,11]) -ax.axis([0.6,11,ymin*0.98,ymax*1.02]) -ax.get_xaxis().set_major_formatter(matplotlib.ticker.ScalarFormatter()) -ax.tick_params(length=10,width=1,labelsize='large',which='major') -show() - - - - -##########################TESTING, PLAYING TO GAIN INTUITION###################################### -#------twiddeling--testing cloud param values---can test other params to by forcing x index or gas_scale values to new value -# 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 #19 gases (PH3, N2, and H are not opacity sources) -#H2O CH4 CO CO2 NH3 N2 HCN H2S PH3 C2H2 C2H6 Na K TiO VO FeH H H2 He -gas_scale=np.array([1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.])*0. #all gases zeroed out--just want to test effect of cloud - - -pdb.set_trace() -for i in range(5): #looping trhough 10 param values - x[5]=-5+i*1.5 - gas_scale=gas_scale*0. - y_mod2,wnocrop2,junk=fx(x,gas_scale) - ax.plot(1E4/wnocrop2, y_mod2*100.,color='blue',alpha=0.2+0.2*i) - gas_scale=gas_scale*0.+1. - y_mod2,wnocrop2,junk=fx(x,gas_scale) - ax.plot(1E4/wnocrop2, y_mod2*100.,color='blue',alpha=0.2+0.2*i) - -show() -pdb.set_trace() - -#truth vector -#pickle.dump(x,open('truth.pic','wb')) - - -pdb.set_trace() - - - - - +def spectrum(Rp= 1.93, Rstar=1.2, M = 0.5, Tirr=1700, logKir=-1.5, logg1=-1, logMet=0., logCtoO=-0.26, logPQCarbon=-5, logPQNitrogen=-5, RayAmp=-5., RaySlp=4., logPc=1.5): + # #the parameters + # #planet/star system params--typically not free parameters in retrieval + # Rp= 1.93# Planet radius in Jupiter Radii--this will be forced to be 10 bar radius--arbitrary (scaling to this is free par) + # Rstar=1.2 #Stellar Radius in Solar Radii + # M = 0.5 #Mass in Jupiter Masses + # #TP profile params (3--Guillot 2010, Parmentier & Guillot 2013--see Line et al. 2013a for implementation) + # Tirr=1700 #terminator **isothermal** temperature--if full redistribution this is equilibrium temp + # logKir=-1.5 #TP profile IR opacity controlls the "vertical" location of the gradient + # logg1=-1 #single channel Vis/IR opacity. Controls the delta T between deep T and TOA T + # #Composition parameters---assumes "chemically consistnat model" described in Kreidberg et al. 2015 + # logMet=0. #. #Metallicity relative to solar log--solar is 0, 10x=1, 0.1x = -1 used -1.01*log10(M)+0.6 + # logCtoO=-0.26 #log C-to-O ratio: log solar is -0.26 + # logPQCarbon=-5 #CH4, CO, H2O Qunech pressure--forces CH4, CO, and H2O to constant value at quench pressure value + # logPQNitrogen=-5 #N2, NH3 Quench pressure--forces N2 and NH3 to "" --ad hoc for chemical kinetics--reasonable assumption + # #generic cloud parameters--Rayleigh Haze(amplitude, slope) and Hard Cloud Top Pressure + # RayAmp=-5. #log Rayleigh Haze Amplitude (relative to H2) + # RaySlp=4. #haze slope--4 is Rayeigh, 0 is "gray" or flat. + # logPc=1.5 #log of the hard gray cloud top pressure + + + + #seting up input state vector. Must be in this order as indicies are hard wired in fx inside fm + # 0 1 2 3 4 5 6 7 8 9 10 11 12 + x=np.array([Tirr, logKir,logg1, logMet, logCtoO, logPQCarbon,logPQNitrogen, Rp, Rstar, M, RayAmp, RaySlp,logPc]) + + + #calling forward model + #thermochemical gas profile scaling factors + # 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 + #H2O CH4 CO CO2 NH3 N2 HCN H2S PH3 C2H2 C2H6 Na K TiO VO FeH H H2 He + gas_scale=np.array([1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.]) #can be made free params if desired (won't affect mmw) + y_mod,wnocrop,atm=fx(x,gas_scale) #returns model spectrum, wavenumber grid, and vertical abundance profiles from chemistry + + + #read in from noise model here + err=np.zeros(len(wnocrop)) + err[:]=50*1E-6 + y_meas=np.zeros(len(err)) + + #adding gaussian noise + for i in range(y_meas.shape[0]): y_meas[i]=y_mod[i]+np.random.randn(1)*err[i] + + #computing chi-square of random noise instance + print np.sum((y_meas-y_mod)**2/err**2)/len(y_meas) + + + #dumping pickles--model then noised up data data + output=[1E4/wnocrop, y_mod] + pickle.dump(output,open("Model.pic","wb")) #spectral model to be noised up by instrument noise model + + output=[1E4/wnocrop, y_meas,err] #noised up "synthetic" spectrum + pickle.dump(output,open("Data.pic","wb")) + + + #plotting------------------------------------------------------------------- + ymin=0.98*np.min(y_mod)*100 + ymax=1.02*np.max(y_mod)*100 + #ymin=2.9 + #ymax=3.05 + + + #plotting structure (TP, mixing ratios, cloud)----------------------- + # atm=np.array([P,T, H2Oarr, CH4arr, COarr,CO2arr, NH3arr, N2arr, H2Sarr, HCNarr, C2H2arr, C2H6arr, H2arr, Naarr, Karr]) + P=atm[0] + T=atm[1] + H2O=atm[2] + CH4=atm[3] + CO=atm[4] + CO2=atm[5] + NH3=atm[6] + HCN=atm[9] + C2H2=atm[10] + Cld=atm[-1] + + #mixing ratios + fig2, ax1=subplots() + ax1.semilogx(H2O,P,'b',ls='--',lw=2,label='H2O') + ax1.semilogx(CH4,P,'black',ls='--',lw=2,label='CH4') + ax1.semilogx(CO,P,'g',ls='--',lw=2,label='CO') + ax1.semilogx(CO2,P,'orange',ls='--',lw=2,label='CO2') + ax1.semilogx(NH3,P,'darkblue',ls='--',lw=2,label='NH3') + ax1.semilogx(HCN,P,'r',ls='--',lw=2,label='HCN') + ax1.semilogx(C2H2,P,'c',ls='--',lw=2,label='C2H2') + ax1.set_xlabel('Mixing Ratio',fontsize=20) + ax1.axis([1E-9,1,100,1E-6]) + + #TP + ax2=ax1.twiny() + ax2.semilogy(T,P,'r-',lw='2',label='TP') + ax2.set_xlabel('Temperature [K]',color='r',fontsize=20) + ax2.axis([0,3000,100,1E-6]) + for tl in ax2.get_xticklabels(): tl.set_color('r') + ax1.set_ylabel('Pressure [bar]',fontsize=20) + ax1.legend(loc=4,frameon=False) + + + #plotting spectrum------------------------------------- + fig1, ax=subplots() + ax.plot(1E4/wnocrop, y_mod*100.) + xlabel('$\lambda$ ($\mu$m)',fontsize=18) + ylabel('(R$_{p}$/R$_{*}$)$^{2} \%$',fontsize=18) + minorticks_on() + ax.errorbar(1E4/wnocrop, y_meas*100, yerr=err*100,xerr=None,fmt='D',color='black',alpha=0.5) + ax.set_xscale('log') + ax.set_xticks([0.6,0.8,1,3,5,7,9,11]) + ax.axis([0.6,11,ymin*0.98,ymax*1.02]) + ax.get_xaxis().set_major_formatter(matplotlib.ticker.ScalarFormatter()) + ax.tick_params(length=10,width=1,labelsize='large',which='major') + show() + + + + + ##########################TESTING, PLAYING TO GAIN INTUITION###################################### + #------twiddeling--testing cloud param values---can test other params to by forcing x index or gas_scale values to new value + # 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 #19 gases (PH3, N2, and H are not opacity sources) + #H2O CH4 CO CO2 NH3 N2 HCN H2S PH3 C2H2 C2H6 Na K TiO VO FeH H H2 He + # gas_scale=np.array([1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.])*0. #all gases zeroed out--just want to test effect of cloud + gas_scale=np.array([1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.])*0. #all gases zeroed out--just want to test effect of cloud + + + ##plotting spectrum------------------------------------- + #fig1, ax=subplots() + #ax.plot(1E4/wnocrop, y_mod*100.) + # + ## pdb.set_trace() + # + ## Clean this plot up! Choose opacity sources or just plot top 5 + #for n,i in enumerate(gas_scale): #looping trhough 10 param values + # gas_scale[n]=1 + # y_mod2,wnocrop2,junk=fx(x,gas_scale) + # ax.plot(1E4/wnocrop2, y_mod2*100.,color='blue',alpha=0.2+0.2*i) + # gas_scale[n] = 0. + # # y_mod2,wnocrop2,junk=fx(x,gas_scale) + # # ax.plot(1E4/wnocrop2, y_mod2*100.,color='blue',alpha=0.2+0.2*i) + # + #xlabel('$\lambda$ ($\mu$m)',fontsize=18) + #ylabel('(R$_{p}$/R$_{*}$)$^{2} \%$',fontsize=18) + #minorticks_on() + #ax.errorbar(1E4/wnocrop, y_meas*100, yerr=err*100,xerr=None,fmt='D',color='black',alpha=0.5) + #ax.set_xscale('log') + #ax.set_xticks([0.6,0.8,1,3,5,7,9,11]) + #ax.axis([0.6,11,ymin*0.98,ymax*1.02]) + #ax.get_xaxis().set_major_formatter(matplotlib.ticker.ScalarFormatter()) + #ax.tick_params(length=10,width=1,labelsize='large',which='major') + #show() + #plt.savefig('./test.png') + # + #pdb.set_trace() + + #truth vector + #pickle.dump(x,open('truth.pic','wb')) + + + + + + + + diff --git a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/thermo.pyc b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/thermo.pyc old mode 100755 new mode 100644 index 0fc2060..5c32056 Binary files a/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/thermo.pyc and b/OPERATIONAL_HOT_JUPITER_SIMPLE_CLOUD_CLEAR/thermo.pyc differ