909,910d908
<                 double	dummy[2], xMinMax[2], yMinMax[2];
<                 int	node_offset = 0;
913a912
>                 int node_offset = 0;
916,946c915,916
<                 if (!nc_get_att_double (cdfid, nVarDimXID, "actual_range", dummy)) {
<                     xMinMax[0] = dummy[0];		
<                     xMinMax[1] = dummy[1];
<                 }
<                 else {
<                     xMinMax[0] = pdfXCoord[0];
<                     xMinMax[1] = pdfXCoord[xdim-1];
<                     node_offset = 0;
<                 }
< 
<                 if (!nc_get_att_double (cdfid, nVarDimYID, "actual_range", dummy)) {
<                     yMinMax[0] = dummy[0];		
<                     yMinMax[1] = dummy[1];
<                 }
<                 else {
<                     yMinMax[0] = pdfYCoord[0];	
<                     yMinMax[1] = pdfYCoord[ydim-1];
<                     node_offset = 0;
<                 }
< 
< #ifdef notdef
<                 // Check for reverse order of y-coordinate
<                 if ( yMinMax[0] > yMinMax[1] ) {
<                     dummy[0] = yMinMax[1];
<                     dummy[1] = yMinMax[0];
<                     yMinMax[0] = dummy[0];
<                     yMinMax[1] = dummy[1];
<                 }
< #endif
< 
<                 poDS->adfGeoTransform[0] = xMinMax[0];
---
>                 poDS->adfGeoTransform[0] = pdfXCoord[0];
>                 poDS->adfGeoTransform[3] = pdfYCoord[0];
948d917
<                 poDS->adfGeoTransform[3] = yMinMax[1];
950,953c919,925
<                 poDS->adfGeoTransform[1] = ( xMinMax[1] - xMinMax[0] ) / 
<                     ( poDS->nRasterXSize + (node_offset - 1) );
<                 poDS->adfGeoTransform[5] = ( yMinMax[0] - yMinMax[1] ) / 
<                     ( poDS->nRasterYSize + (node_offset - 1) );
---
>                 poDS->adfGeoTransform[1] = (( pdfXCoord[xdim-1] - 
>                                               pdfXCoord[0] ) / 
>                                             ( poDS->nRasterXSize + (node_offset - 1) ));
> 
>                 poDS->adfGeoTransform[5] = (( pdfYCoord[ydim-1] - 
>                                               pdfYCoord[0] ) / 
>                                             ( poDS->nRasterYSize + (node_offset - 1) ));