OV-GSO-DC / dustem-wrapper_idl

Commit 14ae4f23bea361531c406c0335f75774f4ba34f2

Authored by Jean-Philippe Bernard
1 parent daa89ea5
Exists in master

improved

Showing 8 changed files with 317 additions and 80 deletions   Show diff stats
Inline Side-by-side
Docs/developers/Code_Structure_Flowcharts/dustem-wrapper_flow_chart1.drawio
  Diff comments   View file @14ae4f2
1   -<mxfile host="Electron" modified="2024-02-29T16:40:49.500Z" agent="5.0 (Macintosh; Intel Mac OS X 10_15_7) AppleWebKit/537.36 (KHTML, like Gecko) draw.io/17.2.1 Chrome/96.0.4664.174 Electron/16.1.0 Safari/537.36" etag="b5mL6BM0y7VNfkiFUNQ-" version="17.2.1" type="device" pages="4"><diagram id="IKMg62xr7mvPSjzlxw8D" name="Dustem_wrap_flow">7V1bk6M2Fv4t+9CP02Xu8Jjpnp5kK8mmth+SeXLRINtUMBDAfdlfv+Jmgy4g04BkW1OVjsGyEPq+c3R0zpF0pz3s37+nbrL7LfZBeKeu/Pc77fFOVZWVo8L/FXc+qju25lQ3tmng14VON56D/4Hml/XdQ+CDrFMwj+MwD5LuTS+OIuDlnXtumsZv3WKbOOw+NXG3ALvx7LkhfvfPwM939Vuo1un+zyDY7ponK2b9fnu3KVy/SbZz/fitdUv7dqc9pHGcV5/27w8gLDqv6Zfqd0+Ub48NS0GUs/zg56/6Iflrtd8mf/7+YH57Be6/rS91La9ueKhf2D9kOdivN0G+zoC/ToHr78F9ksb1W+QfTdfAF0qKj/A7D2Swn7++7YIcPCeuV9x+g4SA99ww2Ebw0oPNBCm88QrSPIAd/FP9xT7w/aLKr1kF/up+pdm6qRinvzr8dhNH+XP9cKUonafx30dMtKJEEIYPcRinZfO0p6cH+A/e94MUMiOIi4cBN4Od9RXvu6YjYOPAe+tW3ZffQbwHefoBixypXeNaE1trrt9ONFHN+t6uRRGtvufWzNweqz6BBz/U+J2BpUrDMoiCXCKII6giCOq8EdT7EIQi+Q4FEr6yu88knINw2ipvOA0anBnI176buxLFQRRNmzeKJoZiMSau3z1XKlUG/DTe+Fk0KSxhhL0D+1/qUwKSdhdJg/vwaNON1QJEaedQgDTsLpDcLVWnVyS3qRtEEsdBgdRN3jg2k3QKkG4Yrt/Ai5bYEk0cTRNBk7u5qlCdAQ2aEkYcRh2ZRHK3VxXcD4BhBiL/p8I5VgASulkWeLA3dvk+rHs1jQ+RD4qnrOAV7J3046/2xY8SIKO5fHxvf/n4UV9Ruzd30y3Ie96hbjHwO845HIRWJxuEPm7upSB08+C169IjdXz9hD/iALaY7upBscviQ+qB+ldtFxxSkY4axxZSUdUxWEUQK/ejVSwpCmT0Blvk9tKahRY3DISE1fNPlDwi8AmW4r4OTiwF70H+V+tz9SvLqC9PvyouBqldUWGY2oMi0GhjQWQA4YiljhQBzEujGLOIACazptbfLrS8s4QQUD1E5YgbuvsXXzqJGKam/Add3Et0cYOuKZTCwQRytMYxUE0wj8ZBH4OG4/qLm/oS+gb3hd3GoKuwykAzoxVFCJDZhb4aKwTI+K1r6v1MA6+GNlnpbxn6itoCgtDoEv6CMF5dO2IzVddGMhX1cq1m4inqf9EHDETKCzKXt5bgNdWFU0aQIc2CaBNLi3IwbGVxD1s1GokcPd4cIg9i6ksoBz1yFvfJgapRoXRfQZlmtZGRKwapVFcKdyypOTrllL3+HL5KNDE0VUQwFYe/ku13wBQBrBTI5AAClohkKrbDHUvcBVNj+ZZCYNZZIqPKJCQRexvaxdyRpKbsFKaPTIBklUmNv0xSU3Y8qGChUBZ2bNEREk8CnogHVDkGYfjhSc3ckXgOW7OI/eNwn5g0FdOGTJnywZK543BPjNSobh85YJ4jkNwtH006faaB0uaeG6lRnT77pFhcJ9d/MLl8lBV3o0ejunxKJKFUShwHBVJRuLvuNKqzp8RRgjho6SjHXAp+IFK9PJVaTaU0MiRMKQp314BGdfLUU0kP9mEQHQ57iSeDa4C7K10bcvXAPtikrifRHEZT4z9W0pdoSQXLBKHucIawMcDkXPJcLFHTVefu4NGpDh7pqDsHSO4eV53q4DkCKZdMUvFEpyMG9zWTOtXLU+FZrmWWiFIRReclBne/nU719lSIert0K5FkSCPgP2hS/T11QgjsYJmpxQYm/4FzILvnH/clk2AygWnxHzWpvp9ax7qwcySabGjyHzGpnp/aBnJl8JkBSJt79FmnOn2ayQnscQklA5TcbZ8mTkNJbPaC1Fu/SjuWJUmde/zZoDp/kjDO1+A9D6qkOwnn8AoSlXvgy6C6gEo4o0MQHYAMerFgyT0abVDdPyWW9SbaEsvB9BBV5R7yMqiOHy/eJ4ccSCjZoFT4j5g9Oy/7pWhKHIfVq8I9p8CgOn3khugjBHPF3/Y5bwOaKI6KbvfdbFdutFH0KnUHDmqPibJDBpIiMH4bQawmE6mIskXGVLtMGLjD53ZQRLXkeBTRmkxrWRRxb8/toNjkN30eRc25b3ao4QRkI/03CaQxGZAGBFJ3Tv9QUO1lQcWdPrcDqjUZqNYAqM6yoOKun9sB1ZkMVEcsScV9QDcDKsRiIlBhTf2gLjymnrc99XWBOplxpGoDoJrLgoq7hW4HVHOqeQusqRdUY+Ex9bydl68LVMueClTL7gd14TG159CuLI9TuWneGX7cFfcwmXmeq4Fhu13RJRPbvVwZa+yiR4Y482yHix5jU7eXepJJt3iTiM9Y3FzicIUmCnXDpBu9Wzh6SIeOVDQR6TRye2nNQoobA8faIMXtLqVnIh3utin2eP2SgZw6THkfYQBZlmqUgarFx5eKj7++HG+43t/bkqX/OeSwmuMgVA1ZioEPP6b5VX16mmasMbGxhpD5Rhpq7LmGGkslSP1Tvk/g3/hwnSCgy6JI6YfLgkBLi7nfNCHIyzLXSqjmMdcw9ExClFYhjRrznZhLcn7UMlRsW3SnmmHx9i8p/LSt+qH+eudBsOhfxxDj5Bpl0EIGG0VvclC4yeB5zg6WIy7qM1uUu/aZLf1Htow+H6a4+gOkAeyMQtrLm5hgQgS1owi3vnl4cJwS2a5wq31YD54h0yz1HDyYoxEfQaxC3UGHBwfhHKtZaKC50ax24WTW1fSnbN02p5s9sy+N0wZyepZi2Mi0mJXTFsppg/G8mck4Pf2ZXM3ZWkce/2gxfIjTJxr/6LCYzOnz6Vs8FhWCKSndMHWY0pqk9EyUPi8/7OIpjdOXqOo/wWlm00MVitOWhTJRH83pbkBCsVDf1tycntyFTT4CkdH0UM4zPSbhNEFxf4LSKiOlmyNaRKE0RsSxITcL3TLMNpaldDP1vExrWkA1rTKrabHiBiaqpkfPELG5psN4vPdknD4vbXKsmqYy2lBPxL0OPa0wcroJ2QvCaQ1horqyRpoepobWhErH3JxWJ+f0yTS2zrKNxyj48dxr8tAubSpnolrweHDP2TaCgtY0U/jfRh+kmf3BVewHZjcaO0941cYDS8/fHu/fPfcqI3sGmvzBO7Bnn5cWK2dBAxrOZrUYTaEUnK508w5Hz4HMVbcic+kpEB4kgySByGZ4EBN+AbJr1DLoIGMQtmJYVsvMFua5US1jsNpRQmkZ1EFiGav7VeufOk7noDMDSzX6qp1bA+HhnyCChDjsIcxrH2ReGiTl5k7XauRYiJGjEfYaXVb94OGLTZzmqRvd0KBgIuEoTeGcz2JLB/y0g4LNOCgoYjkrUXVhmJOMCrqBVstzVHAYXD7ZW7AP3XI1Uwabkz/XaY842x4e6gUnE+gFxFYk5CrqKgF8pU8v1A/7L/ByN9rC9zuN1QjYX46VD5ANTUxyC6M+cnPwtVAB2SyY4eGUXx5/hVV9L/6unuNN/gZ7njpYCJz+yr5aqZ/NY1nGfcNqBw8rXAu27KnN82DLfftqRSFZFpWhlyWl0VeaetlHkcFelksD96W0+6pS8LHtghgLumbHgCFYgki38Tabjep5LdKEYJOTDAnffDENc0EgFcL+1cctyPp083RSinuHH8tlB1/+hN2cQNFSV08hNA/VlbeDIyYGFHzxvItGt1vrxcNtSOpb7CJMgr9LEFSMW4xQrYnmXIP6lTiKz6dfScnqbRFsoWT+c4grMWt9KgXvX83eqyl4DeJDtoaCQhPR00SOqepPtMaPxWhHBCA/hGkN2ckhQMM2QS5AK9rLbHg0ozxK19sF6vqtkn/+XRJEySGnj7nLNaTacJRnQ9JD1BynKkB/BH54OqxXgPYgERyOLcl28Vu53bAIjYHa/yXOgAAtiWLPzb0dtSWCm8+fDuOTtjrXCOaWiiYaTWdvXfRCKs5rpu7Y/KaOWInLaFKShk7DmD2lyLxPX3h1iUNaMIVP1xlH62vXNTpB2Sw8KyetBqLjNTR6Xz1ghNMTFgZstkjbFIODKvbowBpVa6RCkNEBTYrQzLELtdBVLQsvADh6IcVkr9jkZc23dsRKFNLVbhahoY4kr4Ec5KcvnOmvrEjLV+hjZeElus8gUbG8lE6JJA6HCxWu78FqiIVuZqg2EcvXUvGcGGXRpBilyUuQ2u58bceaFXmUSkHUnYEsSIIzu06qykjlN5QXufQ4TtpWnaAKUU81WR9SdBTVz90uRNWfrJXek2o9R5UyN1Pq5zaDbdLc115WP8+2YEbqZ0xVCKKfMUVq9SpS5rUy6lnVzq6fh0L2mBsuiDbxzeoih3t2jLKSTv/RuqgpyKCMxHL7L6SMbLPXBp1dGZFiAnLH0l4Fhew7SzrnfuG57Gzbhl2Kfpp4P5pmnsigtMRy6KFKy+6fio5VWo7GV2mx5R23EtPuiykg+6yx9cO0/5fU5x3SgnZVEHVdNpSwMpr+CwBllbSW+mbMQMNCtgtbkfLdSXbgfHq2sUxvV89+QqWy7i56lG5BVKqOOA0dZxqVimSSOzZXlaqcFz8pD5AqFdXmENFTKa9NKaFzUzgPJUxOrUUnp4oqldJYpaQz23mCreW9oj2PFeW8Y2AFI7CAe9PrzBt5C743vTk2HUEAUl90+EJEUjNvxSM6qUcniPEntTwZZ2JSX8spItbYaZAApL7o4I6IpLauhdTjz1zgTurZzsa5VVKzpqcLT+oLnihedJhLRFIzbxMuOKnty50oqiTHK0LqYO8W/d7mcX3rMdhv4XPD4AX+9cIgWZe7szx58T45VKvBnn4H+Usc/71WVPsd/nefRNth9+rT0zfDcabxmmLb3VSXLa5MlFoIL9O4WDh9Qih1k91vsQ+KEv8H</diagram><diagram id="-igyU8uB8o-wUyM342J-" name="plugin_managment">7ZtRj5s4EIB/yz1EujvpVhiwgccm23Slu+oq7Z1u25eIgpNYSzAFZ5Pcr68JJgEbSJaFhKrpwxYPY2N7xt/YAxkZk9X2Q+xGy4/Ux8FI1/ztyLgf6TrQHJ3/l0p2mcQ2nEywiIkvlI6CR/I/zmsK6Zr4OCkpMkoDRqKy0KNhiD1WkrlxTDdltTkNyk+N3AVWBI+eG6jS/4jPlmIUunWUP2CyWOZPBkiMb+XmymIkydL16aYgMt6PjElMKcuuVtsJDtLJy+clqzetuXvoWIxDdlaFP/9irvNsew+fyLfxQg9NNP5DtPLiBmsxYH+dMLyakZCwWRSsFyS8i2IqhsB2+bzw0UTpJb/n4YRP8nizJAw/Rq6XijfcG7jMDcgi5EWP9xHHXPCCY0b47L4TN1bE99Mmx0lmee1OM2wTAXj8a/K7cxqyR/FwkGqzmD4fDGKkGiQIJjSg8b57xnQ64f+43CcxdwtC04dhN+EzNVYnLp8F3jm8LYjERH7AdIVZvOMq4u7BqMKrgSnKm6OP6EjIlgX/MITMFW65ODR9tBy/EMZ7hSH1OkPerKdaD5atd7DK1axnKNZTTIZD/13Ks9QegZskxOOTsWSrQExqTNehj9OHaLyEt4Q9Fa4/760DRel+W7h1v8sLIR/KU7FQqJUWj9X2pbxerUUSuo493DDsHOhuvMCsQU/MBvZLrFbtW7AfrDBfLotx4DLyUiZ8lU3FEz5Rwkd2cB8Ayu5jAsktsnGLWkUiyw0Z5YYcKDWUTYzSEPcDd1dQi1KFpL7DkrsDCBu7Ja8OWIor/CJ7/tHbDxZovwDMOnzF61CEoeRGspNxCKFrkwzWGdLlo39xGb5Zs9aawClbE8JrWxNdJi5ZQwtM8MzAZA4rMNmS/6CWgUmXIhzqKTABKdRAW2vsl6yPwCv1rQuEMruOgNwRXd44jmf8ROvFJEoBorNdhG8sPLlHt/Rrs9D5SVlonslCe1As1DVp6TttWSht0m3NuoP90BApXb5r3qnLNSykl2v0Q7g8OlTu1n9Els3nGHleNcuyygVl33K+ao1Lqj3lbO3alANqRvCtmOuaWCdJlPvnUFAkH9LstiiSY6LZE4jM6g6fm8ewnAtss0BtxvOWMnhNysC+esoAdJ/9vDxyBpailDMJVke7H2T3hBxpu3bqJCiNz7YvQZzaJKVHV9Ga4VmC/Rtx1LSWFLbMq78sA2qW8q3E6fB9i97jC5dzD3NwWDiT6aC1xJmcIHP6oZlM3xNpKlkdXQJmamr3tn1qsX1ycpZcD2bWULZPbSDYnmXn5uhzTx8Ky6TVjnSJQW23ZpYpeVhHMLOqt5J13ZLUHeMSMKvPue9BNttutzeWnc6xX/8oeEuyn4DZwLLs8vHQ0eXk+Nk4s2Sc9XTSlN8JOs08k/Wtjo+a3nNg0Idv5ItvPY2THfGM9azis8xff/+lhLTfRnA8gvd3SYQ9ZYmUl0AV0grrZU+ljFsAVuXG5/o+N55TMMDz6rw4+oog6oZMSI5tUN1lAavCr+UY2IZMDpj/M5tGH588Ohn/G/6tuZMvFRbpCEyggxPjKS4pYQWhsVFlw8nEcaZTJTTpbzt0aircmrx+IGwzdbvsg203alCB5HkbtRYwqXRdvbeYenDXz0VvbXZdMFDXbfLI02F5WB8CQTnH39Z1kfJ+sp+gbMmf1Hb8zWnlulBfO5R++5B+rXPLgFSkc5FZ5uIFDw2VdlTT8l3xbTDJ3Cb/PQ2nQbFJPjIA1PbAIH38osN+2IRkBDYTUxqd+SaQ8eLxB2KZ+vFndsb77w==</diagram><diagram id="VhkHcouieJL9I0J5Rexb" name="Phangs_fiting_SEDs">7V1Zd6LKFv4t9yEvWSssqpgfO6aTHk6POUOnX1gopdIR8QJGc379KQQUqopBZNImZ63TWiIq3x6/vWtzJYzs7YNrrOafHBMtriBvbq+EuysIBUlS8D/Bymu4AnlNDVdmrmWGa+Cw8Gj9i6JFPlpdWybyUgf6jrPwrVV6ceIsl2jip9YM13U26cOmziL9qStjhqiFx4mxoFf/sUx/Hq6qUDmsv0PWbB5/MpC18BXbiA+Ofok3N0xnk1gS3l4JI9dx/PCRvR2hRXD14usSvu8+49X9F3PR0i/zhq/vn8Tbm5E4+bL+YOvffKSDLzfRWV6MxTr6weba85GtW0vLj763/xpfDPwTVsHDletMkIev7O1mbvnocWVMguUNlgG8Ziys2RI/neAvhly88IJc38KX9E30gm2ZZnDKWy+Em+d4QRVlIB3+L+JXp87Sf4w+HARH+67zvEdBCI6wFouRs3Dc3dcT7u9H+A+vm5aLZcFygg9Dhocvzy19teKfjr8c2iaWoqv3gBwb+e4rPiR69UYWIygjYQZC9HxzkAwoR2vzhFTExxmRMM725z7ghR9EkB0BH6Tgs41npHvI1FdzYznzdN8YL7D+DEimkVTkFJJK10AKFJAUZGhpvgkMWoDHwvA8a4Ivxty3F9FFdZ310kTBh/D4Gdpa/o/E46cdOlL07G6beOnuNX6yxD/lR/JJ+C4oxc8P79s9i9+YCYnnrN0JKhZg33BnyC+2U8hMWWsa4ASAEgO/eM1FC8O3XtI2ngVq9AlfHQv/soP8KAJhCURCMMIfHr0taZSpM0npM2GfmT5TeGmoM2FRMF4Th62CA7y8r8ynP0iTCr4YTxo7lZDy8CscZH4PQ3U1kLPc0d6s7ezZYM7yHRMUu7ZnSl5coeNfa9iDUyJQlLQ0ioLcNYpqFooTx16t/Z1GDijmo6jEBr0zFLUsFD3k66bhGwOEBYoIu1bEOBtOBvobd2r5gxEtUj+lRfV7uhvd3yuff326uxlrj97Pj776i5FiXwdUxzWXC9/kdWHhcN4VMgBMBP7jMPD/Y7xfMCbPs1068GXt49PswQmhBBINiyzfwvv7mjAACoGBSusPCwK1KQjoNJm67tWyK1BDdlWUXFEKhLES9qqWeGU00rQdhmklhPEZvyLXwtczUPcTs7bYHCXTtjzZ70vaJoK0YMoykWyVTtsATJ9JJc1HRtpWIVFiXth2+IJK8iwIQoFE1yyNZTkEuV/CSHEIgBCh0sIoEz4PAsg1xSIQdIXK5381it4AKdrhZBaBqRwipRwDqV06YIIMZ91YwHQz5f8C4pe3c/vjB/77+7/ttz+/5dck9Km1Ddjts2QQplMkTyZsLCl3biramM+1jMegLKdRljpHmY7JElnpecLbuKoSQYcMW8xtmCDSYUgE4sbFuOhGcIYBwySGokpg2GZ+ysQw11sGmmgtp84AY7493WcWncEo5cO4WM+s5dlZ1P44TLFNMp4JMF0bu/5fBHGWkqZTvAIaaQdONkM0nU7hDopYGBZoyoZBHsuSXBMMkKhCChJtLvd0UxIGsiBbGwx0ZYu67M1zSbVxQhX4qExYk4l6nqdJJup5kt6XRD2qlJ9e66dsw9FZ+Qk00j1yx2Dzh+Her0ZgbY/XgvqYnWqF/kLfbrccdg1XcORNnNVe3s/IgzQeCBAsjsS36CeYmGYmVhSmuyzLN2YDrMWwttlR8e3/1o+n24cpEL7/9fILWg/gu9tcDYPnVElK+B6Q6XluRE6QTvActXgt/GUqUcd59YmkR2JefKVXHonMIcE+SjrWJZHMMaDiprqIY4IHBnHCVBcRzISNpieirtixG3SvTB18kfQpTo82lj/XLc+d6ru29MEYJo0hAKSQtNnSwgSW5iyu9V8bzw8akjzdXntIX1lbtPA4z3jJhPOMa+wSqU+gxRo7ExKaf6jLPx2RF+09Ezgjz5Qn4YWeCfSrxC6Tckk5lLKuSST7jyEsV9Osy3lk9iLjD9U3xgtaoOXMn58dp9a0t5BBGjdB6jp0zmxGJsKAAcl8JMU2c1smknRDMvbvw3aAAtxa5SSYuNEtyPrd7VeN1z/88/in/sDrT1/fvNOf/n580J3lzYvjOkvH2kXlN+8fv9/vogXkXV9qNCeTrUCs8nyr0Ryj4bimcK5ayyQHxH0M95R8lRHQRSxFt32TeearMKjTehXTUY1qAmknSsd0AulRWo7pAE14u8jzHXfwIUVRXJsFbjZ2NAN6Hab6v0vST7kJyAit23UTLCLGCITXM34fMobeKssYDNAuLu2wMZnOGypcqlBw9nRMbHzOjY8h21gAZZ8r+26h5E712nw3TcgMMn2STJ8rx0iOYYGqQu53OB+pztxjbeMPtpY4wPwXDVUpRoRKli7FznlGSHvdqN6Itr5rTMKIaMAxHTyROAJB6xrHdjxNNFfn2GJWjp9hBbp94D1ixSj0NGKvHM0FBU+wsbbVpEhntw4RbJ4sK2Vluvst0KeJdJw59ESmqY4fXqscPMmQ1A++Zammg6e6pJoDODJLSDbHK6W56qifrg9cdb0ZQ3mh75chp4R+v33/aJFXcVgFCamXiZM1LfV0NS0KM3dT0nZ9UJZtzIaxj0SoqalpLyx3XheNT5wsjCbwu1TCVINpIBTGRtZW+dI4POtJ91qCTTrX8Egr6SliDeiJpxBlbMsSf2lTr1Z2G+S4GKldlyHAXJexq9FgKcLIe4PrYPdAE6Sj0r3voLvb00AGo/1CDzK6SCdCTBlTGVt02nUidDE0BOSynTngibK01nX1U6B52BCIz+8uFgTSPmkM+9QuCB2X6y4soIo9eHFAJfQqoNI0IuciO77KhlAa0ZGktMygCrUzqHVVhOk0o4CqakMK+0UAAUD4KEWpGMkDngw7yLF9Tcth7ZxnvcMkjrO0maTlcaWFNiRa6pdEk5NRTpBoFXIaPPynpgOJkjuCa5PvfHZz31w4pKmsMBCS1qnNOX/ski6D4tyDeEhRLzEkp0ptXbfrirXTnGfhO4p41xO8R1maM9aCvngPQOzTVapu0z1sPI/9RcvxkEhTmyfKdLIGLKXEiA9HWzc7lotRsc3eBwVydCQ1rqWDInDPqP3aMlGFkHiVnL7StMQ3N/ydA6qQlnhRHCS+vMT3K+tVocZpiWIW0YMGRPyyVk0HNEHjZO3wl84ZcE7MyRL1wW3pB83D16YfOB1K6MdN4BI61o9UxFS4C6Gzxs5++QJVyNcM6QTNEHM1A7akGRP+M699+JNf2J/R88sIOj+/fmSMe9THzzPdW4937enI1O2NPnHcIDVFZq9uUGRKSDVFluaocIzztr0YUwkaQ7KzczZyxLrSYsrGhIyObsM5Zba11FcuMq0didDjHPokPCRyWnqLcxGYcNCh146Vihgqz1t5ITwHaDJhOVuG6iREIaFgrc4mYUKaOf8+mpcaY+vjcxkuNo9L31qu1/YAbBpYglGQGJO62wW29j3QDU7dzgnbkvtOs9iuY4aqMAOxwvCPHU/wdPjHPFDuMtgTyf2iasXEXyDvC9pcaST7KjJigWg8ku4ske7iK+VcbDhAALBvAesqHKB7AiZzNHke4oHybgNkkGmduQ26vp4M8ZIKN0CZGwGobW7IZkJJl5IHyHIhY3UNtgtZmcFzA2YpzDpXM8ZUtSsoL4LfO3bxo1nw6PeKVMhRXfudI11FKozpadR1v8xsqGDKwlV96RBjvGOOvnSVDpEUp6ZU3f5Njj3jeWoredOEdu2l0FqlmhgRpRRWZw51WPHIW7Jk7vOuJ9EHJUVb65Nkq3L10QYiGcu2Lto0L3mZs65rDYbk7oOhC6EdQW9JR8bYuuzS2WCL6rBFte/lGtxsTrG3ULKVTiWbSG3kqnQ6yeYqzU0PYks1zd4eOtPjexAciFt9il+81HyVmIyyh6KzdLXnu5zSFqeswUl5zExjI0GujZBeOQtjo6mckmjUEtJiKkucLB56reRqhoi8n6IqwLYdbAmmuq27e7bR0ot/b0qBeACLVCi7cfFYdcht3S1UB6FLdVAVmYvbWw9tvCpHnKj05g1yEyEAJXvZj73/p0p8UDwjIrPHngifRZh/PJRA3vFXp95elB2u0bWK3NahgALPDB6G/D0MMrrO3yGrmNGjyCPis0tvkqs0477h0IM2tTm767qytQCSO9vI/KS8mSVv3MI3NzKSLdSQEurMG7DZxso7q1uvnWSAqHtpyy1OL2Jj1dj2rjRhV36ebZ9m59QS7JUt1MFOKURqaLgERS6xs0SrzCeqOTtUgKQFicAh5RLbtVQNbt46mu07egB/k1u+anTAckkH3HGh+izH77Olurn7cXHlZLqlzbqqIiR15QYrnXDCnsSmWPUeSrvM51plXpI4ufpNjlJDbUWJ0AZyrG1NSb9CzkzVCpJ+hQyTC1gCkbQQ8RTEimk/fuo6jp88HAew80+OiYIj/gM=</diagram><diagram id="ls9p14zF9e6sMFj8i8uX" name="Phangs_pluggins">7Vpbc5s6EP41fjkz9QDiYh4TJ2460/bkxKfT5MkjgwxquB0hfOmvrzDiJiC+HJu4HefBQctKSPt9u9oVDMDYX38kMHK/hDbyBopkrwfgbqAosmQq7F8q2WSSETAzgUOwzZVKwRT/RHlPLk2wjeKaIg1Dj+KoLrTCIEAWrckgIeGqrrYIvfpTI+ighmBqQa8p/Y5t6vJVKEYpf0DYcfMnyzpfnw9zZb6S2IV2uKqIwP0AjEkY0uzKX4+Rlxovt0vWb9Jxt5gYQQHdp8Noktj/Pgdfp98/0Zunb/Ldw/zbBz7KEnoJX7CdxBT5s8hLHBzMIhcGTjxjIs+DZIZjsuCroZvcRGxhUXoZkdBCMbP37crFFE0jaKXiFSMGk0EPOwFrWmy6iDDBEhGKmaFv+A0f23Y65G2ckUAaSmCk6rJW/qrs7iIM6JQ/XE61KQlfC2xAqoE9bxx6IdlOD0wmY/bH5DYmjCE4TB+GYMyMdstXz2aC1p1mlQuwGMtR6CNKNkwl76ByfDnBNd5clWxRdC5zK0wBXAY5QZ1i5BJDdsFhPABSpQFpAzEU2Depb6RweDCOscVs4VLf4zYlYRLYKH2IxFpojelz5fplC47GW3fryq27Td4I2FKeq42sl6Ll7bLftpV3zOaK7IZXCniw9YQJsdBublNIHETf0APt+FYA1Frwy2UEeZDiZX26baDyJzyGmC2kpI9Up49qCLzIlsl7Vd1bHMisD6RL5tA062NlpmiMxbgANxW1KFWI95/zSHp7ZoK6LNf02UU2gZLyBQrHewE4LLCx7YPiIEn8a3SrRzdDwE5qCW+jPsOb2gUsQdCeIR97qISVySPmnii+wlqDVRFdsm3X6hVW7U1YrblsXkE9ENTWVORcoD7h108v6Mf4o//of/n633rx88d9S3ZZYDfDwSIc0jXtBNHaeJhlIQR0wFjJV+ZZvvJ5Xgig9epss5i/E8qGKSDKAJW1Jji6fqtMJqdBQlfqSIAmEG0p4ehcODRTwnvLlYdA+mu4wLTbj/4gCOT3xqAzIWEZ2WwFl8hDgUPda1ATghoQglqfCcg/yycmWT8A8uHGsRYvxubzU0tQa0B2WH3VMCjzA1CYvnJnPDbNrX/UQVGOrdFY4xERzMySMqZet8n1uu08Zdtbm8bOsk1t500/ZZsmlm0i3fYt23RhIE0VBuqo2Y6ok1rZfPLTgiub27bey2azLgRZcRPcl8yqcAYhG1qvZG7usdZckWTzD05zxO1R7S/NaYWgWZ5f48mR8eQtiu+MJ9olxRP5dw0onacS16T9kKQdvHvSrp86LLWFFPniX4nsGz3eNRspDiM5e8wjg0cRdYr3DsZwv/Bx8AsR4R1gftjQOTNRXzvtG5FWFzD6cAFpaFz+i8HfYg9tsFcsDPd1g+JTinwg/WxuIG774hcJO/TzTPWsbjC6aDfo0QvUphd0H1xevWB/LxCnbOzwAiHHZbvB+b1APvkp5im/EjF6dAOt6QbdpcB7uUGDvSNFZO/ejiDkHmyOvbiBsssNBH2g/K/NgDXLT+sy9fIDRXD/Cw==</diagram></mxfile>
2 1 \ No newline at end of file
  2 +<mxfile host="Electron" modified="2024-03-02T08:15:08.468Z" agent="5.0 (Macintosh; Intel Mac OS X 10_15_7) AppleWebKit/537.36 (KHTML, like Gecko) draw.io/17.2.1 Chrome/96.0.4664.174 Electron/16.1.0 Safari/537.36" etag="psekgSISGgzYhJFoVu_r" version="17.2.1" type="device" pages="4"><diagram id="IKMg62xr7mvPSjzlxw8D" name="Dustem_wrap_flow">7V1bk6M2Fv4t+9CP02Xu8Jjpnp5kK8mmth+SeXLRINtUMBDAfdlfv+Jmgy4g04BkW1OVjsGyEPq+c3R0zpF0pz3s37+nbrL7LfZBeKeu/Pc77fFOVZWVo8L/FXc+qju25lQ3tmng14VON56D/4Hml/XdQ+CDrFMwj+MwD5LuTS+OIuDlnXtumsZv3WKbOOw+NXG3ALvx7LkhfvfPwM939Vuo1un+zyDY7ponK2b9fnu3KVy/SbZz/fitdUv7dqc9pHGcV5/27w8gLDqv6Zfqd0+Ub48NS0GUs/zg56/6Iflrtd8mf/7+YH57Be6/rS91La9ueKhf2D9kOdivN0G+zoC/ToHr78F9ksb1W+QfTdfAF0qKj/A7D2Swn7++7YIcPCeuV9x+g4SA99ww2Ebw0oPNBCm88QrSPIAd/FP9xT7w/aLKr1kF/up+pdm6qRinvzr8dhNH+XP9cKUonafx30dMtKJEEIYPcRinZfO0p6cH+A/e94MUMiOIi4cBN4Od9RXvu6YjYOPAe+tW3ZffQbwHefoBixypXeNaE1trrt9ONFHN+t6uRRGtvufWzNweqz6BBz/U+J2BpUrDMoiCXCKII6giCOq8EdT7EIQi+Q4FEr6yu88knINw2ipvOA0anBnI176buxLFQRRNmzeKJoZiMSau3z1XKlUG/DTe+Fk0KSxhhL0D+1/qUwKSdhdJg/vwaNON1QJEaedQgDTsLpDcLVWnVyS3qRtEEsdBgdRN3jg2k3QKkG4Yrt/Ai5bYEk0cTRNBk7u5qlCdAQ2aEkYcRh2ZRHK3VxXcD4BhBiL/p8I5VgASulkWeLA3dvk+rHs1jQ+RD4qnrOAV7J3046/2xY8SIKO5fHxvf/n4UV9Ruzd30y3Ie96hbjHwO845HIRWJxuEPm7upSB08+C169IjdXz9hD/iALaY7upBscviQ+qB+ldtFxxSkY4axxZSUdUxWEUQK/ejVSwpCmT0Blvk9tKahRY3DISE1fNPlDwi8AmW4r4OTiwF70H+V+tz9SvLqC9PvyouBqldUWGY2oMi0GhjQWQA4YiljhQBzEujGLOIACazptbfLrS8s4QQUD1E5YgbuvsXXzqJGKam/Add3Et0cYOuKZTCwQRytMYxUE0wj8ZBH4OG4/qLm/oS+gb3hd3GoKuwykAzoxVFCJDZhb4aKwTI+K1r6v1MA6+GNlnpbxn6itoCgtDoEv6CMF5dO2IzVddGMhX1cq1m4inqf9EHDETKCzKXt5bgNdWFU0aQIc2CaBNLi3IwbGVxD1s1GokcPd4cIg9i6ksoBz1yFvfJgapRoXRfQZlmtZGRKwapVFcKdyypOTrllL3+HL5KNDE0VUQwFYe/ku13wBQBrBTI5AAClohkKrbDHUvcBVNj+ZZCYNZZIqPKJCQRexvaxdyRpKbsFKaPTIBklUmNv0xSU3Y8qGChUBZ2bNEREk8CnogHVDkGYfjhSc3ckXgOW7OI/eNwn5g0FdOGTJnywZK543BPjNSobh85YJ4jkNwtH006faaB0uaeG6lRnT77pFhcJ9d/MLl8lBV3o0ejunxKJKFUShwHBVJRuLvuNKqzp8RRgjho6SjHXAp+IFK9PJVaTaU0MiRMKQp314BGdfLUU0kP9mEQHQ57iSeDa4C7K10bcvXAPtikrifRHEZT4z9W0pdoSQXLBKHucIawMcDkXPJcLFHTVefu4NGpDh7pqDsHSO4eV53q4DkCKZdMUvFEpyMG9zWTOtXLU+FZrmWWiFIRReclBne/nU719lSIert0K5FkSCPgP2hS/T11QgjsYJmpxQYm/4FzILvnH/clk2AygWnxHzWpvp9ax7qwcySabGjyHzGpnp/aBnJl8JkBSJt79FmnOn2ayQnscQklA5TcbZ8mTkNJbPaC1Fu/SjuWJUmde/zZoDp/kjDO1+A9D6qkOwnn8AoSlXvgy6C6gEo4o0MQHYAMerFgyT0abVDdPyWW9SbaEsvB9BBV5R7yMqiOHy/eJ4ccSCjZoFT4j5g9Oy/7pWhKHIfVq8I9p8CgOn3khugjBHPF3/Y5bwOaKI6KbvfdbFdutFH0KnUHDmqPibJDBpIiMH4bQawmE6mIskXGVLtMGLjD53ZQRLXkeBTRmkxrWRRxb8/toNjkN30eRc25b3ao4QRkI/03CaQxGZAGBFJ3Tv9QUO1lQcWdPrcDqjUZqNYAqM6yoOKun9sB1ZkMVEcsScV9QDcDKsRiIlBhTf2gLjymnrc99XWBOplxpGoDoJrLgoq7hW4HVHOqeQusqRdUY+Ex9bydl68LVMueClTL7gd14TG159CuLI9TuWneGX7cFfcwmXmeq4Fhu13RJRPbvVwZa+yiR4Y482yHix5jU7eXepJJt3iTiM9Y3FzicIUmCnXDpBu9Wzh6SIeOVDQR6TRye2nNQoobA8faIMXtLqVnIh3utin2eP2SgZw6THkfYQBZlmqUgarFx5eKj7++HG+43t/bkqX/OeSwmuMgVA1ZioEPP6b5VX16mmasMbGxhpD5Rhpq7LmGGkslSP1Tvk/g3/hwnSCgy6JI6YfLgkBLi7nfNCHIyzLXSqjmMdcw9ExClFYhjRrznZhLcn7UMlRsW3SnmmHx9i8p/LSt+qH+eudBsOhfxxDj5Bpl0EIGG0VvclC4yeB5zg6WIy7qM1uUu/aZLf1Htow+H6a4+gOkAeyMQtrLm5hgQgS1owi3vnl4cJwS2a5wq31YD54h0yz1HDyYoxEfQaxC3UGHBwfhHKtZaKC50ax24WTW1fSnbN02p5s9sy+N0wZyepZi2Mi0mJXTFsppg/G8mck4Pf2ZXM3ZWkce/2gxfIjTJxr/6LCYzOnz6Vs8FhWCKSndMHWY0pqk9EyUPi8/7OIpjdOXqOo/wWlm00MVitOWhTJRH83pbkBCsVDf1tycntyFTT4CkdH0UM4zPSbhNEFxf4LSKiOlmyNaRKE0RsSxITcL3TLMNpaldDP1vExrWkA1rTKrabHiBiaqpkfPELG5psN4vPdknD4vbXKsmqYy2lBPxL0OPa0wcroJ2QvCaQ1horqyRpoepobWhErH3JxWJ+f0yTS2zrKNxyj48dxr8tAubSpnolrweHDP2TaCgtY0U/jfRh+kmf3BVewHZjcaO0941cYDS8/fHu/fPfcqI3sGmvzBO7Bnn5cWK2dBAxrOZrUYTaEUnK508w5Hz4HMVbcic+kpEB4kgySByGZ4EBN+AbJr1DLoIGMQtmJYVsvMFua5US1jsNpRQmkZ1EFiGav7VeufOk7noDMDSzX6qp1bA+HhnyCChDjsIcxrH2ReGiTl5k7XauRYiJGjEfYaXVb94OGLTZzmqRvd0KBgIuEoTeGcz2JLB/y0g4LNOCgoYjkrUXVhmJOMCrqBVstzVHAYXD7ZW7AP3XI1Uwabkz/XaY842x4e6gUnE+gFxFYk5CrqKgF8pU8v1A/7L/ByN9rC9zuN1QjYX46VD5ANTUxyC6M+cnPwtVAB2SyY4eGUXx5/hVV9L/6unuNN/gZ7njpYCJz+yr5aqZ/NY1nGfcNqBw8rXAu27KnN82DLfftqRSFZFpWhlyWl0VeaetlHkcFelksD96W0+6pS8LHtghgLumbHgCFYgki38Tabjep5LdKEYJOTDAnffDENc0EgFcL+1cctyPp083RSinuHH8tlB1/+hN2cQNFSV08hNA/VlbeDIyYGFHzxvItGt1vrxcNtSOpb7CJMgr9LEFSMW4xQrYnmXIP6lTiKz6dfScnqbRFsoWT+c4grMWt9KgXvX83eqyl4DeJDtoaCQhPR00SOqepPtMaPxWhHBCA/hGkN2ckhQMM2QS5AK9rLbHg0ozxK19sF6vqtkn/+XRJEySGnj7nLNaTacJRnQ9JD1BynKkB/BH54OqxXgPYgERyOLcl28Vu53bAIjYHa/yXOgAAtiWLPzb0dtSWCm8+fDuOTtjrXCOaWiiYaTWdvXfRCKs5rpu7Y/KaOWInLaFKShk7DmD2lyLxPX3h1iUNaMIVP1xlH62vXNTpB2Sw8KyetBqLjNTR6Xz1ghNMTFgZstkjbFIODKvbowBpVa6RCkNEBTYrQzLELtdBVLQsvADh6IcVkr9jkZc23dsRKFNLVbhahoY4kr4Ec5KcvnOmvrEjLV+hjZeElus8gUbG8lE6JJA6HCxWu78FqiIVuZqg2EcvXUvGcGGXRpBilyUuQ2u58bceaFXmUSkHUnYEsSIIzu06qykjlN5QXufQ4TtpWnaAKUU81WR9SdBTVz90uRNWfrJXek2o9R5UyN1Pq5zaDbdLc115WP8+2YEbqZ0xVCKKfMUVq9SpS5rUy6lnVzq6fh0L2mBsuiDbxzeoih3t2jLKSTv/RuqgpyKCMxHL7L6SMbLPXBp1dGZFiAnLH0l4Fhew7SzrnfuG57Gzbhl2Kfpp4P5pmnsigtMRy6KFKy+6fio5VWo7GV2mx5R23EtPuiykg+6yx9cO0/5fU5x3SgnZVEHVdNpSwMpr+CwBllbSW+mbMQMNCtgtbkfLdSXbgfHq2sUxvV89+QqWy7i56lG5BVKqOOA0dZxqVimSSOzZXlaqcFz8pD5AqFdXmENFTKa9NKaFzUzgPJUxOrUUnp4oqldJYpaQz23mCreW9oj2PFeW8Y2AFI7CAe9PrzBt5C743vTk2HUEAUl90+EJEUjNvxSM6qUcniPEntTwZZ2JSX8spItbYaZAApL7o4I6IpLauhdTjz1zgTurZzsa5VVKzpqcLT+oLnihedJhLRFIzbxMuOKnty50oqiTHK0LqYO8W/d7mcX3rMdhv4XPD4AX+9cIgWZe7szx58T45VKvBnn4H+Usc/71WVPsd/nefRNth9+rT0zfDcabxmmLb3VSXLa5MlFoIL9O4WDh9Qih1k91vsQ+KEv8H</diagram><diagram id="-igyU8uB8o-wUyM342J-" name="plugin_managment">7ZtRj5s4EIB/yz1EujvpVhiwgccm23Slu+oq7Z1u25eIgpNYSzAFZ5Pcr68JJgEbSJaFhKrpwxYPY2N7xt/YAxkZk9X2Q+xGy4/Ux8FI1/ztyLgf6TrQHJ3/l0p2mcQ2nEywiIkvlI6CR/I/zmsK6Zr4OCkpMkoDRqKy0KNhiD1WkrlxTDdltTkNyk+N3AVWBI+eG6jS/4jPlmIUunWUP2CyWOZPBkiMb+XmymIkydL16aYgMt6PjElMKcuuVtsJDtLJy+clqzetuXvoWIxDdlaFP/9irvNsew+fyLfxQg9NNP5DtPLiBmsxYH+dMLyakZCwWRSsFyS8i2IqhsB2+bzw0UTpJb/n4YRP8nizJAw/Rq6XijfcG7jMDcgi5EWP9xHHXPCCY0b47L4TN1bE99Mmx0lmee1OM2wTAXj8a/K7cxqyR/FwkGqzmD4fDGKkGiQIJjSg8b57xnQ64f+43CcxdwtC04dhN+EzNVYnLp8F3jm8LYjERH7AdIVZvOMq4u7BqMKrgSnKm6OP6EjIlgX/MITMFW65ODR9tBy/EMZ7hSH1OkPerKdaD5atd7DK1axnKNZTTIZD/13Ks9QegZskxOOTsWSrQExqTNehj9OHaLyEt4Q9Fa4/760DRel+W7h1v8sLIR/KU7FQqJUWj9X2pbxerUUSuo493DDsHOhuvMCsQU/MBvZLrFbtW7AfrDBfLotx4DLyUiZ8lU3FEz5Rwkd2cB8Ayu5jAsktsnGLWkUiyw0Z5YYcKDWUTYzSEPcDd1dQi1KFpL7DkrsDCBu7Ja8OWIor/CJ7/tHbDxZovwDMOnzF61CEoeRGspNxCKFrkwzWGdLlo39xGb5Zs9aawClbE8JrWxNdJi5ZQwtM8MzAZA4rMNmS/6CWgUmXIhzqKTABKdRAW2vsl6yPwCv1rQuEMruOgNwRXd44jmf8ROvFJEoBorNdhG8sPLlHt/Rrs9D5SVlonslCe1As1DVp6TttWSht0m3NuoP90BApXb5r3qnLNSykl2v0Q7g8OlTu1n9Els3nGHleNcuyygVl33K+ao1Lqj3lbO3alANqRvCtmOuaWCdJlPvnUFAkH9LstiiSY6LZE4jM6g6fm8ewnAtss0BtxvOWMnhNysC+esoAdJ/9vDxyBpailDMJVke7H2T3hBxpu3bqJCiNz7YvQZzaJKVHV9Ga4VmC/Rtx1LSWFLbMq78sA2qW8q3E6fB9i97jC5dzD3NwWDiT6aC1xJmcIHP6oZlM3xNpKlkdXQJmamr3tn1qsX1ycpZcD2bWULZPbSDYnmXn5uhzTx8Ky6TVjnSJQW23ZpYpeVhHMLOqt5J13ZLUHeMSMKvPue9BNttutzeWnc6xX/8oeEuyn4DZwLLs8vHQ0eXk+Nk4s2Sc9XTSlN8JOs08k/Wtjo+a3nNg0Idv5ItvPY2THfGM9azis8xff/+lhLTfRnA8gvd3SYQ9ZYmUl0AV0grrZU+ljFsAVuXG5/o+N55TMMDz6rw4+oog6oZMSI5tUN1lAavCr+UY2IZMDpj/M5tGH588Ohn/G/6tuZMvFRbpCEyggxPjKS4pYQWhsVFlw8nEcaZTJTTpbzt0aircmrx+IGwzdbvsg203alCB5HkbtRYwqXRdvbeYenDXz0VvbXZdMFDXbfLI02F5WB8CQTnH39Z1kfJ+sp+gbMmf1Hb8zWnlulBfO5R++5B+rXPLgFSkc5FZ5uIFDw2VdlTT8l3xbTDJ3Cb/PQ2nQbFJPjIA1PbAIH38osN+2IRkBDYTUxqd+SaQ8eLxB2KZ+vFndsb77w==</diagram><diagram id="VhkHcouieJL9I0J5Rexb" name="Phangs_fiting_SEDs">7V1Zd6O4Ev4t9yEvOSccJLE+TpxOLzPTW2ZLv3AwxjbdtvEFnKTn14+wwQZJLMYgYYe+58xNsOOFr1T11Vel0hUaLV/eBvZ6/rs/cRdXUJ68XKG7KwihASD+v/jKz90VoEJ1d2UWeJPk2uHCg/evm1yUk6sbb+KGuSdGvr+IvHX+ouOvVq4T5a7ZQeA/55829Rf5d13bM5e68ODYC/rq394kmu+uGlA/XH/nerN5+s5AM3ePLO30yck3Cef2xH/OXEJvrtAo8P1o99PyZeQu4ruX3pfd390XPLr/YIG7iur8wef3j8rtzUhxPm0+LK0vkWuBTzfJqzzZi03yhSebMHKXlrfyouRzRz/Tm4G/wjr+cR34jhviO3v7PPci92FtO/HlZ2wD+Jq98GYr/KuDP5gb4AtPbhB5+Jb+kjyw9CaT+CVvwx3csiQjQ9GAevivgh+d+qvoIXlzED87CvwfexRQ/AxvsRj5Cz/Yfjx0fz/C//D1iRdgW/D8+M1cO8S355a+W+lXxx/OfclcSu7eW9dfulHwEz8lefRGUxIoU2NGye/PB8uAWnJtnrGK9Hl2Yoyz/Wsf8MI/JJAdAR+k4FvaP1wrdCfWem6vZqEV2eMFXj8DknkkdS2HpC4aSEQBSUHmria/xA4txmNhh6Hn4Jsxj5aL5KYG/mY1ceM3kfFv7osX/ZP5+XGLjpr8dveSeejuZ/rLCn+Vf7K/7P4Kqunvh7/b/pb+YSEkob8JHLfagCM7mLlRtZ9yJzlvTQOcAVBl4JdeC9yFHXlPeR/PAjV5h8++h7/ZwX50RHgChTCM3RdP/izrlKlXUvOvBJGaf6XdraFeCZuC/TPztHX8hLDsI8v5NzLVig8mk87OIKx89xEONr+Hofky0IrC0d6tbf3Z4M7KAxNURPszvYxXWPjb2sshKBEoqmYeRaSJRtEoQtHxl+tNtF2RA4rlKOqpQxeGolmEYuhG1sSO7AHCioUIRS/ENBvOEv3nYOpFgxOtWn46x+X3eDe6v9c/fv/97mZsPoTffo2M74wU+zqWOq6lUvicnwsP0/kAFQCYIf7jHfH/bby/YDs/Ztt04NMmwi+zB2cHJVBpWDTtFt7ft4QB0AkMDHr9sCAwuoKATpOp+94suwItZFdVyRW1gDBWaL/UMo+MRqa5xTC/CGH6ip/dwMP3M17uJ2ZtqTvKpm1ltt+XtE0BecPUNCLZqp22AZh/JYN0HwVpW4NEiXlj+egFjewZIVRh0S1bY10NQeuXMVIaAiBMqLYxakTMgwBKXakIhFxhyOUfjZI3QE52OFlFYC4OhVocg6hdmzBBRrDujDDdTOU/gfLpzXz56wf56/u/lm++fSmvSVhT7yVWt89SQZhOXc1x2FhS4Xyim2O51DMeg7KWR1kVjjLNyTJZ6XnC2/lSJUiHBjnmNkwQaRqSgPgcYFwsO36FAcMshopBYMgzP2ViWBot45Xorab+AGO5P91nFsJgVMthXGxm3ursPGp/AqbCU4xnAkzXxq7/l0BctEjzKV6FjLQFp1ghmk6ncAtFagwLd8qGQRtrqtYSDJCoQiKVdpd7uSkLA1mQbQ0GurJF3fbutaTWNKEGelQhrNlEvSzSZBP1MkvvS6KeVMpPr/VTvuHorPwEGeneDcbg+Tc7uF+PwGY53iDjoTjV2sUL6+XlRcKh4QqOQsdf7+39jCJI50SAUHFUmWOcYGJamFhRmG6zrMieDbBWw8qzo+LL/71/Hm/fTgH6+ufTd+i9BV+D7moYsmSoaib2gMLIc6NISD0hcrQStfCHaSQdl9UnshGJefP1XkUkMocEe5Z0bEgilWNA8aa2hGNCBwZpwtSWEMyEjZYnkq7YcRB3r0x9fJOsKU6Pnr1obnlhMLW2bemDM8w6QwBII+HZ0sIEltYsrq3vz2EUNySF1nITutbae3EXoRTaT4VwnnGNXYcEKWXJ9l3V2JmQ0PpDW/HpiLxoH5nAGUWmMguvjEygXyV2jfTzVECpG5oUsv8Ywno1zbaCR2EvMn5T69l+chfuahbNz05T6zpaaCCPG1JFU+fCZmSCBgxIliOp8MxtmUjSDck4vg/bASpw46pJMHGjW5Ctu9vPpmx9+PvhD+utbD1+/uWd9fjXw1vLX908+YG/8r0tK795//D1fssW3PD6UtmcRrYCscrzXNkco+G4JTrXrGVSAsqewz1mH2UQukSlENs3Wea+Kkmd2StORzWqIdJP1OZ0iIwonDkdoAXvwA0jPxhiSBWL41ngZmNHK6DXu1T/tST9hkwk/SxqzTdMsIQYOzbe0H49YoyhkxVCxmAAvrjwUWMKgzfUpVyh4OzlmNT5nJseQ7axAMo/N47dqOZO9dZiNy3IDDZ9kk2fq8ZIjmHBDpjc73A+Vl24x3qJ39hbYYL5rztUpRgMlSxdKsJ1RkhH3aTe6L5Ege3sGNGAYw5HUydXs3gc+USaZK7OscWskjjDIrp90D3ShVEZaZReBZoLIk+ws7bVrEkXtw4Rap6m6XVtWvwW6NNMOs0cemLTVMePbDYmTxok14fM2app8tSWVUsAM7OMZUuyXlurTvrp+qBVt5sx1Df6fjlyyuj32/ePNnkD0ypIWL1GvFjXVk9X0xKauZ2Stu2D8pb2bBj7SDWyEZUMTbionb5wtjKaAfBSFVMgEz2FOmPEElfFNCVoPelfy+hJ50qQzJqxIl0CPYkVioa9WeZf3tkbjQMHOTBG5Rs0ECwNGtsqDbYijHw4BA+my0Kky9JEd0Ejur89D2Q83G8XQkYXGUXIznTWsDe+UYSuh+4QufRwThArU3QBFNFS7A6Ij+8uFgRSETcZ7JYvCIIrdhfGqNIQXs2oUK8YlWnmxU+NbPqqy6FMoilJ5yyiotZF1LaKwnSeUaFW8bDCfmlAVPKv6w2pPJCJkpVBTu7r2g5blz3bnSdxnKct1C2Pqy7wsGi1XxZNDkc5waINKJnw8D8jTyRqbgpuzb7LBc59f+GQp7JoIDm6gEUDOe/WZYicexAPOeoroORAFt2xq7Suc55F7KgSXk+IHnV1znQV9CV6AGKrrt50py4ARoHH4RQvFFrbPNGmD2VglDMis4JWd9XkmVkfsqllV8iNLIEmi+QEg5drGjzsF11qLQ01yPErMmd+pHQ3/Z20d8MQbe9AVXIRAbMwxNng6/ZJ98zgDWhKZqaSRbSgAQU/bDZbAiY5nEHm7fJpob2tJVBLWmyv07Lc+hWF9PbczV85U39voHLzV5ubP5DJzlC5ux4gR/4omx/+kBfLj+6PpxH0v33+lTG40Rr/mFnhZrxtNHcn1vLZcvwgzjDdSa+OGpqorjFRWMvEgGOcfu2NlcqzGPZbnHqRUVrnWA1hQkaT1N3EsaW3staBO/G2WkCPU+GT8FDJueccJxww4aA51FZcSoSmMFyHO3gO0BTCcrZC00mIQqKjheuUESakhZPsk8mnKbYRfi07wO5xFXmrzWY5AJsHlhAGVMbMbb7Atr6bucP52SXkLLuDtEi0OmY8CpNuVZI8Np9gkDzmEzWRlE4hd34aDVN4RJ7w2V0GX3wXGVwgGXRk+SvXCvCd8i+VDiByvojMsQGYCQld23fmrvNjIAT14wbRjWcIJwR0nTzL8bIrboCylAIYPPfkMqGkS8IDZKWQmcJXX50ZcgNmOcyELzPGgLQrqC3i7zsO8E+z+KfXRVUIzmnogokKYw4addsFZkPHDSs4ooBfMS3hqr1kiDGmsWSxiEqGyNMgTb3pNm6qmCNTW8K7lrNbr2iKPCGr0ajTwo3a7eT3oKZNm30yaUNrPpuAdNsqd5um5cjLHFbdKgXSxFOgC1Ebq+OrKK2R0U9RXDEbfFEbvqj1nVgdcsZk70qNALsfC17/6LJuwyysadm6UMsmp2A0VdFVQkXXuxv/w7ZqWrM99JWnhwgc5Fprih+81CzVAGwohGWpPd+j1AKjL3Q2KpQ4+JrU+nvubExD0jNdWChvppoqaYq5/6c1c0TkgYgGgrwDbA19msNwYh79uPir5vtxAaxaPcUNiceuhNoNusU9P8L6EXVNSttWDz24hkS8UO0eRHL3HwDkKaAtnd1pEG+UDnco+mQ6wZwVWP58qIKy51+dejQom6nRxYnSZqFY8y7kDUPqvuMXolN3yKpe9Ih0JBp27d1tjebTd8w6aFdbsi1OlK89HIe6r4Y2zHEOM39FtXqnNlzn8LSlvQ7P6ti00ypz5EYsjePYITZWne3Mymt19WfR9mnoTStkr25xDgpVD6mB3ypUJDPzr7GUaEha5mUM4l3MWOc6ZFsKX0/V3aas44W+o+vRXe5nbDEAazUDsODi9FmOzmdbdXdnaUn1bJrTTltDR9m1Eu8sR5XCWHvevq6g3kNr1+RSryyrqqQ1dPjY/HPjaBXyCEByIG1LSb9O1Kx0syLpp2hyhUqgkB4CKDzS/n6PcEuPDqhM9krKWU017m5XLft5htAYBfOTnvMrS2vaTqUTWYlac4W2ZuJ05SU7PSvZw+BPNs4FHnh/Uj6pk4oqY08KZ0Gr9bJCm+4qn5PWPc3q2NJb1wyjbkK5A1kYxVCABDMUQ887qzjvy9TRms41ILfwIHKiYMe+i3UcRuywrns1wWBqOK7jsBbN2FDjYwjbmWBAFXgY+jrfLfO0vH79cPfwYHmOv476hRGnKRM6wWJZm+D5YkSrxX2IGI3UmOOGwN4ADuUPVLv+IVR+pIaA6s2PyAWIyDkNk3NQoEX1a8teOXM/sOxoZm/CEFy/Mr+jEum5wvHkCzZIrWu/QueK9tDzoJqeB4qtvKoEhdSaCr9AIU936W7GHPuOs84Q8ZfyqyQ6StrhloDBOqWZr8OhW9qzMsfcu9jx4KcBCZQ8kPtRKMJEDsaxGta79/0/MulEHAiGphiiIzitG1qBO7Xmro1v8aXiYKZbY/rSOM84/OBEItWk3JghTqZq5BU7VDmBVMBxxq1wLKVudie23Cj8QKWjO4iJjULpMMi6JZT9dKhOa4OsYypeLd8jD+iFosenMg6usEafLpwlmAQKCoOs8UWBlhdzR/Mslrbl+AP1Zm7RIzQb8cy7w4MEMkyjcqPcGXTdpYGxukoodEDmvonmoL40VX1lpEuo8BxsHZUWHDtWZjqc/n8oY3Pqycv036m1Ncpsk2uXVs9oNi1J4YVZvQLNkk4eBKjdoPVrH7IpKQbmo7v/EGV3ZJq5pm6D7yrosDH12HZrDu2pyCTGRMhAdO8I62iMEtokbH1gPli2PuTm6wNoctn6MISuj+6On+bduN3m8Xht2j+sa/9C9RIg63mLV8iWzPqsSIElrEjB9i6QFbXe29z69pvMBtjz3YGDGIM3eqgStpcLmHppKgDM3J4z4oO2JCGSSohWsb2AHFqwP6SvoYSIfw18P8o+PbDX89/9iRs/4z8=</diagram><diagram id="ls9p14zF9e6sMFj8i8uX" name="Phangs_pluggins">7Vpbc5s6EP41fjkz9QDiYh4TJ2460/bkxKfT5MkjgwxquB0hfOmvrzDiJiC+HJu4HefBQctKSPt9u9oVDMDYX38kMHK/hDbyBopkrwfgbqAosmQq7F8q2WSSETAzgUOwzZVKwRT/RHlPLk2wjeKaIg1Dj+KoLrTCIEAWrckgIeGqrrYIvfpTI+ighmBqQa8p/Y5t6vJVKEYpf0DYcfMnyzpfnw9zZb6S2IV2uKqIwP0AjEkY0uzKX4+Rlxovt0vWb9Jxt5gYQQHdp8Noktj/Pgdfp98/0Zunb/Ldw/zbBz7KEnoJX7CdxBT5s8hLHBzMIhcGTjxjIs+DZIZjsuCroZvcRGxhUXoZkdBCMbP37crFFE0jaKXiFSMGk0EPOwFrWmy6iDDBEhGKmaFv+A0f23Y65G2ckUAaSmCk6rJW/qrs7iIM6JQ/XE61KQlfC2xAqoE9bxx6IdlOD0wmY/bH5DYmjCE4TB+GYMyMdstXz2aC1p1mlQuwGMtR6CNKNkwl76ByfDnBNd5clWxRdC5zK0wBXAY5QZ1i5BJDdsFhPABSpQFpAzEU2Depb6RweDCOscVs4VLf4zYlYRLYKH2IxFpojelz5fplC47GW3fryq27Td4I2FKeq42sl6Ll7bLftpV3zOaK7IZXCniw9YQJsdBublNIHETf0APt+FYA1Frwy2UEeZDiZX26baDyJzyGmC2kpI9Up49qCLzIlsl7Vd1bHMisD6RL5tA062NlpmiMxbgANxW1KFWI95/zSHp7ZoK6LNf02UU2gZLyBQrHewE4LLCx7YPiIEn8a3SrRzdDwE5qCW+jPsOb2gUsQdCeIR97qISVySPmnii+wlqDVRFdsm3X6hVW7U1YrblsXkE9ENTWVORcoD7h108v6Mf4o//of/n633rx88d9S3ZZYDfDwSIc0jXtBNHaeJhlIQR0wFjJV+ZZvvJ5Xgig9epss5i/E8qGKSDKAJW1Jji6fqtMJqdBQlfqSIAmEG0p4ehcODRTwnvLlYdA+mu4wLTbj/4gCOT3xqAzIWEZ2WwFl8hDgUPda1ATghoQglqfCcg/yycmWT8A8uHGsRYvxubzU0tQa0B2WH3VMCjzA1CYvnJnPDbNrX/UQVGOrdFY4xERzMySMqZet8n1uu08Zdtbm8bOsk1t500/ZZsmlm0i3fYt23RhIE0VBuqo2Y6ok1rZfPLTgiub27bey2azLgRZcRPcl8yqcAYhG1qvZG7usdZckWTzD05zxO1R7S/NaYWgWZ5f48mR8eQtiu+MJ9olxRP5dw0onacS16T9kKQdvHvSrp86LLWFFPniX4nsGz3eNRspDiM5e8wjg0cRdYr3DsZwv/Bx8AsR4R1gftjQOTNRXzvtG5FWFzD6cAFpaFz+i8HfYg9tsFcsDPd1g+JTinwg/WxuIG774hcJO/TzTPWsbjC6aDfo0QvUphd0H1xevWB/LxCnbOzwAiHHZbvB+b1APvkp5im/EjF6dAOt6QbdpcB7uUGDvSNFZO/ejiDkHmyOvbiBsssNBH2g/K/NgDXLT+sy9fIDRXD/Cw==</diagram></mxfile>
3 3 \ No newline at end of file
... ...
LabTools/IRAP/JPB/make_phangs_ssps_isrf_prediction.pro
  Diff comments   View file @14ae4f2
1 1 PRO make_phangs_ssps_isrf_prediction,source_name=source_name,save=save,help=help
2 2  
3   -;make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save
4   -
5 3 ;+
6 4 ; NAME:
7 5 ; make_phangs_ssps_isrf_prediction
8 6 ; CALLING SEQUENCE:
9   -; make_phangs_ssps_isrf_prediction[,source_name=][,/save]
  7 +; make_phangs_ssps_isrf_prediction[,source_name=][,/save][,/help]
10 8 ; PURPOSE:
11 9 ; predicts the minimum ISRF in voronoi bins
12 10 ; INPUTS:
... ... @@ -21,11 +19,13 @@ PRO make_phangs_ssps_isrf_prediction,source_name=source_name,save=save,help=help
21 19 ; OPTIONAL OUTPUT:
22 20 ; None
23 21 ; PROCEDURE AND SUBROUTINE USED
24   -; None
  22 +; ISRF saved is the one predicted by dustem_plugin_phangs_stellar_isrf for each Muse Voronoi bin,
  23 +; scaled in amplitude to match the Muse filter data for that Voronoi bin.
  24 +; The G0 value computed is wrt the Mathis field at 1 mic.
