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