improved

Jean-Philippe Bernard
1 parent 4cd2e9bc
Showing 6 changed files with 52 additions and 34 deletions Show diff stats
Docs/developers/Code_Structure_Flowcharts/dustem-wrapper_flow_chart1.drawio
LabTools/IRAP/JPB/check_phangs_ssps_isrf_prediction.pro
LabTools/IRAP/JPB/make_phangs_isrf_classes.pro
LabTools/IRAP/JPB/phangs_brute_force_fit_with_isrf_grid.pro
LabTools/IRAP/JPB/phangs_make_astrosat_images.pro
LabTools/IRAP/JPB/srun/make_phangs_grids.pro
-<mxfile host="Electron" modified="2024-02-27T10:03:36.573Z" 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="QFmzlgI8t46DTHHOp26z" version="17.2.1" type="device" pages="4"><diagram id="IKMg62xr7mvPSjzlxw8D" name="Dustem_wrap_flow">7V1bk6M2Fv4t+9CP02Xu8Jjpnp5kK8mmth+SeXLRINtUMBDAfdlfv+Jmgy4g04BkW1OVji3LQtb3naOjc46kO+1h//49dZPdb7EPwjt15b/faY93qqooKxP+ryj5qEosw6kKtmng15VOBc/B/0BduKpLD4EPsk7FPI7DPEi6hV4cRcDLO2VumsZv3WqbOOw+NXG3ACt49twQL/0z8PNdVWqr1qn8ZxBsd82TFbP+fXu3qVz/kmzn+vFbq0j7dqc9pHGcV6/27w8gLAavGZfqe0+UT48dS0GUs3zh56/6Iflrtd8mf/7+YH57Be6/rS91K69ueKh/sH/IcrBfb4J8nQF/nQLX34P7JI3rX5F/NEMDf1BSvISfeSCD4/z1bRfk4DlxvaL4DRIClrlhsI3gWw92E6Sw4BWkeQAH+Kf6g33g+0WTX7MK/NX9SrN1UzFOf3X46SaO8uf64UpRO0/jv4+YaEWNIAwf4jBOy+5pT08P8B8s94MUMiOIi4cBN4OD9RUfu2YgYOfAe6uoHsvvIN6DPP2AVZpPVzWuNbG15v3biSaqWZftWhTR6jK3Zub22PQJPPiixu8MLFUalkEU5BJBHEEVQVDnjaDehyAUyXcokPAnu/tMwjkIp63yhtOgwZmBfO27uStRHETRtHmjaGIoFnPi+t1zpVJlwE/jjZ9Fk8ISRjg6cPylPiUgaXeRNLhPjzbdWC1AlHYOBUjD7gLJ3VJ1ekVym7pBJHEcFEjd5I1js0inAOmG4foNvGiJLdHE0TQRNLmbqwrVGdCgKWHEYdSRRSR3e1XB/QAYZiDyfyqcYwUgoZtlgQdHY5fvw3pU0/gQ+aB4ygq+g6OTfvzVfvOjBMho3j6+tz98/KjfUYc3d9MtyHt+Q91j4HecczgIrUE2CGPclKUgdPPgtevSIw18/YQ/4gD2mO7qQbHL4kPqgfpbbRcc0pCOGscW0lA1MFhDECv3o1UtKSpk9A5b5P7SuoVWNwyEhNXzT5Q8IvAJluK+Dk4sBe9B/lfrdfUty6jfnr5VvBmkdkWFYWoPikCjjQWRAYQjljpSBDAvjWLMIgKYzJpaf7/Q+s4SQkD1EJUzbujuX3zpJGJYmvKfdHEv0cVNuqZQCgcTyNEax0A1wTwaB30MGo7rr27qS+gb3Bd2G5OuwioDzYpWFCFAVhf6aqwQIPO3rqn3M028Gtplpb9n6E/UFhCERpfwF4Tx6toRm6m6NpKpqJdrNRNPUf+LPmAgUn4gc31rCV5TXThlBBnSLIg2sbQoB8NWFvewVaORyNHjzSHyIKa+hHLQI2dxXxyoGhVK9xWUaVYbGblikEp1pXDHkpqjUy7Z69fhq0QTQ1NFBFNx+CvZfgdMEcBKgUwOIGCJSKZiO9yxxF0wNZZvKQRmnSUyqkxCErG3oV3MHUlqyk5h+sgESFaZ1PjLJDVlx4MKFgplYccWAyHxJOCJeECVYxCGH57UzB2J57A1i9g/DveFSdMwbcqUKR8smTsO98RIjer2kRPmOQLJ3fLRpNNnGiht7rmRGtXps0+KzXVy/weTy0dZcTd6NKrLp0QSSqXEcVAgFYW7606jOntKHCWIg5aOcsyl4Aci1ctTqdVUSiNDwpSicHcNaFQnT72U9OAYBtHhsJd4MrgGuLvStSFXDxyDTep6Es1hNDX+cyV9i5ZUsEwQ6g5nCBsDTK4lz8USNV117g4energkY66c4Dk7nHVqQ6eI5ByyyQVT3Q5YnDfM6lTvTwVnuVeZokoFVF0XWJw99vpVG9Phai3S7cSSYY0Av6TJtXfUyeEwAGWmVpsYPKfOAeye/5xXzIJJhOYFv9Zk+r7qXWsCwdHosmGJv8Zk+r5qW0gVwafGYC0uUefdarTp1mcwBGXUDJAyd32aeI0lMRmL0i99au0Y1mS1LnHnw2q8ycJ43wN3vOgSrqTcA7vIFG5B74MqguohDM6BNEByKAXC5bco9EG1f1TYlkfoi2xHEwPUVXuIS+D6vjx4n1yyIGEkg1Khf+M2XPysl+KpsRxWL0q3HMKDKrTRx6IPkIwV/xtn/MOoIniqBh238125UEbxahST+CgjpgoJ2QgKQLjjxHEWjKRhihHZEx1yoSBO3xuB0VUS45HEW3JtJZFEff23A6KTX7T51HUnPvmhBpOQDbSf5NAGpMBaUAgdef0DwXVXhZU3OlzO6Bak4FqDYDqLAsq7vq5HVCdyUB1xJJU3Ad0M6BCLCYCFbbUD+rCc+p5x1NfF6iTGUeqNgCquSyouFvodkA1p1q3wJZ6QTUWnlPPO3n5ukC17KlAtex+UBeeU3su7cryOJWH5p3hx11xD5OZ57kaGI7bFV0ysdPLlbHGLnpliDPPcbjoNTZ1f6k3mXSrN4n4jNXNJS5XaKJQN0y60aeFo5d06EhDE5FOI/eX1i2kujFwrQ1S3e5SeibS4W6b4ozXLxnIqdOU9xEGkGWpRpmoWnx8qfj468uxwPX+3pYs/c8hh80cJ6FqylIMfPoxza/q09M0c42JzTWEzDfSVGPPNdVYKkHqn/J9Av/Gh+sEAd0WRUo/XBYEWlrM/aYJQV6WuVZCNY+5hqFnEqK0CmnWmO/GXJLzo5ah4tiiO9UMi1//ksJX22oc6o93HgSL/nEMMU6uUQYtZLJR9CYHhZsMnufsYLnior6zRblr39nSf2XL6Pthind/gDSAg1FIe1mICSZEUDuKcOuThwfHKZHtCrfah/XgHTLNVs/Bizka8RHEKtQddHpwEM6xmoUGmhvNahdOZl1Nf8vWbXO6OTP70jhtILdnKYaNLItZOW2hnDYY75uZjNPT38nV3K115PGPFsOHOH2i8Y8Oi8mcPp++xWNRIZiS0g1ThymtSUrPROnz8sMuntI4fYmq/hOcZjY9VKE4bVkoE/XRnO4GJBQL9W3NzenJXdjkKxAZTQ/lPNNjEk4TFPcnKK0yUrq5okUUSmNEHBtys9Ajw2xjWUo3S8/LtKYFVNMqs5oWK25gomp69AoRW2s6jNd7T8bp89Imx6ppKqMN9UTc69DTCiOnm5C9IJzWECaqK2uk6WFqaEuodMzNaXVyTp9MY+ss23iMgh/PvSYP7dKWciaqBY8X95xtIyhoSzOF/230QZrZH1zFvmB2o7HzhFdtPLD0/O3x/t1zrzKyZ6DJH7wDe/Z5abFyFTSg4WxWi9EUSsHpSjfvcPQayFx1GzKXXgLhQTJIEohshgcx4Qcgu0Ytg04yBuEohmW1zGxhnhvVMgarHSWUlkEdJJaxul+1/qnjdA66MrBUo6/ZuTUQHv4JIkiIwx7CvPZB5qVBUh7udK1GjoUYORrhrNFl1Q8evtjEaZ660Q1NCiYSjtIUzvkstnTATzsp2IyTgiKWsxJVF4Y5yaygG2izPGcFh8Hlk70F+9AtdzNlsDv5c532iLPt4aHecDKBXkBsRUKuoq4SwFf69EL9sP8CL3ejLfx9p7kaAfvLsfEBsqGJSW5h1EduDr4WKiCbBTM8nPLL46+wqe/F39VzvMnf4MhTJwuB01/Zdyv1s3ksy7gfWO3gYYVrwZY9tXkebLkfX60oJMuiMvSypDT6SlMv+ygy2Mt6aeC+lHZfVQs+tl0RY0HX7BgwBEsQ6TbeZrNRPa9FmhBscpIh4ZsvpmEuCKRCOL/6eARZn26eTkpx7/Bjue3gy59wmBMoWurqKYTmobrydnDGxICCPzzvotEd1nrzcBuSuohdhEnwdwmCinGLEao10ZprUL8SZ/H59CspWb0tgi2UzH8OcSVmrVel4P2rOXs1Ba9BfMjWUFBoInpayDE1/Yne+LEY/YgA5IcwvSE7OQTo2CbIBehFe5sNj26UV+l6u0Bdv1Xyz39Igig55PQ5d7mOVAeO8uxIeoia61QFGI/AD0+X9QrQHySCw7En2S5+K48bFqEzUPu/xBkQoCdR7Lm5t6P2RHDz+dNhfNJR5xrB3FLRRKPp7K2L3kjFec/UHZvf1BErcRlNStLQZRizpxRZ9+kL7y5xSBum8OU642x97bpGJyibhVflpN1AdLyGZu+rB4xwe8LCgM0WaZticlDFnh1Yo2qNVAgyO6BJEZo5dqMWuqtl4Q0ARy+kmOwVm7ys+daOWIlCutrNIjTUkeQ1kIv89IUz/ZUVafsKfa4svET3GSQqlpfSqZHE4XClwvU92Ayx0s1M1SZi+VoqnhOjLJoUozR5CVLbna/tWLMij1IpiLozkA1JcGXXSVUZqfyG8iKXnsdJx6oTVCHqqSbrQ4qOovq525Wo+pO10XtSq+eoUuZuSv3cZrBNWvvay+rn2TbMSP2MqQpB9DOmSK1eRcq8V0Y9q9nZ9fNQyB5zwwXRJr5ZXeRwz45RVtLpP1oXNRUZlJFYbv+FlJFt9tqgsysjUkxAnljaq6CQc2dJ99wvvJad7diwS9FPE59H06wTGZSWWA49VGnZ/UvRsUrL0fgqLba841Zi2n2xBGRfNba+mPZ/k/q8Q1rQrgqirsuOEnZG078BoKyS9lLfjBloWMhxYStSvjvJDpxPzzaW6e3q2U+oVNbTRY/SLYhK1RGnoeNMo1KRTHLH5qpSlfPiJ+UFUqWi2hwieirltSkldG0K16GExam16OJUUaVSGquUdGY7T7C9vFd05rGinHcNrGAEFvBsep35IG/Bz6Y3x6YjCEDqiw5fiEhq5qN4RCf16AQx/qSWN+NMTOpruUXEGrsMEoDUFx3cEZHU1rWQevydC9xJPdvdOLdKatb0dOFJfcELxYsOc4lIauZjwgUntX25C0WV5HhFSB3s3WLc2zyuix6D/RY+Nwxe4F8vDJJ1eTrLkxfvk0O1G+zpd5C/xPHfa0W13+F/90m0HXavPj19MxxnGq8pdtxN9bbFlYlSC+HbNC42Tp8QSt1k91vsg6LG/wE=</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">7V1pc6NIEv0t+8ERu46wgqri/NiWW+6eqw/HzLbnC4FQSaIbhBaQJc+v3+IUVBUISxySGk/EtCkhjnpZmS9fJvgGjZ3do2esl7+7M2zfQGG2u0EPNxBCAQDyTzjyGo8ApArxyMKzZsnYfuDJ+gcng+luG2uG/cKOgevagbUuDpruaoXNoDBmeJ67Le42d+3iWdfGAjMDT6Zhs6P/tWbBMh5VobIf/4CtxTI9M5C1+BPHSHdO7sRfGjN3mxtC72/Q2HPdIP7N2Y2xHc5eOi/x9yYln2YX5uFVUOcLnz8+i/d3Y9H8tPnF0b8EWAef7pKjvBj2Jrnh2cYPsKNbKytIrjt4TSeD3MI6/HXtuSb2yczeb5dWgJ/WhhkOb4kNkDHDthYrsmmSC8MeGXjBXmCRKX2XfOBYs1l4yHs/hlsYCUgVZSDt/y+ST+fuKnhKTg7CvQPP/ZGhgMI9LNseu7brRZeHJpMx+SHjM8sjtmC54cmw4ZPpuWdnK711cnF4lxtKZu8Ruw4OvFeyS/qpmCCZ2rKUbG/3hgHlZGyZMwqUjBmJLS6yQ+/hIr8kiL0BPcig5xg/sO7jmb5eGquFrwfG1CbLZwCyAKQgFoDU+sYRMTgyiOHV7F3ozkI4bMP3LZPMxTJw7GROPXezmuHwJALZwjsr+Jb7/TkCR0q2Hna5jx5e040VuZVv+Y34W1BKt/ffi7bSL5Yi4rsbz8SH7TcwvAUODnspPCv4ahbfHIASB790zMO2EVgvRQ/PAzU5w2fXIneWmY9EuQGZMov4tpMv5R0ydRxUPA6UpOJx4mlhjkPMwHjN7bYOd/DLL1elLhdIlVdF7y6plHnH598bezb/x9u/yNj/v2//lQSitb1ZWKv/3Ej3N9LDyF9jk1kbRdvnebLcQomcUeyugMS6nvl8Dk0z5/xsPA8yn5XbcSZPZUluxiFJ2uHIAhSORYtteSSpNY+UeaHnvBPieiQmKsjyPeJhMR5r2mTCRBZ4knMSz8rpiErRRESV8hZ1vY5EsxjKhBryOrJMu50CO23Hj8iM1SZOJKNFER8a6FAlr824T298SCnDMcxKdHKzhjNwWgpEKoRkeXZvIKplIJqus94E0XocQKwEUUkpYW8gamUg+jjQZ0ZgDAhWL0Ox72WYxve8SLD15lYweNADi0/rcPE9P4wnE+WP778/3E21J//vXwP1O0ecuw1F0ttRJXrmq20R4u2hw7nQNKbov02zAcP8sYiI+6dNQA6TYVOeMBFODiPu3QAEMihCoArs6uEhoLaFACuwMdN+XBoEGhBmDukyTaRP0RE/Y88i8xku9hMFnxTfvOJTZfpnknwhKvmSj02+6KxIleplX0ckSNxp7UZoPMqaEUIH7LlhW6wrPspnZYq0+EgHn7qWSAU7iOCoHfmRUjmBIFReFy2KooJaebJswF0VrPw4FMFqE1xOiG6NJd3NhT+B+On90vn1F+Hrx7+c939/qS5h6nNrF1bDLlIymM+xHGnRHChZPVrRpkKlR6wPclazTuNd7yCzRCyXh14mui0vVEhRDVnsMJ3hYsiSjwTDrUdg0Y3wCAOEOQgz35pC2GVGyoWwMlKG69Bazd0BxUpnKqK+nSlb3CugGBVcL86dnk2wFLtU3rn4smWwfTW9ZIlefgmd6ltACusrWyuhc1Fgi1jMrLevHTWmAR2hP5Wimk/Nq8JMPjWvMvQzSc0zVejUxiCJSvIhdZw3Z+YnyEgT7E3B9jfDm6zHYONMN0h9Kk+54tih73a7EQkTN3Dsm+46s/4LiiZtcwIKYQl2GDO4kJYmWAykUbYVGIsB1YOodtlI8eV/1rfn+8c5QF//fPkOrUfw1bvIChIz3xxUSiFQqJQJiJxWurZKSFwI2IXVEA04qrcXvC2EN0IfyIXwVXsqbnNXX54DcOdXYTlF+VLorXuP7qqFwpHcQKZbRRkC21wJiTuPrIqTPGow9cKmnrlLbkMn7kXfWsFSt3xvrkeP+lxdqDjJTwFAw4g61Ha4wLLazq3+fesHYZ+WrzsbH+tra4dtf+QbL9cYO2jaDQAnk+80drTXhf1TxI6DMYGTZ/JjR795JW2XjMuvGzxEmh2Ra+w0eJQ2aJOT6lvjBdt4tQiWFyc+th0tZFDEDXX56GEV7WKRpGjAgGQ1kmKXpVUukmyfNonvV/iMRLO4Sb3jxrZmPwr68+d3H/Tnv54e9Y9PXycRI8D+1Wb7Mt1cJPXN2Djd1g1RtuM6RkdAzHjac/5TDmkTR4iDV9Nto40QN7UmcdP65G0KbZ3o2GY9EdFRo2PeBlgh0cN+4HpDnDjE1EDfcQKwCuRtnM7/LIk9EyZg36Iw4IktRmi8vvHzCC6MWI84rTHd4tKN4lIavKGyl1WuU3JJvdGZay4S/fg645+Pjt0IdRy7WdFlsOlmbfpCdET6zQ9QVeiHPi7HqksfL3fIia0VIZj/4KHydPgBDm6FvFuGCtmom9QU8S7wDDNmRAOORfJE4wiQ1jeO3USa5I1kby1YVcQZHtE9S90jXSkHI404kKdGwgxsrYs3b9Kgnj0LI1lW6tp0B0+Ad2vScebQG3mi3a2gHU2eZEivD6Fjq2bJU1NWPQKEmeUseyQotbVqVTobrbrljKG+0ffqyBmjB3Qrem2TVwmtgpTV0/3xbVs9WzFLaGb0erio18lyjMU1vi/3JKqpqcUo3OnjwPwmRLbMlsfvWgVTDRaBUCDbNNipXprSszPpUMupSVdDj7SakSJeEr01OMvEl+V+iq5ePTpsUEypxffl8CcVVoaMqEZDrIgg7w+hg9/nTImOSv+xg+1gLwIZvtYwjiDjqwwi1MusVM4rCrsNImwxNAbkuoM5EKiytNZ39ROxOmwMxB8frhYE2j9pHb5ykA9Cz+W6aydUaUg/TKhQn4RK06ici+74qv2iN6ojSelYQUWNK6hNVYTZNOOAVNWLFfYqAAFAxShFOZLJA4GmHTX/8EBjdti45tnsuzXe5mlLRcu3lRZ6sWipV4uGQmMWrcKRBvf/qUUi0fVzudXqZtZcOKSpPBoIae/U5Uuc+BVcjsSZgbhPUa+RkjOltr7bdcXGZc6LiB2HdNcmo0ddmVPsVeYEgHoWVzn2Udz9w+VpvOiYD4mstHmiTedrwIU/8kViAKo0pLZ6PMufgwIVa+QuKS6fTedDv9J+Y5moQlm8Sr8DpW2Lb+/t9yOgoqLFi+Jg8SdYfK9Zrwq1kZYrZlE9aEAkH2vHrQENaSNZ2/8UcwaSE49kiTlxV+uD1eEbWx8kHcqtj7swJPS8PgqM6eBTCGezMvrlPyqqXhnSCStDrFwZsJ2VQTb3f4E93n3/h+zR+/8D</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>
 \ No newline at end of file
