in cdDriv the following:
if( vary.lgVaryOn )
{
	if( lgTraceInput )
		fprintf( ioQQQ, "cdDriv: calling optimize_do\n");
	/* option to drive optimizer set if OPTIMIZE was in input stream */
	optimize_do();
}
else
{
	if( lgTraceInput )
		fprintf( ioQQQ, "cdDriv: calling cloudy\n");
	/* optimize did not occur, only compute one model, call cloudy */
	cloudy();
}

optimize_do chooses which optimization routine to call
