Spherical Terrain in T3D

for ( S32 y = 0; y < blockSize; y++ )
{
	for ( S32 x = 0; x < blockSize; x++ )
	{

		F32 x1 = x+1;
		F32 y1 = y+1;

		x1 /= 256;
		y1 /= 256;

		height =  sqrt((F32)1 - ((x1*x1)/2) - ((y1*y1)/2) + (((x1*x1)*(y1*y1))/3) );


		F32 z1 = height;

		widthY =	sqrt((F32)1 - ((z1*z1)/2) - ((x1*x1)/2) + (((z1*z1)*(x1*x1))/3) );//x*2;
		lengthX =	sqrt((F32)1 - ((y1*y1)/2) - ((z1*z1)/2) + (((y1*y1)*(z1*z1))/3) );//x*2;

		lengthX *= 256;

		widthY *= 256;;

		height *= 256;

		Con::printf("setterraindata: %d %d %f %d" , x , y ,height, floatToFixed( height ));

		file->setWidthY( x, y, floatToFixed( widthY ) );
		file->setLengthX( x, y, floatToFixed( lengthX ) );

		file->setHeight( x, y, floatToFixed( height ) );


		count++;

	}
}

terrain->updateGrid( Point2I::Zero, Point2I( blockSize, blockSize ) );
terrain->updateGridMaterials( Point2I::Zero, Point2I( blockSize, blockSize ) );
terrain->registerObject( terrainName );

SimGroup *missionGroup;
if( Sim::findObject( "MissionGroup", missionGroup ) ){
	missionGroup->addObject( terrain );
}

inline void TerrainFile::setWidthY( U32 x, U32 y, U16 width )
{
   x %= mSize;
   y %= mSize;
   mWidthMap[ x + ( y * mSize ) ] = width;
}
inline void TerrainFile::setLengthX( U32 x, U32 y, U16 length )
{
   x %= mSize;
   y %= mSize;
   mLengthMap[ x + ( y * mSize ) ] = length;
}
inline void TerrainFile::setHeight( U32 x, U32 y, U16 height )
{
   x %= mSize;
   y %= mSize;
   mHeightMap[ x + ( y * mSize ) ] = height;
}

for ( U32 y = 0; y < smVBStride; y++ )
{
	for ( U32 x = 0; x < smVBStride; x++ )
	{

		// We clamp here to keep the geometry from reading across
		// one side of the height map to the other causing walls
		// around the edges of the terrain.
		//   terpoint.x = mClamp( mPoint.x + x * stepSize, 0, blockSize - 1 );
		//    terpoint.y = mClamp( mPoint.y + y * stepSize, 0, blockSize - 1 );
		gridPt.x = mClamp( mPoint.x + x * stepSize, 0, blockSize - 1 );
		gridPt.y = mClamp( mPoint.y + y * stepSize, 0, blockSize - 1 );

		// Setup this point.
		point.x = (F32)gridPt.x * squareSize;
		point.y = (F32)gridPt.y * squareSize;

		lengthX = fixedToFloat( file->getLengthX( gridPt.x, gridPt.y ) );
		widthY = fixedToFloat( file->getWidthY( gridPt.x, gridPt.y ) );
		height = fixedToFloat( file->getHeight( gridPt.x, gridPt.y ) );

		Point3F terdata = file->getterdata(count);

		vert->point.x = lengthX;// point.x;//terdata.x;//point.x*2;
		vert->point.y = widthY;//point.y;//terdata.y;//point.y;
		vert->point.z = height;//50.0755;//terdata.z;//height;

		count++;

		Con::printf("updatevertexbuffer: %f %f %f" ,widthY, lengthX,height);
		// Get the normal.
		mTerrain->getSmoothNormal( point, &normal, true, false );
		vert->normal = normal;

		// Get the tangent z.
		vert->tangentZ = fixedToFloat( file->getHeight( gridPt.x + 1, gridPt.y ) ) - height;

		// Test the empty state for this vert.
		if ( file->isEmptyAt( gridPt.x, gridPt.y ) )
		{
			mHasEmpty = true;
			mEmptyVertexList.push_back( vbcounter );
		}

		vbcounter++;
		++vert;
	}
}