25 25 ; SIDE EFFECTS:
26   -; None
  26 +; produces file _isrf_min_prediction.sav
27 27 ; EXAMPLE:
28   -; make_phangs_ssps_isrf_prediction,/save
  28 +; make_phangs_ssps_isrf_prediction,source_name='ngc0628',/save
29 29 ; MODIFICATION HISTORY:
30 30 ; written by Jean-Philippe Bernard
31 31 ;-
... ... @@ -36,12 +36,12 @@ IF keyword_set(help) THEN BEGIN
36 36 ENDIF
37 37  
38 38 ;===== predicts the minimum ISRF in voronoi bins
  39 +;===== The prediction is scaled to the Muse filter data observations
39 40  
40 41 win=0L
41 42 ;window,win,xsize=900,ysize=1000 & win=win+1
42 43  
43 44 ;=== This is where the data is read from and the ISRFs will be stored
44   -;data_dir='/Volumes/PILOT_FLIGHT1/PHANGS-JWST/DR1/'
45 45 data_dir=!phangs_data_dir+'/ISRF/WORK/'
46 46  
47 47 use_model='DBP90' ;Example with default keywords uses the DBP90 model
... ... @@ -50,9 +50,6 @@ use_polarization=0 ; initialize Dustemwrap in no polarization mode
50 50 ;== INITIALISE DUSTEM
51 51 dustem_init,model=use_model,polarization=use_polarization,show_plots=show_plots
52 52  
53   -;IF keyword_set(from_restore) THEN goto,from_restore
54   -;IF keyword_set(from_classes_restore) THEN goto,from_classes_restore
55   -
56 53 use_source_name='ngc0628'
57 54 IF keyword_set(source_name) THEN use_source_name=source_name
58 55  
... ... @@ -69,8 +66,27 @@ obp=[1.1,0,1.15,1]
69 66 ;dustem_init,show_plots=show_plots
70 67 ;needed only for NHCO
71 68 restore,data_dir+use_source_name+'_jwst_images.sav',/verb
  69 +;% RESTORE: Restored variable: JWST_IMAGES.
  70 +;% RESTORE: Restored variable: FILTERS.
  71 +;% RESTORE: Restored variable: HREF.
  72 +;% RESTORE: Restored variable: NHCO.
