Create image example

Example for generating a rest-frame physical scale image. This example will: - Build a parametric galaxy (see make_sfzh and make_sed) - Define its morphology - Calculate rest-frame luminosities for the UVJ bands - Make an image of the galaxy, including an RGB image.

llllllllll!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!llllllllll
ll!!!!!!!!!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!!!!!!!!!ll
!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!
!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!
!!iiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiii!!
iiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiii
iiiiii>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>iiiiii
>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>
>>>>>>>>>><<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<>>>>>>>>>>
>>>><<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++____________________________________________________________________++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<>>>>
<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~++++++++++++++++++++__________________________________--------------------------------__________________________________++++++++++++++++++++~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<
<<<<<<~~~~~~~~~~~~~~~~++++++++++++++++________________________----------------------------------------------------------------------------________________________++++++++++++++++~~~~~~~~~~~~~~~~<<<<<<
~~~~~~~~~~~~~~++++++++++++++++__________________------------------------????????????????????????????????????????????????????????------------------------__________________++++++++++++++++~~~~~~~~~~~~~~
~~~~~~~~++++++++++++++________________------------------??????????????????????????????????]]]]]]]]]]]]]]]]]]]]??????????????????????????????????------------------________________++++++++++++++~~~~~~~~
~~++++++++++++______________----------------????????????????????]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]????????????????????----------------______________++++++++++++~~
++++++++______________--------------????????????????]]]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]]]]]]]]]????????????????--------------______________++++++++
++____________--------------??????????????]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]]]??????????????--------------____________++
________------------????????????]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]????????????------------________
____----------????????????]]]]]]]]]]]][[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[]]]]]]]]]]]]????????????----------____
--------????????????]]]]]]]]]][[[[[[[[[[[[}}}}}}}}}}}}}}{{{{{{{{{{{{{{{{{{1111111111111111111111111111111111111111111111111111{{{{{{{{{{{{{{{{{{}}}}}}}}}}}}}}[[[[[[[[[[[[]]]]]]]]]]????????????--------
----??????????]]]]]]]]]][[[[[[[[[[}}}}}}}}}}}}{{{{{{{{{{{{{{11111111111111111111))))))))))))))))))))))))))))))))))))))))11111111111111111111{{{{{{{{{{{{{{}}}}}}}}}}}}[[[[[[[[[[]]]]]]]]]]??????????----
????????]]]]]]]]]][[[[[[[[[[}}}}}}}}}}{{{{{{{{{{{{11111111111111))))))))))))))))))))(((((((((((((((((((((((((((((((())))))))))))))))))))11111111111111{{{{{{{{{{{{}}}}}}}}}}[[[[[[[[[[]]]]]]]]]]????????
????]]]]]]]][[[[[[[[[[}}}}}}}}}}{{{{{{{{{{1111111111))))))))))))))((((((((((((((((((((((((||||||||||||||||||||(((((((((((((((((((((((())))))))))))))1111111111{{{{{{{{{{}}}}}}}}}}[[[[[[[[[[]]]]]]]]????
]]]]]]]][[[[[[[[}}}}}}}}}}{{{{{{{{1111111111))))))))))))((((((((((((||||||||||||||||||||||||||tttttttttttt||||||||||||||||||||||||||(((((((((((())))))))))))1111111111{{{{{{{{}}}}}}}}}}[[[[[[[[]]]]]]]]
]]]][[[[[[[[}}}}}}}}{{{{{{{{1111111111))))))))((((((((((((||||||||||||ttttttttttttttttttttttttttttfffftttttttttttttttttttttttttttt||||||||||||(((((((((((())))))))1111111111{{{{{{{{}}}}}}}}[[[[[[[[]]]]
[[[[[[[[}}}}}}}}{{{{{{{{11111111))))))))((((((((||||||||||||ttttttttttttfffffffffffffffffffffffffffffffffffffffffffffffffffffffftttttttttttt||||||||||||(((((((())))))))11111111{{{{{{{{}}}}}}}}[[[[[[[[
[[[[}}}}}}{{{{{{{{11111111))))))))((((((((||||||||ttttttttttffffffffffffjjjjjjjjjjjjjjjjjjjjjjjjrrrrrrrrjjjjjjjjjjjjjjjjjjjjjjjjfffffffffffftttttttttt||||||||(((((((())))))))11111111{{{{{{{{}}}}}}[[[[
}}}}}}{{{{{{{{11111111))))))((((((((||||||||ttttttttffffffffffjjjjjjjjjjrrrrrrrrrrrrrrrrrrrrxxxxxxxxxxxxxxxxrrrrrrrrrrrrrrrrrrrrjjjjjjjjjjfffffffffftttttttt||||||||(((((((())))))11111111{{{{{{{{}}}}}}
}}{{{{{{{{111111))))))))((((((||||||||ttttttttffffffffjjjjjjjjrrrrrrrrrrxxxxxxxxxxxxxxxxxxnnnnnnnnnnnnnnnnnnnnxxxxxxxxxxxxxxxxxxrrrrrrrrrrjjjjjjjjfffffffftttttttt||||||||(((((())))))))111111{{{{{{{{}}
{{{{{{111111))))))))((((((||||||ttttttttffffffjjjjjjjjrrrrrrrrxxxxxxxxxxnnnnnnnnnnnnnnuuuuuuuuuuuuuuuuuuuuuuuuuuuunnnnnnnnnnnnnnxxxxxxxxxxrrrrrrrrjjjjjjjjfffffftttttttt||||||(((((())))))))111111{{{{{{
{{11111111))))))((((((||||||ttttttffffffjjjjjjjjrrrrrrxxxxxxxxnnnnnnnnnnuuuuuuuuuuuuvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvuuuuuuuuuuuunnnnnnnnnnxxxxxxxxrrrrrrjjjjjjjjfffffftttttt||||||(((((())))))11111111{{
111111))))))((((((||||||ttttttffffffjjjjjjrrrrrrxxxxxxnnnnnnnnuuuuuuuuvvvvvvvvvvvvccccccccccccccccccccccccccccccccccccvvvvvvvvvvvvuuuuuuuunnnnnnnnxxxxxxrrrrrrjjjjjjfffffftttttt||||||(((((())))))111111
11))))))((((((||||||ttttttffffffjjjjjjrrrrxxxxxxnnnnnnnnuuuuuuvvvvvvvvcccccccccczzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzccccccccccvvvvvvvvuuuuuunnnnnnnnxxxxxxrrrrjjjjjjfffffftttttt||||||(((((())))))11
))))))((((((||||ttttttffffffjjjjrrrrrrxxxxxxnnnnnnuuuuuuvvvvvvcccccczzzzzzzzzzXXXXXXXXXXXXYYYYYYYYYYYYYYYYYYYYXXXXXXXXXXXXzzzzzzzzzzccccccvvvvvvuuuuuunnnnnnxxxxxxrrrrrrjjjjfffffftttttt||||(((((())))))
))))((((||||||ttttffffffjjjjjjrrrrxxxxxxnnnnuuuuuuvvvvvvcccccczzzzzzXXXXXXXXYYYYYYYYYYUUUUUUUUUUUUUUUUUUUUUUUUUUUUYYYYYYYYYYXXXXXXXXzzzzzzccccccvvvvvvuuuuuunnnnxxxxxxrrrrjjjjjjfffffftttt||||||(((())))
((((((||||||ttttffffffjjjjrrrrxxxxxxnnnnuuuuuuvvvvcccccczzzzzzXXXXXXYYYYYYUUUUUUUUJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJUUUUUUUUYYYYYYXXXXXXzzzzzzccccccvvvvuuuuuunnnnxxxxxxrrrrjjjjfffffftttt||||||((((((
((((||||ttttttffffjjjjrrrrrrxxxxnnnnuuuuuuvvvvcccczzzzzzXXXXYYYYYYUUUUUUJJJJJJJJCCCCCCCCCCLLLLLLLLLLLLLLLLLLLLCCCCCCCCCCJJJJJJJJUUUUUUYYYYYYXXXXzzzzzzccccvvvvuuuuuunnnnxxxxrrrrrrjjjjfffftttttt||||((((
((||||ttttttffffjjjjrrrrxxxxxxnnnnuuuuvvvvcccczzzzzzXXXXYYYYUUUUUUJJJJCCCCCCLLLLLLLLQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQLLLLLLLLCCCCCCJJJJUUUUUUYYYYXXXXzzzzzzccccvvvvuuuunnnnxxxxxxrrrrjjjjfffftttttt||||((
||||ttttffffffjjjjrrrrxxxxnnnnuuuuuuvvvvcccczzzzXXXXYYYYUUUUJJJJCCCCCCLLLLLLQQQQQQ00000000OOOOOOOOOOOOOOOOOOOO00000000QQQQQQLLLLLLCCCCCCJJJJUUUUYYYYXXXXzzzzccccvvvvuuuuuunnnnxxxxrrrrjjjjfffffftttt||||
||ttttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXXXYYYYUUUUJJJJCCCCLLLLQQQQQQ0000OOOOOOOOZZZZZZZZZZZZZZZZZZZZZZZZZZZZOOOOOOOO0000QQQQQQLLLLCCCCJJJJUUUUYYYYXXXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttttt||
ttttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXXXYYUUUUJJJJCCCCLLLLQQQQ0000OOOOZZZZZZmmmmmmmmwwwwwwwwwwwwwwwwwwwwmmmmmmmmZZZZZZOOOO0000QQQQLLLLCCCCJJJJUUUUYYXXXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttttt
ttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXYYYYUUUUJJJJCCLLLLQQQQ0000OOOOZZZZmmmmwwwwwwqqqqqqppppppppppppppppqqqqqqwwwwwwmmmmZZZZOOOO0000QQQQLLLLCCJJJJUUUUYYYYXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttt
ttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzXXXXYYYYUUUUJJCCCCLLLLQQ0000OOOOZZmmmmwwwwqqqqppppppddddddddbbbbbbbbddddddddppppppqqqqwwwwmmmmZZOOOO0000QQLLLLCCCCJJUUUUYYYYXXXXzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftt
ttffffjjjjrrrrxxxxnnnnuuvvvvcccczzzzXXYYYYUUUUJJCCCCLLQQQQ00OOOOZZZZmmwwwwqqppppddddbbbbbbkkkkkkkkkkkkkkkkkkkkbbbbbbddddppppqqwwwwmmZZZZOOOO00QQQQLLCCCCJJUUUUYYYYXXzzzzccccvvvvuunnnnxxxxrrrrjjjjfffftt
ffffjjjjrrrrxxxxnnnnuuuuvvcccczzzzXXXXYYUUUUJJCCCCLLQQQQ00OOOOZZmmmmwwqqqqppddddbbkkkkhhhhhhaaaaaaaaaaaaaaaahhhhhhkkkkbbddddppqqqqwwmmmmZZOOOO00QQQQLLCCCCJJUUUUYYXXXXzzzzccccvvuuuunnnnxxxxrrrrjjjjffff
ffffjjjjrrrrxxxxnnuuuuvvvvcccczzXXXXYYUUUUJJCCCCLLQQQQ00OOOOZZmmwwwwqqppppddbbkkkkhhaaaaoooo****************ooooaaaahhkkkkbbddppppqqwwwwmmZZOOOO00QQQQLLCCCCJJUUUUYYXXXXzzccccvvvvuuuunnxxxxrrrrjjjjffff
ffffjjrrrrxxxxnnnnuuuuvvvvcczzzzXXYYYYUUJJJJCCLLLLQQ0000OOZZmmmmwwqqppddddbbkkhhaaaaoo****####MMMMMMMMMMMM####****ooaaaahhkkbbddddppqqwwmmmmZZOO0000QQLLLLCCJJJJUUYYYYXXzzzzccvvvvuuuunnnnxxxxrrrrjjffff
ffjjjjrrrrxxxxnnnnuuuuvvcccczzXXXXYYYYUUJJCCCCLLQQQQ00OOZZZZmmwwqqppppddbbkkhhaaoooo**##MMMMWW&&&&&&&&&&&&WWMMMM##**ooooaahhkkbbddppppqqwwmmZZZZOO00QQQQLLCCCCJJUUYYYYXXXXzzccccvvuuuunnnnxxxxrrrrjjjjff
ffjjjjrrrrxxxxnnnnuuvvvvcccczzXXXXYYUUUUJJCCCCLLQQ0000OOZZmmwwwwqqppddbbkkhhhhaaoo**##MMWW&&8888%%%%%%%%8888&&WWMM##**ooaahhhhkkbbddppqqwwwwmmZZOO0000QQLLCCCCJJUUUUYYXXXXzzccccvvvvuunnnnxxxxrrrrjjjjff
ffjjjjrrrrxxxxnnuuuuvvvvcccczzXXXXYYUUUUJJCCLLLLQQ00OOOOZZmmwwqqqqppddbbkkhhaaoo**##MMWW&&88%%BB@@$$$$@@BB%%88&&WWMM##**ooaahhkkbbddppqqqqwwmmZZOOOO00QQLLLLCCJJUUUUYYXXXXzzccccvvvvuuuunnxxxxrrrrjjjjff
ffjjjjrrrrxxxxnnuuuuvvvvcccczzXXXXYYUUUUJJCCLLLLQQ00OOOOZZmmwwqqqqppddbbkkhhaaoo**##MMWW&&88%%BB@@$$$$@@BB%%88&&WWMM##**ooaahhkkbbddppqqqqwwmmZZOOOO00QQLLLLCCJJUUUUYYXXXXzzccccvvvvuuuunnxxxxrrrrjjjjff
ffjjjjrrrrxxxxnnnnuuvvvvcccczzXXXXYYUUUUJJCCCCLLQQ0000OOZZmmwwwwqqppddbbkkhhhhaaoo**##MMWW&&8888%%%%%%%%8888&&WWMM##**ooaahhhhkkbbddppqqwwwwmmZZOO0000QQLLCCCCJJUUUUYYXXXXzzccccvvvvuunnnnxxxxrrrrjjjjff
ffjjjjrrrrxxxxnnnnuuuuvvcccczzXXXXYYYYUUJJCCCCLLQQQQ00OOZZZZmmwwqqppppddbbkkhhaaoooo**##MMMMWW&&&&&&&&&&&&WWMMMM##**ooooaahhkkbbddppppqqwwmmZZZZOO00QQQQLLCCCCJJUUYYYYXXXXzzccccvvuuuunnnnxxxxrrrrjjjjff
ffffjjrrrrxxxxnnnnuuuuvvvvcczzzzXXYYYYUUJJJJCCLLLLQQ0000OOZZmmmmwwqqppddddbbkkhhaaaaoo****####MMMMMMMMMMMM####****ooaaaahhkkbbddddppqqwwmmmmZZOO0000QQLLLLCCJJJJUUYYYYXXzzzzccvvvvuuuunnnnxxxxrrrrjjffff
ffffjjjjrrrrxxxxnnuuuuvvvvcccczzXXXXYYUUUUJJCCCCLLQQQQ00OOOOZZmmwwwwqqppppddbbkkkkhhaaaaoooo****************ooooaaaahhkkkkbbddppppqqwwwwmmZZOOOO00QQQQLLCCCCJJUUUUYYXXXXzzccccvvvvuuuunnxxxxrrrrjjjjffff
ffffjjjjrrrrxxxxnnnnuuuuvvcccczzzzXXXXYYUUUUJJCCCCLLQQQQ00OOOOZZmmmmwwqqqqppddddbbkkkkhhhhhhaaaaaaaaaaaaaaaahhhhhhkkkkbbddddppqqqqwwmmmmZZOOOO00QQQQLLCCCCJJUUUUYYXXXXzzzzccccvvuuuunnnnxxxxrrrrjjjjffff
ttffffjjjjrrrrxxxxnnnnuuvvvvcccczzzzXXYYYYUUUUJJCCCCLLQQQQ00OOOOZZZZmmwwwwqqppppddddbbbbbbkkkkkkkkkkkkkkkkkkkkbbbbbbddddppppqqwwwwmmZZZZOOOO00QQQQLLCCCCJJUUUUYYYYXXzzzzccccvvvvuunnnnxxxxrrrrjjjjfffftt
ttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzXXXXYYYYUUUUJJCCCCLLLLQQ0000OOOOZZmmmmwwwwqqqqppppppddddddddbbbbbbbbddddddddppppppqqqqwwwwmmmmZZOOOO0000QQLLLLCCCCJJUUUUYYYYXXXXzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftt
ttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXYYYYUUUUJJJJCCLLLLQQQQ0000OOOOZZZZmmmmwwwwwwqqqqqqppppppppppppppppqqqqqqwwwwwwmmmmZZZZOOOO0000QQQQLLLLCCJJJJUUUUYYYYXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttt
ttttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXXXYYUUUUJJJJCCCCLLLLQQQQ0000OOOOZZZZZZmmmmmmmmwwwwwwwwwwwwwwwwwwwwmmmmmmmmZZZZZZOOOO0000QQQQLLLLCCCCJJJJUUUUYYXXXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttttt
||ttttttffffjjjjrrrrxxxxnnnnuuuuvvvvcccczzzzXXXXYYYYUUUUJJJJCCCCLLLLQQQQQQ0000OOOOOOOOZZZZZZZZZZZZZZZZZZZZZZZZZZZZOOOOOOOO0000QQQQQQLLLLCCCCJJJJUUUUYYYYXXXXzzzzccccvvvvuuuunnnnxxxxrrrrjjjjfffftttttt||
||||ttttffffffjjjjrrrrxxxxnnnnuuuuuuvvvvcccczzzzXXXXYYYYUUUUJJJJCCCCCCLLLLLLQQQQQQ00000000OOOOOOOOOOOOOOOOOOOO00000000QQQQQQLLLLLLCCCCCCJJJJUUUUYYYYXXXXzzzzccccvvvvuuuuuunnnnxxxxrrrrjjjjfffffftttt||||
((||||ttttttffffjjjjrrrrxxxxxxnnnnuuuuvvvvcccczzzzzzXXXXYYYYUUUUUUJJJJCCCCCCLLLLLLLLQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQQLLLLLLLLCCCCCCJJJJUUUUUUYYYYXXXXzzzzzzccccvvvvuuuunnnnxxxxxxrrrrjjjjfffftttttt||||((
((((||||ttttttffffjjjjrrrrrrxxxxnnnnuuuuuuvvvvcccczzzzzzXXXXYYYYYYUUUUUUJJJJJJJJCCCCCCCCCCLLLLLLLLLLLLLLLLLLLLCCCCCCCCCCJJJJJJJJUUUUUUYYYYYYXXXXzzzzzzccccvvvvuuuuuunnnnxxxxrrrrrrjjjjfffftttttt||||((((
((((((||||||ttttffffffjjjjrrrrxxxxxxnnnnuuuuuuvvvvcccccczzzzzzXXXXXXYYYYYYUUUUUUUUJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJUUUUUUUUYYYYYYXXXXXXzzzzzzccccccvvvvuuuuuunnnnxxxxxxrrrrjjjjfffffftttt||||||((((((
))))((((||||||ttttffffffjjjjjjrrrrxxxxxxnnnnuuuuuuvvvvvvcccccczzzzzzXXXXXXXXYYYYYYYYYYUUUUUUUUUUUUUUUUUUUUUUUUUUUUYYYYYYYYYYXXXXXXXXzzzzzzccccccvvvvvvuuuuuunnnnxxxxxxrrrrjjjjjjfffffftttt||||||(((())))
))))))((((((||||ttttttffffffjjjjrrrrrrxxxxxxnnnnnnuuuuuuvvvvvvcccccczzzzzzzzzzXXXXXXXXXXXXYYYYYYYYYYYYYYYYYYYYXXXXXXXXXXXXzzzzzzzzzzccccccvvvvvvuuuuuunnnnnnxxxxxxrrrrrrjjjjfffffftttttt||||(((((())))))
11))))))((((((||||||ttttttffffffjjjjjjrrrrxxxxxxnnnnnnnnuuuuuuvvvvvvvvcccccccccczzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzccccccccccvvvvvvvvuuuuuunnnnnnnnxxxxxxrrrrjjjjjjfffffftttttt||||||(((((())))))11
111111))))))((((((||||||ttttttffffffjjjjjjrrrrrrxxxxxxnnnnnnnnuuuuuuuuvvvvvvvvvvvvccccccccccccccccccccccccccccccccccccvvvvvvvvvvvvuuuuuuuunnnnnnnnxxxxxxrrrrrrjjjjjjfffffftttttt||||||(((((())))))111111
{{11111111))))))((((((||||||ttttttffffffjjjjjjjjrrrrrrxxxxxxxxnnnnnnnnnnuuuuuuuuuuuuvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvuuuuuuuuuuuunnnnnnnnnnxxxxxxxxrrrrrrjjjjjjjjfffffftttttt||||||(((((())))))11111111{{
{{{{{{111111))))))))((((((||||||ttttttttffffffjjjjjjjjrrrrrrrrxxxxxxxxxxnnnnnnnnnnnnnnuuuuuuuuuuuuuuuuuuuuuuuuuuuunnnnnnnnnnnnnnxxxxxxxxxxrrrrrrrrjjjjjjjjfffffftttttttt||||||(((((())))))))111111{{{{{{
}}{{{{{{{{111111))))))))((((((||||||||ttttttttffffffffjjjjjjjjrrrrrrrrrrxxxxxxxxxxxxxxxxxxnnnnnnnnnnnnnnnnnnnnxxxxxxxxxxxxxxxxxxrrrrrrrrrrjjjjjjjjfffffffftttttttt||||||||(((((())))))))111111{{{{{{{{}}
}}}}}}{{{{{{{{11111111))))))((((((((||||||||ttttttttffffffffffjjjjjjjjjjrrrrrrrrrrrrrrrrrrrrxxxxxxxxxxxxxxxxrrrrrrrrrrrrrrrrrrrrjjjjjjjjjjfffffffffftttttttt||||||||(((((((())))))11111111{{{{{{{{}}}}}}
[[[[}}}}}}{{{{{{{{11111111))))))))((((((((||||||||ttttttttttffffffffffffjjjjjjjjjjjjjjjjjjjjjjjjrrrrrrrrjjjjjjjjjjjjjjjjjjjjjjjjfffffffffffftttttttttt||||||||(((((((())))))))11111111{{{{{{{{}}}}}}[[[[
[[[[[[[[}}}}}}}}{{{{{{{{11111111))))))))((((((((||||||||||||ttttttttttttfffffffffffffffffffffffffffffffffffffffffffffffffffffffftttttttttttt||||||||||||(((((((())))))))11111111{{{{{{{{}}}}}}}}[[[[[[[[
]]]][[[[[[[[}}}}}}}}{{{{{{{{1111111111))))))))((((((((((((||||||||||||ttttttttttttttttttttttttttttfffftttttttttttttttttttttttttttt||||||||||||(((((((((((())))))))1111111111{{{{{{{{}}}}}}}}[[[[[[[[]]]]
]]]]]]]][[[[[[[[}}}}}}}}}}{{{{{{{{1111111111))))))))))))((((((((((((||||||||||||||||||||||||||tttttttttttt||||||||||||||||||||||||||(((((((((((())))))))))))1111111111{{{{{{{{}}}}}}}}}}[[[[[[[[]]]]]]]]
????]]]]]]]][[[[[[[[[[}}}}}}}}}}{{{{{{{{{{1111111111))))))))))))))((((((((((((((((((((((((||||||||||||||||||||(((((((((((((((((((((((())))))))))))))1111111111{{{{{{{{{{}}}}}}}}}}[[[[[[[[[[]]]]]]]]????
????????]]]]]]]]]][[[[[[[[[[}}}}}}}}}}{{{{{{{{{{{{11111111111111))))))))))))))))))))(((((((((((((((((((((((((((((((())))))))))))))))))))11111111111111{{{{{{{{{{{{}}}}}}}}}}[[[[[[[[[[]]]]]]]]]]????????
----??????????]]]]]]]]]][[[[[[[[[[}}}}}}}}}}}}{{{{{{{{{{{{{{11111111111111111111))))))))))))))))))))))))))))))))))))))))11111111111111111111{{{{{{{{{{{{{{}}}}}}}}}}}}[[[[[[[[[[]]]]]]]]]]??????????----
--------????????????]]]]]]]]]][[[[[[[[[[[[}}}}}}}}}}}}}}{{{{{{{{{{{{{{{{{{1111111111111111111111111111111111111111111111111111{{{{{{{{{{{{{{{{{{}}}}}}}}}}}}}}[[[[[[[[[[[[]]]]]]]]]]????????????--------
____----------????????????]]]]]]]]]]]][[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{{}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[]]]]]]]]]]]]????????????----------____
________------------????????????]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]????????????------------________
++____________--------------??????????????]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}[[[[[[[[[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]]]??????????????--------------____________++
++++++++______________--------------????????????????]]]]]]]]]]]]]]]]]]]]]][[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[[]]]]]]]]]]]]]]]]]]]]]]????????????????--------------______________++++++++
~~++++++++++++______________----------------????????????????????]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]]????????????????????----------------______________++++++++++++~~
~~~~~~~~++++++++++++++________________------------------??????????????????????????????????]]]]]]]]]]]]]]]]]]]]??????????????????????????????????------------------________________++++++++++++++~~~~~~~~
~~~~~~~~~~~~~~++++++++++++++++__________________------------------------????????????????????????????????????????????????????????------------------------__________________++++++++++++++++~~~~~~~~~~~~~~
<<<<<<~~~~~~~~~~~~~~~~++++++++++++++++________________________----------------------------------------------------------------------------________________________++++++++++++++++~~~~~~~~~~~~~~~~<<<<<<
<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~++++++++++++++++++++__________________________________--------------------------------__________________________________++++++++++++++++++++~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<
>>>><<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++____________________________________________________________________++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<>>>>
>>>>>>>>>><<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<>>>>>>>>>>
>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>
iiiiii>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>iiiiii
iiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiii
!!iiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiii!!
!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!
!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!
ll!!!!!!!!!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!!!!!!!!!ll
llllllllll!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!llllllllll

import matplotlib.pyplot as plt
from unyt import Msun, Myr, kpc

from synthesizer.emission_models import ReprocessedEmission
from synthesizer.grid import Grid
from synthesizer.instruments import UVJ, Instrument
from synthesizer.parametric import SFH, Stars, ZDist
from synthesizer.parametric.galaxy import Galaxy
from synthesizer.parametric.morphology import Sersic2D

if __name__ == "__main__":
    # Define the morphology using a simple effective radius and slope
    morph = Sersic2D(r_eff=1 * kpc, sersic_index=1.0, ellipticity=0.5)

    # Get the location of this script, __file__ is the absolute path of this
    # script, however we just want to directory
    # script_path = os.path.abspath(os.path.dirname(__file__))

    # Define the grid
    grid_name = "test_grid"

    grid = Grid(grid_name)

    # Create the emission model
    model = ReprocessedEmission(grid)

    # Define the SFZDist
    Z_p = {"metallicity": 0.01}
    metal_dist = ZDist.DeltaConstant(**Z_p)
    sfh_p = {"max_age": 100 * Myr}
    sfh = SFH.Constant(**sfh_p)  # constant star formation
    sfzh = Stars(
        grid.log10ages,
        grid.metallicities,
        sf_hist=sfh,
        metal_dist=metal_dist,
        initial_mass=10**9 * Msun,
        morphology=morph,
    )

    # Initialise a parametric Galaxy
    galaxy = Galaxy(sfzh)

    # Generate stellar spectra
    galaxy.stars.get_spectra(model)

    # Get a UVJ filter set
    filters = UVJ()

    # Get photometry
    galaxy.stars.get_photo_lnu(filters)

    # Define geometry of the images
    resolution = 0.01 * kpc  # resolution in kpc
    npix = 100
    fov = resolution.value * npix * kpc

    # Generate images using the instrument-owned imaging path
    instrument = Instrument(
        "DemoInstrument",
        resolution=resolution,
        filters=filters,
    )
    img = galaxy.get_images_luminosity(
        "reprocessed",
        fov=fov,
        instrument=instrument,
    )

    # Make and plot an rgb image
    img.make_rgb_image(rgb_filters={"R": "J", "G": "V", "B": "U"})
    fig, ax, _ = img.plot_rgb_image()

    plt.show()

    # Print an ASCII representation of one filter image
    img["J"].print_ascii()