diff --git a/utils/gr-sim/dots/8086_emulator.c b/utils/gr-sim/dots/8086_emulator.c
new file mode 100644
index 00000000..4df1aa78
--- /dev/null
+++ b/utils/gr-sim/dots/8086_emulator.c
@@ -0,0 +1,247 @@
+#include <stdio.h>
+#include <unistd.h>
+
+#include "8086_emulator.h"
+
+unsigned short stack[4096];
+unsigned short ax,bx,cx,dx,si,di,bp,cs,ds,es,fs;
+int cf=0,of=0,zf=0,sf=0;
+//int sp=0;
+
+/* unsigned multiply */
+void mul_16(unsigned short value) {
+	unsigned int result;
+
+	result=ax*value;
+
+//	printf("%x*%x=%x ",value,ax,result);
+
+	ax=result&0xffff;
+	dx=result>>16;
+
+//	printf("%x:%x x=%d y=%d\n",dx,ax,
+//		dx&0xff,(dx>>8)&0xff);
+
+
+
+
+	if (dx==0) {
+		of=0;
+		cf=0;
+	}
+	else {
+		of=1;
+		cf=1;
+	}
+
+}
+
+/* 
+static void imul(short value) {
+	int result;
+
+	result=ax*value;
+
+	ax=result&0xffff;
+	dx=result>>16;
+
+}
+*/
+
+
+
+/* signed multiply */
+void imul_8(char value) {
+
+	short result;
+	char src;
+
+	src=ax;
+
+	result=src*value;
+
+//	printf("imul: %d*%d=%d ",src,value,result);
+
+	ax=result;
+
+	if (ax==(ax&0xff)) {
+		cf=0;
+		of=0;
+	}
+	else {
+		cf=1;
+		of=1;
+	}
+
+}
+
+/* signed multiply */
+/* DX:AX = AX * value */
+void imul_16(short value) {
+
+	int result;
+	short src;
+
+	src=ax;
+
+	result=src*value;
+
+//	printf("imul: %d*%d=%d ",src,value,result);
+
+	ax=(result&0xffff);
+	dx=((result>>16)&0xffff);
+
+	if (dx==0) {
+		cf=0;
+		of=0;
+	}
+	else {
+		cf=1;
+		of=1;
+	}
+
+}
+
+/* signed multiply */
+void imul_16_bx(short value) {
+
+	int result;
+	short src;
+
+	src=bx;
+
+	result=src*value;
+
+//	printf("imul: %d*%d=%d ",src,value,result);
+
+	bx=(result&0xffff);
+
+	if (bx==result) {
+		cf=0;
+		of=0;
+	}
+	else {
+		cf=1;
+		of=1;
+	}
+
+}
+
+/* signed multiply */
+void imul_16_dx(short value) {
+
+	int result;
+	short src;
+
+	src=dx;
+
+	result=src*value;
+
+//	printf("imul: %d*%d=%d ",src,value,result);
+
+	dx=(result&0xffff);
+
+	if (dx==result) {
+		cf=0;
+		of=0;
+	}
+	else {
+		cf=1;
+		of=1;
+	}
+
+}
+
+
+
+
+
+/* unsigned divide */
+void div_8(unsigned char value) {
+
+	unsigned char r,q;
+	unsigned int result,remainder;
+
+//	printf("Dividing %d (%x) by %d (%x): ",ax,ax,value,value);
+
+	if (value==0) {
+		printf("Divide by zero!\n");
+		return;
+	}
+
+	result=ax/value;
+	remainder=ax%value;
+
+	q=result;
+	r=remainder;
+
+//	printf("Result: q=%d r=%d\n",q,r);
+
+	ax=(r<<8)|(q&0xff);
+
+}
+
+
+/* signed divide */
+/* DX:AX / word.  AX=quotient, DX=remainder */
+void idiv_16(unsigned short value) {
+
+	short r,q,divisor;
+	int temp32,result,remainder;
+
+	divisor=value;
+
+	temp32=ax;
+	temp32&=0xffff;
+	temp32|=(dx<<16);
+
+//	printf("Dividing %d (%x = %04x:%04x) by %d (%x): ",
+//		temp32,temp32,dx,ax,value,value);
+
+	if (divisor==0) {
+		printf("Divide by zero!\n");
+		return;
+	}
+
+	result=temp32/divisor;
+	remainder=temp32%divisor;
+
+	q=result;
+	r=remainder;
+
+//	printf("q=%d r=%d\n",q,r);
+
+	ax=q;
+	dx=r;
+}
+
+unsigned short sar(unsigned short value,int shift) {
+
+	short temp;
+
+	temp=value;
+	temp>>=shift;
+
+	return temp;
+
+}
+
+void push(int value) {
+	//printf("Pushing %x\n",value);
+	stack[sp]=value;
+	sp++;
+}
+
+short pop(void) {
+	if (sp==0) {
+		printf("Stack underflow!\n");
+		return 0;
+	}
+
+	sp--;
+
+	//printf("Popping %x\n",stack[sp]);
+
+	return stack[sp];
+
+}
diff --git a/utils/gr-sim/dots/8086_emulator.h b/utils/gr-sim/dots/8086_emulator.h
new file mode 100644
index 00000000..9b572e45
--- /dev/null
+++ b/utils/gr-sim/dots/8086_emulator.h
@@ -0,0 +1,17 @@
+extern unsigned short stack[4096];
+extern unsigned short ax,bx,cx,dx,si,di,bp,cs,ds,es,fs;
+extern int cf,of,zf,sf;
+extern unsigned short sp;
+
+void mul_16(unsigned short value);
+void imul_8(char value);
+void imul_16(short value);
+void imul_16_bx(short value);
+void imul_16_dx(short value);
+void div_8(unsigned char value);
+void idiv_16(unsigned short value);
+
+unsigned short sar(unsigned short value,int shift);
+
+void push(int value);
+short pop(void);
diff --git a/utils/gr-sim/dots/Makefile b/utils/gr-sim/dots/Makefile
new file mode 100644
index 00000000..ca91365f
--- /dev/null
+++ b/utils/gr-sim/dots/Makefile
@@ -0,0 +1,30 @@
+CC = gcc
+CFLAGS = -O2 -Wall -g
+
+SDL_LIBS= `sdl-config --libs`
+SDL_INCLUDE= `sdl-config --cflags`
+
+all:	dots
+
+###
+
+dots:	dots.o vga_emulator.o 8086_emulator.o ../gr-sim.a
+	$(CC) -o dots dots.o vga_emulator.o 8086_emulator.o ../gr-sim.a $(LFLAGS) $(SDL_LIBS)
+
+dots.o:	dots.c vga_emulator.h
+	$(CC) $(CFLAGS) $(SDL_INCLUDE) -c dots.c
+
+###
+
+8086_emulator.o:	8086_emulator.c 8086_emulator.h
+	$(CC) $(CFLAGS) -c 8086_emulator.c
+
+###
+
+vga_emulator.o:	vga_emulator.c vga_emulator.h
+	$(CC) $(CFLAGS) $(SDL_INCLUDE) -c vga_emulator.c
+
+###
+
+clean:	
+	rm -f *~ *.o dots hellmood_memories
diff --git a/utils/gr-sim/dots/dots.c b/utils/gr-sim/dots/dots.c
new file mode 100644
index 00000000..2e63355b
--- /dev/null
+++ b/utils/gr-sim/dots/dots.c
@@ -0,0 +1,916 @@
+#include <stdio.h>
+#include <string.h>
+#include <stdlib.h>
+#include <unistd.h>
+#include <malloc.h>
+#include <math.h>
+
+#include "8086_emulator.h"
+#include "vga_emulator.h"
+
+#include "../gr-sim.h"
+#include "../tfv_zp.h"
+
+#include "sin1024.h"
+
+#define	MAXDOTS	1024
+
+#define BOTTOM 8000
+
+static short gravitybottom=BOTTOM;
+
+static short bpmin=30000;
+static short bpmax=-30000;
+static short gravity=0;
+static short dotnum=0;
+static short gravityd=16;
+
+//???,-1280,-960,-640,-320};
+
+static short rows[200];
+
+//dot dw	MAXDOTS dup(0,0,0,0,0,0,0,0) ;x,y,z,oldposshadow,oldpos,-,-,-
+
+static struct {
+	short	x;	// 0
+	short	y;	// 2
+	short	z;	// 4
+	short	old1;	// 6 oldpos shadow
+	short	old2;	// 8 oldpos
+	short	old3;	// 10
+	short	old4;	// 12
+	short	yadd;	// 14
+} dot[MAXDOTS];
+
+static short rotsin=0;
+static short rotcos=0;
+
+static char *bgpic;
+
+static int depthtable1[128];
+static int depthtable2[128];
+static int depthtable3[128];
+//static int depthtable4[128];
+
+static unsigned char depthtable1_bytes[512];
+static unsigned char depthtable2_bytes[512];
+static unsigned char depthtable3_bytes[512];
+//static unsigned char depthtable4_bytes[512];
+
+
+static void drawdots(void) {
+	int temp32;
+	int yy;
+
+				//	CBEG
+	ax=0xa000;			// mov	ax,0a000h
+	es=ax;				// mov	es,ax
+	ax=cs;				// mov	ax,cs
+	ds=ax;				// mov	ds,ax
+
+	/* why [2]? */
+	fs=bgpic[2];			// mov	fs,cs:_bgpic[2]
+	cx=dotnum;			// mov	cx,cs:_dotnum
+	si=0;				// mov	si,OFFSET dot
+
+label1:
+	push(cx);			// push	cx
+	ax=dot[si].x;			// mov	ax,ds:[si+0] ;X
+	imul_16(rotsin);		// imul	ds:_rotsin
+	ax=ax;				// mov	ax,ax
+	cx=dx;				// mov	cx,dx
+	ax=dot[si].z;			// mov	ax,ds:[si+4] ;Z
+	imul_16(rotcos);		// imul	ds:_rotcos
+	ax=ax-bx;			// sub	ax,bx
+	dx=dx-cx;			// sub	dx,cx
+	bp=dx;				// mov	bp,dx
+	bp=bp+9000;			// add	bp,9000
+
+	ax=dot[si].x;			// mov	ax,ds:[si+0] ;X
+	imul_16(rotcos);		// imul	ds:_rotcos
+	bx=ax;				// mov	bx,ax
+	cx=dx;				// mov	cx,dx
+	ax=dot[si].z;			// mov	ax,ds:[si+4] ;Z
+	imul_16(rotsin);		// imul	ds:_rotsin
+
+	temp32=ax+bx;			// add	ax,bx
+	ax=ax+bx;			//
+	dx=dx+cx;			// adc	dx,cx
+	if (temp32&(1<<16)) dx=dx+1;
+
+
+//	printf("Before: ax=0x%04X dx=%04X\n",ax,dx);
+	ax=(ax>>8)|(dx<<8);		// shrd	ax,dx,8
+	dx=sar(dx,8);			// sar	dx,8
+//	printf("After: ax=0x%04X dx=%04X\n",ax,dx);
+
+	bx=ax;				// mov	bx,ax
+	cx=dx;				// mov	cx,dx
+	ax=(ax>>3)|(dx<<13);		// shrd	ax,dx,3
+
+	dx=sar(dx,3);			// sar	dx,3
+	temp32=ax+bx;			// add	ax,bx
+	ax=ax+bx;
+	dx=dx+cx;			// adc	dx,cx
+	if (temp32&(1<<16)) dx=dx+1;
+
+	temp32=(dx<<16)|(ax&0xfffff);
+	idiv_16(bp);			// idiv bp
+	ax=ax+160;			// add	ax,160
+	push(ax);			// push	ax
+
+	/* if off end of screen, no need for shadow */
+
+	if (ax>319) goto label2;	// cmp	ax,319
+
+					// ja	@@2
+	/**********/
+	/* shadow */
+	/**********/
+
+	ax=0;				// xor	ax,ax
+	dx=8;				// mov	dx,8
+	idiv_16(bp);			// idiv	bp
+	ax=ax+100;			// add	ax,100
+
+	/* if shadow off screen, don't draw */
+	if (ax>199) goto label2;	// cmp	ax,199
+					// ja	@@2
+	bx=ax;				// mov	bx,ax
+
+	// not needed, it's a C array
+	//bx=bx<<1;			// shl	bx,1
+	bx=rows[bx];			// mov	bx,ds:_rows[bx]
+	ax=pop();			// pop	ax
+	bx=bx+ax;			// add	bx,ax
+	push(ax);			// push	ax
+
+//	printf("Drawing shadow at %d,%d\n",bx%320,bx/320);
+
+
+
+	/* erase old shadow (?)*/
+	di=dot[si].old1;		// mov	di,ds:[si+6]
+
+	yy=((di/320)*48)/200;
+	if (yy>23) color_equals(5);
+	else color_equals(0);
+	plot( (di%320)/8,yy);
+
+	//ax=bgpic[di];			// mov	ax,fs:[di]
+	//framebuffer[di]=ax;		// mov	es:[di],ax
+
+//	framebuffer[di]=bgpic[di];
+//	framebuffer[di+1]=bgpic[di+1];
+	framebuffer_write(di,bgpic[di]);
+	framebuffer_write(di+1,bgpic[di+1]);
+
+
+
+	/* draw new shadow (?) */
+//	ax=87+87*256;			// mov	ax,87+87*256
+//	framebuffer[bx]=ax;		// mov	word ptr es:[bx],ax
+
+//	framebuffer[bx]=87;
+//	framebuffer[bx+1]=87;
+	framebuffer_write(bx,87);
+	framebuffer_write(bx+1,87);
+
+//	bx/320 -> 200     200->48  *48/200
+
+	yy=((bx/320)*48)/200;
+	color_equals(0);
+	plot( (bx%320)/8,yy);
+
+//	printf("Plotting at %d,%d\n",(bx%320)/8,(bx/320)/5);
+
+	/* save this to erase next time */
+	dot[si].old1=bx;		// mov	ds:[si+6],bx
+
+	/********/
+	/* ball */
+	/********/
+//	ax=gravity;			// mov	ax,ds:_gravity
+//	dot[si].yadd+=ax;		// add	ds:[si+14],ax
+
+//	if (si==100) printf("Gravity: %hd (%04x)  Yadd: %hd (%04x)\n",
+//			gravity,gravity,dot[si].yadd,dot[si].yadd);
+
+	dot[si].yadd+=gravity;
+
+//	if (si==100) printf("\tyadd after yadd+=gravity: %hd (%04x)\n",
+//			dot[si].yadd,dot[si].yadd);
+
+//	ax=dot[si].y;			// mov	ax,ds:[si+2] ;Y
+//	ax+=dot[si].yadd;		// add	ax,ds:[si+14]
+	ax=dot[si].y+dot[si].yadd;
+
+//	if (si==100) printf("\tax=y+yadd: %hu (%04x) = "
+//			"%hd (%04x) + %hd (%04x)\n",
+//			ax,ax,dot[si].y,dot[si].y,dot[si].yadd,dot[si].yadd);
+
+//	if (si==100) printf("\tcomparing if (ax<gravitybottom): "
+//			"%hu < %hd\n",ax,gravitybottom);
+
+
+	temp32=ax;
+	if (temp32&0x8000) temp32|=0xffff0000;
+	if (temp32<gravitybottom) goto label4; //cmp	ax,ds:_gravitybottom
+					// jl	@@4
+
+//	if (si==100) printf("\twas greater than (not less)\n");
+
+	push(ax);			// push	ax
+
+//	ax=dot[si].yadd;		// mov	ax,ds:[si+14]
+//	ax=-ax;				// neg	ax
+	ax=-dot[si].yadd;
+
+//	if (si==100) printf("\tax is -yadd: %hu %x\n",ax,ax);
+
+//	if (si==100) printf("\tabout to multiply gravityd*ax: "
+//			"%hd (%x) * %hu (%x)\n",gravityd,gravityd,
+//			ax,ax);
+	imul_16(gravityd);		// imul	cs:_gravityd
+
+//	if (si==100) printf("\tresult dx:ax is %x:%x (%hu)\n",
+//			dx,ax,ax);
+
+	ax=sar(ax,4);			// sar	ax,4
+
+//	if (si==100) printf("\tyadd=(ax>>4 is %x (%hu))\n",
+//			ax,ax);
+	dot[si].yadd=ax;		// mov	ds:[si+14],ax
+	ax=pop();			// pop	ax
+//	if (si==100) printf("\trestoring ax=%x, adding yadd %x\n",
+//		ax,dot[si].yadd);
+
+	ax+=dot[si].yadd;		// add	ax,ds:[si+14]
+
+//	if (si==100) printf("\tax=%x\n",ax);
+
+label4:
+	dot[si].y=ax;			// mov	ds:[si+2],ax
+
+//	if (si==100) printf("\tdot[si].y=%x\n",dot[si].y);
+
+	if (ax&0x8000) {		// cwd
+		dx=0xffff;
+	}
+	else {
+		dx=0;
+	}
+
+//	if (si==100) printf("\tdx:ax = %04hx:%04hx\n",dx,ax);
+
+	dx=(dx<<6)|(ax>>10);		// shld	dx,ax,6
+	ax=ax<<6;			// shl	ax,6
+
+//	if (si==100) printf("\tdx:ax <<6 = %04hx:%04hx, bp=%04hx\n",dx,ax,bp);
+
+	idiv_16(bp);			// idiv	bp
+//	if (si==100) printf("\tY ax=%d\n",ax);
+	ax=ax+100;			// add	ax,100
+	if (ax>199) goto label3;	// cmp	ax,199
+					// ja	@@3
+//	if (si==100) printf("\tdraw Y ax=%d\n",ax);
+	bx=ax;				// mov	bx,ax
+
+	// not needed, C array
+	//bx=bx<<1;			// shl	bx,1
+	bx=rows[bx];			// mov	bx,ds:_rows[bx]
+
+	ax=pop();			// pop	ax
+	bx=bx+ax;			// add	bx,ax
+
+	di=dot[si].old2;		// mov	di,ds:[si+8]
+//	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+//	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+//	framebuffer[di+2]=bgpic[di+2];
+//	framebuffer[di+3]=bgpic[di+3];
+
+	yy=((di/320)*48)/200;
+	if (yy>23) color_equals(5);
+	else color_equals(0);
+	plot( (di%320)/8,yy);
+
+	framebuffer_write(di,bgpic[di]);
+	framebuffer_write(di+1,bgpic[di+1]);
+	framebuffer_write(di+2,bgpic[di+2]);
+	framebuffer_write(di+3,bgpic[di+3]);
+
+
+
+
+	di=di+320;			// add	di,320
+//	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+//	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+//	framebuffer[di+2]=bgpic[di+2];
+//	framebuffer[di+3]=bgpic[di+3];
+
+	framebuffer_write(di,bgpic[di]);
+	framebuffer_write(di+1,bgpic[di+1]);
+	framebuffer_write(di+2,bgpic[di+2]);
+	framebuffer_write(di+3,bgpic[di+3]);
+
+
+	di=di+320;			// add	di,320
+//	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+//	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+//	framebuffer[di+2]=bgpic[di+2];
+//	framebuffer[di+3]=bgpic[di+3];
+
+	framebuffer_write(di,bgpic[di]);
+	framebuffer_write(di+1,bgpic[di+1]);
+	framebuffer_write(di+2,bgpic[di+2]);
+	framebuffer_write(di+3,bgpic[di+3]);
+
+				//;;	add	di,320
+				//;;	mov	eax,fs:[di]
+				//;;	mov	es:[di],eax
+
+	bp=bp>>6;		// shr	bp,6
+	bp=bp&(~3L);		// and	bp,not 3
+
+	temp32=bp;
+	if (temp32&0x8000) temp32|=0xffff0000;
+	if (temp32>=bpmin) goto label_t1;	// cmp	bp,cs:_bpmin
+					// jge	@@t1
+	bpmin=bp;			// mov	cs:_bpmin,bp
+label_t1:
+	temp32=bp;
+	if (temp32&0x8000) temp32|=0xffff0000;
+	if (temp32<=bpmax) goto label_t2;	// cmp	bp,cs:_bpmax
+					// jle	@@t2
+	bpmax=bp;			// mov	cs:_bpmax,bp
+label_t2:
+//	eax=depthtable1[bp];		// mov	ax,word ptr ds:_depthtable1[bp]
+					// mov	word ptr es:[bx+1],ax
+
+
+	yy=((bx/320)*48)/200;
+	color_equals(6);
+	plot( (bx%320)/8,yy);
+
+	framebuffer[bx+1]=depthtable1_bytes[bp];
+	framebuffer[bx+2]=depthtable1_bytes[bp+1];
+
+	//eax=depthtable2[bp];		// mov	eax,ds:_depthtable2[bp]
+					// mov	dword ptr es:[bx+320],eax
+	framebuffer[bx+320]=depthtable2_bytes[bp];
+	framebuffer[bx+321]=depthtable2_bytes[bp+1];
+	framebuffer[bx+322]=depthtable2_bytes[bp+2];
+	framebuffer[bx+323]=depthtable2_bytes[bp+3];
+
+//	eax=depthtable3[bp];	// mov	ax,word ptr ds:_depthtable3[bp]
+				// mov	word ptr es:[bx+641],ax
+	framebuffer[bx+641]=depthtable3_bytes[bp];
+	framebuffer[bx+642]=depthtable3_bytes[bp+1];
+	dot[si].old2=bx;	// mov	ds:[si+8],bx
+
+
+//labelz:
+	cx=pop();		// pop	cx
+	si=si+1;		// add	si,16	point to next dot
+	cx=cx-1;
+	if (cx!=0) goto label1;	// loop	@@1
+label0:
+	return;
+				// @@0:	CEND
+
+label2:
+	/* This is called when we are off the screen */
+	/* erases old but didn't draw new */
+
+	/* erase old dot */
+	di=dot[si].old2;		// mov	di,ds:[si+8]
+
+
+	yy=((di/320)*48)/200;
+	if (yy>23) color_equals(5);
+	else color_equals(0);
+	plot( (di%320)/8,yy);
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	di=di+320;			// add	di,320
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	di=di+320;			// add  di,320
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	ax=(framebuffer[di]|(framebuffer[di+1]<<8));
+
+	/* doing something to shadow here? */
+
+	di=dot[si].old1;		// mov	di,ds:[si+6]
+	dot[si].old1=ax;		// mov	ds:[si+6],ax
+
+	framebuffer[di]=bgpic[di];	// mov	ax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],ax
+
+	bx=pop();			// pop	bx
+	cx=pop();			// pop	cx
+	si=si+1;			// add	si,16
+	cx=cx-1;			// loop	@@1
+	if (cx!=0) goto label1;
+	goto label0;			// jmp	@@0
+
+label3:
+	/* erase old dot */
+	di=dot[si].old2;		// mov	di,ds:[si+8]
+
+	yy=((di/320)*48)/200;
+	if (yy>23) color_equals(5);
+	else color_equals(0);
+	plot( (di%320)/8,yy);
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	di=di+320;			// add	di,320
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	di=di+320;			// add	di,320
+
+	framebuffer[di]=bgpic[di];	// mov	eax,fs:[di]
+	framebuffer[di+1]=bgpic[di+1];	// mov	es:[di],eax
+	framebuffer[di+2]=bgpic[di+2];
+	framebuffer[di+3]=bgpic[di+3];
+
+	bx=pop();			// pop	bx
+	cx=pop();			// pop	cx
+	si=si+1;			// add	si,16
+	cx=cx-1;			// loop	@@1
+	if (cx!=0) goto label1;
+	goto label0;			// jmp	@@0
+
+}
+
+static void setpalette(char *pal) {
+
+	int c;
+
+	// push	bp
+	// mov	bp,sp
+	// push	si
+	// push	di
+	// push	ds
+	// mov	si,[bp+6]
+	// mov	ds,[bp+8]
+					// mov	dx,3c8h
+					// mov	al,0
+	outp(0x3c8,0);			//out	dx,al
+
+	for(c=0;c<768;c++) outp(0x3c9,pal[c]);
+	grsim_update();
+					// inc	dx
+					//mov	cx,768
+					//rep	outsb
+					//	pop	ds
+					//	pop	di
+					//	pop	si
+					//	pop	bp
+					//	ret
+
+}
+
+//short face[]={
+//	2248,-2306,0,		// from face.inc
+//	30000,30000,30000
+//};
+
+/* wait for VGA border start */
+static short dis_waitb(void) {
+
+// descr: Waits for border start
+// waitb_1 PROC NEAR
+//        call    checkkeys
+
+//        IFDEF INDEMO
+//        sti
+//        mov     ax,cs:copperframecount
+//@@v:    cmp     cs:copperframecount,ax
+//        je      @@v
+//@@q:    mov     ax,cs:copperframecount
+//       mov     cs:copperframecount,0
+//        ELSE
+
+//        mov     dx,3dah
+//@@1:    in      al,dx
+//        test    al,8
+//        jnz     @@1
+//@@2:    in      al,dx
+//        test    al,8
+ //       jz      @@2
+
+  //      mov     ax,1 ;number of frames taken            ;TEMP!
+
+//        ENDIF
+//        ret
+//waitb_1 ENDP
+
+
+	/* approximate by 70Hz sleep */
+	usleep(14286);
+
+	return 1;
+}
+
+static short dis_exit(void) {
+	return 0;
+}
+
+static short dis_indemo(void) {
+	return 0;
+}
+
+//char far *vram=(char far *)0xa0000000L;
+static unsigned char *vram=framebuffer;
+
+static char pal[768];
+static char pal2[768];
+
+static short isin(short deg) {
+	return(sin1024[deg&1023]);
+}
+
+static short icos(short deg) {
+	return(sin1024[(deg+256)&1023]);
+}
+
+static void setborder(short color) {
+
+	//printf("Setting border to %d\n",color);
+
+	// to write attribute register:
+	//	read/write address to $3c0
+	//	data written to $3c0, read from $3c1
+	// flip flop tracks if it's index/data, you reset that
+	//	by reading $3da
+
+	//mov	dx,3dah			// input status reg #1
+	//in	al,dx			// resets index/addr flip-flop
+
+	//mov	dx,3c0h			// attribute access
+	//mov	al,11h+32		// $11=overscan (border color)
+					// 32 is PAS bit
+	//out	dx,al
+
+	//mov	al,color
+	//out	dx,al
+}
+
+static short cols[]={
+0,0,0,
+4,25,30,
+8,40,45,
+16,55,60};
+
+static short dottaul[1024];
+
+int main(int argc,char **argv) {
+
+//	short timer=30000;
+	short dropper,repeat;
+	short frame=0;
+	short rota=-1*64;
+//	short fb=0;
+	short rot=0,rots=0;
+	short a,b,c,d,i,j=0;//,mode;
+	short grav,gravd;
+	short f=0;
+	int ch;
+
+	//dis_partstart();
+
+	dotnum=512;
+	for(a=0;a<dotnum;a++) {
+		dottaul[a]=a;
+	}
+
+	for(a=0;a<500;a++) {
+		b=rand()%dotnum;
+		c=rand()%dotnum;
+		d=dottaul[b];
+		dottaul[b]=dottaul[c];
+		dottaul[c]=d;
+	}
+
+	dropper=22000;
+
+	for(a=0;a<dotnum;a++) {
+		dot[a].x=0;
+		dot[a].y=2560-dropper;
+		dot[a].z=0;
+		dot[a].yadd=0;
+	}
+	//mode=7;
+	grav=3;
+	gravd=13;
+	gravitybottom=8105;
+	i=-1;
+
+	for(a=0;a<500;a++) { // scramble
+		b=rand()%dotnum;
+		c=rand()%dotnum;
+		d=dot[b].x; dot[b].x=dot[c].x; dot[c].x=d;
+		d=dot[b].y; dot[b].y=dot[c].y; dot[c].y=d;
+		d=dot[b].z; dot[b].z=dot[c].z; dot[c].z=d;
+	}
+
+	/* setup rows lookup table */
+	for(a=0;a<200;a++) {
+		rows[a]=a*320;
+	}
+
+//	set_default_pal();
+
+//	mode13h_graphics_init("dots",2);
+
+	grsim_init();
+
+	gr();
+
+	clear_screens();
+	soft_switch(MIXCLR);
+
+	ram[DRAW_PAGE]=0;
+
+	color_equals(5);
+	for(a=24;a<48;a++) hlin(0,0,40,a);
+
+
+//	set mode 13h
+//	_asm mov ax,13h
+//	_asm int 10h
+
+	/* set palette address to 0 */
+	outp(0x3c8,0);
+
+	/* set up colors for first 64 colors */
+	for(a=0;a<16;a++) {
+		for(b=0;b<4;b++) {
+			c=100+a*9;
+			outp(0x3c9,cols[b*3+0]);
+			outp(0x3c9,cols[b*3+1]*c/256);
+			outp(0x3c9,cols[b*3+2]*c/256);
+		}
+	}
+
+	/* set palette for color 255 */
+	outp(0x3c8,255);
+
+	/* some sort of purplish color? */
+	outp(0x3c9,31);
+	outp(0x3c9,0);
+	outp(0x3c9,15);
+
+
+	/* set colors starting from 64 ... 164 */
+	/* looks like a grey gradient of some sort */
+	outp(0x3c8,64);
+
+	for(a=0;a<100;a++) {
+		c=64-256/(a+4);
+		c=c*c/64;
+		outp(0x3c9,c/4);
+		outp(0x3c9,c/4);
+		outp(0x3c9,c/4);
+	}
+
+	/* read out the VGA card's idea of palette (?) */
+	outp(0x3c7,0);
+	for(a=0;a<768;a++) pal[a]=inp(0x3c9);
+
+	/* clear palette to all 0 */
+	/* this lets up setup background while not visible */
+	outp(0x3c8,0);
+	for(a=0;a<768;a++) outp(0x3c9,0);
+
+	/* put grey gradient on bottom half of screen? */
+	for(a=0;a<100;a++) {
+		memset(vram+(100+a)*320,a+64,320);
+	}
+
+	/* set up depth table? */
+	/* this has further away balls a darker color */
+	for(a=0;a<128;a++) {
+		c=a-(43+20)/2;
+		c=c*3/4;
+		c+=8;
+		if(c<0) c=0; else if(c>15) c=15;
+		c=15-c;
+		depthtable1[a]=0x202+0x04040404*c;
+		depthtable2[a]=0x02030302+0x04040404*c;
+		depthtable3[a]=0x202+0x04040404*c;
+		//depthtable4[a]=0x02020302+0x04040404*c;
+	}
+
+	/* make a byte-wise copy of this */
+	/* the original code just indexes byte-wise into integer data */
+	/* which is a pain */
+	for(a=0;a<128;a++) {
+		depthtable1_bytes[(a*4)+0]=(depthtable1[a]>>0)&0xff;
+		depthtable1_bytes[(a*4)+1]=(depthtable1[a]>>8)&0xff;
+		depthtable1_bytes[(a*4)+2]=(depthtable1[a]>>16)&0xff;
+		depthtable1_bytes[(a*4)+3]=(depthtable1[a]>>24)&0xff;
+
+		depthtable2_bytes[(a*4)+0]=(depthtable2[a]>>0)&0xff;
+		depthtable2_bytes[(a*4)+1]=(depthtable2[a]>>8)&0xff;
+		depthtable2_bytes[(a*4)+2]=(depthtable2[a]>>16)&0xff;
+		depthtable2_bytes[(a*4)+3]=(depthtable2[a]>>24)&0xff;
+
+		depthtable3_bytes[(a*4)+0]=(depthtable3[a]>>0)&0xff;
+		depthtable3_bytes[(a*4)+1]=(depthtable3[a]>>8)&0xff;
+		depthtable3_bytes[(a*4)+2]=(depthtable3[a]>>16)&0xff;
+		depthtable3_bytes[(a*4)+3]=(depthtable3[a]>>24)&0xff;
+
+
+
+	}
+
+
+	/* allocate space for background */
+	//bgpic=halloc(64000L,1L);
+	bgpic=calloc(65536L,1L);
+
+	/* backup background */
+	memcpy(bgpic,vram,64000);
+
+	grsim_update();
+
+	/* Fade back in from black to palette */
+	a=0;
+	for(b=64;b>=0;b--) {
+		for(c=0;c<768;c++) {
+			a=pal[c]-b;
+			if(a<0) a=0;
+			pal2[c]=a;
+		}
+		/* wait for retrace (delay) */
+		dis_waitb();
+		dis_waitb();
+		outp(0x3c8,0);
+		for(c=0;c<768;c++) outp(0x3c9,pal2[c]);
+		grsim_update();
+	}
+
+	while(!dis_exit() && frame<2450) {
+
+		/* code sets border color */
+		/* then waits for it to end, as a timing thing? */
+		setborder(0);
+
+		/* when not in demo this defaults to 1? */
+		repeat=dis_waitb();
+
+		if(frame>2300) setpalette(pal2);
+
+		setborder(1);
+
+		if(dis_indemo()) {
+			/* ?? music synchronization? */
+//			a=dis_musplus();
+//			if(a>-4 && a<0) break;
+		}
+
+		repeat=1;
+
+		while(repeat--) {
+
+		frame++;
+		if(frame==500) f=0;
+		i=dottaul[j];
+		j++; j%=dotnum;
+
+		/* initial spin */
+		if(frame<500) {
+			dot[i].x=isin(f*11)*40;
+			dot[i].y=icos(f*13)*10-dropper;
+			dot[i].z=isin(f*17)*40;
+			dot[i].yadd=0;
+//			printf("%d: %d,%d,%d,%d\n",i,
+//				dot[i].x,dot[i].y,dot[i].z,dot[i].yadd);
+		}
+		/* bouncing ring */
+		else if(frame<900) {
+			dot[i].x=icos(f*15)*55;
+			dot[i].y=dropper;
+			dot[i].z=isin(f*15)*55;
+			dot[i].yadd=-260;
+		}
+		/* fountain */
+		else if(frame<1700) {
+			a=sin1024[frame&1023]/8;
+			dot[i].x=icos(f*66)*a;
+			dot[i].y=8000;
+			dot[i].z=isin(f*66)*a;
+			dot[i].yadd=-300;
+		}
+		/* swirling */
+		else if(frame<2360) {
+				/*
+				a=rand()/128+32;
+				dot[i].y=8000-a*80;
+				b=rand()&1023;
+				a+=rand()&31;
+				dot[i].x=sin1024[b]*a/3+(a-50)*7;
+				dot[i].z=sin1024[(b+256)&1023]*a/3+(a-40)*7;
+				dot[i].yadd=300;
+				if(frame>1640 && !(frame&31) && grav>-2) grav--;
+				*/
+			dot[i].x=rand()-16384;
+			dot[i].y=8000-rand()/2;
+			dot[i].z=rand()-16384;
+			dot[i].yadd=0;
+			if(frame>1900 && !(frame&31) && grav>0) grav--;
+		}
+		/* palette to white */
+		else if(frame<2400) {
+			a=frame-2360;
+			for(b=0;b<768;b+=3) {
+				c=pal[b+0]+a*3;
+				if(c>63) c=63;
+				pal2[b+0]=c;
+				c=pal[b+1]+a*3;
+				if(c>63) c=63;
+				pal2[b+1]=c;
+				c=pal[b+2]+a*4;
+				if(c>63) c=63;
+				pal2[b+2]=c;
+			}
+		}
+		/* palette to black */
+		else if(frame<2440) {
+			a=frame-2400;
+			for(b=0;b<768;b+=3) {
+				c=63-a*2;
+				if(c<0) c=0;
+				pal2[b+0]=c;
+				pal2[b+1]=c;
+				pal2[b+2]=c;
+			}
+		}
+		if(dropper>4000) dropper-=100;
+
+		rotcos=icos(rot)*64; rotsin=isin(rot)*64;
+		rots+=2;
+
+		if(frame>1900) {
+			rot+=rota/64;
+			rota--;
+		}
+		else rot=isin(rots);
+
+		f++;
+		gravity=grav;
+		gravityd=gravd;
+
+		}
+
+
+		drawdots();
+
+		grsim_update();
+
+#if 0
+		for(i=0;i<dotnum;i++) {
+			printf("%d: %d,%d,%d  %d\n",i,
+				dot[i].x,dot[i].y,dot[i].z,dot[i].yadd);
+		}
+#endif
+
+again:
+		ch=grsim_input();
+		if (ch==27) {
+			return 0;
+		}
+//		else if (ch==0) goto again;
+
+	}
+
+//	restores 80x25 color text mode
+//	if (!dis_indemo()) {
+//		_asm mov ax,3h
+//		_asm int 10h
+//	}
+	return 0;
+}
diff --git a/utils/gr-sim/dots/sin1024.h b/utils/gr-sim/dots/sin1024.h
new file mode 100644
index 00000000..5f0d8c83
--- /dev/null
+++ b/utils/gr-sim/dots/sin1024.h
@@ -0,0 +1,59 @@
+int16_t sin1024[]={
+   0,  1,  3,  4,  6,  7,  9, 10, 12, 14, 15, 17, 18, 20, 21, 23,
+  25, 26, 28, 29, 31, 32, 34, 36, 37, 39, 40, 42, 43, 45, 46, 48,
+  49, 51, 53, 54, 56, 57, 59, 60, 62, 63, 65, 66, 68, 69, 71, 72,
+  74, 75, 77, 78, 80, 81, 83, 84, 86, 87, 89, 90, 92, 93, 95, 96,
+  97, 99,100,102,103,105,106,108,109,110,112,113,115,116,117,119,
+ 120,122,123,124,126,127,128,130,131,132,134,135,136,138,139,140,
+ 142,143,144,146,147,148,149,151,152,153,155,156,157,158,159,161,
+ 162,163,164,166,167,168,169,170,171,173,174,175,176,177,178,179,
+ 181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,
+ 197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,211,
+ 212,213,214,215,216,217,217,218,219,220,221,221,222,223,224,225,
+ 225,226,227,227,228,229,230,230,231,232,232,233,234,234,235,235,
+ 236,237,237,238,238,239,239,240,241,241,242,242,243,243,244,244,
+ 244,245,245,246,246,247,247,247,248,248,249,249,249,250,250,250,
+ 251,251,251,251,252,252,252,252,253,253,253,253,254,254,254,254,
+ 254,254,255,255,255,255,255,255,255,255,255,255,255,255,255,255,
+ 256,255,255,255,255,255,255,255,255,255,255,255,255,255,255,254,
+ 254,254,254,254,254,253,253,253,253,252,252,252,252,251,251,251,
+ 251,250,250,250,249,249,249,248,248,247,247,247,246,246,245,245,
+ 244,244,244,243,243,242,242,241,241,240,239,239,238,238,237,237,
+ 236,235,235,234,234,233,232,232,231,230,230,229,228,227,227,226,
+ 225,225,224,223,222,221,221,220,219,218,217,217,216,215,214,213,
+ 212,211,211,210,209,208,207,206,205,204,203,202,201,200,199,198,
+ 197,196,195,194,193,192,191,190,189,188,187,186,185,184,183,182,
+ 181,179,178,177,176,175,174,173,171,170,169,168,167,166,164,163,
+ 162,161,159,158,157,156,155,153,152,151,149,148,147,146,144,143,
+ 142,140,139,138,136,135,134,132,131,130,128,127,126,124,123,122,
+ 120,119,117,116,115,113,112,110,109,108,106,105,103,102,100, 99,
+  97, 96, 95, 93, 92, 90, 89, 87, 86, 84, 83, 81, 80, 78, 77, 75,
+  74, 72, 71, 69, 68, 66, 65, 63, 62, 60, 59, 57, 56, 54, 53, 51,
+  49, 48, 46, 45, 43, 42, 40, 39, 37, 36, 34, 32, 31, 29, 28, 26,
+  25, 23, 21, 20, 18, 17, 15, 14, 12, 10,  9,  7,  6,  4 , 3 , 1,
+   0, -1, -3, -4, -6, -7, -9,-10,-12,-14,-15,-17,-18,-20,-21,-23,
+ -25,-26,-28,-29,-31,-32,-34,-36,-37,-39,-40,-42,-43,-45,-46,-48,
+ -49,-51,-53,-54,-56,-57,-59,-60,-62,-63,-65,-66,-68,-69,-71,-72,
+ -74,-75,-77,-78,-80,-81,-83,-84,-86,-87,-89,-90,-92,-93,-95,-96,
+ -97, -99,-100,-102,-103,-105,-106,-108,-109,-110,-112,-113,-115,-116,-117,-119,
+-120,-122,-123,-124,-126,-127,-128,-130,-131,-132,-134,-135,-136,-138,-139,-140,
+-142,-143,-144,-146,-147,-148,-149,-151,-152,-153,-155,-156,-157,-158,-159,-161,
+-162,-163,-164,-166,-167,-168,-169,-170,-171,-173,-174,-175,-176,-177,-178,-179,
+-181,-182,-183,-184,-185,-186,-187,-188,-189,-190,-191,-192,-193,-194,-195,-196,
+-197,-198,-199,-200,-201,-202,-203,-204,-205,-206,-207,-208,-209,-210,-211,-211,
+-212,-213,-214,-215,-216,-217,-217,-218,-219,-220,-221,-221,-222,-223,-224,-225,
+-225,-226,-227,-227,-228,-229,-230,-230,-231,-232,-232,-233,-234,-234,-235,-235,
+-236,-237,-237,-238,-238,-239,-239,-240,-241,-241,-242,-242,-243,-243,-244,-244,
+-244,-245,-245,-246,-246,-247,-247,-247,-248,-248,-249,-249,-249,-250,-250,-250,
+-251,-251,-251,-251,-252,-252,-252,-252,-253,-253,-253,-253,-254,-254,-254,-254,
+-254,-254,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,
+-256,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-255,-254,
+-254,-254,-254,-254,-254,-253,-253,-253,-253,-252,-252,-252,-252,-251,-251,-251,
+-251,-250,-250,-250,-249,-249,-249,-248,-248,-247,-247,-247,-246,-246,-245,-245,-244,-244,-244,-243,-243,-242,-242,-241,-241,-240,-239,-239,-238,-238,-237,-237,
+-236,-235,-235,-234,-234,-233,-232,-232,-231,-230,-230,-229,-228,-227,-227,-226,-225,-225,-224,-223,-222,-221,-221,-220,-219,-218,-217,-217,-216,-215,-214,-213,
+-212,-211,-211,-210,-209,-208,-207,-206,-205,-204,-203,-202,-201,-200,-199,-198,-197,-196,-195,-194,-193,-192,-191,-190,-189,-188,-187,-186,-185,-184,-183,-182,
+-181,-179,-178,-177,-176,-175,-174,-173,-171,-170,-169,-168,-167,-166,-164,-163,-162,-161,-159,-158,-157,-156,-155,-153,-152,-151,-149,-148,-147,-146,-144,-143,
+-142,-140,-139,-138,-136,-135,-134,-132,-131,-130,-128,-127,-126,-124,-123,-122,-120,-119,-117,-116,-115,-113,-112,-110,-109,-108,-106,-105,-103,-102,-100,-99,
+-97,-96,-95,-93,-92,-90,-89,-87,-86,-84,-83,-81,-80,-78,-77,-75,-74,-72,-71,-69,-68,-66,-65,-63,-62,-60,-59,-57,-56,-54,-53,-51,
+-49,-48,-46,-45,-43,-42,-40,-39,-37,-36,-34,-32,-31,-29,-28,-26,-25,-23,-21,-20,-18,-17,-15,-14,-12,-10,-9,-7,-6,-4,-3,-1
+};
diff --git a/utils/gr-sim/dots/vga_emulator.c b/utils/gr-sim/dots/vga_emulator.c
new file mode 100644
index 00000000..f2ac239f
--- /dev/null
+++ b/utils/gr-sim/dots/vga_emulator.c
@@ -0,0 +1,319 @@
+#include <stdio.h>
+#include <unistd.h>
+
+#include <SDL.h>
+
+#include "vga_emulator.h"
+
+/* technically bigger than 320x200 but some code writes off edge... */
+/* in theory 16-bit pointer can't write beyond this */
+unsigned char framebuffer[65536];
+
+static SDL_Surface *sdl_screen=NULL;
+
+static struct palette pal;
+
+static int debug=0;
+
+static int xsize=320,ysize=200,scale=1;
+
+int mode13h_graphics_init(char *name, int new_scale) {
+
+	int mode;
+
+	scale=new_scale;
+	xsize=320*scale;
+	ysize=200*scale;
+
+	mode=SDL_SWSURFACE|SDL_HWPALETTE|SDL_HWSURFACE;
+
+	if ( SDL_Init(SDL_INIT_VIDEO) < 0 ) {
+		fprintf(stderr,
+			"Couldn't initialize SDL: %s\n", SDL_GetError());
+		return -1;
+	}
+
+	/* Clean up on exit */
+	atexit(SDL_Quit);
+
+	/* assume 32-bit color */
+	sdl_screen = SDL_SetVideoMode(xsize, ysize, 32, mode);
+
+	if ( sdl_screen == NULL ) {
+		fprintf(stderr, "ERROR!  Couldn't set %dx%d video mode: %s\n",
+			xsize,ysize,SDL_GetError());
+		return -1;
+	}
+	SDL_EnableKeyRepeat(SDL_DEFAULT_REPEAT_DELAY, SDL_DEFAULT_REPEAT_INTERVAL);
+
+	if (name==NULL) {
+		SDL_WM_SetCaption("Linux/C/SDL","unknown");
+	}
+	else {
+		SDL_WM_SetCaption(name,name);
+	}
+
+	return 0;
+}
+
+int mode13h_graphics_update(void) {
+
+        unsigned int *t_pointer;
+
+        int x,y,temp,in_ptr,out_ptr,i,j;
+
+        /* point to SDL output pixels */
+        t_pointer=((Uint32 *)sdl_screen->pixels);
+
+	out_ptr=0;
+	for(y=0;y<200;y++) {
+		for(j=0;j<scale;j++) {
+	        for(x=0;x<320;x++) {
+			in_ptr=(y*320)+x;
+	                temp=(pal.red[framebuffer[in_ptr]]<<16)|
+        	                (pal.green[framebuffer[in_ptr]]<<8)|
+                	        (pal.blue[framebuffer[in_ptr]]<<0)|0;
+			for(i=0;i<scale;i++) {
+				t_pointer[out_ptr]=temp;
+				out_ptr++;
+			}
+		}
+		}
+        }
+
+        SDL_UpdateRect(sdl_screen, 0, 0, xsize, ysize);
+
+        return 0;
+}
+
+
+/* output of vgapal https://github.com/canidlogic/vgapal */
+static unsigned char default_pal[3*256]=
+{ 0,   0,   0,    0,   0, 170,    0, 170,   0,    0, 170, 170,
+170,   0,   0,  170,   0, 170,  170,  85,   0,  170, 170, 170,
+ 85,  85,  85,   85,  85, 255,   85, 255,  85,   85, 255, 255,
+255,  85,  85,  255,  85, 255,  255, 255,  85,  255, 255, 255,
+  0,   0,   0,   20,  20,  20,   32,  32,  32,   44,  44,  44,
+ 56,  56,  56,   69,  69,  69,   81,  81,  81,   97,  97,  97,
+113, 113, 113,  130, 130, 130,  146, 146, 146,  162, 162, 162,
+182, 182, 182,  203, 203, 203,  227, 227, 227,  255, 255, 255,
+  0,   0, 255,   65,   0, 255,  125,   0, 255,  190,   0, 255,
+255,   0, 255,  255,   0, 190,  255,   0, 125,  255,   0,  65,
+255,   0,   0,  255,  65,   0,  255, 125,   0,  255, 190,   0,
+255, 255,   0,  190, 255,   0,  125, 255,   0,   65, 255,   0,
+  0, 255,   0,    0, 255,  65,    0, 255, 125,    0, 255, 190,
+  0, 255, 255,    0, 190, 255,    0, 125, 255,    0,  65, 255,
+125, 125, 255,  158, 125, 255,  190, 125, 255,  223, 125, 255,
+255, 125, 255,  255, 125, 223,  255, 125, 190,  255, 125, 158,
+255, 125, 125,  255, 158, 125,  255, 190, 125,  255, 223, 125,
+255, 255, 125,  223, 255, 125,  190, 255, 125,  158, 255, 125,
+125, 255, 125,  125, 255, 158,  125, 255, 190,  125, 255, 223,
+125, 255, 255,  125, 223, 255,  125, 190, 255,  125, 158, 255,
+182, 182, 255,  199, 182, 255,  219, 182, 255,  235, 182, 255,
+255, 182, 255,  255, 182, 235,  255, 182, 219,  255, 182, 199,
+255, 182, 182,  255, 199, 182,  255, 219, 182,  255, 235, 182,
+255, 255, 182,  235, 255, 182,  219, 255, 182,  199, 255, 182,
+182, 255, 182,  182, 255, 199,  182, 255, 219,  182, 255, 235,
+182, 255, 255,  182, 235, 255,  182, 219, 255,  182, 199, 255,
+  0,   0, 113,   28,   0, 113,   56,   0, 113,   85,   0, 113,
+113,   0, 113,  113,   0,  85,  113,   0,  56,  113,   0,  28,
+113,   0,   0,  113,  28,   0,  113,  56,   0,  113,  85,   0,
+113, 113,   0,   85, 113,   0,   56, 113,   0,   28, 113,   0,
+  0, 113,   0,    0, 113,  28,    0, 113,  56,    0, 113,  85,
+  0, 113, 113,    0,  85, 113,    0,  56, 113,    0,  28, 113,
+ 56,  56, 113,   69,  56, 113,   85,  56, 113,   97,  56, 113,
+113,  56, 113,  113,  56,  97,  113,  56,  85,  113,  56,  69,
+113,  56,  56,  113,  69,  56,  113,  85,  56,  113,  97,  56,
+113, 113,  56,   97, 113,  56,   85, 113,  56,   69, 113,  56,
+ 56, 113,  56,   56, 113,  69,   56, 113,  85,   56, 113,  97,
+ 56, 113, 113,   56,  97, 113,   56,  85, 113,   56,  69, 113,
+ 81,  81, 113,   89,  81, 113,   97,  81, 113,  105,  81, 113,
+113,  81, 113,  113,  81, 105,  113,  81,  97,  113,  81,  89,
+113,  81,  81,  113,  89,  81,  113,  97,  81,  113, 105,  81,
+113, 113,  81,  105, 113,  81,   97, 113,  81,   89, 113,  81,
+ 81, 113,  81,   81, 113,  89,   81, 113,  97,   81, 113, 105,
+ 81, 113, 113,   81, 105, 113,   81,  97, 113,   81,  89, 113,
+  0,   0,  65,   16,   0,  65,   32,   0,  65,   48,   0,  65,
+ 65,   0,  65,   65,   0,  48,   65,   0,  32,   65,   0,  16,
+ 65,   0,   0,   65,  16,   0,   65,  32,   0,   65,  48,   0,
+ 65,  65,   0,   48,  65,   0,   32,  65,   0,   16,  65,   0,
+  0,  65,   0,    0,  65,  16,    0,  65,  32,    0,  65,  48,
+  0,  65,  65,    0,  48,  65,    0,  32,  65,    0,  16,  65,
+ 32,  32,  65,   40,  32,  65,   48,  32,  65,   56,  32,  65,
+ 65,  32,  65,   65,  32,  56,   65,  32,  48,   65,  32,  40,
+ 65,  32,  32,   65,  40,  32,   65,  48,  32,   65,  56,  32,
+ 65,  65,  32,   56,  65,  32,   48,  65,  32,   40,  65,  32,
+ 32,  65,  32,   32,  65,  40,   32,  65,  48,   32,  65,  56,
+ 32,  65,  65,   32,  56,  65,   32,  48,  65,   32,  40,  65,
+ 44,  44,  65,   48,  44,  65,   52,  44,  65,   60,  44,  65,
+ 65,  44,  65,   65,  44,  60,   65,  44,  52,   65,  44,  48,
+ 65,  44,  44,   65,  48,  44,   65,  52,  44,   65,  60,  44,
+ 65,  65,  44,   60,  65,  44,   52,  65,  44,   48,  65,  44,
+ 44,  65,  44,   44,  65,  48,   44,  65,  52,   44,  65,  60,
+ 44,  65,  65,   44,  60,  65,   44,  52,  65,   44,  48,  65,
+  0,   0,   0,    0,   0,   0,    0,   0,   0,    0,   0,   0,
+  0,   0,   0,    0,   0,   0,    0,   0,   0,    0,   0,   0
+};
+
+void set_default_pal(void) {
+
+	int i;
+
+	for(i=0;i<256;i++) {
+		pal.red[i]=default_pal[i*3];
+		pal.green[i]=default_pal[(i*3)+1];
+		pal.blue[i]=default_pal[(i*3)+2];
+
+	}
+	return;
+
+}
+
+int graphics_input(void) {
+
+	SDL_Event event;
+	int keypressed;
+
+	while ( SDL_PollEvent(&event)) {
+
+		switch(event.type) {
+
+		case SDL_KEYDOWN:
+			keypressed=event.key.keysym.sym;
+			switch (keypressed) {
+
+			case SDLK_ESCAPE:
+				return 27;
+
+			default:
+				return keypressed;
+			}
+			break;
+		}
+	}
+	return 0;
+
+}
+
+void framebuffer_write_20bit(int address, int value) {
+	int real_addr;
+
+	real_addr=address-0xa0000;
+	if ((real_addr<0) || (real_addr>320*200)) return;
+
+	framebuffer[real_addr]=value;
+
+}
+
+void framebuffer_write(int address, int value) {
+
+	if (debug) {
+	if ((address<0) || (address>320*200)) {
+		fprintf(stderr,
+			"Error!  Framebuffer write of 0x%x (%d,%d) "
+			"out of bounds\n",
+			address,address%320,address/320);
+		return;
+	}
+	}
+
+	framebuffer[address]=value;
+
+}
+
+
+
+
+static int dac_write_address=0,dac_write_which=0;
+static int dac_read_address=0,dac_read_which=0;
+
+void outp(short address, int value) {
+
+	switch(address) {
+
+	/* DAC read address */
+	case 0x3c7:
+		if (debug) fprintf(stderr,"Setting read address to %d\n",value);
+		dac_read_address=value;
+		dac_read_which=0;
+		break;
+
+	/* DAC write address */
+	case 0x3c8:
+		if (debug) fprintf(stderr,"Setting write address to %d\n",value);
+		dac_write_address=value;
+		dac_write_which=0;
+		break;
+
+	/* PEL/DAC data */
+	/* Note colors are 0..63 */
+	case 0x3c9:
+		if (dac_write_which==0) {
+			pal.red[dac_write_address]=value*4;
+			if (debug) fprintf(stderr,"Color %d R=0x%x ",
+				dac_write_address,value);
+		}
+		if (dac_write_which==1) {
+			pal.green[dac_write_address]=value*4;
+			if (debug) fprintf(stderr,"G=0x%x ",value);
+		}
+		if (dac_write_which==2) {
+			pal.blue[dac_write_address]=value*4;
+			if (debug) fprintf(stderr,"B=0x%x\n",value);
+		}
+		dac_write_which++;
+		if (dac_write_which==3) {
+			dac_write_address++;
+			dac_write_which=0;
+			if (dac_write_address>255) {
+				if (debug) fprintf(stderr,"Palette overflow!\n");
+				dac_write_address=0;	/* FIXME: is this right? */
+			}
+		}
+		break;
+
+
+	default:
+		printf("outp to unknown port 0x%x\n",address);
+	}
+}
+
+int inp(short address) {
+
+	int retval=0;
+
+	switch(address) {
+
+	/* DAC read address */
+	case 0x3c9:
+		if (dac_read_which==0) {
+			retval=pal.red[dac_read_address]/4;
+			if (debug) fprintf(stderr,"Color %d R=0x%x ",
+				dac_read_address,retval);
+		}
+		if (dac_read_which==1) {
+			retval=pal.green[dac_read_address]/4;
+			if (debug) fprintf(stderr,"G=0x%x ",retval);
+		}
+		if (dac_read_which==2) {
+			retval=pal.blue[dac_read_address]/4;
+			if (debug) fprintf(stderr,"B=0x%x\n",retval);
+		}
+		dac_read_which++;
+		if (dac_read_which==3) {
+			dac_read_address++;
+			dac_read_which=0;
+			if (dac_read_address>255) {
+				if (debug) fprintf(stderr,"Palette overflow!\n");
+				dac_read_address=0;	/* FIXME: is this right? */
+			}
+		}
+		break;
+
+
+		default:
+			printf("inp from unknown port 0x%x\n",address);
+	}
+
+        return retval;
+}
diff --git a/utils/gr-sim/dots/vga_emulator.h b/utils/gr-sim/dots/vga_emulator.h
new file mode 100644
index 00000000..2e06468a
--- /dev/null
+++ b/utils/gr-sim/dots/vga_emulator.h
@@ -0,0 +1,18 @@
+extern unsigned char framebuffer[65536];
+
+struct palette {
+        unsigned char red[256];
+        unsigned char green[256];
+        unsigned char blue[256];
+};
+
+int mode13h_graphics_init(char *name, int scale);
+int mode13h_graphics_update(void);
+void set_default_pal(void);
+int graphics_input(void);
+
+void framebuffer_write_20bit(int address, int value);
+void framebuffer_write(int address, int value);
+
+void outp(short address, int value);
+int inp(short addr);