72 73 ;needed for stellar parameters and voronoi bin info
73 74 restore,data_dir+use_source_name+'_muse_images.sav',/verb
  75 +;% RESTORE: Restored variable: ST_TEMPLATES.
  76 +;% RESTORE: Restored variable: ST_MUSE_WEIGHTS.
  77 +;% RESTORE: Restored variable: VORONOI_ID.
  78 +;% RESTORE: Restored variable: AGE_VALUES.
  79 +;% RESTORE: Restored variable: METALICITY_VALUES.
  80 +;% RESTORE: Restored variable: BINS.
  81 +;% RESTORE: Restored variable: HREF.
  82 +;=== needed for absolute scaling of ISRF
  83 +restore,data_dir+use_source_name+'_muse_filters_data.sav',/verb
  84 +;% RESTORE: Restored variable: MUSE_IMAGES.
  85 +;% RESTORE: Restored variable: MUSE_FILTERS.
  86 +;% RESTORE: Restored variable: HREF.
  87 +restore,data_dir+use_source_name+'_muse_seds_muse_pixels.sav',/verb
  88 +;% RESTORE: Restored variable: ALL_SEDS.
  89 +;% RESTORE: Restored variable: ALL_SEDS_INDICES.
74 90  
75 91 Nvor=max(voronoi_id)
76 92  
... ... @@ -86,7 +102,15 @@ lambir=dustem_get_wavelengths(isrf_wavelengths=isrf_wavelengths)
86 102 show_each=100
87 103 Nlamb=n_elements(isrf_wavelengths)
88 104  
89   -ISRFS=fltarr(Nlamb,Nvor)
  105 +ISRFS=dblarr(Nlamb,Nvor)
  106 +G0s=dblarr(Nvor)
  107 +
  108 +;==== get Mathis field for G0 calculations
  109 +file=!dustem_soft_dir+'/data/ISRF_MATHIS.DAT'
  110 +readcol,file,Mathis_wavs,Mathis_ISRF
  111 +Mathis_ISRF=interpol(Mathis_ISRF,Mathis_wavs,isrf_wavelengths)
  112 +Mathis_1mic=interpol(Mathis_ISRF,isrf_wavelengths,1.0)
  113 +