+<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>
 \ No newline at end of file
@@ -8,8 +8,8 @@ PRO check_phangs_ssps_isrf_prediction,show_map=show_map,from_restore=from_restor
  
 from_restore:
  
-;make_phangs_isrf_classes,/save
-make_phangs_isrf_classes
+;make_phangs_isrf_classes,bidon,/save
+make_phangs_isrf_classes,bidon
  
 stop
  
-PRO make_phangs_isrf_classes,save=save,help=help
+PRO make_phangs_isrf_classes,bidon,save=save,help=help
  
 ;+
 ; NAME:
-PRO phangs_brute_force_fit_with_isrf_grid,source_name=source_name,normalize=normalize,from_restore=from_restore
+PRO phangs_brute_force_fit_with_isrf_grid,source_name=source_name,normalize=normalize,from_restore=from_restore,save=save
  
-;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628'
+;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/normalize,/save
+;phangs_brute_force_fit_with_isrf_grid,source_name='ngc0628',/save
  
 use_source_name='ngc0628'
 IF keyword_set(source_name) THEN use_source_name=source_name
@@ -18,6 +19,7 @@ IF keyword_set(from_restore) THEN GOTO,from_restore
  
 dustem_init,show_plots=show_plots
  
+;==== needed for NHCO
 restore,data_dir+use_source_name+'_jwst_images.sav',/verb
 ;% RESTORE: Restored variable: JWST_IMAGES.
 ;% RESTORE: Restored variable: FILTERS.
@@ -31,9 +33,9 @@ restore,data_dir+use_source_name+&#39;_muse_images.sav&#39;,/verb
 ;% RESTORE: Restored variable: METALICITY_VALUES.
 ;% RESTORE: Restored variable: BINS.
 ;% RESTORE: Restored variable: HREF.
-restore,data_dir+use_source_name+'_muse_NH.sav',/verb
+;restore,data_dir+use_source_name+'_muse_NH.sav',/verb
 ;% RESTORE: Restored variable: NH_MAP.
-restore,data_dir+use_source_name+'_astrosat_data.sav',/verb
+restore,data_dir+use_source_name+'_astrosat_data.sav',/verb   ;should really be called _astrosat_images.sav
 ;% RESTORE: Restored variable: ASTROSAT_IMAGES.
 ;% RESTORE: Restored variable: ASTROSAT_FILTERS.
 ;% RESTORE: Restored variable: HREF.
@@ -41,11 +43,11 @@ restore,data_dir+use_source_name+&#39;_astrosat_data.sav&#39;,/verb
 file=data_dir+use_source_name+'_isrf_classes_one_ratio.sav'
 restore,file,/verb
 ;% RESTORE: Restored variable: CLASSES.
-;% RESTORE: Restored variable: VOR_CLASS.
+;% RESTORE: Restored variable: VOR_CLASS.      ;This is used
 ;% RESTORE: Restored variable: MAP_CLASSES.
-ind=where(vor_class NE 0)
-class_min=min(vor_class[ind])
-class_max=max(vor_class[ind])
+;ind=where(vor_class NE 0)
+;class_min=min(vor_class[ind])
+;class_max=max(vor_class[ind])
 restore,data_dir+use_source_name+'_isrf_min_prediction.sav',/verb
 ;% RESTORE: Restored variable: ISRFS.
 ;% RESTORE: Restored variable: OBJECT_DISTANCE.
@@ -58,8 +60,9 @@ restore,data_dir+use_source_name+&#39;_isrf_min_prediction.sav&#39;,/verb
 ;class_min=13
 ;class_max=13
  
-class_min=8
-class_max=17
+;class_min=8 & class_max=17
+;class_min=15 & class_max=15
+class_min=8 & class_max=15
  
 restore,data_dir+use_source_name+'_jwst_seds_muse_pixels.sav',/verb
 ;% RESTORE: Restored variable: ALL_SEDS.
@@ -102,11 +105,8 @@ FOR ii=0L,Nvor-1 DO BEGIN
 ENDFOR
  
 ;==== Brute force fit using ISRF classes
-grid_brute_force:
  
-;stop
-
-Nvor=max(voronoi_id)
+;Nvor=max(voronoi_id)
 G0s=fltarr(Nvor)+la_undef()
 Ypahs=fltarr(Nvor)+la_undef()
 YVSGs=fltarr(Nvor)+la_undef()
@@ -247,10 +247,10 @@ ENDFOR
 ;save,GOs,Ypahs,Yvsgs,file=dir+source_name+'DBP90_JWST_G0_YPAH_YVSG.fits'
 IF keyword_set(save) THEN BEGIN
 	file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_on-voronoi_with-ISRFclasses'+norm_str+'.sav'
-	save,G0s,Ypahs,Yvsgs,facts,chi2s,rchi2s,dG0s,dYpahs,dYvsgs,file=file_save
+	save,G0s,Ypahs,Yvsgs,facts,chi2s,rchi2s,dG0s,dYpahs,dYvsgs,G0_hit,Ypah_hit,Yvsg_hit,file=file_save
 	message,'Wrote '+file_save,/continue
 	file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+'.sav'
-	save,G0_map,Ypah_map,Yvsg_map,fact_map,chi2_map,rchi2_map,dG0_map,dYpah_map,dYvsg_map,file=file_save
+	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
 	message,'Wrote '+file_save,/continue
 ENDIF
  
@@ -261,19 +261,36 @@ restore,file_save,/verb
 file_save=data_dir+use_source_name+'_DBP90_JWST_G0_YPAH_YVSG_maps_with-ISRFclasses'+norm_str+'.sav'
 restore,file_save,/verb
  
+stop
  
-;G0 histogram
+;===== histograms
 win=0L
-window,win,xsize=800,ysize=800 & win=win+1
-
+window,win,xsize=700,ysize=900 & win=win+1
+!p.multi=[0,1,3]
 ind=where(G0s NE la_undef())
 res=histogram(G0s[ind],locations=xv)
 cgplot,xv,res,psym=10,title='G0s histogram',xtit='G0',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
  
 ind=where(Ypahs NE la_undef())
-res=histogram(Ypahs[ind],locations=xv)
-cgplot,xv,res,psym=10,title='Ypahs histogram',xtit='Ypah',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
+res=histogram(alog10(Ypahs[ind]),locations=xv,Nbins=100)
+cgplot,xv,res,psym=10,title='Ypahs histogram',xtit='log(Ypah)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
  
+ind=where(Yvsgs NE la_undef())
+res=histogram(alog10(Yvsgs[ind]),locations=xv,Nbins=100)
+cgplot,xv,res,psym=10,title='Yvsgs histogram',xtit='log(Yvsg)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
+
+Ypahs_f=la_div(Ypahs,facts)
+ind=where(Ypahs_f NE la_undef())
+res=histogram(alog10(Ypahs_f[ind]),locations=xv,Nbins=100)
+cgplot,xv,res,psym=10,title='Ypahs histogram',xtit='log(Ypah_f)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
+
+Yvsgs_f=la_div(Yvsgs,facts)
+ind=where(Yvsgs_f NE la_undef())
+res=histogram(alog10(Yvsgs_f[ind]),locations=xv,Nbins=100)
+cgplot,xv,res,psym=10,title='Yvsgs histogram',xtit='log(Yvsg_f)',ytit='Number',/ylog,yrange=[1,max(res)],/ysty
+
+;==== show maps
+win=0L
 window,win,xsize=800,ysize=900 & win=win+1
  
 obp=[1.1,0.,1.15,1]
@@ -281,22 +298,21 @@ obp=[1.1,0.,1.15,1]
 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)'
  
 window,win,xsize=800,ysize=900 & win=win+1
-imrange=[-3.,2]
+imrange=[-4.,-3.]   ;normalize=1
+imrange=[-4.,-3.5]+2   ;normalize=0
 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)'
  
 window,win,xsize=800,ysize=900 & win=win+1
-imrange=[-2.,6]
-;imrange=[-3.,0.]-3.  ;why ??
-;image_cont20,la_log10(YVSG_map),Href,/square,imrange=[-2.,6],axis_color_table=1,image_color_table='jpbloadct',/silent,off_bar=obp,title='log10(Yvsg)'
+imrange=[-4.,-3.]+0.5   ;normalize=1
+imrange=[-2.2,-1.6]   ;normalize=0
 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)'
  
 window,win,xsize=800,ysize=900 & win=win+1
-imrange=[3,7]
-;imrange=[3,7]+6.
+imrange=[2,7]+1
 image_cont20,la_log10(chi2_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='log10(Chi2)',off_bar=obp
  
 window,win,xsize=800,ysize=900 & win=win+1
-imrange=[-3,2]
+imrange=[-3,2]   ;normalize=0
 image_cont20,la_log10(fact_map),Href,/square,imrange=imrange,axis_color_table=1,image_color_table='jpbloadct',/silent,title='log10(fact)',off_bar=obp
  
 window,win,xsize=800,ysize=900 & win=win+1
@@ -318,6 +334,7 @@ image_cont20,la_log10(la_mul(Ypah_map,fact_map)),Href,/square,imrange=imrange,ax
  
  
 stop
+
 jwst_4_coutours=fltarr((size(jwst_images))[1],(size(jwst_images))[2],2)
 jwst_4_coutours[*,*,0]=jwst_images[*,*,0,5]
 jwst_4_coutours[*,*,1]=jwst_images[*,*,0,5]
@@ -61,7 +61,7 @@ astrosat_data:
 ;in fact Angstrom-1 cm-2 erg s-1 PER PIXEL ?!
  
 ;This is to get the reference header
-restore,save_data_dir+'ngc0628_jwst_images.sav',/verb
+restore,save_data_dir+use_source_name+'_jwst_images.sav',/verb
 ;% RESTORE: Restored variable: JWST_IMAGES.
 ;% RESTORE: Restored variable: FILTERS.
 ;% RESTORE: Restored variable: HREF.
@@ -21,7 +21,8 @@ bidon=1
  
 ;make grid 
 ;make_sed_phangs_tables,/test,grid_type=2,/show_seds,isrf_class=15
-make_sed_phangs_tables,grid_type=2,isrf_class=18
+make_sed_phangs_tables,grid_type=2,isrf_class=19
+;make_sed_phangs_tables,grid_type=2,isrf_class=18
 ;make_sed_phangs_tables,grid_type=2,isrf_class=17
 ;make_sed_phangs_tables,grid_type=2,isrf_class=16
 ;make_sed_phangs_tables,grid_type=2,isrf_class=15