IRAF command or task to convert 2D_fits to 3D_fits

Please let me know the IRAF command or task which convert 2D_fits(bayered monochrome single lay data) to 3D_fits(debayerd with 3-lay color data).
I have a color CMOS integrated smart telescope(SeestarS50).
Stacked images by S50 are already 3D, but I need to check all subframes(2D) for photometry.
e.g.)
I can get G-channel data by imcopy, or imslice from 3D_fits. I need to convert 2D to 3D before it.
I know other program can handle this kind of conversion, but I prefer list processing in IRAF operation, thanks.
If this is not proper place for my topic, I will move to the best room.

Best regards,
Tatsuji