Qualcuno usa NCL?

[ugo1]

Ciao,

rieccomi

dopo esser riuscito a creare una mappa dai shapefile ho cominciato fare i primi grafici da dati reali ed ho creto il seguente script

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


begin


url = "http://nomad1.ncep.noaa.gov:9090/dods/gfs_master/gfs20110417/"
filename = url + "gfs_master_00z"

exists = isfilepresent(filename)
if(.not.exists) then
print("OPeNDAP isfilepresent test unsuccessful.")
print("Either file doesn't exist, or NCL does not have OPeNDAP capabilities on this system")
else
print("OPeNDAP isfilepresent test successful.")
a = addfile(filename,"r")
vars = getfilevarnames(a)
end if

temp = a->tmp2m(:,48,9) ; temperature
temp@units = "(C)"
temp = temp - 273.15
print(temp)

;--- Open workstation
wks_type = "png" ; use "newpdf" instead of "pdf" for smaller files
wks = gsn_open_wks(wks_type,"naturalearth")

;---Set some resources for the map.
res = True
res@cnMonoLineColor = True ; Turn off the drawing of contours lines in one color.
res@tiMainString = "Temperature (C)" ; Create a title.
field1 = create "field1" coordArraysClass defaultapp
"caXArray": (temp!0)
"caYArray": (temp)
end create
plot = gsn_contour(wks,field1,True) ; Draw a contour plot.
end

il problema che non mi crea i dati dove sbaglio le variabili sono correttamente lette (verificato la funzione print).

le caratteristiche della variabile "tmp2m" si trovano qui http://nomad1.ncep.noaa.gov:9090/dod...master_00z.dds

Grazie

Ugo

[ugo1]

Ciao,

Ho un altro problema. Le linee di livello di NCL sono poco arrotondate rispetto a quelle di Grads (ovvamente i dati sono stati presi dalla stessa fonte)

italia grads http://i103.photobucket.com/albums/m136/ugo73/t2.png


italia NCL http://i103.photobucket.com/albums/m...learth_con.png

Sugerimenti sistemare NCL?

Ugo

[Danyastoria]

Scusate, approfitto della discussione già aperta per chiedervi una cosa: Qualcuno ha una vaga idea di come si usi ncl ncar?
Io non sono in grado di usarlo. ho provato a trascrivere qualche esempio già prestampato (aprire un file .txt) File Handling: Reading ASCII data
ma non funziona comunque. Qualcuno mi potrebbe aiutare??In che cartella va messo lo script dell'esempio? e il file di testo che deve leggere? e come lancio il programma? ng4ex nome_programma -clean ??

[ugo1]

Scusate, approfitto della discussione già aperta per chiedervi una cosa: Qualcuno ha una vaga idea di come si usi ncl ncar?
Io non sono in grado di usarlo. ho provato a trascrivere qualche esempio già prestampato (aprire un file .txt) File Handling: Reading ASCII data
ma non funziona comunque. Qualcuno mi potrebbe aiutare??In che cartella va messo lo script dell'esempio? e il file di testo che deve leggere? e come lancio il programma? ng4ex nome_programma -clean ??

Ciao,

ti passo in pm un link da cui scaricare uno script con materiale extra (shapefiles)

Ugo