90 114 FOR vid=first_vid,Nvor-1 DO BEGIN
91 115 IF vid mod show_each EQ 0 THEN BEGIN
92 116 t2=systime(0,/sec)
... ... @@ -97,22 +121,44 @@ FOR vid=first_vid,Nvor-1 DO BEGIN
97 121 ENDIF
98 122 ;===== get the SSP weights for the given voronoi bins
99 123 weights=phangs_binid2weights(st_muse_weights,vid,st_templates,age_values,metalicity_values,reddening=reddening)
100   - amplitude=1.
101   - reddening=0. ;This sets Muse reddening to 0. Will be contsrained later
102   - fpd=phangs_stellar_continuum_plugin_weight2params(weights,parameter_values=val,redenning=reddening,/force_include_reddening,amplitude=amplitude,/force_include_amplitude)
  124 + amplitude=1. ;amplitude is set to 1, and will be corrected by G0 below
  125 + use_reddening=0. ;This sets Muse reddening to 0. Will be contsrained later
  126 + fpd=phangs_stellar_continuum_plugin_weight2params(weights,parameter_values=val,redenning=use_reddening,/force_include_reddening,amplitude=amplitude,/force_include_amplitude)
103 127 Nparams=n_elements(val)
104 128 key=intarr(Nparams)
105 129 FOR i=0L,Nparams-1 DO BEGIN
106 130 toto=dustem_parameter_description2type(fpd[i],string_name=string_name,key=one_key)
107 131 key[i]=one_key
108 132 ENDFOR
109   - ISRFS[*,vid]=dustem_plugin_phangs_stellar_isrf(key=key,val=val,object_distance=object_distance,object_thickness=object_thickness)
  133 + ;==== compute the prediction stellar continuum for comparison to the Muse filter data.
  134 + ;==== Here we use the reddening from MUSE
  135 + val_sed=val
  136 + val_sed[1]=reddening ;for the SED, we use E(B-V)
  137 + spectrum=dustem_plugin_phangs_stellar_continuum(key=key,val=val_sed)
  138 + sed=interpol(spectrum[*,0],lambir,dustem_filter2wav(muse_filters[0:2]))
  139 + muse_sed=(*all_seds[vid]).stokesI
  140 + facts=muse_sed/sed
  141 + muse_factor=la_mean(facts) ;This takes the average value over the MUSE bands
  142 + ;print,G0,reddening
  143 + ;stop
  144 + val_isrf=val ;for the ISRF, we use E(B-V)=0
  145 + ISRFS[*,vid]=dustem_plugin_phangs_stellar_isrf(key=key,val=val_isrf,object_distance=object_distance,object_thickness=object_thickness)
  146 + ISRFS[*,vid]=ISRFS[*,vid]*muse_factor ;=== scale the amplitude of ISRFs to the muse data
  147 + ;=== compute corresponding G0
  148 + G0=interpol(ISRFS[*,vid],isrf_wavelengths,1.)/Mathis_1mic
  149 + G0s[vid]=G0
110 150 ENDFOR
111 151  
112 152 IF keyword_set(save) THEN BEGIN
113   - save,ISRFS,object_distance,object_thickness,source_name,file=data_dir+'ngc0628_isrf_min_prediction.sav',/verb
  153 + file_save=data_dir+use_source_name+'_isrf_min_prediction.sav'
  154 + save,ISRFS,G0s,object_distance,object_thickness,use_source_name,file=file_save,/verb
  155 + phangs_make_isrf_product,use_source_name,voronoi_id,ISRFS,G0s,object_distance,object_thickness
114 156 ENDIF
115 157  
  158 +;cgplot,lindgen(Nvor),G0s,/ylog,yrange=[1.e-5,1],xtit='Voronoi ID',ytit='G0'
  159 +
  160 +stop
  161 +
116 162 the_end:
117 163  
118 164 END
119 165 \ No newline at end of file
... ...
LabTools/IRAP/JPB/make_sed_phangs_tables.pro
  Diff comments   View file @14ae4f2
... ... @@ -2,6 +2,7 @@ PRO make_sed_phangs_tables,test=test,show_seds=show_seds,grid_type=grid_type,isr
2 2  
3 3 ;make_sed_phangs_tables
4 4 ;make_sed_phangs_tables,/test,grid_type=2,/show_seds,isrf_class=15
  5 +;make_sed_phangs_tables,/test,grid_type=3,/show_seds
5 6  
6 7 ;=== This is needed for Phangs work
7 8 ;defsysv,'!phangs_data_dir','/Volumes/PILOT_FLIGHT1/PHANGS/'
... ... @@ -18,24 +19,20 @@ IF keyword_set(grid_type) THEN BEGIN
18 19 ENDIF
19 20  
20 21 CASE use_grid_type OF
21   - 1: BEGIN ;======= This is to produce a grid with the DBP90 model for G0 and PAH abundance
  22 + 1: BEGIN ;======= This is to produce a grid with the DBP90 model + Mathis field for G0 and PAH abundance
22 23  
23 24 !quiet=1
24 25 model='DBP90'
25   -
26 26 ;===== This is to do a small test_table
27 27 pd = [ $
28 28 '(*!dustem_params).G0', $ ;G0
29 29 '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
30 30 '(*!dustem_params).grains(1).mdust_o_mh'] ;VSG mass fraction
31   -
32 31 iv_min = [0.1,1.e-4,1.e-4]
33 32 iv_max = [100,1.e-1,1.e-1]
34 33 log=[1,1,1]
35   -
36 34 fpd=[]
37 35 fiv=[]
38   -
39 36 IF keyword_set(test) THEN BEGIN
40 37 iv_Nvalues=[5,3,3]
41 38 table_name=dir+'TEST_'+model+'_JWST_G0_YPAH_YVSG_4Phangs.fits'
... ... @@ -43,17 +40,14 @@ CASE use_grid_type OF
43 40 iv_Nvalues=[50,20,20]
44 41 table_name=dir+model+'_JWST_G0_YPAH_YVSG_4Phangs.fits'
45 42 ENDELSE
46   -
47 43 dustem_init,model=model
48 44 ;=== select filters to be used for the grid
49 45 filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names]
50   -
51 46 dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,log=log,show_seds=show_seds
52 47 END
53   - 2: BEGIN
  48 + 2: BEGIN ;This is DPB90 model with PHANGS ISRF classes
54 49 !quiet=1
55 50 model='DBP90' ;This is to use the DBP90 dust model
56   -
57 51 ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
58 52 ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
59 53 pd = [ $
... ... @@ -63,7 +57,6 @@ CASE use_grid_type OF
63 57 iv_min = [0.1,1.e-4,1.e-4]
64 58 iv_max = [100,1.e-1,1.e-1]
65 59 log=[1,1,1]
66   -
67 60 Nclass=31L
68 61 fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0','dustem_plugin_phangs_class_isrf_1'] ;ISRF class to be used
69 62 dustem_init,model=model
... ... @@ -73,7 +66,6 @@ CASE use_grid_type OF
73 66 ;=== select filters to be used for the grid
74 67 filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names, $
75 68 (*!dustem_filters).iras.filter_names,(*!dustem_filters).pacs.filter_names,(*!dustem_filters).spire.filter_names]
76   -
77 69 IF keyword_set(isrf_class) THEN BEGIN ;do just one table
78 70 use_isrf_class=isrf_class
79 71 ENDIF ELSE BEGIN
... ... @@ -91,6 +83,47 @@ CASE use_grid_type OF
91 83 ENDELSE
92 84 dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,log=log,show_seds=show_seds
93 85 END
  86 + 3: BEGIN ;This is DPB90 model with DL07 ISRF prescription (a la Chastanet)
  87 + !quiet=1
  88 + model='DBP90' ;This is to use the DBP90 dust model
  89 + ;======= This is to produce a grid with the above model for G0 and PAH abundance and ISRF from the Muse data
  90 + ;======= Note that g0 is fixed parameters is actually not used, replaced by dustem_plugin_phangs_class_isrf_2
  91 + pd = [ $
  92 + 'dustem_plugin_dl07_isrf_model_1', $ ;gamma factor
  93 + 'dustem_plugin_dl07_isrf_model_3', $ ;Umin factor
  94 + '(*!dustem_params).grains(0).mdust_o_mh',$ ;PAH0 mass fraction
  95 + '(*!dustem_params).grains(1).mdust_o_mh' $ ;VSG mass fraction
  96 + ]
  97 + ;umin: 0.1-100
  98 + ;gamma: 0.0,0.02
  99 + ;umax: 1.e7
  100 + ;alpha: 2.
  101 + iv_min = [0.1,0. ,1.e-4,1.e-4]
  102 + iv_max = [100,0.02,1.e-1,1.e-1]
  103 + log=[1,1,1,1]
  104 + fpd=['(*!dustem_params).gas.g0','(*!dustem_params).g0', $
  105 + 'dustem_plugin_dl07_isrf_model_2', $ ;alpha
  106 + 'dustem_plugin_dl07_isrf_model_4' $ ;Umax
  107 + ]
  108 + dustem_init,model=model
  109 + fortran_user=dustem_set_up_fortran(/random_name) ;use a random fortran number
  110 + !dustem_verbose=0
  111 + (*!dustem_params).KEYWORDS='quiet '+(*!dustem_params).KEYWORDS ;This makes Fortran be quiet too
  112 + ;=== select filters to be used for the grid
  113 + filters=[(*!dustem_filters).nircam.filter_names,(*!dustem_filters).miri.filter_names, $
  114 + (*!dustem_filters).iras.filter_names,(*!dustem_filters).pacs.filter_names,(*!dustem_filters).spire.filter_names]
  115 + fiv=[-1.,1.,2.,1.e7] ;This sets the fixed parameters
  116 + ;define the number of free parameters in the grid and the grid fits table name
  117 + IF keyword_set(test) THEN BEGIN
  118 + iv_Nvalues=[2,2,3,3]
  119 + table_name=dir+'TEST_'+model+'_DL07ISRF_G0_YPAH_YVSG_4Phangs'+'.fits'
  120 + ENDIF ELSE BEGIN
  121 + iv_Nvalues=[6,5,20,20]
  122 + table_name=dir+model+'_DL07ISRF_G0_YPAH_YVSG_4Phangs'+'.fits'
  123 + ENDELSE
  124 + dustem_make_sed_table,model,pd,iv_min,iv_max,iv_Nvalues,fpd=fpd,fiv=fiv,filename=table_name,filters=filters,log=log,show_seds=show_seds
  125 + END
  126 +
