diff --git a/src/osgTerrain/DataSet.cpp b/src/osgTerrain/DataSet.cpp
index ff595def4..a2ff1c8ed 100644
--- a/src/osgTerrain/DataSet.cpp
+++ b/src/osgTerrain/DataSet.cpp
@@ -64,14 +64,14 @@ CoordinateSystemType getCoordinateSystemType(const osg::CoordinateSystemNode* lh
     
      
     
-    std::cout<<"getCoordinateSystemType("<<projection_string<<")"<<std::endl;
-    std::cout<<"    lhsSR.IsGeographic()="<<lhsSR.IsGeographic()<<std::endl;
-    std::cout<<"    lhsSR.IsProjected()="<<lhsSR.IsProjected()<<std::endl;
-    std::cout<<"    lhsSR.IsLocal()="<<lhsSR.IsLocal()<<std::endl;
+    osg::notify(osg::INFO)<<"getCoordinateSystemType("<<projection_string<<")"<<std::endl;
+    osg::notify(osg::INFO)<<"    lhsSR.IsGeographic()="<<lhsSR.IsGeographic()<<std::endl;
+    osg::notify(osg::INFO)<<"    lhsSR.IsProjected()="<<lhsSR.IsProjected()<<std::endl;
+    osg::notify(osg::INFO)<<"    lhsSR.IsLocal()="<<lhsSR.IsLocal()<<std::endl;
 
     free(projection_string);
 
-    if (strcmp(lhsSR.GetRoot()->GetValue(),"GEOCCS")==0) std::cout<<"    lhsSR. is GEOCENTRIC "<<std::endl;
+    if (strcmp(lhsSR.GetRoot()->GetValue(),"GEOCCS")==0) osg::notify(osg::INFO)<<"    lhsSR. is GEOCENTRIC "<<std::endl;
     
 
     if (strcmp(lhsSR.GetRoot()->GetValue(),"GEOCCS")==0) return GEOCENTRIC;    
@@ -94,7 +94,7 @@ double getAngularUnits(const osg::CoordinateSystemNode* lhs)
 
     char* str;
     double result = lhsSR.GetAngularUnits(&str);
-    std::cout<<"lhsSR.GetAngularUnits("<<str<<") "<<result<<std::endl;
+    osg::notify(osg::INFO)<<"lhsSR.GetAngularUnits("<<str<<") "<<result<<std::endl;
 
     return result;
 }
@@ -112,11 +112,11 @@ double getLinearUnits(const osg::CoordinateSystemNode* lhs)
 
     char* str;
     double result = lhsSR.GetLinearUnits(&str);
-    std::cout<<"lhsSR.GetLinearUnits("<<str<<") "<<result<<std::endl;
+    osg::notify(osg::INFO)<<"lhsSR.GetLinearUnits("<<str<<") "<<result<<std::endl;
 
-    std::cout<<"lhsSR.IsGeographic() "<<lhsSR.IsGeographic()<<std::endl;
-    std::cout<<"lhsSR.IsProjected() "<<lhsSR.IsProjected()<<std::endl;
-    std::cout<<"lhsSR.IsLocal() "<<lhsSR.IsLocal()<<std::endl;
+    osg::notify(osg::INFO)<<"lhsSR.IsGeographic() "<<lhsSR.IsGeographic()<<std::endl;
+    osg::notify(osg::INFO)<<"lhsSR.IsProjected() "<<lhsSR.IsProjected()<<std::endl;
+    osg::notify(osg::INFO)<<"lhsSR.IsLocal() "<<lhsSR.IsLocal()<<std::endl;
     
     return result;
 }
@@ -129,10 +129,12 @@ bool areCoordinateSystemEquivalent(const osg::CoordinateSystemNode* lhs,const os
     // if one CS is NULL then true false
     if (!lhs || !rhs)
     {
-        //std::cout<<"areCoordinateSystemEquivalent lhs="<<lhs<<"  rhs="<<rhs<<std::endl;
+        osg::notify(osg::INFO)<<"areCoordinateSystemEquivalent lhs="<<lhs<<"  rhs="<<rhs<<" return true"<<std::endl;
         return false;
     }
     
+    osg::notify(osg::INFO)<<"areCoordinateSystemEquivalent lhs="<<lhs->getCoordinateSystem()<<"  rhs="<<rhs->getCoordinateSystem()<<std::endl;
+
     // use compare on ProjectionRef strings.
     if (lhs->getCoordinateSystem() == rhs->getCoordinateSystem()) return true;
     
@@ -159,7 +161,7 @@ bool areCoordinateSystemEquivalent(const osg::CoordinateSystemNode* lhs,const os
 #if 0
     int result2 = lhsSR.IsSameGeogCS(&rhsSR);
 
-     std::cout<<"areCoordinateSystemEquivalent "<<std::endl
+     osg::notify(osg::INFO)<<"areCoordinateSystemEquivalent "<<std::endl
               <<"LHS = "<<lhs->getCoordinateSystem()<<std::endl
               <<"RHS = "<<rhs->getCoordinateSystem()<<std::endl
               <<"result = "<<result<<"  result2 = "<<result2<<std::endl;
@@ -206,7 +208,7 @@ DataSet::SourceData* DataSet::SourceData::readData(Source* source)
                 }
                 else if (gdalDataSet->GetGCPCount()>0 && gdalDataSet->GetGCPProjection())
                 {
-                    std::cout << "    Using GCP's"<< std::endl;
+                    osg::notify(osg::INFO) << "    Using GCP's"<< std::endl;
 
 
                     /* -------------------------------------------------------------------- */
@@ -220,7 +222,7 @@ DataSet::SourceData* DataSet::SourceData::readData(Source* source)
 
                     if ( hTransformArg == NULL )
                     {
-                        std::cout<<" failed to create transformer"<<std::endl;
+                        osg::notify(osg::INFO)<<" failed to create transformer"<<std::endl;
                         return NULL;
                     }
 
@@ -234,7 +236,7 @@ DataSet::SourceData* DataSet::SourceData::readData(Source* source)
                                                  adfDstGeoTransform, &nPixels, &nLines )
                         != CE_None )
                     {
-                        std::cout<<" failed to create warp"<<std::endl;
+                        osg::notify(osg::INFO)<<" failed to create warp"<<std::endl;
                         return NULL;
                     }
 
@@ -251,7 +253,7 @@ DataSet::SourceData* DataSet::SourceData::readData(Source* source)
                 }
                 else
                 {
-                    std::cout << "    No GeoTransform or GCP's - unable to compute position in space"<< std::endl;
+                    osg::notify(osg::INFO) << "    No GeoTransform or GCP's - unable to compute position in space"<< std::endl;
                     
                     data->_geoTransform.set( 1.0,    0.0,    0.0,    0.0,
                                              0.0,    1.0,    0.0,    0.0,
@@ -306,7 +308,7 @@ const DataSet::SpatialProperties& DataSet::SourceData::computeSpatialProperties(
         if (_gdalDataSet)
         {
 
-            //std::cout<<"Projecting bounding volume for "<<_source->getFileName()<<std::endl;
+            //osg::notify(osg::INFO)<<"Projecting bounding volume for "<<_source->getFileName()<<std::endl;
 
             
             // insert into the _spatialPropertiesMap for future reuse.
@@ -324,7 +326,7 @@ const DataSet::SpatialProperties& DataSet::SourceData::computeSpatialProperties(
 
             if (!hTransformArg)
             {
-                std::cout<<" failed to create transformer"<<std::endl;
+                osg::notify(osg::INFO)<<" failed to create transformer"<<std::endl;
                 return sp;
             }
         
@@ -335,7 +337,7 @@ const DataSet::SpatialProperties& DataSet::SourceData::computeSpatialProperties(
                                          adfDstGeoTransform, &nPixels, &nLines )
                 != CE_None )
             {
-                std::cout<<" failed to create warp"<<std::endl;
+                osg::notify(osg::INFO)<<" failed to create warp"<<std::endl;
                 return sp;
             }
 
@@ -355,7 +357,7 @@ const DataSet::SpatialProperties& DataSet::SourceData::computeSpatialProperties(
         }
 
     }
-    std::cout<<"DataSet::DataSource::assuming compatible coordinates."<<std::endl;
+    osg::notify(osg::INFO)<<"DataSet::DataSource::assuming compatible coordinates."<<std::endl;
     return *this;
 }
 
@@ -366,28 +368,28 @@ bool DataSet::SourceData::intersects(const SpatialProperties& sp) const
 
 void DataSet::SourceData::read(DestinationData& destination)
 {
-    std::cout<<"A"<<std::endl;
+    osg::notify(osg::INFO)<<"A"<<std::endl;
 
     if (!_source) return;
     
-    std::cout<<"B"<<std::endl;
+    osg::notify(osg::INFO)<<"B"<<std::endl;
 
     switch (_source->getType())
     {
     case(Source::IMAGE):
-        std::cout<<"B.1"<<std::endl;
+        osg::notify(osg::INFO)<<"B.1"<<std::endl;
         readImage(destination);
         break;
     case(Source::HEIGHT_FIELD):
-        std::cout<<"B.2"<<std::endl;
+        osg::notify(osg::INFO)<<"B.2"<<std::endl;
         readHeightField(destination);
         break;
     case(Source::MODEL):
-        std::cout<<"B.3"<<std::endl;
+        osg::notify(osg::INFO)<<"B.3"<<std::endl;
         readModels(destination);
         break;
     }
-    std::cout<<"C"<<std::endl;
+    osg::notify(osg::INFO)<<"C"<<std::endl;
 }
 
 void DataSet::SourceData::readImage(DestinationData& destination)
@@ -401,7 +403,7 @@ void DataSet::SourceData::readImage(DestinationData& destination)
         osg::BoundingBox intersect_bb(s_bb.intersect(d_bb));
         if (!intersect_bb.valid())
         {
-            std::cout<<"Reading image but it does not intesection destination - ignoring"<<std::endl;
+            osg::notify(osg::INFO)<<"Reading image but it does not intesection destination - ignoring"<<std::endl;
             return;
         }
         
@@ -416,8 +418,8 @@ void DataSet::SourceData::readImage(DestinationData& destination)
        int destWidth = osg::minimum((int)ceilf((float)destination._image->s()*(intersect_bb.xMax()-d_bb.xMin())/(d_bb.xMax()-d_bb.xMin())),(int)destination._image->s())-destX;
        int destHeight = osg::minimum((int)ceilf((float)destination._image->t()*(intersect_bb.yMax()-d_bb.yMin())/(d_bb.yMax()-d_bb.yMin())),(int)destination._image->t())-destY;
 
-        std::cout<<"   copying from "<<windowX<<"\t"<<windowY<<"\t"<<windowWidth<<"\t"<<windowHeight<<std::endl;
-        std::cout<<"             to "<<destX<<"\t"<<destY<<"\t"<<destWidth<<"\t"<<destHeight<<std::endl;
+        osg::notify(osg::INFO)<<"   copying from "<<windowX<<"\t"<<windowY<<"\t"<<windowWidth<<"\t"<<windowHeight<<std::endl;
+        osg::notify(osg::INFO)<<"             to "<<destX<<"\t"<<destY<<"\t"<<destWidth<<"\t"<<destHeight<<std::endl;
         
         // which band do we want to read from...        
         int numBands = _gdalDataSet->GetRasterCount();
@@ -468,7 +470,7 @@ void DataSet::SourceData::readImage(DestinationData& destination)
             int lineSpace=-(int)(destination._image->getRowSizeInBytes());
 
             unsigned char* imageData = destination._image->data(destX,destY+destHeight-1);
-            std::cout << "reading RGB"<<std::endl;
+            osg::notify(osg::INFO) << "reading RGB"<<std::endl;
 
             bool directCopy = false;
             if (directCopy)
@@ -531,7 +533,7 @@ void DataSet::SourceData::readImage(DestinationData& destination)
         }
         else
         {
-            std::cout<<"Warning image not read as Red, Blue and Green bands not present"<<std::endl;
+            osg::notify(osg::INFO)<<"Warning image not read as Red, Blue and Green bands not present"<<std::endl;
         }
 
 
@@ -540,10 +542,10 @@ void DataSet::SourceData::readImage(DestinationData& destination)
 
 void DataSet::SourceData::readHeightField(DestinationData& destination)
 {
-    std::cout<<"In DataSet::SourceData::readHeightField"<<std::endl;
+    osg::notify(osg::INFO)<<"In DataSet::SourceData::readHeightField"<<std::endl;
     if (destination._heightField.valid())
     {
-        std::cout<<"Reading height field"<<std::endl;
+        osg::notify(osg::INFO)<<"Reading height field"<<std::endl;
 
         osg::BoundingBox s_bb = getExtents(destination._cs.get());
         
@@ -553,7 +555,7 @@ void DataSet::SourceData::readHeightField(DestinationData& destination)
 
         if (!intersect_bb.valid())
         {
-            std::cout<<"Reading height field but it does not intesection destination - ignoring"<<std::endl;
+            osg::notify(osg::INFO)<<"Reading height field but it does not intesection destination - ignoring"<<std::endl;
             return;
         }
         
@@ -568,8 +570,8 @@ void DataSet::SourceData::readHeightField(DestinationData& destination)
        int destWidth = osg::minimum((int)ceilf((float)destination._heightField->getNumColumns()*(intersect_bb.xMax()-d_bb.xMin())/(d_bb.xMax()-d_bb.xMin())),(int)destination._heightField->getNumColumns())-destX;
        int destHeight = osg::minimum((int)ceilf((float)destination._heightField->getNumRows()*(intersect_bb.yMax()-d_bb.yMin())/(d_bb.yMax()-d_bb.yMin())),(int)destination._heightField->getNumRows())-destY;
 
-        std::cout<<"   copying from "<<windowX<<"\t"<<windowY<<"\t"<<windowWidth<<"\t"<<windowHeight<<std::endl;
-        std::cout<<"             to "<<destX<<"\t"<<destY<<"\t"<<destWidth<<"\t"<<destHeight<<std::endl;
+        osg::notify(osg::INFO)<<"   copying from "<<windowX<<"\t"<<windowY<<"\t"<<windowWidth<<"\t"<<windowHeight<<std::endl;
+        osg::notify(osg::INFO)<<"             to "<<destX<<"\t"<<destY<<"\t"<<destWidth<<"\t"<<destHeight<<std::endl;
 
 
         // which band do we want to read from...        
@@ -611,47 +613,47 @@ void DataSet::SourceData::readHeightField(DestinationData& destination)
             }
 
 
-            if (bandSelected->GetUnitType()) std::cout << "bandSelected->GetUnitType()=" << bandSelected->GetUnitType()<<std::endl;
-            else std::cout << "bandSelected->GetUnitType()= null" <<std::endl;
+            if (bandSelected->GetUnitType()) osg::notify(osg::INFO) << "bandSelected->GetUnitType()=" << bandSelected->GetUnitType()<<std::endl;
+            else osg::notify(osg::INFO) << "bandSelected->GetUnitType()= null" <<std::endl;
             
 
             int success = 0;
             float noDataValue = bandSelected->GetNoDataValue(&success);
             if (success)
             {
-                std::cout<<"We have NoDataValue = "<<noDataValue<<std::endl;
+                osg::notify(osg::INFO)<<"We have NoDataValue = "<<noDataValue<<std::endl;
             }
             else
             {
-                std::cout<<"We have no NoDataValue"<<std::endl;
+                osg::notify(osg::INFO)<<"We have no NoDataValue"<<std::endl;
                 noDataValue = 0.0f;
             }
 
             float offset = bandSelected->GetOffset(&success);
             if (success)
             {
-                std::cout<<"We have Offset = "<<offset<<std::endl;
+                osg::notify(osg::INFO)<<"We have Offset = "<<offset<<std::endl;
             }
             else
             {
-                std::cout<<"We have no Offset"<<std::endl;
+                osg::notify(osg::INFO)<<"We have no Offset"<<std::endl;
                 offset = 0.0f;
             }
 
             float scale = bandSelected->GetScale(&success);
             if (success)
             {
-                std::cout<<"We have Scale = "<<scale<<std::endl;
+                osg::notify(osg::INFO)<<"We have Scale = "<<scale<<std::endl;
             }
             else
             {
                 scale = destination._dataSet->getVerticalScale();
-                std::cout<<"We have no Scale from file so use DataSet vertical scale of "<<scale<<std::endl;
+                osg::notify(osg::INFO)<<"We have no Scale from file so use DataSet vertical scale of "<<scale<<std::endl;
                 //scale = (xyInDegrees /*&& !destination._dataSet->getConvertFromGeographicToGeocentric()*/) ? 1.0f/111319.0f : 1.0f;
 
             }
             
-            std::cout<<"********* getLinearUnits = "<<getLinearUnits(_cs.get())<<std::endl;
+            osg::notify(osg::INFO)<<"********* getLinearUnits = "<<getLinearUnits(_cs.get())<<std::endl;
 
             // read data into temporary array
             float* heightData = new float [ destWidth*destHeight ];
@@ -659,11 +661,11 @@ void DataSet::SourceData::readHeightField(DestinationData& destination)
             // raad the data.
             osg::HeightField* hf = destination._heightField.get();
 
-            std::cout<<"reading height field"<<std::endl;
+            osg::notify(osg::INFO)<<"reading height field"<<std::endl;
             //bandSelected->RasterIO(GF_Read,windowX,_numValuesY-(windowY+windowHeight),windowWidth,windowHeight,floatdata,destWidth,destHeight,GDT_Float32,numBytesPerZvalue,lineSpace);
             bandSelected->RasterIO(GF_Read,windowX,_numValuesY-(windowY+windowHeight),windowWidth,windowHeight,heightData,destWidth,destHeight,GDT_Float32,0,0);
 
-            std::cout<<"    scaling height field"<<std::endl;
+            osg::notify(osg::INFO)<<"    scaling height field"<<std::endl;
 
             float noDataValueFill = 0.0f;
             bool ignoreNoDataValue = true;
@@ -693,7 +695,7 @@ void DataSet::SourceData::readModels(DestinationData& destination)
 {
     if (_model.valid())
     {
-        std::cout<<"Raading model"<<std::endl;
+        osg::notify(osg::INFO)<<"Raading model"<<std::endl;
         destination._models.push_back(_model);
     }
 }
@@ -711,7 +713,7 @@ void DataSet::Source::setSortValueFromSourceDataResolution()
 
 void DataSet::Source::loadSourceData()
 {
-    std::cout<<"DataSet::Source::loadSourceData() "<<_filename<<std::endl;
+    osg::notify(osg::INFO)<<"DataSet::Source::loadSourceData() "<<_filename<<std::endl;
     
     _sourceData = SourceData::readData(this);
     
@@ -724,20 +726,20 @@ void DataSet::Source::assignCoordinateSystemAndGeoTransformAccordingToParameterP
     {
         _sourceData->_cs = _cs;
         
-        std::cout<<"assigning CS from Source to Data."<<std::endl;
+        osg::notify(osg::INFO)<<"assigning CS from Source to Data."<<std::endl;
         
     }
     else
     {
         _cs = _sourceData->_cs;
-        std::cout<<"assigning CS from Data to Source."<<std::endl;
+        osg::notify(osg::INFO)<<"assigning CS from Data to Source."<<std::endl;
     }
     
     if (getGeoTransformPolicy()==PREFER_CONFIG_SETTINGS)
     {
         _sourceData->_geoTransform = _geoTransform;
 
-        std::cout<<"assigning GeoTransform from Source to Data."<<_geoTransform<<std::endl;
+        osg::notify(osg::INFO)<<"assigning GeoTransform from Source to Data."<<_geoTransform<<std::endl;
 
     }
     else if (getGeoTransformPolicy()==PREFER_CONFIG_SETTINGS_BUT_SCALE_BY_FILE_RESOLUTION)
@@ -757,13 +759,13 @@ void DataSet::Source::assignCoordinateSystemAndGeoTransformAccordingToParameterP
 
         _sourceData->_geoTransform = _geoTransform;
 
-        std::cout<<"assigning GeoTransform from Source to Data."<<_geoTransform<<std::endl;
+        osg::notify(osg::INFO)<<"assigning GeoTransform from Source to Data."<<_geoTransform<<std::endl;
 
     }
     else
     {
         _geoTransform = _sourceData->_geoTransform;
-        std::cout<<"assigning GeoTransform from Data to Source."<<_geoTransform<<std::endl;
+        osg::notify(osg::INFO)<<"assigning GeoTransform from Data to Source."<<_geoTransform<<std::endl;
     }
     
     _sourceData->computeExtents();
@@ -786,13 +788,13 @@ bool DataSet::Source::needReproject(const osg::CoordinateSystemNode* cs, double
     // always need to reproject imagery with GCP's.
     if (_sourceData->_hasGCPs)
     {
-        std::cout<<"Need to to reproject due to presence of GCP's"<<std::endl;
+        osg::notify(osg::INFO)<<"Need to to reproject due to presence of GCP's"<<std::endl;
         return true;
     }
 
     if (!areCoordinateSystemEquivalent(_cs.get(),cs))
     {
-        std::cout<<"Need to do reproject !areCoordinateSystemEquivalent(_cs.get(),cs)"<<std::endl;
+        osg::notify(osg::INFO)<<"Need to do reproject !areCoordinateSystemEquivalent(_cs.get(),cs)"<<std::endl;
 
         return true;
     }
@@ -816,19 +818,19 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
     if (!_sourceData) return 0;
     if (_type==MODEL) return 0;
     
-    std::cout<<"reprojecting to file "<<filename<<std::endl;
+    osg::notify(osg::INFO)<<"reprojecting to file "<<filename<<std::endl;
 
     GDALDriverH hDriver = GDALGetDriverByName( "GTiff" );
         
     if (hDriver == NULL)
     {       
-    std::cout<<"Unable to load driver for "<<"GTiff"<<std::endl;
+    osg::notify(osg::INFO)<<"Unable to load driver for "<<"GTiff"<<std::endl;
         return 0;
     }
     
     if (GDALGetMetadataItem( hDriver, GDAL_DCAP_CREATE, NULL ) == NULL )
     {
-        std::cout<<"GDAL driver does not support create for "<<osgDB::getFileExtension(filename)<<std::endl;
+        osg::notify(osg::INFO)<<"GDAL driver does not support create for "<<osgDB::getFileExtension(filename)<<std::endl;
         return 0;
     }
 
@@ -843,7 +845,7 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
 
     if (!hTransformArg)
     {
-        std::cout<<" failed to create transformer"<<std::endl;
+        osg::notify(osg::INFO)<<" failed to create transformer"<<std::endl;
         return 0;
     }
 
@@ -854,18 +856,18 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
                                  adfDstGeoTransform, &nPixels, &nLines )
         != CE_None )
     {
-        std::cout<<" failed to create warp"<<std::endl;
+        osg::notify(osg::INFO)<<" failed to create warp"<<std::endl;
         return 0;
     }
     
     if (targetResolution>0.0f)
     {
-        std::cout<<"recomputing the target transform size"<<std::endl;
+        osg::notify(osg::INFO)<<"recomputing the target transform size"<<std::endl;
         
         double currentResolution = sqrt(osg::square(adfDstGeoTransform[1])+osg::square(adfDstGeoTransform[2])+
                                         osg::square(adfDstGeoTransform[4])+osg::square(adfDstGeoTransform[5]));
 
-        std::cout<<"        default computed resolution "<<currentResolution<<" nPixels="<<nPixels<<" nLines="<<nLines<<std::endl;
+        osg::notify(osg::INFO)<<"        default computed resolution "<<currentResolution<<" nPixels="<<nPixels<<" nLines="<<nLines<<std::endl;
 
         double extentsPixels = sqrt(osg::square(adfDstGeoTransform[1])+osg::square(adfDstGeoTransform[2]))*(double)(nPixels-1);
         double extentsLines = sqrt(osg::square(adfDstGeoTransform[4])+osg::square(adfDstGeoTransform[5]))*(double)(nLines-1);
@@ -876,14 +878,14 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
         adfDstGeoTransform[4] *= ratio;
         adfDstGeoTransform[5] *= ratio;
         
-        std::cout<<"    extentsPixels="<<extentsPixels<<std::endl;
-        std::cout<<"    extentsLines="<<extentsLines<<std::endl;
-        std::cout<<"    targetResolution="<<targetResolution<<std::endl;
+        osg::notify(osg::INFO)<<"    extentsPixels="<<extentsPixels<<std::endl;
+        osg::notify(osg::INFO)<<"    extentsLines="<<extentsLines<<std::endl;
+        osg::notify(osg::INFO)<<"    targetResolution="<<targetResolution<<std::endl;
         
         nPixels = (int)ceil(extentsPixels/sqrt(osg::square(adfDstGeoTransform[1])+osg::square(adfDstGeoTransform[2])))+1;
         nLines = (int)ceil(extentsLines/sqrt(osg::square(adfDstGeoTransform[4])+osg::square(adfDstGeoTransform[5])))+1;
 
-        std::cout<<"        target computed resolution "<<targetResolution<<" nPixels="<<nPixels<<" nLines="<<nLines<<std::endl;
+        osg::notify(osg::INFO)<<"        target computed resolution "<<targetResolution<<" nPixels="<<nPixels<<" nLines="<<nLines<<std::endl;
         
     }
 
@@ -921,7 +923,7 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
     GDALTransformerFunc pfnTransformer = GDALGenImgProjTransform;
 
     
-    std::cout<<"Setting projection "<<cs->getCoordinateSystem()<<std::endl;
+    osg::notify(osg::INFO)<<"Setting projection "<<cs->getCoordinateSystem()<<std::endl;
 
 /* -------------------------------------------------------------------- */
 /*      Copy the color table, if required.                              */
@@ -994,7 +996,7 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
             double new_noDataValue = 0;
             if (success)
             {
-                std::cout<<"\tassinging no data value "<<noDataValue<<" to band "<<i+1<<std::endl;
+                osg::notify(osg::INFO)<<"\tassinging no data value "<<noDataValue<<" to band "<<i+1<<std::endl;
             
                 psWO->padfSrcNoDataReal[i] = noDataValue;
                 psWO->padfSrcNoDataImag[i] = 0.0;
@@ -1028,7 +1030,7 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
             double new_noDataValue = 0.0;
             if (success)
             {
-                std::cout<<"\tassinging no data value "<<noDataValue<<" to band "<<i+1<<std::endl;
+                osg::notify(osg::INFO)<<"\tassinging no data value "<<noDataValue<<" to band "<<i+1<<std::endl;
             
                 psWO->padfSrcNoDataReal[i] = noDataValue;
                 psWO->padfSrcNoDataImag[i] = 0.0;
@@ -1069,7 +1071,7 @@ DataSet::Source* DataSet::Source::doReproject(const std::string& filename, osg::
                                GDALGetRasterYSize( hDstDS ) );
     }
 
-    std::cout<<"new projection is "<<GDALGetProjectionRef(hDstDS)<<std::endl;
+    osg::notify(osg::INFO)<<"new projection is "<<GDALGetProjectionRef(hDstDS)<<std::endl;
 
 /* -------------------------------------------------------------------- */
 /*      Cleanup.                                                        */
@@ -1184,15 +1186,15 @@ void DataSet::DestinationTile::computeMaximumSourceResolution(CompositeSource* s
 {
     for(CompositeSource::source_iterator itr(sourceGraph);itr.valid();++itr)
     {
-    
-        if ((*itr)->getMaxLevel()<_level)
+        Source* source = itr->get();
+        if (!source || source->getMaxLevel()<_level)
         {
             // skip the contribution of this source since this destination tile exceeds its contribution level.
             continue;
         }
 
-        SourceData* data = (*itr)->getSourceData();
-        if (data && (*itr)->getType()!=Source::MODEL)
+        SourceData* data = source->getSourceData();
+        if (data && source->getType()!=Source::MODEL)
         {
 
             SpatialProperties sp = data->computeSpatialProperties(_cs.get());
@@ -1321,10 +1323,10 @@ void DataSet::DestinationTile::allocate()
         _terrain->_heightField->setYInterval(dem_dy);
 
         //float xMax = _terrain->_heightField->getOrigin().x()+_terrain->_heightField->getXInterval()*(float)(dem_numColumns-1);
-        //std::cout<<"ErrorX = "<<xMax-_extents.xMax()<<std::endl;
+        //osg::notify(osg::INFO)<<"ErrorX = "<<xMax-_extents.xMax()<<std::endl;
 
         //float yMax = _terrain->_heightField->getOrigin().y()+_terrain->_heightField->getYInterval()*(float)(dem_numRows-1);
-        //std::cout<<"ErrorY = "<<yMax-_extents.yMax()<<std::endl;
+        //osg::notify(osg::INFO)<<"ErrorY = "<<yMax-_extents.yMax()<<std::endl;
 
     }
 
@@ -1356,14 +1358,14 @@ void DataSet::DestinationTile::setNeighbours(DestinationTile* left, DestinationT
     _neighbour[ABOVE_LEFT] = above_left;
     
     
-//     std::cout<<"LEFT="<<_neighbour[LEFT]<<std::endl;
-//     std::cout<<"LEFT_BELOW="<<_neighbour[LEFT_BELOW]<<std::endl;
-//     std::cout<<"BELOW="<<_neighbour[BELOW]<<std::endl;
-//     std::cout<<"BELOW_RIGHT="<<_neighbour[BELOW_RIGHT]<<std::endl;
-//     std::cout<<"RIGHT="<<_neighbour[RIGHT]<<std::endl;
-//     std::cout<<"RIGHT_ABOVE="<<_neighbour[RIGHT_ABOVE]<<std::endl;
-//     std::cout<<"ABOVE="<<_neighbour[ABOVE]<<std::endl;
-//     std::cout<<"ABOVE_LEFT="<<_neighbour[ABOVE_LEFT]<<std::endl;
+//     osg::notify(osg::INFO)<<"LEFT="<<_neighbour[LEFT]<<std::endl;
+//     osg::notify(osg::INFO)<<"LEFT_BELOW="<<_neighbour[LEFT_BELOW]<<std::endl;
+//     osg::notify(osg::INFO)<<"BELOW="<<_neighbour[BELOW]<<std::endl;
+//     osg::notify(osg::INFO)<<"BELOW_RIGHT="<<_neighbour[BELOW_RIGHT]<<std::endl;
+//     osg::notify(osg::INFO)<<"RIGHT="<<_neighbour[RIGHT]<<std::endl;
+//     osg::notify(osg::INFO)<<"RIGHT_ABOVE="<<_neighbour[RIGHT_ABOVE]<<std::endl;
+//     osg::notify(osg::INFO)<<"ABOVE="<<_neighbour[ABOVE]<<std::endl;
+//     osg::notify(osg::INFO)<<"ABOVE_LEFT="<<_neighbour[ABOVE_LEFT]<<std::endl;
     
     for(int i=0;i<NUMBER_OF_POSITIONS;++i)
     {
@@ -1377,7 +1379,7 @@ void DataSet::DestinationTile::checkNeighbouringTiles()
     {
         if (_neighbour[i] && _neighbour[i]->_neighbour[(i+4)%NUMBER_OF_POSITIONS]!=this)
         {
-            std::cout<<"Error:: Tile "<<this<<"'s _neighbour["<<i<<"] does not point back to it."<<std::endl;
+            osg::notify(osg::INFO)<<"Error:: Tile "<<this<<"'s _neighbour["<<i<<"] does not point back to it."<<std::endl;
         }
     }
 }
@@ -1587,7 +1589,7 @@ const char* edgeString(DataSet::DestinationTile::Position position)
 void DataSet::DestinationTile::setTileComplete(bool complete)
 {
     _complete = complete;
-    std::cout<<"setTileComplete("<<complete<<") for "<<_level<<"\t"<<_tileX<<"\t"<<_tileY<<std::endl;
+    osg::notify(osg::INFO)<<"setTileComplete("<<complete<<") for "<<_level<<"\t"<<_tileX<<"\t"<<_tileY<<std::endl;
 }
 
 void DataSet::DestinationTile::equalizeEdge(Position position)
@@ -1608,7 +1610,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
     osg::Image* image1 = _imagery.valid()?_imagery->_image.get() : 0;
     osg::Image* image2 = tile2->_imagery.valid()?tile2->_imagery->_image.get() : 0;
 
-    //std::cout<<"Equalizing edge "<<edgeString(position)<<" of \t"<<_level<<"\t"<<_tileX<<"\t"<<_tileY
+    //osg::notify(osg::INFO)<<"Equalizing edge "<<edgeString(position)<<" of \t"<<_level<<"\t"<<_tileX<<"\t"<<_tileY
     //         <<"  neighbour "<<tile2->_level<<"\t"<<tile2->_tileX<<"\t"<<tile2->_tileY<<std::endl;
              
              
@@ -1621,8 +1623,8 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
         image1->getDataType()==GL_UNSIGNED_BYTE)
     {
 
-        //std::cout<<"   Equalizing image1= "<<image1<<                         " with image2 = "<<image2<<std::endl;
-        //std::cout<<"              data1 = 0x"<<std::hex<<(int)image1->data()<<" with data2  = 0x"<<(int)image2->data()<<std::endl;
+        //osg::notify(osg::INFO)<<"   Equalizing image1= "<<image1<<                         " with image2 = "<<image2<<std::endl;
+        //osg::notify(osg::INFO)<<"              data1 = 0x"<<std::hex<<(int)image1->data()<<" with data2  = 0x"<<(int)image2->data()<<std::endl;
 
         unsigned char* data1 = 0;
         unsigned char* data2 = 0;
@@ -1638,7 +1640,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
             data2 = image2->data(image2->s()-1,1); // RIGHT hand side
             delta2 = image2->getRowSizeInBytes();
             num = (image1->t()==image2->t())?image2->t()-2:0; // note miss out corners.
-            //std::cout<<"       left "<<num<<std::endl;
+            //osg::notify(osg::INFO)<<"       left "<<num<<std::endl;
             break;
         case BELOW:
             data1 = image1->data(1,0); // BELOW hand side
@@ -1646,7 +1648,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
             data2 = image2->data(1,image2->t()-1); // ABOVE hand side
             delta2 = 3;
             num = (image1->s()==image2->s())?image2->s()-2:0; // note miss out corners.
-            //std::cout<<"       below "<<num<<std::endl;
+            //osg::notify(osg::INFO)<<"       below "<<num<<std::endl;
             break;
         case RIGHT:
             data1 = image1->data(image1->s()-1,1); // LEFT hand side
@@ -1654,7 +1656,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
             data2 = image2->data(0,1); // RIGHT hand side
             delta2 = image2->getRowSizeInBytes();
             num = (image1->t()==image2->t())?image2->t()-2:0; // note miss out corners.
-            //std::cout<<"       right "<<num<<std::endl;
+            //osg::notify(osg::INFO)<<"       right "<<num<<std::endl;
             break;
         case ABOVE:
             data1 = image1->data(1,image1->t()-1); // ABOVE hand side
@@ -1662,10 +1664,10 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
             data2 = image2->data(1,0); // BELOW hand side
             delta2 = 3;
             num = (image1->s()==image2->s())?image2->s()-2:0; // note miss out corners.
-            //std::cout<<"       above "<<num<<std::endl;
+            //osg::notify(osg::INFO)<<"       above "<<num<<std::endl;
             break;
         default :
-            //std::cout<<"       default "<<num<<std::endl;
+            //osg::notify(osg::INFO)<<"       default "<<num<<std::endl;
             break;
         }
 
@@ -1696,7 +1698,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
             data1 += delta1;
             data2 += delta2;
 
-            //std::cout<<"    equalizing colour "<<(int)data1<<"  "<<(int)data2<<std::endl;
+            //osg::notify(osg::INFO)<<"    equalizing colour "<<(int)data1<<"  "<<(int)data2<<std::endl;
 
         }
 
@@ -1707,7 +1709,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
 
     if (heightField1 && heightField2)
     {
-        //std::cout<<"   Equalizing heightfield"<<std::endl;
+        //osg::notify(osg::INFO)<<"   Equalizing heightfield"<<std::endl;
 
         float* data1 = 0;
         float* data2 = 0;
@@ -1753,7 +1755,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
         {
             float z = (*data1 + *data2)/2.0f;
 
-            //std::cout<<"    equalizing height"<<std::endl;
+            //osg::notify(osg::INFO)<<"    equalizing height"<<std::endl;
 
             *data1 = z;
             *data2 = z;
@@ -1768,7 +1770,7 @@ void DataSet::DestinationTile::equalizeEdge(Position position)
 
 void DataSet::DestinationTile::equalizeBoundaries()
 {
-    std::cout<<"DataSet::DestinationTile::equalizeBoundaries()"<<std::endl;
+    osg::notify(osg::INFO)<<"DataSet::DestinationTile::equalizeBoundaries()"<<std::endl;
 
     equalizeCorner(LEFT_BELOW);
     equalizeCorner(BELOW_RIGHT);
@@ -1803,7 +1805,7 @@ void DataSet::DestinationTile::optimizeResolution()
 
         if (minHeight==maxHeight)
         {
-            std::cout<<"******* We have a flat tile ******* "<<std::endl;
+            osg::notify(osg::INFO)<<"******* We have a flat tile ******* "<<std::endl;
 
             unsigned int minimumSize = 8;
 
@@ -1873,7 +1875,6 @@ osg::StateSet* DataSet::DestinationTile::createStateSet()
         break;
       case(DataSet::MIP_MAPPING_IMAGERY):
         {
-            osg::notify(osg::NOTICE)<<"Mip mapped imagery not currently supported, falling back to hardware mip mapping."<<std::endl;
             texture->setFilter(osg::Texture::MIN_FILTER,osg::Texture::LINEAR_MIPMAP_LINEAR);
             texture->setFilter(osg::Texture::MAG_FILTER,osg::Texture::LINEAR);
         }
@@ -1887,8 +1888,6 @@ osg::StateSet* DataSet::DestinationTile::createStateSet()
     bool compressedImageSupported = inlineImageFile;
     bool mipmapImageSupported = inlineImageFile;
 
-    std::cout<<"__________________image->getPixelFormat()="<<image->getPixelFormat()<<std::endl;
-
     if (compressedImageSupported && 
         _dataSet->getTextureType()==COMPRESSED_TEXTURE &&
         (image->getPixelFormat()==GL_RGB || image->getPixelFormat()==GL_RGBA))
@@ -1912,7 +1911,7 @@ osg::StateSet* DataSet::DestinationTile::createStateSet()
 
         texture->setInternalFormatMode(osg::Texture::USE_IMAGE_DATA_FORMAT);
 
-        std::cout<<">>>>>>>>>>>>>>>compress image.<<<<<<<<<<<<<<"<<std::endl;
+        osg::notify(osg::INFO)<<">>>>>>>>>>>>>>>compressed image.<<<<<<<<<<<<<<"<<std::endl;
 
     }
     else
@@ -1930,7 +1929,7 @@ osg::StateSet* DataSet::DestinationTile::createStateSet()
 
     if (!inlineImageFile)
     {
-        std::cout<<"Writing out imagery to "<<imageName<<std::endl;
+        osg::notify(osg::INFO)<<"Writing out imagery to "<<imageName<<std::endl;
         osgDB::writeImageFile(*_imagery->_image,_imagery->_image->getFileName().c_str());
     }
 
@@ -1943,7 +1942,7 @@ osg::Node* DataSet::DestinationTile::createHeightField()
 
     if (_terrain.valid() && _terrain->_heightField.valid())
     {
-        std::cout<<"--- Have terrain build tile ----"<<std::endl;
+        osg::notify(osg::INFO)<<"--- Have terrain build tile ----"<<std::endl;
 
         osg::HeightField* hf = _terrain->_heightField.get();
         
@@ -1953,7 +1952,7 @@ osg::Node* DataSet::DestinationTile::createHeightField()
     }
     else 
     {
-        std::cout<<"**** No terrain to build tile from use flat terrain fallback ****"<<std::endl;
+        osg::notify(osg::INFO)<<"**** No terrain to build tile from use flat terrain fallback ****"<<std::endl;
         // create a dummy height field to file in the gap
         osg::HeightField* hf = new osg::HeightField;
         hf->allocate(2,2);
@@ -2012,7 +2011,7 @@ static inline osg::Vec3 computeLocalSkirtVector(const osg::EllipsoidModel* et, c
 
 osg::Node* DataSet::DestinationTile::createPolygonal()
 {
-    std::cout<<"--------- DataSet::DestinationTile::createDrawableGeometry() ------------- "<<std::endl;
+    osg::notify(osg::INFO)<<"--------- DataSet::DestinationTile::createDrawableGeometry() ------------- "<<std::endl;
 
     const osg::EllipsoidModel* et = _dataSet->getEllipsoidModel();
     bool mapLatLongsToXYZ = _dataSet->getConvertFromGeographicToGeocentric() && et;
@@ -2022,7 +2021,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
     
     if (_terrain.valid() && _terrain->_heightField.valid())
     {
-        std::cout<<"--- Have terrain build tile ----"<<std::endl;
+        osg::notify(osg::INFO)<<"--- Have terrain build tile ----"<<std::endl;
         grid = _terrain->_heightField.get();
     }
     else
@@ -2039,7 +2038,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
             if (longitude_range>45.0) numColumns = (unsigned int)ceilf((float)numColumns*sqrtf(longitude_range/45.0));
             if (latitude_range>45.0) numRows = (unsigned int)ceilf((float)numRows*sqrtf(latitude_range/45.0));
             
-            std::cout<<"numColumns = "<<numColumns<<"  numRows="<<numRows<<std::endl;
+            osg::notify(osg::INFO)<<"numColumns = "<<numColumns<<"  numRows="<<numRows<<std::endl;
         }
         else
         {
@@ -2057,7 +2056,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
 
     if (!grid)
     {
-        std::cout<<"**** No terrain to build tile from use flat terrain fallback ****"<<std::endl;
+        osg::notify(osg::INFO)<<"**** No terrain to build tile from use flat terrain fallback ****"<<std::endl;
         
         return 0;
     }
@@ -2208,7 +2207,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
                 double phi = 2.0 * asin (d*0.5/globe_radius); // d/globe_radius;
                 double beta = theta+phi;
                 double cutoff = osg::PI_2 - 0.1;
-                //std::cout<<"theta="<<theta<<"\tphi="<<phi<<" beta "<<beta<<std::endl;
+                //osg::notify(osg::INFO)<<"theta="<<theta<<"\tphi="<<phi<<" beta "<<beta<<std::endl;
                 if (phi<cutoff && beta<cutoff)
                 {
 
@@ -2219,7 +2218,7 @@ osg::Node* DataSet::DestinationTile::createPolygonal()
                 }
                 else
                 {
-                    //std::cout<<"Turning off cluster culling for wrap around tile."<<std::endl;
+                    //osg::notify(osg::INFO)<<"Turning off cluster culling for wrap around tile."<<std::endl;
                     useClusterCullingCallback = false;
                 }
             }
@@ -2477,19 +2476,18 @@ void DataSet::DestinationTile::readFrom(CompositeSource* sourceGraph)
 {
     allocate();
 
-    std::cout<<"DestinationTile::readFrom() "<<std::endl;
+    osg::notify(osg::INFO)<<"DestinationTile::readFrom() "<<std::endl;
     for(CompositeSource::source_iterator itr(sourceGraph);itr.valid();++itr)
     {
     
         Source* source = itr->get();
-        
-        if (_level>=source->getMinLevel() && _level<=source->getMaxLevel()) 
+        if (source && _level>=source->getMinLevel() && _level<=source->getMaxLevel()) 
         {
             SourceData* data = (*itr)->getSourceData();
             if (data)
             {
-                std::cout<<"DataSet::DestinationTile::readFrom -> SourceData::read() "<<std::endl;
-                std::cout<<"    destination._level="<<_level<<"\t"<<source->getMinLevel()<<"\t"<<source->getMaxLevel()<<std::endl;
+                osg::notify(osg::INFO)<<"DataSet::DestinationTile::readFrom -> SourceData::read() "<<std::endl;
+                osg::notify(osg::INFO)<<"    destination._level="<<_level<<"\t"<<source->getMinLevel()<<"\t"<<source->getMaxLevel()<<std::endl;
                 if (_imagery.valid()) data->read(*_imagery);
                 if (_terrain.valid()) data->read(*_terrain);
             }
@@ -2512,8 +2510,7 @@ void DataSet::DestinationTile::addRequiredResolutions(CompositeSource* sourceGra
     for(CompositeSource::source_iterator itr(sourceGraph);itr.valid();++itr)
     {
         Source* source = itr->get();
-
-        if (source->intersects(*this))
+        if (source && source->intersects(*this))
         {
             if (source->getType()==Source::IMAGE)
             {
@@ -2579,7 +2576,7 @@ void DataSet::CompositeDestination::addRequiredResolutions(CompositeSource* sour
 
 void DataSet::CompositeDestination::readFrom(CompositeSource* sourceGraph)
 {
-    std::cout<<"CompositeDestination::readFrom() "<<std::endl;
+    osg::notify(osg::INFO)<<"CompositeDestination::readFrom() "<<std::endl;
 
     // handle leaves
     for(TileList::iterator titr=_tiles.begin();
@@ -2885,7 +2882,7 @@ void DataSet::CompositeDestination::unrefLocalData()
         ++titr)
     {
         DestinationTile* tile = titr->get();
-        std::cout<<"   unref tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
+        osg::notify(osg::INFO)<<"   unref tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
         tile->unrefData();
     }
 }
@@ -3047,7 +3044,7 @@ DataSet::CompositeDestination* DataSet::createDestinationGraph(CompositeDestinat
 
         if (needToDivideX && needToDivideY)
         {
-            std::cout<<"Need to Divide X + Y for level "<<currentLevel<<std::endl;
+            osg::notify(osg::INFO)<<"Need to Divide X + Y for level "<<currentLevel<<std::endl;
             // create four tiles.
             osg::BoundingBox bottom_left(extents.xMin(),extents.yMin(),extents.zMin(),xCenter,yCenter,extents.zMax());
             osg::BoundingBox bottom_right(xCenter,extents.yMin(),extents.zMin(),extents.xMax(),yCenter,extents.zMax());
@@ -3107,7 +3104,7 @@ DataSet::CompositeDestination* DataSet::createDestinationGraph(CompositeDestinat
         }
         else if (needToDivideX)
         {
-            std::cout<<"Need to Divide X only"<<std::endl;
+            osg::notify(osg::INFO)<<"Need to Divide X only"<<std::endl;
 
             // create two tiles.
             osg::BoundingBox left(extents.xMin(),extents.yMin(),extents.zMin(),xCenter,extents.yMax(),extents.zMax());
@@ -3147,7 +3144,7 @@ DataSet::CompositeDestination* DataSet::createDestinationGraph(CompositeDestinat
         }
         else if (needToDivideY)
         {
-            std::cout<<"Need to Divide Y only"<<std::endl;
+            osg::notify(osg::INFO)<<"Need to Divide Y only"<<std::endl;
 
             // create two tiles.
             osg::BoundingBox top(extents.xMin(),yCenter,extents.zMin(),extents.xMax(),extents.yMax(),extents.zMax());
@@ -3186,7 +3183,7 @@ DataSet::CompositeDestination* DataSet::createDestinationGraph(CompositeDestinat
         }
         else
         {
-            std::cout<<"No Need to Divide"<<std::endl;
+            osg::notify(osg::INFO)<<"No Need to Divide"<<std::endl;
         }
     }
     
@@ -3207,7 +3204,7 @@ void DataSet::computeDestinationGraphFromSources(unsigned int numLevels)
                 if (sd->_cs.valid())
                 {
                     _destinationCoordinateSystem = sd->_cs;
-                    std::cout<<"Setting coordinate system to "<<_destinationCoordinateSystem->getCoordinateSystem()<<std::endl;
+                    osg::notify(osg::INFO)<<"Setting coordinate system to "<<_destinationCoordinateSystem->getCoordinateSystem()<<std::endl;
                     break;
                 }
             }
@@ -3218,7 +3215,7 @@ void DataSet::computeDestinationGraphFromSources(unsigned int numLevels)
     {
         CoordinateSystemType cst = getCoordinateSystemType(_destinationCoordinateSystem.get());
 
-        std::cout<<"new DataSet::createDestination()"<<std::endl;
+        osg::notify(osg::INFO)<<"new DataSet::createDestination()"<<std::endl;
         if (cst!=GEOGRAPHIC && getConvertFromGeographicToGeocentric())
         {
             // need to use the geocentric coordinate system as a base for creating an geographic intermediate
@@ -3255,15 +3252,15 @@ void DataSet::computeDestinationGraphFromSources(unsigned int numLevels)
             if (sd)
             {
                 osg::BoundingBox local_extents(sd->getExtents(_intermediateCoordinateSystem.get()));
-                std::cout<<"local_extents = xMin()"<<local_extents.xMin()<<" "<<local_extents.xMax()<<std::endl;
-                std::cout<<"                yMin()"<<local_extents.yMin()<<" "<<local_extents.yMax()<<std::endl;
+                osg::notify(osg::INFO)<<"local_extents = xMin()"<<local_extents.xMin()<<" "<<local_extents.xMax()<<std::endl;
+                osg::notify(osg::INFO)<<"                yMin()"<<local_extents.yMin()<<" "<<local_extents.yMax()<<std::endl;
                 extents.expandBy(local_extents);
             }
         }
     }
 
-    std::cout<<"extents = xMin()"<<extents.xMin()<<" "<<extents.xMax()<<std::endl;
-    std::cout<<"          yMin()"<<extents.yMin()<<" "<<extents.yMax()<<std::endl;
+    osg::notify(osg::INFO)<<"extents = xMin()"<<extents.xMin()<<" "<<extents.xMax()<<std::endl;
+    osg::notify(osg::INFO)<<"          yMin()"<<extents.yMin()<<" "<<extents.yMax()<<std::endl;
 
     // then create the destinate graph accordingly.
 
@@ -3290,7 +3287,7 @@ struct DerefLessFunctor
 {
     bool operator () (const T& lhs, const T& rhs)
     {
-        return lhs->getSortValue() > rhs->getSortValue();
+        return (lhs.valid() && rhs.valid()) && (lhs->getSortValue() > rhs->getSortValue());
     }
 };
 
@@ -3316,7 +3313,7 @@ void DataSet::CompositeSource::sort()
     // sort the composite sources internal data
     for(ChildList::iterator itr=_children.begin();itr!=_children.end();++itr)
     {
-        (*itr)->sort();
+        if (itr->valid()) (*itr)->sort();
     }
 }
 
@@ -3333,33 +3330,35 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
         for(CompositeSource::source_iterator sitr(_sourceGraph.get());sitr.valid();++sitr)
         {
             Source* source = sitr->get();
-            std::cout<<"Source File "<<source->getFileName()<<std::endl;
-            
-
-            const Source::ResolutionList& resolutions = source->getRequiredResolutions();
-            std::cout<<"    resolutions.size() "<<resolutions.size()<<std::endl;
-            std::cout<<"    { "<<std::endl;
-            Source::ResolutionList::const_iterator itr;
-            for(itr=resolutions.begin();
-                itr!=resolutions.end();
-                ++itr)
+            if (source) 
             {
-                std::cout<<"        resX="<<itr->_resX<<" resY="<<itr->_resY<<std::endl;
-            }
-            std::cout<<"    } "<<std::endl;
+                osg::notify(osg::INFO)<<"Source File "<<source->getFileName()<<std::endl;
 
-            source->consolodateRequiredResolutions();
-            
-            std::cout<<"    consolodated resolutions.size() "<<resolutions.size()<<std::endl;
-            std::cout<<"    consolodated { "<<std::endl;
-            for(itr=resolutions.begin();
-                itr!=resolutions.end();
-                ++itr)
-            {
-                std::cout<<"        resX="<<itr->_resX<<" resY="<<itr->_resY<<std::endl;
-            }
-            std::cout<<"    } "<<std::endl;
 
+                const Source::ResolutionList& resolutions = source->getRequiredResolutions();
+                osg::notify(osg::INFO)<<"    resolutions.size() "<<resolutions.size()<<std::endl;
+                osg::notify(osg::INFO)<<"    { "<<std::endl;
+                Source::ResolutionList::const_iterator itr;
+                for(itr=resolutions.begin();
+                    itr!=resolutions.end();
+                    ++itr)
+                {
+                    osg::notify(osg::INFO)<<"        resX="<<itr->_resX<<" resY="<<itr->_resY<<std::endl;
+                }
+                osg::notify(osg::INFO)<<"    } "<<std::endl;
+
+                source->consolodateRequiredResolutions();
+
+                osg::notify(osg::INFO)<<"    consolodated resolutions.size() "<<resolutions.size()<<std::endl;
+                osg::notify(osg::INFO)<<"    consolodated { "<<std::endl;
+                for(itr=resolutions.begin();
+                    itr!=resolutions.end();
+                    ++itr)
+                {
+                    osg::notify(osg::INFO)<<"        resX="<<itr->_resX<<" resY="<<itr->_resY<<std::endl;
+                }
+                osg::notify(osg::INFO)<<"    } "<<std::endl;
+            }
 
         }
 
@@ -3372,7 +3371,7 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
         for(CompositeSource::source_iterator itr(_sourceGraph.get());itr.valid();++itr)
         {
             Source* source = itr->get();
-            if (source->needReproject(_intermediateCoordinateSystem.get()))
+            if (source && source->needReproject(_intermediateCoordinateSystem.get()))
             {
                 // do the reprojection to a tempory file.
                 std::string newFileName = temporyFilePrefix + osgDB::getStrippedName(source->getFileName()) + ".tif";
@@ -3381,6 +3380,11 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
                 
                 // replace old source by new one.
                 if (newSource) *itr = newSource;
+                else
+                {
+                    osg::notify(osg::WARN)<<"Failed to reproject"<<source->getFileName()<<std::endl;
+                    *itr = 0;
+                }
             }
         }
     }
@@ -3390,7 +3394,7 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
         for(CompositeSource::source_iterator itr(_sourceGraph.get());itr.valid();++itr)
         {
             Source* source = itr->get();
-            source->buildOverviews();
+            if (source) source->buildOverviews();
         }
     }
 
@@ -3399,8 +3403,11 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
         for(CompositeSource::source_iterator itr(_sourceGraph.get());itr.valid();++itr)
         {
             Source* source = itr->get();
-            source->setSortValueFromSourceDataResolution();
-            std::cout<<"sort "<<source->getFileName()<<" value "<<source->getSortValue()<<std::endl;
+            if (source)
+            {
+                source->setSortValueFromSourceDataResolution();
+                osg::notify(osg::INFO)<<"sort "<<source->getFileName()<<" value "<<source->getSortValue()<<std::endl;
+            }
             
         }
         
@@ -3409,17 +3416,23 @@ void DataSet::updateSourcesForDestinationGraphNeeds()
         _sourceGraph->sort();
     }
     
-    std::cout<<"Using source_lod_iterator itr"<<std::endl;
+    osg::notify(osg::INFO)<<"Using source_lod_iterator itr"<<std::endl;
     for(CompositeSource::source_lod_iterator csitr(_sourceGraph.get(),CompositeSource::LODSourceAdvancer(0.0));csitr.valid();++csitr)
     {
-        std::cout<<"  LOD "<<(*csitr)->getFileName()<<std::endl;
+        Source* source = csitr->get();
+        if (source)
+        {
+            osg::notify(osg::INFO)<<"  LOD "<<(*csitr)->getFileName()<<std::endl;
+        }
     }
-    std::cout<<"End of Using Source Iterator itr"<<std::endl;
+    osg::notify(osg::INFO)<<"End of Using Source Iterator itr"<<std::endl;
     
 }
 
 void DataSet::populateDestinationGraphFromSources()
 {
+    osg::notify(osg::NOTICE)<<std::endl<<"started DataSet::populateDestinationGraphFromSources)"<<std::endl;
+
     if (_databaseType==LOD_DATABASE)
     {
 
@@ -3437,12 +3450,13 @@ void DataSet::populateDestinationGraphFromSources()
         //  compute x and y range
         //  from top row down to bottom row equalize boundairies a write out
     }
+    osg::notify(osg::NOTICE)<<"completed DataSet::populateDestinationGraphFromSources)"<<std::endl;
 }
 
 
 void DataSet::_readRow(Row& row)
 {
-    std::cout<<"_readRow "<<row.size()<<std::endl;
+    osg::notify(osg::NOTICE)<<"_readRow "<<row.size()<<std::endl;
     for(Row::iterator citr=row.begin();
         citr!=row.end();
         ++citr)
@@ -3453,7 +3467,7 @@ void DataSet::_readRow(Row& row)
             ++titr)
         {
             DestinationTile* tile = titr->get();
-            std::cout<<"   reading tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
+            osg::notify(osg::NOTICE)<<"   reading tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
             tile->readFrom(_sourceGraph.get());
         }
     }
@@ -3461,7 +3475,7 @@ void DataSet::_readRow(Row& row)
 
 void DataSet::_equalizeRow(Row& row)
 {
-    std::cout<<"_equalizeRow "<<row.size()<<std::endl;
+    osg::notify(osg::NOTICE)<<"_equalizeRow "<<row.size()<<std::endl;
     for(Row::iterator citr=row.begin();
         citr!=row.end();
         ++citr)
@@ -3472,7 +3486,7 @@ void DataSet::_equalizeRow(Row& row)
             ++titr)
         {
             DestinationTile* tile = titr->get();
-            std::cout<<"   equalizing tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
+            osg::notify(osg::NOTICE)<<"   equalizing tile level="<<tile->_level<<" X="<<tile->_tileX<<" Y="<<tile->_tileY<<std::endl;
             tile->equalizeBoundaries();
             tile->setTileComplete(true);
         }
@@ -3481,7 +3495,7 @@ void DataSet::_equalizeRow(Row& row)
 
 void DataSet::_writeRow(Row& row)
 {
-    std::cout<<"_writeRow "<<row.size()<<std::endl;
+    osg::notify(osg::NOTICE)<<"_writeRow "<<row.size()<<std::endl;
     for(Row::iterator citr=row.begin();
         citr!=row.end();
         ++citr)
@@ -3497,14 +3511,14 @@ void DataSet::_writeRow(Row& row)
                 std::string filename = parent->getSubTileName();
                 if (node.valid())
                 {
-                    std::cout<<"   writeSubTile filename="<<filename<<std::endl;
+                    osg::notify(osg::NOTICE)<<"   writeSubTile filename="<<filename<<std::endl;
                     osgDB::writeNodeFile(*node,filename);
                     parent->setSubTilesGenerated(true);
                     parent->unrefSubTileData();
                 }
                 else
                 {
-                    std::cout<<"   failed to writeSubTile node for tile, filename="<<filename<<std::endl;
+                    osg::notify(osg::WARN)<<"   failed to writeSubTile node for tile, filename="<<filename<<std::endl;
                 }
             }
         }
@@ -3522,12 +3536,12 @@ void DataSet::_writeRow(Row& row)
 
             if (node.valid())
             {
-                std::cout<<"   writeNodeFile = "<<cd->_level<<" X="<<cd->_tileX<<" Y="<<cd->_tileY<<" filename="<<filename<<std::endl;
+                osg::notify(osg::NOTICE)<<"   writeNodeFile = "<<cd->_level<<" X="<<cd->_tileX<<" Y="<<cd->_tileY<<" filename="<<filename<<std::endl;
                 osgDB::writeNodeFile(*node,filename);
             }
             else
             {
-                std::cout<<"   faild to write node for tile = "<<cd->_level<<" X="<<cd->_tileX<<" Y="<<cd->_tileY<<" filename="<<filename<<std::endl;
+                osg::notify(osg::WARN)<<"   faild to write node for tile = "<<cd->_level<<" X="<<cd->_tileX<<" Y="<<cd->_tileY<<" filename="<<filename<<std::endl;
             }
         }
     }
@@ -3535,11 +3549,14 @@ void DataSet::_writeRow(Row& row)
 
 void DataSet::createDestination(unsigned int numLevels)
 {
+    osg::notify(osg::NOTICE)<<std::endl<<"started DataSet::createDestination("<<numLevels<<")"<<std::endl;
 
     computeDestinationGraphFromSources(numLevels);
     
     updateSourcesForDestinationGraphNeeds();
 
+    osg::notify(osg::NOTICE)<<"completed DataSet::createDestination("<<numLevels<<")"<<std::endl;
+
 }
 
 osg::Node* DataSet::decorateWithCoordinateSystemNode(osg::Node* subgraph)
@@ -3567,7 +3584,7 @@ void DataSet::writeDestination()
     if (_destinationGraph.valid())
     {
         std::string filename = _tileBasename+_tileExtension;
-        std::cout<<"DataSet::writeDestination("<<filename<<")"<<std::endl;
+        osg::notify(osg::NOTICE)<<std::endl<<"started DataSet::writeDestination("<<filename<<")"<<std::endl;
 
         if (_databaseType==LOD_DATABASE)
         {
@@ -3593,7 +3610,7 @@ void DataSet::writeDestination()
                 // skip is level is empty.
                 if (level.empty()) continue;
                 
-                std::cout<<"New level"<<std::endl;
+                osg::notify(osg::INFO)<<"New level"<<std::endl;
 
                 Level::iterator prev_itr = level.begin();
                 _readRow(prev_itr->second);
@@ -3615,10 +3632,12 @@ void DataSet::writeDestination()
                 _writeRow(prev_itr->second);
             }
         }
+        osg::notify(osg::NOTICE)<<"completed DataSet::writeDestination("<<filename<<")"<<std::endl;
+
     }
     else
     {
-        std::cout<<"Error: no scene graph to output, no file written."<<std::endl;
+        osg::notify(osg::WARN)<<"Error: no scene graph to output, no file written."<<std::endl;
     }
 
 }