Why do you want a first order approximation instead of a proper second order one? And why do you want to compute unconditional welfare and the consumption equivalent based on simulated moments?
You may want to look at DSGE_mod/Born_Pfeifer_2018/Welfare at master · JohannesPfeifer/DSGE_mod · GitHub