Writing ascii file

Hello,

I want to write a 256*128 dataset in asii file. I want my format to be:
lon,lat,slp_mon
where lon and lat to be the actual longitude and latitude values from the
netcdf file that I am using.
Here is my code that I am trying to use:

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

;***********************************************
;read in netcdf data
;***********************************************

 in = addfile("./ctrl.cam2.h0.1997-11-01-00000.nc","r")
 slp = in->PSL(123:246,:,:)
 lat = in ->lat
 lon = in->lon
 slp=slp/100 ;convert Pa to hPa

 printVarSummary(lon)
;***********************************************
;Calculate monthly mean
;***********************************************

slp_mon = dim_avg( slp(lat|:, lon|:, time|:) ) ;to reorder the array so
that time is the rightmost dimension

printVarSummary(slp_mon)
printVarSummary(lat)
printVarSummary(lon)

;************************************************
;Write the monthly mean to ascii file
;************************************************
slp_1d = new ((/32768/) , "float")
lon_1d = new ((/32768/) , "float")
lat_1d = new ((/32768/) , "float")

data = new ( (/ 3, 32768 /) , "float")

data (0,:) = slp_1d
data (1,:) = lon_1d
data (2,:) = lat_1d

index =0
do j =1,255
  do i =1,127
  slp_1d(index)= slp_mon(i,j)
  lon_1d(index) = j
  lat_1d(index) = i
  end do
end do

fmtx = "f3.4,6x,f3.4,6x,f4.5"

write_matrix(data,fmtx, False)

;************************************************
;Check
;************************************************
print("OK")

end

Thanks,

-- 
*************************************
Soumik Basu
Graduate Student, Research Assistant
International Arctic Research Center
University of Alaska Fairbanks
Fairbanks,Alaska,USA
**************************************
-----------------------------------------
This email was sent using SquirrelMail.
   "Webmail for nuts!"
http://squirrelmail.org/
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk