Raccord-MLS

Fichiers de données MLS_Aura:explication des relevés
On dispose des données MLS  (données satellites) de temperature $$T(lon,lat,k)$$  fonction de :
 * longitude
 * latitude
 * niveau de pression k

Il y a 55 niveaux de pressions k:

level_hdf=[1000.0,825.4042,681.29205,562.3413,464.15887,383.11868,316.22775,261.01572,215.44347,177.82794,146.77992,121.15276,100.0,82.54042,68.1292,56.23413,46.41589,38.311867,31.622776,26.101572,21.544348,17.782795,14.677993,12.115276,10.0,8.254042,6.8129206,5.623413,4.6415887,3.8311868,3.1622777,2.6101573,2.1544347,1.7782794,1.4677993,1.2115277,1.0,0.68129206,0.4641589,0.31622776,0.21544346,0.14677992,0.1,0.046415888,0.021544347,0.01,0.004641589,0.0021544348,0.0010,4.641589E-4,2.1544346E-4,1.0E-4,4.6415887E-5,2.1544347E-5,1.0E-5]';

La température est donc relevée pendant une  journée (24 heures) lors de la trajectoire du satellite

temp_hdf: est une matrice ou chaque vecteur colonne représente les relevées de Température à k (niveau de pression) fixé. Il y a donc T(t,longitude(t),latitude(t) )) pris à chaque instant (à une latitude et longitude données par la trajectoire du satellitz à chaque instant t).

temp_hdf: est une matrice ou chaque vecteur colonne représente les relevées de Température à: k (niveau de pression) fixé. Il y a donc T(t,longitude(t),latitude(t) )) pris à chaque instant (à une latitude et longitude données par la trajectoire du satellitz à chaque instant t).



Nomenclature des données MLS
Pour plus d' infos sur la campagne MLS aura: http://mirador.gsfc.nasa.gov/collections/ML2T__003.shtml

On a un fichier MLS par jour:

Ex: MLS-Aura_L2GP-Temperature_v03-30-c01_2009d001.he5'

- _L - _v - _c_ d.

For example: MLS-Aura_L2GP-H2O_v03-30-c01_2007d059.he5 Where:
 * instrument = name of instrument = MLS.
 * satellite = name of satellite = Aura.
 * level = data processing level = 2GP (Level-2 geophysical product).
 * product = H2O, BRO, O3, Temperature, etc.
 * major-minor = algorithm version, made up of major and minor number 03-30
 * xx = file cycle number
 * yyyy = four digit calendar year
 * ddd = day number in the year (001=January 1)
 * ext = file extension, i.e. format = he5 (HDF-EOS5)

Rmq:


 * -Fevrier: jours 34 35 et 39 renomés  (C02 ou C03 remplacés par C01) pour faciliter la lecture matlab.

Carte de température MLS-Aura
On souhaite disposer d'une carte de température sur une grille régulière (en longitude et latitude). A l' aide d' une interpolation linéaire on peut obtenir une estimation de la température aux positions géographique ((longitude,latitude))souhaitées.

Données era interim
Via le serveur Climserv on peut acceder aux données erai : ce qui nous permet d'obtenir les données de température: Mis dans un dossier:
 * mensuelle (dans 1xmonthly>>...)
 * journalière s par tranches de 6 h00, ce qui nous donne des gros fivhiers (ta.2009.apheiGlobal_075.nc) qui contiennent les relevées des température toutes les 6 heures pendant 30 jours du mois.

/net/nfs/tmp15/pincemaille/donnees_erai/ Les données sont lues  dans les fichiers correspondants:

geopt.200901.aphei.GLOBAL_075.nc, ta.200901.aphei.GLOBAL_075.nc ...

Le fichier 'namelist.ERAI' va chercher ces fichiers: fname = '/net/nfs/tmp15/pincemaille/donnees_erai/{v}.{y}{m}.aphei.GLOBAL_075.nc', Remarque:  il s' agit d'un fichier mensuelle

On va donc enregistrer tous les fichiers Erai qui njous interessent (ceux de janvier/fevrier/juillet/aout) dans '/net/nfs/tmp15/pincemaille/donnees_erai'. Pour ça il faut aller dans le navigateur firefox puis faire:

>>Edit>>Preference>>General>>Save files to>>/net/nfs/tmp15/pincemaille/donnees_erai (ou autre chose: c'est l' emplacement ou l' on veut sauver nos fichiers)

Faire tourner le modèle
Par la suite à l' aide de ces données sources ERAI on peut faire tourner Raccord en spécifiant bien les paramètres dans le fichier param.dat (datemin et date max, nombre de niveaux de pression...). Ainsi on peut obtenir les données de sortiie chaque jour (ex:/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/2009-01-01/) (01 janvier: dans le bon repertpire on va chercher:tracco.200901.nc)

A- Dans le fichier param.dat:

1) remplir là ou iront les fichiers en sortie du modèle: repOutput /tmp15/pincemaille/output_raccord_v2_20121001/Janvier_2009/2009-01-30/ 2)Modifier les dates pour qu'elles conincident

B- Dans le fichier namelist.ERAI (attention

1) Mettre la date cohérente avec la date dans le fichier de paramètres

C- Faire tourner le modèle sous linux, là ou est enregistré le programme raccord.exe, lancer la comande shell avec le fichier de paramètre

(ps le modèle ira tous seul chercher le fichier namelist.ERAI donc ne pas modifier ce nom comme namelista.ERAI mais plutot copier dans un autre dossier)

%Données Raccord: tableau des chemins raccord Tab_rep_o={}; for i=1:9 Tab_rep_o={ Tab_rep_o{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ao09/','2009-08-0',num2str(i),'/']}; end for i=10:30 Tab_rep_o={ Tab_rep_o{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ao09/','2009-08-',num2str(i),'/']}; end Tab_datei={}; for i=1:9 Tab_datei={ Tab_datei{:} ['2009-08-0',num2str(i)]}; end for i=10:30 Tab_datei={ Tab_datei{:} ['2009-08-',num2str(i)]}; end Tab_datef={Tab_datei{2:end} ['2009-08-',num2str(i+1)]} % Trouver les positions grâce à: ftell(fid) puis lire fscanf(fid, '%s' ,1) cd ~/Bureau/CNRS/Logiciels/Raccord_deux_versions/raccord_v2_20121001/source_aug/ % 1) Modification du fichier param30.dat (attention c'est lui qui est % spécifique)  i=1;   for i=21:length(Tab_rep_o) %Ouverture du fichier fid=fopen('parama.dat','r+') rep_o=Tab_rep_o{i}    %rep_o='/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/2009-02-03/' fseek(fid,474,'bof') %Placement du curseur fprintf(fid,'%s',rep_o) %%Ecriture par dessus  fseek(fid,769,'bof')  date_deb=Tab_datei{i}   %date_deb='2009-01-08' %fscanf(fid, '%s' ,1)  fprintf(fid,'%s',date_deb) fseek(fid,911,'bof')  date_fin=Tab_datef{i}  % date_fin='2009-01-09' fprintf(fid,'%s',date_fin) fclose(fid); % 2) Modif du fichier namelist.ERAI fid=fopen('namelist.ERAI','r+') fseek(fid,122,'bof') fprintf(fid,'%s', date_deb) fseek(fid,159,'bof') %Placement du curseur fprintf(fid,'%s', date_fin)% fclose(fid); % B) Se placer dans le répertoire ou l' on souhaite que les fichiers soient % compilés !./raccord.exe <para % fseek(fid,122,'bof') % fscanf(fid, '%s' ,1) % ftell(fid)   end

Dans Matlab se placer là ou se trouve le programme avec:
 * Juillet 2009

-la bonne namelist et

-le bon fichier de paramètre et

-le bon fichier par /net/nfs/home/pincemai/Bureau/CNRS/Logiciels/Raccord_deux_versions/raccord_v2_20121001/source_ju

%Données Raccord: tableau des chemins raccord Tab_rep_o={}; for i=1:9 Tab_rep_o={ Tab_rep_o{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ju09/','2009-07-0',num2str(i),'/']}; end for i=10:30 Tab_rep_o={ Tab_rep_o{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ju09/','2009-07-',num2str(i),'/']}; end Tab_datei={}; for i=1:9 Tab_datei={ Tab_datei{:} ['2009-07-0',num2str(i)]}; end for i=10:30 Tab_datei={ Tab_datei{:} ['2009-07-',num2str(i)]}; end Tab_datef={Tab_datei{2:end} ['2009-07-',num2str(i+1)]} % Trouver les positions grâce à: ftell(fid) puis lire fscanf(fid, '%s' ,1) cd ~/Bureau/CNRS/Logiciels/Raccord_deux_versions/raccord_v2_20121001/source_ju % 1) Modification du fichier param30.dat (attention c'est lui qui est % spécifique)  i=1;   for i=1:length(Tab_rep_o) %Ouverture du fichier fid=fopen('paramj.dat','r+') rep_o=Tab_rep_o{i}    %rep_o='/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/2009-02-03/' fseek(fid,474,'bof') %Placement du curseur fprintf(fid,'%s',rep_o) %%Ecriture par dessus  fseek(fid,769,'bof')  date_deb=Tab_datei{i}   %date_deb='2009-01-08' %fscanf(fid, '%s' ,1)  fprintf(fid,'%s',date_deb) fseek(fid,911,'bof')  date_fin=Tab_datef{i}  % date_fin='2009-01-09' fprintf(fid,'%s',date_fin) fclose(fid); % 2) Modif du fichier namelist.ERAI fid=fopen('namelist.ERAI','r+') fseek(fid,122,'bof') fprintf(fid,'%s', date_deb) fseek(fid,159,'bof') %Placement du curseur fprintf(fid,'%s', date_fin)% fclose(fid); % B) Se placer dans le répertoire ou l' on souhaite que les fichiers soient % compilés !./raccord.exe<par % fseek(fid,122,'bof') % fscanf(fid, '%s' ,1) % ftell(fid)   end

Tracé graphique des cartes de température

 * Code source pour le tracé graphique:

%A modifier: 1)Fichier Raccord 2)Fichier MLS-Aura 3) Chemin (pour %enregistrement 4) Titre Fichier for k=46:55 %Choix desz paramètres (grilles d'interpolation...) t=1 %choix de la période (premier quart de journée (6 premieres heures ,...) vect_longitude=[-180:1:180]'; % Afin d' interpoler T(lon_hdf,lat_hdf) sur une grille régulière vect_latitude=[-90:1:90]'; % vect_longitude=lon_ncs ;vect_latitude=lat_nc %Fichiers Raccord: tableau des chemins raccord Tab_path_Raccord={}; for i=1:9 Tab_path_Raccord={ Tab_path_Raccord{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/','2009-02-0',num2str(i)]}  %Tab_path_Raccord={  Tab_path_Raccord{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/','2009-01-0',num2str(i)]}; end for i=10:27 Tab_path_Raccord={ Tab_path_Raccord{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/','2009-02-',num2str(i)]}; end file_Raccord='tracco.200902.nc' % file_Raccord='tracco.200901.nc' %Fichiers MLS-Aura Path_hdf='/net/nfs/tmp15/pincemaille/MLS_Aura/fevrier_09/' %Path_hdf='/net/nfs/tmp15/pincemaille/MLS_Aura/' Tab_file_MLS={}; % Tab_file_MLS={'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d001.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d002.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d003.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d004.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d005.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d006.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d007.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d008.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d009.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d010.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d011.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d012.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d013.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d014.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d015.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d016.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d017.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d018.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d019.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d020.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d021.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d022.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d023.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d024.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d025.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d026.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d027.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d028.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d029.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d030.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d031.he5'} for i=32:59 Tab_file_MLS= {Tab_file_MLS{:} ['MLS-Aura_L2GP-Temperature_v03-30-c01_2009d0',num2str(i),'.he5']} end %Donnees MLS et Raccord [Moyenne_hdf,Variance_hdf]=moy_mls(k,vect_longitude,vect_latitude,Path_hdf,Tab_file_MLS); [Moyenne_nc,Variance_nc]=moy_raccord(k,vect_longitude,vect_latitude,Tab_path_Raccord,file_Raccord); %Tracé de la carte Matrice_moy_hdf=Moyenne_hdf;Matrice_moy_nc=Moyenne_nc;Matrice_var_hdf=Variance_hdf;Matrice_var_nc=Variance_nc; vecteur_abcisse=vect_longitude;vecteur_ordonne=vect_latitude;level_hdf=[1000.0,825.4042,681.29205,562.3413,464.15887,383.11868,316.22775,261.01572,215.44347,177.82794,146.77992,121.15276,100.0,82.54042,68.1292,56.23413,46.41589,38.311867,31.622776,26.101572,21.544348,17.782795,14.677993,12.115276,10.0,8.254042,6.8129206,5.623413,4.6415887,3.8311868,3.1622777,2.6101573,2.1544347,1.7782794,1.4677993,1.2115277,1.0,0.68129206,0.4641589,0.31622776,0.21544346,0.14677992,0.1,0.046415888,0.021544347,0.01,0.004641589,0.0021544348,0.0010,4.641589E-4,2.1544346E-4,1.0E-4,4.6415887E-5,2.1544347E-5,1.0E-5]'; titre_fenetre_graph=['Moyenne des Temperatures (Fevrier 2009) - Niveau de Pression:',num2str(level_hdf(k)),'(HPA)'];%titre_fenetre_graph=['Moyenne des Temperatures (Janvier 2009) - Niveau de Pression:',num2str(level_hdf(k)),'(HPA)']; [f]=carte_temperature(Matrice_moy_hdf,Matrice_moy_nc,Matrice_var_hdf,Matrice_var_nc,vecteur_abcisse,vecteur_ordonne,titre_fenetre_graph) %set(f,'visible','on') %Enregistrement titre_fichier=[num2str(k),'_Pression_',num2str(level_hdf(k)),'.m.jpg'] % saveas(f,strcat(titre_dossier,'/',titre_fichier)); chemin=('/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_Fevrier/')% chemin=('/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_Janvier/') saveas(f,[chemin,titre_fichier]) end


 * Essai d' automatisation:

%A modifier: 1)Fichier Raccord (t_path_rac) 2)Fichier MLS-Aura (path_hdf, 3) Chemin (pour %enregistrement 4) Titre Fichier (>> cd ~/Bureau/CNRS/Programme_sources/) %Paramètres d' entrée: t_path_rac='/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ao09/2009-08-0'%t_path_rac='/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Ju09/2009-07-0'% t_path_rac='/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/2009-02-0' file_Raccord='tracco.200908.nc'%file_Raccord='tracco.200907.nc'% indice_fin_raccord=30 Path_hdf='/net/nfs/tmp15/pincemaille/MLS_Aura/aout_09'%Path_hdf='/net/nfs/tmp15/pincemaille/MLS_Aura/juillet_09' prem_file_hdf='MLS-Aura_L2GP-Temperature_v03-30-c01_2009d213' % prem_file_hdf='MLS-Aura_L2GP-Temperature_v03-30-c01_2009d182'%'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d032' indice_debut_jour_mls=213 indice_fin_mls= 243 titre_graphe='Moyenne des Temperatures (Août 2009) - Niveau de Pression:' titre_fenetre_graph= titre_graphe chemin_sauvegarde='/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_aout/' %'/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_Fevrier/' cd ~/Bureau/CNRS/Programme_sources/ %Début du rogramme for k=1:55 %Choix desz paramètres (grilles d'interpolation...) t=1 %choix de la période (premier quart de journée (6 premieres heures ,...) vect_longitude=[-180:1:180]'; % Afin d' interpoler T(lon_hdf,lat_hdf) sur une grille régulière vect_latitude=[-90:1:90]'; % vect_longitude=lon_ncs ;vect_latitude=lat_nc %Fichiers Raccord: tableau des chemins raccord %t_path_rac='/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/Fe09/2009-02-0'%exemple Tab_path_Raccord={};  %indice_fin_raccord=27  for i=1:9   Tab_path_Raccord={  Tab_path_Raccord{:} [t_path_rac,num2str(i)]}  %Tab_path_Raccord={  Tab_path_Raccord{:} ['/net/nfs/tmp15/pincemaille/output_raccord_v2_20121001/','2009-01-0',num2str(i)]};  end  for i=10:indice_fin_raccord       Tab_path_Raccord={  Tab_path_Raccord{:} [t_path_rac(1:end-1),num2str(i)]};  end   % file_Raccord='tracco.200902.nc' % file_Raccord='tracco.200901.nc'  %Fichiers MLS-Aura  %prem_file_hdf='MLS-Aura_L2GP-Temperature_v03-30-c01_2009d032' %indice_debut_jour_mls=32 %indice_fin_mls= 59 Tab_file_MLS={}; % Tab_file_MLS={'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d001.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d002.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d003.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d004.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d005.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d006.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d007.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d008.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d009.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d010.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d011.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d012.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d013.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d014.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d015.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d016.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d017.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d018.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d019.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d020.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d021.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d022.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d023.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d024.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d025.he5';'MLS-Aura_L2GP-Temperature_v03-30-c02_2009d026.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d027.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d028.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d029.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d030.he5';'MLS-Aura_L2GP-Temperature_v03-30-c01_2009d031.he5'} for i=indice_debut_jour_mls:indice_fin_mls Tab_file_MLS= {Tab_file_MLS{:} [prem_file_hdf(1:end-3),num2str(i),'.he5']} end %Donnees MLS et Raccord [Moyenne_hdf,Variance_hdf]=moy_mls(k,vect_longitude,vect_latitude,Path_hdf,Tab_file_MLS); [Moyenne_nc,Variance_nc]=moy_raccord(k,vect_longitude,vect_latitude,Tab_path_Raccord,file_Raccord); %Tracé de la carte Matrice_moy_hdf=Moyenne_hdf;Matrice_moy_nc=Moyenne_nc;Matrice_var_hdf=Variance_hdf;Matrice_var_nc=Variance_nc; vecteur_abcisse=vect_longitude;vecteur_ordonne=vect_latitude;level_hdf=[1000.0,825.4042,681.29205,562.3413,464.15887,383.11868,316.22775,261.01572,215.44347,177.82794,146.77992,121.15276,100.0,82.54042,68.1292,56.23413,46.41589,38.311867,31.622776,26.101572,21.544348,17.782795,14.677993,12.115276,10.0,8.254042,6.8129206,5.623413,4.6415887,3.8311868,3.1622777,2.6101573,2.1544347,1.7782794,1.4677993,1.2115277,1.0,0.68129206,0.4641589,0.31622776,0.21544346,0.14677992,0.1,0.046415888,0.021544347,0.01,0.004641589,0.0021544348,0.0010,4.641589E-4,2.1544346E-4,1.0E-4,4.6415887E-5,2.1544347E-5,1.0E-5]'; %titre_graphe='Moyenne des Temperatures (Fevrier 2009) - Niveau de Pression:' %titre_fenetre_graph=[titre_graphe,num2str(level_hdf(k)),'(HPA)'];%titre_fenetre_graph=['Moyenne des Temperatures (Janvier 2009) - Niveau de Pression:',num2str(level_hdf(k)),'(HPA)']; [f]=carte_temperature(Matrice_moy_hdf,Matrice_moy_nc,Matrice_var_hdf,Matrice_var_nc,vecteur_abcisse,vecteur_ordonne,titre_fenetre_graph) %set(f,'visible','on') %Enregistrement titre_fichier=[num2str(k),'_Pression_',num2str(level_hdf(k)),'.m.jpg'] % saveas(f,strcat(titre_dossier,'/',titre_fichier)); %chemin_sauvegarde='/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_Fevrier/' chemin=(chemin_sauvegarde)% chemin=('/net/nfs/tmp15/pincemaille/output_matlab/Moyenne_Janvier/') saveas(f,[chemin,titre_fichier]) end

Commandes Shell utilisées
1) Pour copier tous les fichiers jpg du serveur (pincemai@adresse_Ip) sur  mon pc:

Rmq: Pour copier des fichiers qui sont sur un serveur donné (par exemple sur le compte pincemai à jussieu donc sur sirroco) vers son pc perso. En shell, on se place dans le dossier ou l' on veut copier les fichiers  (sur son pc perso ) puis on execute la commande: scp pincemai@194.254.42.104:/net/nfs/...adresse_ou_sont_stockés_fichiers_désirés/fichier.extension. Avec 194.254.42.104 adresse ip du serveur distant. pincemai@eole /tmp15/hauche> hostname -i 194.254.42.104 pincemai@eole /tmp15/hauche> scp * pincemai@194.254.42.104:/net/nfs/home/pincemai/Bureau/CNRS/Mimosa Password: Je suis su eole (à distance donc en ssh) dans le répertoire hauche/ et je veux copier tous les fichiers qui y sont (commande shell:*) sur le eme serveur
 * scp pincemai@134.157.176.57:/net/nfs/tmp15/pincemaille/output_matlab/Essai_paslarge/*.jpg.

2)Pour copier des fichiers sur le serveur (je suis depuis mon pc)


 * scp aire.png pincemai@134.157.176.57:/net/nfs/tmp15/pincemaille/output_matlab/Essai_paslarge

3)Se connecter au serveur (en mode graphique, rajouter -X)
 * ssh pincemai@134.157.176.57
 * ssh -X pincemai@sirocco.aero.jussieu.fr
 * ssh -X foehn
 * ssh -X 3104595@ssh.ufr-info-p6.jussieu.fr

4) Obtenir une adresse ip


 * hostname -i

5) Copier des données flux de sortie

history >> hist.txt

6)  Retrouver ou est installé un programme donné par une certaine commande $ type -p commande (ex type -p okular)

Site utiles
Pour s'inspirer de script matlab:

http://www.mathworks.com/matlabcentral/fileexchange/index?page=2&term=colorbar&utf8=%E2%9C%93

Pour aller chercher données MLS

http://mirador.gsfc.nasa.gov/cgi-bin/mirador/presentNavigation.pl?tree=project&dataset=ML2T.003&project=MLS&dataGroup=L2_V003&version=003

Commandes Matlab utiles
Pour avoir de l' aide:


 * doc

Pour utiliser les commmandes linux
 * help pcolor
 * !gedit param.dat

Tentative d' automatiser raccord:
(voir en code source) ===%Modifier param.dat,par, param2.dat !raccord.exe param30.dat

fid=fopen('test.txt','r') Dir = fscanf(fid, '%s' ,1) Dir2 = fscanf(fid, '%s' ,1) dates = fscanf(fid, '%g' ,6) fclose(fid);

%Ouverture du fichier fid = fopen('test.txt','r+');

%Placement du curseur fseek(fid,1,'bof')  %placé en faisant des teste (50 Octés ?!?)

%écriture par dessus eval='2010165' fprintf(eval) ftell

strmatch fclose(fid);===