Fonction squeeze et fichier netcdf

Bonjour à tous,
Je voudrais plotter un fichier netcdf, ci-après vous trouverez les variables que je veux plotter et leurs dimensions (30 correspond au nombre de lignes de la matrice, 38 et 39 c'est les colonnes et normalement, il y a 42/43 niveaux de profondeurs :

Citation:

---

double std_u(s_rho,eta_u,xi_u), shape = [30 39 42]
:long_name = "standanrd deviation of u-momentum component"
:units = "meter second-1"
double std_v(s_rho,eta_v,xi_v), shape = [30 38 43]
:long_name = "standanrd deviation ofv-momentum component"
:units = "meter second-1"

---

Je voudrais plotter seulement le niveau 5 (correspondant à 5 m de profondeur).
Pour ce faire, j'ai utilisé la fonction squeeze et je ne sais pas si c'est correcte...
ci-après le code :

Code:

---

1 2 3 4 5 6 7 8 9 10 11

x = nc_varget('std_u_v_2006.nc','std_u'); y = nc_varget('std_u_v_2006.nc','std_v'); janx=squeeze(x(5,:,:)); jany=squeeze(y(5,:,:)); janxr=janx(1:30,2:37); janyr=jany(1:30,2:37); janxr(janxr==1e37)=NaN; janyr(janyr==1e37)=NaN;   quiver(janxr,janyr, 6) contourf(janxr,8);

---

Je voudrais savoir si c'est correcte ou est ce que j'ai tout faux :(

Ne serait-ce pas plutôt ceci ?

Code:

---

1 2	janx = x(:,:,5); jany = y(:,:,5);

---

Bonjour DUT,
Tu trouveras, ci-joint le fichier que j’obtiens avec mon code initial (IM1). Si je change, et c'est logique de mettre la troisième dimension en troisième lieu : le résultat ne ressemble plus à la réalité (voir IM2)
PS : J'ai mis les outputs ainsi que le fichier netcdf (pas volumineux) en pièce jointe

Pièce jointe 131829

Pièce jointe 131830

Pièce jointe 131831

Je pense que la réponse se trouve ci-dessous :

Code:

---

1 2 3 4 5

>> x = ncread('std_u_v_2006.nc','std_u'); >> whos Name Size Bytes Class Attributes x 42x39x30 393120 double

---

Merci beaucoup Dut pour ton aide si précieuse!!!
Bonne fin d'après midi :ccool: