--- a/extra/ncArray/inst/private/ncarray_example_file.m	Wed Dec 18 08:51:51 2013 +0000
+++ b/extra/ncArray/inst/private/ncarray_example_file.m	Thu Dec 19 15:06:03 2013 +0000
@@ -4,21 +4,24 @@
       ~isempty(which('ncwrite'))
   % use matlab netcdf high level interface
 
+  %dtype = 'single';
+  dtype = 'double';
+
   % Variables
-  nccreate(filename,'lon','Format','classic','Datatype','single',...
+  nccreate(filename,'lon','Format','classic','Datatype',dtype,...
            'Dimensions',{'x',220, 'y',144});
   ncwriteatt(filename,'lon','long_name','Longitude')
   ncwriteatt(filename,'lon','units','degrees_east')
   
-  nccreate(filename,'lat','Datatype','single','Dimensions',{'x',220, 'y',144});
+  nccreate(filename,'lat','Datatype',dtype,'Dimensions',{'x',220, 'y',144});
   ncwriteatt(filename,'lat','long_name','Latitude')
   ncwriteatt(filename,'lat','units','degrees_north')
   
-  nccreate(filename,'time','Datatype','single','Dimensions',{'time',1});
+  nccreate(filename,'time','Datatype',dtype,'Dimensions',{'time',1});
   ncwriteatt(filename,'time','long_name','Time')
   ncwriteatt(filename,'time','units','days since 1858-11-17 00:00:00 GMT')
   
-  nccreate(filename,'SST','Datatype','single','Dimensions',...
+  nccreate(filename,'SST','Datatype',dtype,'Dimensions',...
            {'x',220, 'y',144, 'time',1});
   ncwriteatt(filename,'SST','missing_value',single(9999))
   ncwriteatt(filename,'SST','_FillValue',single(9999))
@@ -40,21 +43,21 @@
 
   % variables
 
-  nc{'lon'} = ncfloat('y','x');  % 31680 elements 
+  nc{'lon'} = ncdouble('y','x');  % 31680 elements 
   nc{'lon'}.long_name = ncchar('Longitude');
   nc{'lon'}.units = ncchar('degrees_east');
   
-  nc{'lat'} = ncfloat('y','x');  % 31680 elements 
+  nc{'lat'} = ncdouble('y','x');  % 31680 elements 
   nc{'lat'}.long_name = ncchar('Latitude');
   nc{'lat'}.units = ncchar('degrees_north');
 
-  nc{'time'} = ncfloat('time');  % 1 elements 
+  nc{'time'} = ncdouble('time');  % 1 elements 
   nc{'time'}.long_name = ncchar('Time');
   nc{'time'}.units = ncchar('days since 1858-11-17 00:00:00 GMT');
   
-  nc{'SST'} = ncfloat('time','y','x');  % 31680 elements 
-  nc{'SST'}.missing_value = ncfloat(9999);
-  nc{'SST'}.FillValue_ = ncfloat(9999);
+  nc{'SST'} = ncdouble('time','y','x');  % 31680 elements 
+  nc{'SST'}.missing_value = ncdouble(9999);
+  nc{'SST'}.FillValue_ = ncdouble(9999);
   nc{'SST'}.units = ncchar('degC');
   nc{'SST'}.long_name = ncchar('Sea Surface Temperature');
   nc{'SST'}.coordinates = ncchar('lat lon');