94 127 ENDCASE
95 128  
96 129  
... ...
LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
  Diff comments   View file @14ae4f2
1   -PRO phangs_brute_force_fit_with_isrf_grid,source_name=source_name,normalize=normalize,from_restore=from_restore,save=save
  1 +PRO phangs_brute_force_fit_with_isrf_grid,source_name=source_name,normalize=normalize,from_restore=from_restore,save=save,force_mathis=force_mathis,help=help
2 2  
3 3 ;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/normalize,/save
  4 +;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/normalize,/force_mathis,/save
  5 +
4 6 ;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/save
5 7  
  8 +;+
  9 +; NAME:
  10 +; phangs_brute_force_fit_with_isrf_grid
  11 +; CALLING SEQUENCE:
  12 +; phangs_brute_force_fit_with_isrf_grid[,source_name=][,/normalize][,/save][,/help][,/from_restore][,/force_mathis]
  13 +; PURPOSE:
  14 +; Brute-Force fit SEDs of a given phangs Galaxy using ISRF grid
  15 +; INPUTS:
  16 +; None
  17 +; OPTIONAL KEYWORDS:
  18 +; save = if set, save the classes
  19 +; help = if set, print this help
  20 +; OUTPUTS:
  21 +; None
  22 +; OPTIONAL INPUT:
  23 +; source_name = source name (default='ngc0628')
  24 +; OPTIONAL OUTPUT:
  25 +; None
  26 +; PROCEDURE AND SUBROUTINE USED
  27 +;
  28 +; SIDE EFFECTS:
  29 +; produces file *_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses*.sav
  30 +; EXAMPLE:
  31 +; phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/normalize,/save
  32 +; MODIFICATION HISTORY:
  33 +; written by Jean-Philippe Bernard
  34 +;-
  35 +
  36 +IF keyword_set(help) THEN BEGIN
  37 + doc_library,'phangs_brute_force_fit_with_isrf_grid'
  38 + goto,the_end
  39 +ENDIF
  40 +
6 41 use_source_name='ngc0628'
7 42 IF keyword_set(source_name) THEN use_source_name=source_name
8 43  
... ... @@ -10,6 +45,10 @@ norm_str=&#39;&#39;
10 45 IF keyword_set(normalize) THEN BEGIN
11 46 norm_str='_normalized'
12 47 ENDIF
  48 +mathis_str=''
  49 +IF keyword_set(force_mathis) THEN BEGIN
  50 + mathis_str='_mathis'
  51 +ENDIF
13 52 pdp_define_la_common
14 53  
15 54 data_dir=!phangs_data_dir+'/ISRF/WORK/'
... ... @@ -19,12 +58,20 @@ IF keyword_set(from_restore) THEN GOTO,from_restore
19 58  
20 59 dustem_init,show_plots=show_plots
21 60  
22   -;==== needed for NHCO
23   -restore,data_dir+use_source_name+'_jwst_images.sav',/verb
  61 +;==== needed for NHCO,NHHI
  62 +restore,data_dir+use_source_name+'_co_images.sav',/verb
  63 +;% RESTORE: Restored variable: NHCO.
  64 +;% RESTORE: Restored variable: HREF.
  65 +restore,data_dir+use_source_name+'_hi_images.sav',/verb
  66 +;% RESTORE: Restored variable: NHI.
  67 +;% RESTORE: Restored variable: HREF.
  68 +
  69 +;restore,data_dir+use_source_name+'_jwst_images.sav',/verb
24 70 ;% RESTORE: Restored variable: JWST_IMAGES.
25 71 ;% RESTORE: Restored variable: FILTERS.
26 72 ;% RESTORE: Restored variable: HREF.
27 73 ;% RESTORE: Restored variable: NHCO.
  74 +
28 75 restore,data_dir+use_source_name+'_muse_images.sav',/verb
29 76 ;% RESTORE: Restored variable: ST_TEMPLATES.
30 77 ;% RESTORE: Restored variable: ST_MUSE_WEIGHTS.
... ... @@ -55,14 +102,12 @@ restore,data_dir+use_source_name+&#39;_isrf_min_prediction.sav&#39;,/verb
55 102 ;% RESTORE: Restored variable: SOURCE_NAME.
56 103  
57 104 ;For test
58   -;class_min=13
59   -;class_max=13
60   -;class_min=13
61   -;class_max=13
  105 +class_min=8 & class_max=18
62 106  
63   -;class_min=8 & class_max=17
64   -;class_min=15 & class_max=15
65   -class_min=8 & class_max=15
  107 +IF keyword_set(force_mathis) THEN BEGIN
  108 + class_min=0L
  109 + class_max=0L
  110 +ENDIF
66 111  
67 112 restore,data_dir+use_source_name+'_jwst_seds_muse_pixels.sav',/verb
68 113 ;% RESTORE: Restored variable: ALL_SEDS.
... ... @@ -80,7 +125,7 @@ restore,data_dir+use_source_name+&#39;_astrosat_seds_muse_pixels.sav&#39;,/verb
80 125 astrosat_seds=all_seds
81 126 astrosat_filters=(*astrosat_seds[0]).filter
82 127  
83   -Nvor=max(voronoi_id)
  128 +Nvor=long(max(voronoi_id))
84 129 all_filters=[jwst_filters,astrosat_filters,muse_filters]
85 130 Nfilters=n_elements(all_filters)
86 131 all_seds=fltarr(Nfilters,Nvor)
... ... @@ -120,7 +165,8 @@ G0_hit=fltarr(Nvor)+la_undef()
120 165 Ypah_hit=fltarr(Nvor)+la_undef()
121 166 Yvsg_hit=fltarr(Nvor)+la_undef()
122 167  
123   -use_NHmap=NHCO ;used NH map in 1.e21 (from CO)
  168 +;use_NHmap=NHCO ;used NH map in 1.e21 (from CO)
  169 +use_NHmap=la_add(NHCO,NHI) ;used NH map in 1.e21 (from CO+HI)
124 170 ;use_NHmap=NH_map ;used NH map in 1.e21 (from MUSE)
125 171  
126 172 ;==== This is all the filters in a given SED
... ... @@ -143,7 +189,11 @@ model=&#39;DBP90&#39;
143 189  
144 190 bidon=dustem_get_wavelengths(isrf_wavelengths=isrf_wavelengths)
145 191  
146   -fixed_parameters_description=['dustem_plugin_phangs_class_isrf_2'] ;This is the ISRF G0 factor
  192 +IF not keyword_set(force_mathis) THEN BEGIN
  193 + fixed_parameters_description=['dustem_plugin_phangs_class_isrf_2'] ;This is the ISRF G0 factor, fixed for grid
  194 +ENDIF ELSE BEGIN
  195 + ;left blanck on purpose
  196 +ENDELSE
147 197  
148 198 file=!dustem_soft_dir+'/data/ISRF_MATHIS.DAT'
149 199 ;file='/Users/jpb/Soft_Libraries/dustem_fortran/data/ISRF.DAT'
... ... @@ -158,7 +208,13 @@ FOR isrf_class=class_min,class_max DO BEGIN
158 208 isrf_class_str='_isrfclass'+strtrim(isrf_class,2)
159 209 table_name=grids_data_dir+model+'_MuseISRF_JWST_G0_YPAH_YVSG_4Phangs'+isrf_class_str+'.fits'
160 210 ;select Voronoi bins with that ISRF class
161   - ind_class=where(vor_class EQ isrf_class,Nvor_class)
  211 + IF not keyword_set(force_mathis) THEN BEGIN
  212 + ind_class=where(vor_class EQ isrf_class,Nvor_class)
  213 + ENDIF ELSE BEGIN
  214 + ;stop
  215 + ind_class=lindgen(Nvor)
  216 + Nvor_class=Nvor
  217 + ENDELSE
162 218 reset=1
163 219 IF Nvor_class NE 0 THEN BEGIN
164 220 FOR vvid=0LL,Nvor_class-1 DO BEGIN
... ... @@ -167,13 +223,10 @@ FOR isrf_class=class_min,class_max DO BEGIN
167 223 ENDIF
168 224 ;=== Get the G0 value for that vid
169 225 vid=ind_class[vvid]
170   - ;print,vid
171   - ;IF vid EQ 12641 THEN BEGIN
172   - ; message,'stopped at vid =12641',/continue
173   - ; stop
174   - ;ENDIF
175 226 G0=interpol(ISRFs[*,vid],isrf_wavelengths,1.0)/Mathis_1mic
176   - fixed_parameters_values=[G0]
  227 + IF not keyword_set(force_mathis) THEN BEGIN
  228 + fixed_parameters_values=[G0]
  229 + ENDIF
177 230 sed=*seds_ptr[vid]
178 231 ;=== restrict sed to requested filters
179 232 sed=sed[ind_filters]
... ... @@ -193,8 +246,11 @@ FOR isrf_class=class_min,class_max DO BEGIN
193 246 ,fixed_parameters_description=fixed_parameters_description,fixed_parameters_values=fixed_parameters_values $
194 247 ,weighted_params=weighted_params,reset=reset)
195 248 reset=0
196   - G0s[vid]=G0
197   - ;Ypahs[vid]=params[1]/fact ;extenssive quantities must be divided by normalization factor
  249 + IF not keyword_set(force_mathis) THEN BEGIN
  250 + G0s[vid]=G0
  251 + ENDIF ELSE BEGIN
  252 + G0s[vid]=weighted_params[0]
  253 + ENDELSE
198 254 Ypahs[vid]=weighted_params[1];/fact ;Caution, extenssive quantities must be divided by normalization factor. We don't here
199 255 Yvsgs[vid]=weighted_params[2];/fact
200 256 facts[vid]=fact
... ... @@ -246,19 +302,19 @@ ENDFOR
246 302  
247 303 ;save,GOs,Ypahs,Yvsgs,file=dir+source_name+'DBP90_JWST_G0_YPAH_YVSG.fits'
248 304 IF keyword_set(save) THEN BEGIN
249   - file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses'+norm_str+'.sav'
  305 + file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses'+norm_str+mathis_str+'.sav'
250 306 save,G0s,Ypahs,Yvsgs,facts,chi2s,rchi2s,dG0s,dYpahs,dYvsgs,G0_hit,Ypah_hit,Yvsg_hit,file=file_save
251 307 message,'Wrote '+file_save,/continue
252   - file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+'.sav'
  308 + file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+mathis_str+'.sav'
253 309 save,G0_map,Ypah_map,Yvsg_map,fact_map,chi2_map,rchi2_map,dG0_map,dYpah_map,dYvsg_map,Ypah_hit_map,Yvsg_hit_map,file=file_save
254 310 message,'Wrote '+file_save,/continue
255 311 ENDIF
256 312  
257 313 from_restore:
258 314  
259   -file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses'+norm_str+'.sav'
  315 +file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses'+norm_str+mathis_str+'.sav'
260 316 restore,file_save,/verb
261   -file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+'.sav'
  317 +file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+mathis_str+'.sav'
262 318 restore,file_save,/verb
263 319  
264 320 stop
... ... @@ -267,7 +323,7 @@ stop
267 323 win=0L
268 324 window,win,xsize=700,ysize=900 & win=win+1
269 325 !p.multi=[0,1,3]
270   -ind=where(G0s NE la_undef())
  326 +ind=where(G0s NE la_undef() and G0s GT 0)
271 327 res=histogram(G0s[ind],locations=xv)
272 328 cgplot,xv,res,psym=10,title='G0s histogram',xtit='G0',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
273 329  
... ... @@ -279,32 +335,24 @@ ind=where(Yvsgs NE la_undef())
279 335 res=histogram(alog10(Yvsgs[ind]),locations=xv,Nbins=100)
280 336 cgplot,xv,res,psym=10,title='Yvsgs histogram',xtit='log(Yvsg)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
281 337  
282   -Ypahs_f=la_div(Ypahs,facts)
283   -ind=where(Ypahs_f NE la_undef())
284   -res=histogram(alog10(Ypahs_f[ind]),locations=xv,Nbins=100)
285   -cgplot,xv,res,psym=10,title='Ypahs histogram',xtit='log(Ypah_f)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
286   -
287   -Yvsgs_f=la_div(Yvsgs,facts)
288   -ind=where(Yvsgs_f NE la_undef())
289   -res=histogram(alog10(Yvsgs_f[ind]),locations=xv,Nbins=100)
290   -cgplot,xv,res,psym=10,title='Yvsgs histogram',xtit='log(Yvsg_f)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
291   -
292 338 ;==== show maps
  339 +obp=[1.1,0.,1.15,1]
293 340 win=0L
294   -window,win,xsize=800,ysize=900 & win=win+1
295 341  
296   -obp=[1.1,0.,1.15,1]
  342 +window,win,xsize=800,ysize=900 & win=win+1
  343 +imrange=[-4.5,-2.5]
297 344 ;image_cont20,G0_map,Href,/square,imrange=[-2,40],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='G0'
298   -image_cont20,la_log10(G0_map),Href,/square,imrange=[-0.2,2],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(G0)'
  345 +image_cont20,la_log10(G0_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(G0)'
299 346  
300 347 window,win,xsize=800,ysize=900 & win=win+1
301   -imrange=[-4.,-3.] ;normalize=1
302   -imrange=[-4.,-3.5]+2 ;normalize=0
  348 +imrange=[-4.,-2.] ;normalize=1
  349 +;imrange=[-4.,-3.5]+2 ;normalize=0
303 350 image_cont20,la_log10(Ypah_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Ypah)'
304 351  
305 352 window,win,xsize=800,ysize=900 & win=win+1
306   -imrange=[-4.,-3.]+0.5 ;normalize=1
307   -imrange=[-2.2,-1.6] ;normalize=0
  353 +;imrange=[-4.,-3.]+0.5 ;normalize=1
  354 +imrange=[-4.,-2.]+0.5 ;normalize=1
  355 +;imrange=[-2.2,-1.6] ;normalize=0
308 356 image_cont20,la_log10(YVSG_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Yvsg)'
309 357  
310 358 window,win,xsize=800,ysize=900 & win=win+1
... ... @@ -325,13 +373,20 @@ image_cont20,Yvsg_hit_map,Href,/square,imrange=imrange,axis_color_table=1,image_
325 373  
326 374 window,win,xsize=800,ysize=900 & win=win+1
327 375 imrange=[-2,2]
328   -image_cont20,use_NH_map,Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='Yvsg hit',off_bar=obp
  376 +image_cont20,la_log10(use_NHmap),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='NH map',off_bar=obp
329 377  
330 378 window,win,xsize=800,ysize=900 & win=win+1
331 379 imrange=0
332 380 imrange=[-2.0,-1.45]
333 381 image_cont20,la_log10(la_mul(Ypah_map,fact_map)),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='Ypah *fact',off_bar=obp
334 382  
  383 +;=== other plots
  384 +
  385 +
  386 +cgplot,G0s,Ypahs,/xlog,/ylog,xr=[1.e-5,1.e-2],/xsty,psym=3,xtit='G0',ytit='Ypah'
  387 +cgplot,G0s,Yvsgs,/xlog,/ylog,xr=[1.e-5,1.e-2],/xsty,psym=3,xtit='G0',ytit='Yvsg'
  388 +cgplot,Yvsgs,Ypahs,/xlog,/ylog,xr=[1.e-4,1.e-1],yr=[1.e-4,1.e-1],/xsty,psym=3,xtit='Yvsg',ytit='Ypah'
  389 +
335 390  
336 391 stop
337 392  
... ... @@ -359,4 +414,6 @@ wset,3 &amp; image_cont20,alog10(chi2_map),href,/square,imrange=[2,6],image_color_ta
359 414 window,4,xsize=800,ysize=1000
360 415 image_cont20,alog10(fact_map),href,/square,imrange=[-1,5],image_color_table='jpbloadct',/silent,title='fact',off_bar=obp,axis_color_table=1,levels=levs
361 416  
  417 +the_end:
  418 +
362 419 END
... ...
LabTools/IRAP/JPB/phangs_isrf_pipeline.pro
  Diff comments   View file @14ae4f2
... ... @@ -10,6 +10,8 @@ phangs_make_jwst_images,source_name=&#39;ngc0628&#39;,/save,/show,/nostop
10 10 phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
11 11 ;make _astrosat_data.sav
12 12 phangs_make_astrosat_images,source_name='ngc0628',/save,/show,/nostop
  13 +phangs_make_hi_images,source_name='ngc0628',/save,/show,/nostop
  14 +phangs_make_co_alma_images,source_name='ngc0628',/save,/show,/nostop
13 15  
14 16 ;==== extract SEDs on MUSE pixels
15 17 ;produces _jwst_seds_muse_pixels.sav, _muse_seds_muse_pixels.sav', _astrosat_seds_muse_pixels.sav
... ...
LabTools/IRAP/JPB/phangs_make_jwst_images.pro
  Diff comments   View file @14ae4f2
1   -PRO phangs_make_jwst_images,source_name=source_name,save=save,show_images=show_images,nostop=nostop
  1 +PRO phangs_make_jwst_images,source_name=source_name,save=save,show_images=show_images,nostop=nostop,help=help
2 2  
3   -;phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop
  3 +;+
  4 +; NAME:
  5 +; phangs_make_jwst_images
  6 +; PURPOSE:
  7 +; makes JWST images usable by the Phangs ISRF project
  8 +; CATEGORY:
  9 +; Dustem Phangs
  10 +; CALLING SEQUENCE:
  11 +; phangs_make_jwst_images[,source_name=][,/save][,/show_images][,/nostop]
  12 +; INPUTS:
  13 +; None
  14 +; OPTIONAL INPUT PARAMETERS:
  15 +; source_name : source name (default = 'ngc0628')
  16 +; OUTPUTS:
  17 +; None
  18 +; OPTIONAL OUTPUT PARAMETERS:
  19 +; None
  20 +; ACCEPTED KEY-WORDS:
  21 +; help = If set, print this help
  22 +; save = If set, save result
  23 +; show_images= if set, show images
  24 +; nostop = if set, does not stop
  25 +; COMMON BLOCKS:
  26 +; None
  27 +; SIDE EFFECTS:
  28 +; A file is written
  29 +; RESTRICTIONS:
  30 +; None
  31 +; PROCEDURE:
  32 +; input JWST images are asumed to be in units of MJy/sr
  33 +; EXAMPLES
  34 +; phangs_make_jwst_images,source_name='ngc0628',/save,/show,/nostop
  35 +; MODIFICATION HISTORY:
  36 +; Written by J.-Ph. Bernard (2023)
  37 +; Evolution details on the DustEMWrap gitlab.
  38 +; See http://dustemwrap.irap.omp.eu/ for FAQ and help.
  39 +;-
  40 +
  41 +IF keyword_set(help) THEN BEGIN
  42 + doc_library,'phangs_make_jwst_images'
  43 + goto,the_end
  44 +ENDIF
4 45  
5 46 use_source_name='ngc0628'
6 47 IF keyword_set(source_name) THEN use_source_name=source_name
... ... @@ -11,7 +52,6 @@ obp=[1.1,0,1.15,1]
11 52 win=0L
12 53  
13 54 ;data_dir='/Volumes/PILOT_FLIGHT1/PHANGS-JWST/DR1/'
14   -NH_data_dir=!phangs_data_dir+'/phangs_drive/PHANGS-JWST/DR1/'
15 55 data_dir=!phangs_data_dir+'/phangs_drive/PHANGS-JWST/DR1/'
16 56 save_data_dir=!phangs_data_dir+'/ISRF/WORK/'
17 57  
... ... @@ -27,7 +67,17 @@ IF keyword_set(show_images) THEN BEGIN
27 67 window,win & win=win+1
28 68 image_cont20,d,href,/square,imrange=[-0.1,5],image_color_table='jpbloadct',/silent,tit=tit
29 69 ENDIF
30   -;stop
  70 +
  71 +;=== save reference header
  72 +IF keyword_set(save) THEN BEGIN
  73 + save_file=save_data_dir+use_source_name+'_ref_header.sav'
  74 + save,href,file=save_file
  75 + message,'Saved '+save_file,/continue
  76 +ENDIF
  77 +
  78 +stop
  79 +
  80 +
31 81 Nx=sxpar(href,'NAXIS1')
32 82 Ny=sxpar(href,'NAXIS2')
33 83  
... ... @@ -228,5 +278,6 @@ IF keyword_set(save) THEN BEGIN
228 278 message,'Saved '+save_file,/continue
229 279 ENDIF
230 280  
  281 +the_end:
231 282  
232 283 END
233 284 \ No newline at end of file
... ...
LabTools/IRAP/JPB/phangs_make_muse_filters_images.pro
  Diff comments   View file @14ae4f2
1 1 PRO phangs_make_muse_filters_images,source_name=source_name,save=save,show_images=show_images,nostop=nostop
2 2  
3   -;phangs_make_muse_filters_images,source_name='ngc0628',/show,/nostop
4   -;phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
  3 +;+
  4 +; NAME:
  5 +; phangs_make_muse_filters_images
  6 +; PURPOSE:
  7 +; makes Muse filter images usable by the Phangs ISRF project
  8 +; CATEGORY:
  9 +; Dustem Phangs
  10 +; CALLING SEQUENCE:
  11 +; phangs_make_muse_filters_images[,source_name=][,/save][,/show_images][,/nostop]
  12 +; INPUTS:
  13 +; None
  14 +; OPTIONAL INPUT PARAMETERS:
  15 +; source_name : source name (default = 'ngc0628')
  16 +; OUTPUTS:
  17 +; None
  18 +; OPTIONAL OUTPUT PARAMETERS:
  19 +; None
  20 +; ACCEPTED KEY-WORDS:
  21 +; help = If set, print this help
  22 +; save = If set, save result
  23 +; show_images= if set, show images
  24 +; nostop = if set, does not stop
  25 +; COMMON BLOCKS:
  26 +; None
  27 +; SIDE EFFECTS:
  28 +; A file is written
  29 +; RESTRICTIONS:
  30 +; None
  31 +; PROCEDURE:
  32 +; input MUSE images are asumed to be in units of ?
  33 +;
  34 +; EXAMPLES
  35 +; phangs_make_muse_filters_images,source_name='ngc0628',/show,/nostop
  36 +; phangs_make_muse_filters_images,source_name='ngc0628',/save,/show,/nostop
  37 +; MODIFICATION HISTORY:
  38 +; Written by J.-Ph. Bernard (2023)
  39 +; Evolution details on the DustEMWrap gitlab.
  40 +; See http://dustemwrap.irap.omp.eu/ for FAQ and help.
  41 +;-
  42 +
  43 +IF keyword_set(help) THEN BEGIN
  44 + doc_library,'phangs_make_muse_filters_images'
  45 + goto,the_end
  46 +ENDIF
5 47  
6 48 use_source_name='ngc0628'
7 49 IF keyword_set(source_name) THEN use_source_name=source_name
... ... @@ -148,4 +190,6 @@ ENDFOR
148 190  
149 191 save,NH_map,file=save_data_dir+use_source_name+'_muse_NH.sav'
150 192  
  193 +the_end:
  194 +
151 195 END
152 196 \ No newline at end of file
... ...
LabTools/IRAP/JPB/srun/make_phangs_grids.pro
  Diff comments   View file @14ae4f2
... ... @@ -21,7 +21,11 @@ bidon=1
21 21  
22 22 ;make grid
23 23 ;make_sed_phangs_tables,/test,grid_type=2,/show_seds,isrf_class=15
24   -make_sed_phangs_tables,grid_type=2,isrf_class=19
  24 +make_sed_phangs_tables,grid_type=2,isrf_class=23
  25 +;make_sed_phangs_tables,grid_type=2,isrf_class=22
  26 +;make_sed_phangs_tables,grid_type=2,isrf_class=21
  27 +;make_sed_phangs_tables,grid_type=2,isrf_class=20
  28 +;make_sed_phangs_tables,grid_type=2,isrf_class=19
25 29 ;make_sed_phangs_tables,grid_type=2,isrf_class=18
26 30 ;make_sed_phangs_tables,grid_type=2,isrf_class=17
27 31 ;make_sed_phangs_tables,grid_type=2,isrf_class=16
... ...