def plotRealImg(sim, cam, rawdata, t: int, odir: Path=None, fg=None):
ncols = len(cam)<EOL>T = nans(ncols, dtype=datetime)<EOL><INDENT>if asi is not None:<EOL><INDENT>ncols=<NUM_LIT:3><EOL>if isinstance(asi,(tuple,list)):<EOL><INDENT>pass<EOL><DEDENT>elif isinstance(asi,(str,Path)):<EOL><INDENT>asi = Path(asi).expanduser()<EOL>if asi.is_dir():<EOL><INDENT>asi=list(asi.glob('<STR_LIT>'))<EOL><DEDENT><DEDENT><DEDENT><DEDENT>if fg is None:<EOL><INDENT>doclose = True<EOL>fg, axs = subplots(nrows=<NUM_LIT:1>, ncols=ncols, figsize=(<EOL><NUM_LIT:15>, <NUM_LIT:12>), dpi=DPI, facecolor='<STR_LIT>')<EOL>axs = atleast_1d(axs)  <EOL><DEDENT>else:  <EOL><INDENT>doclose = False<EOL>fg.clf()<EOL>axs = [fg.add_subplot(<NUM_LIT:1>, ncols, i + <NUM_LIT:1>) for i in range(ncols)]<EOL><DEDENT>for i, C in enumerate(cam):<EOL><INDENT>if C.usecam:  <EOL><INDENT>T[i] = updateframe(t, rawdata[i], None, cam[i], axs[i], fg)<EOL><DEDENT>elif C.name == '<STR_LIT>':  <EOL><INDENT>dasc = dio.load(C.fn, treq=T[sim.useCamBool][<NUM_LIT:0>])<EOL>C.tKeo = dasc.time<EOL>updateframe(<NUM_LIT:0>, dasc.values, dasc.wavelength, C,<EOL>axs[i], fg)  <EOL>try:<EOL><INDENT>overlayrowcol(axs[i], C.hlrows, C.hlcols)<EOL><DEDENT>except AttributeError:<EOL><INDENT>pass  <EOL><DEDENT><DEDENT>else:<EOL><INDENT>logging.error(f'<STR_LIT>')<EOL><DEDENT>if i == <NUM_LIT:0>:<EOL><INDENT>axs[<NUM_LIT:0>].set_ylabel(datetime.strftime(<EOL>T[<NUM_LIT:0>], '<STR_LIT>')).set_color('<STR_LIT>')<EOL><DEDENT><DEDENT>writeplots(fg, '<STR_LIT>', T[<NUM_LIT:0>], odir=odir,<EOL>dpi=sim.dpi, facecolor='<STR_LIT:k>', doclose=doclose)<EOL>
sim: histfeas/simclass.py
cam: camclass.py
rawdata: nframe x ny x nx ndarray
t: integer index to read
odir: output directory (where to write results)

plots both cameras together,
and magnetic zenith 1-D cut line
and 1 degree radar beam red circle centered on magnetic zenith