ascii data

Hi,

I've just intalled the NCL for x86_64 and I have a problem reading and
writing ascii data.
I just read a file (asciiread) and then write the data in a new file.
So, the new data are not equal to the data on the original file.

I read this:

... ... ...
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
... ... ...

and after write the data in a new file:

... ... ...
-999.000 9.000 -9999.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -999.000 9.000 -9999.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -999.000 9.000 -9999.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -999.000 9.000 -9999.000 -9999.000
-9999.000 -9999.000 -9999.000 -9999.000 -999.000 9.000 -9999.000
... ... ...

Some -9999 data are broken in two data: a -999 and a 9.
The script that I used are:

begin

matrix = asciiread("file_name.asc", (/360, 720/), "float")

opt = True
opt_at_fout = "new_file.asc"

write_matrix (matrix, "720f10.3", opt)

                
end

Help, please!

Thanks,

Varejão

--
######################################
Claudeci Gomes Varejão Júnior
Computer Science
Federal University of Viçosa
Viçosa    Minas Gerais    Brasil
http://www.dpi.ufv.br
http://madeira.dea.ufv.br
######################################
       "Cuidado com o stress....é melhor chegar atrasado nesta vida,
do que adiantado na outra!"
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk