[FR] question à propos des objets pp (pour trouver le champ tracé)

[aymeric-spiga]

(question par S. Fromang)

Ma question est la suivante: je travaille dans ipython. Lorsque j'appelle le script, je me recupère donc un objet (mettons "dpsi") en retour de l'appel, en plus du plot. Si j'appelle la méthode plot() sur cette objet, il me refait le plot. Mais si j'essaie de visualiser ce qui est stocké dans dpsi.f, ca ne correspond pas. Voila la sequence:

dpsi=anomaly.anomaly(var='psi',tmin=0.,tmax=0.,tdiff=15.,zlevel=0)

dpsi.plot()

imshow(dspi.f)

C'est l'appel a imshow(dspi.f) qui pose problème. Pourquoi?

Du coup, je ne comprends pas ou est stockée l'info sur le champ que je veux visualiser. Il s'en sort lors de l'appel a la méthode plot, mais manifestement l'info n'est pas dans dpsi.f.

[aymeric-spiga]

C'est parce que ton usage devient assez avancé : le recours aux objets pp de ppclass suppose que tu veux les fonctionnalités de lecture netCDF + plot, un package intégré en quelque sorte

ton objet dpsi est un objet pp

il est toujours possible de faire ce que tu dis mais pas avec dpsi.f, ce serait plutôt un machin genre

dpsi.p[n].f

pour aller chercher l'objet plot qui est dans ton objet pp, avec n l'index du plot (par exemple si tu ne requiers qu'un seul champ, un seul fichier, etc, c'est n=0)

En fait tu vas te rendre compte que tes capacités sont limitées. A un moment donc, sans perdre de l'utilité de ppclass, tu peux séparer

• la lecture de ton netCDF : utiliser .getf() ou .getfd() les méthodes de ppclass pour choper direct le champ netCDF yeah = pp(file="truc.nc",var="u",z=0).getf()
• le plot : définir un objet pplot que tu peux customiser après pl = ppplot.plot2d() pl.f = yeah pl.makeshow()

Il y a des scripts exemple quelque part, tous mes scripts en fait. Ta question était peut être juste de la curiosité, mais comme ça j'en dis plus !