RE: [ncl-talk] netcdf to ascii conversion

Hi,

Thanks Dennis and Arindam for helping me.

I use the script mentioned in example3 in website
http://www.ncl.ucar.edu/Document/Functions/Built-in/asciiwrite.shtml
but I got a error massage

fatal:new: data type must either be a keyword or string

in the line where new variable data is declared i.e.
data = new(fName, npts, "string")

Thanks
Debasish

-----Original Message-----
From: Dennis Shea [mailto:shea_at_ucar.edu]
Sent: Monday, June 04, 2007 7:51 AM
To: Arindam Chakraborty
Cc: Debasish Pai Mazumder; ncl-talk_at_ucar.edu
Subject: Re: netcdf to ascii conversion

Based on Debasish's question and Arindam's responses, I created
a new Example 3 at:

http://www.ncl.ucar.edu/Document/Functions/Built-in/asciiwrite.shtml

=======================
Arindam Chakraborty wrote:
> Oh, sorry Debasish, I did a last minute editing on that line before
> sending. sprintf takes two arguments and while testing I used lat and
> lon separately. But at the time of sending I though of combining them
> together keeping c in mind, but I was wrong. You can replace that line
> with the followings lines:
> ascval = ascval+sprintf("%8.2f ",lat(j))
> ascval = ascval+sprintf("%8.2f ",lon(i))
>
>
> It should work after this. Thanks for pointing this out.
>
> arindam
>
> Debasish Pai Mazumder wrote:
>> Hi Arindam,
>>
>> Thanks for your reply
>>
>> I used your code but it give some error massage in line
>>
>> ascval = ascval+sprintf("%8.2f %8.2f ",lat(j),lon(i))
>>
>> Error massage: fatal:syntax error: function sprintf expects 2 arguments,
got
>> 3
>> fatal:Syntax Error in block, block not executed
>>
>> Thanks
>>
>> Debasish
>>
>>
>> -----Original Message-----
>> From: Arindam Chakraborty [mailto:arindam_at_met.fsu.edu]
>> Sent: Sunday, June 03, 2007 2:45 PM
>> To: Debasish Pai Mazumder
>> Cc: ncl-talk_at_ucar.edu
>> Subject: Re: netcdf to ascii conversion
>>
>> hi Debasish,
>> there may be an easy way out for this in NCL, but the following should
>> also work:
>>
>> ...
>> ...
>>
>> ; var is the variable name,
>> ; lon is the longitude dimension name,
>> ; lat is the latitude dimension name.
>>
>> lon = var&lon
>> lat = var&lat
>>
>> ; the following assumes the order of dimensions in the variables as
>> (time,lat,lon)
>> dimsz = dimsizes(var)
>> nt = dimsz(0)
>> ny = dimsz(1)
>> nx = dimsz(2)
>>
>> fout = "ascfile.dat" ; output ascii file name
>> system("echo > "+fout) ; create an empty file. Existing file will be
>> overwritten.
>>
>> do j = 0, ny-1
>> do i = 0, nx-1
>> ascval = ""
>> ascval = ascval+sprintf("%8.2f %8.2f ",lat(j),lon(i))
>> do k = 0, nt-1
>> ascval = ascval+sprintf("%6.1f ",var(k,j,i))
>> end do
>> system("echo "+ascval+" >> "+fout)
>> end do
>> end do
>>
>> ...
>> ...
>>
>> Hope this helps,
>>
>> arindam
>>
>>
>> Debasish Pai Mazumder wrote:
>>
>>> Hi all,
>>>
>>> I have a netcdf file for air temperature. The dimension of air
>>> temperature is given below,
>>>
>>> Air_temp(time,lat,lon) where time=124, lat=16, lon=59
>>>
>>> I am trying to convert this fie in to ascii file. I used the code
>>> provided by Dennis Shea in website
>>> http://www.ncl.ucar.edu/Support/talk_archives/2006/0113.html and the
>>> program gave me different ascii files for lat, lon and temperature but
>>> all those values in ascii file are written in a one column. I also
>>> used write_matrix option but I guess it only works for two dimensional
>>> data.
>>>
>>> I am trying to make a ascii file with different rows and columns like
>>>
>>> lat lon time=0 time=1... time=123
>>>
>>> 41.5 25 260 262.6 267.8
>>>
>>> 41.5 27.5 261.2 262.3 264.5
>>>
>>> .
>>>
>>> .
>>>
>>> .
>>>
>>> 41.5 160 254.6 262.9 270.5
>>>
>>> 43 25 260.2 264.4 267.8
>>>
>>> 43 27.5 ..................
>>>
>>> ..........................
>>>
>>> .........................
>>>
>>> Is there any way to get ascii file with the given format from netcdf
>>> file in NCL
>>>
>>> Thanks
>>>
>>> Debasish
>>>
>>>
>>>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 04 2007 - 13:06:59 MDT