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-talkReceived on Sun Nov 20 2005 - 20:03:16 MST
This archive was generated by hypermail 2.2.0 : Mon Nov 21 2005 - 10:28:38 MST