Discussion :

[Maffalda]

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 :

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Jerome Briot]

Ne serait-ce pas plutôt ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Maffalda]

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





std_u_v_2006.7z

[Jerome Briot]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
>> x = ncread('std_u_v_2006.nc','std_u');
>> whos
  Name       Size                Bytes  Class     Attributes

  x         42x39x30            393120  double

[Maffalda]

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