
Dear All,
           I am a research student in environment. I have only little programming knowledge. I am currently doing the last project about rainfall impact on ground water quality in an area. It happens that I have to use R to read rainfall data (3 dimension) from ASC file (*.asc), and then write them into one NCDF file (*.nc).

           I have been working very hard on study R, but I still can not fix the problem. Could someone please as kind as point out that what the problems are in my R script?

         Firstly, this is an example of data in asc file:
                    NCOLS      241

                    NROWS      160

                    XLLCORNER    140.00012207031    

                    YLLCORNER   -39.000000000000    

                    CELLSIZE    0.50000000000000E-01

                    NODATA_VALUE   -99.00000    

                    166.30  160.87  155.23  149.33  143.83  138.52  133.29  128.34  123.76  119.21

                    115.06  110.90  107.22  103.69  100.40   97.29   94.58   92.15   90.00   87.91

                    86.20   84.57   83.22   81.94   81.11   80.38   79.37   78.73   79.70   79.62
          -----------------------------------------------------------------------------------------------------------------------      



          And then this is the script I wrote:

 setwd("E:/grid")

#defining dimension

x=dim.def.ncdf("Lon","degreesE",140.0251:146.6751)

y=dim.def.ncdf("Lat","degreesN",(-31.025):(-38.975))

t=dim.def.ncdf("Time","1968-01",1:12,unlim=TRUE)

 
 
#setup variable

varmr=var.def.ncdf("mr","mm",list(x,y,t),-99.00,

longname="monthly rainfall")

 
#create ncdf file

ncnew=create.ncdf("rainfall.nc", varmr)

 
#read input

files=list.files(pattern=".asc")

mrain=matrix(0:0,0,3)

 
for(i in files)

{rainfall=data.frame(readAsciiGrid(i))

  mrain=cbind(mrain,rainfall)

}

put.var.ncdf(ncnew, mrain)
 

close.ncdf(ncnew)

-------------------------------------------------------------------------------

[[elided Hotmail spam]]





Many thanks,
Jun



 		 	   		  
	[[alternative HTML version deleted]]

