begin

	dirasc = "./"
        filasc = systemfunc("cd "+dirasc+" ; ls ascii.txt")    ;wrong with others months
        pthasc = dirasc+filasc
        nfil   = dimsizes(filasc)
        print(filasc)
        print("----")

;----- Read the values in as 1D array of strings to get rows and columns.
        data   = asciiread(pthasc,-1,"string")
        ncols  = dimsizes(str_split(data(0),","))
        nrows  = dimsizes(data)
	
;----- Open pc variable (for each latitude/longitude point)
        lat    = stringtodouble( str_get_field(data(1:), 8, ",") )    
        lon    = stringtodouble( str_get_field(data(1:), 9, ",") )  
       ;lon    = 360.0+(lon)     ; -180 180 
        pc     = stringtodouble( str_get_field(data(1:), 10, ",") )

	var1       = lat
	var1       = 1
        var        = doubletofloat(var1)
        var@lat1d = lat
	var@lon1d = lon


;------ Create coordinates  
	
	
        lonW  = (/-119.995/) 
        lonE  = (/-33.005/) 
        latS  = (/-55.995/)  
        latN  = (/32.995/)   

	lat0  = latS
	lat1  = latN
	lon0  = lonW
	lon1  = lonE
	dgrad = 0.01        

	npointslat = toint(abs(lat1-lat0)/dgrad)+1
	npointslon = toint(abs(lon1-lon0)/dgrad)+1

	lat2 = fspan(lat0,lat1,npointslat)*1.d
	lon2 = fspan(lon0,lon1,npointslon)*1.d

	lat2!0         = "lat" 
	lat2&lat       =  lat2
	lat2@long_name = "lat" 
	lat2@units     = "degrees_north" 
	
	lon2!0         = "lon" 
	lon2&lon       =  lon2
	lon2@long_name = "lon" 
	lon2@units     = "degrees_east" 
	
	nlats = dimsizes(lat2)
	nlons = dimsizes(lon2)
	
;----- Create a new var to pc's values
	new_data           =  new((/npointslat,npointslon/),"float")
	new_data           =  0.     ; start
	new_data!0         = "lat" 
	new_data!1         = "lon" 
	new_data&lat       =  lat2
	new_data&lon       =  lon2

;------ 	PROBLEM
	do i=0,nrows-1
	   laat                  = lat(i)
	   loon                  = lon(i)
	   new_data({la},{lo}) = new_data({laat},{loon})+(/var(i)/) 
	end do


	new_data = where(new_data .gt. 0, new_data, new_data@_FillValue)