Changeset 640
- Timestamp:
- 04/16/14 16:03:12 (10 years ago)
- Location:
- trunk
- Files:
-
- 144 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/adm/Doxyfile_varamma
r112 r640 30 30 # 31 31 # (remove the last "\") 32 #33 32 # 34 33 # TODO -
trunk/adm/makefile
r518 r640 338 338 @mkdir -p $(DIRTMP)/test/ 339 339 340 341 340 install_all : \ 342 341 install_html \ -
trunk/docs/docs_dev/Makefile
r628 r640 24 24 # moins bourrin dans l'usage de makefile_non_sphinx et la génération du 25 25 # logo 26 #27 26 # 28 27 # EVOLUTIONS -
trunk/docs/docs_dev/source/conf.py
r639 r640 224 224 htmlhelp_basename = 'VARAMMAdoc' 225 225 226 227 226 # -- Options for LaTeX output -------------------------------------------------- 228 227 … … 256 255 # If false, no module index is generated. 257 256 #latex_domain_indices = True 258 259 257 260 258 # -- Options for manual page output -------------------------------------------- -
trunk/docs/docs_dev/source/guides/contents.rst
r638 r640 40 40 contributing 41 41 42 43 42 Indices and tables 44 43 ================== -
trunk/docs/docs_dev/source/guides/data_content.rst
r636 r640 243 243 :file:`ftp://trmmopen.gsfc.nasa.gov/pub/merged/combinedMicro/V5/{yyyy}/` 244 244 245 246 247 245 Terminology 248 246 +++++++++++ 249 247 250 248 :file:`3B40RT.{yyyy}{mm}{dd}{hh}.6.bin.gz` 251 252 249 253 250 Size … … 670 667 taille : 158Mo pour classif et height, 313Mo pour tb108 671 668 672 673 669 terminologie 674 670 :file:`${PROJECT_ID}/MSG/{yyyy}/{mm}/{dd}/{yyyymmddhhmn}_msg-{product}-map_15min.nc` … … 711 707 712 708 .. _data_epsat: 713 714 709 715 710 EPSAT data -
trunk/docs/docs_dev/source/guides/glossary.rst
r638 r640 7 7 .. 8 8 .. $URL$ 9 ..10 9 .. 11 10 .. - fplod 20110706T123647Z cratos.locean-ipsl.upmc.fr (Linux) -
trunk/docs/docs_dev/source/guides/html.xsl
r30 r640 21 21 22 22 <xsl:import href="http://docbook.sourceforge.net/release/xsl/1.75.2/xhtml/docbook.xsl"/> 23 24 23 25 24 <xsl:output method="xml" -
trunk/docs/docs_dev/source/guides/main_many01_html.xsl
r32 r640 29 29 $Id$ 30 30 31 - fplod 2009-05-15T10:42:45Z aedon.locean-ipsl.upmc.fr (Darwin) 31 32 32 - fplod 2009-05-15T10:42:45Z aedon.locean-ipsl.upmc.fr (Darwin)33 33 * creation 34 34 -
trunk/docs/docs_dev/source/guides/makefile
r89 r640 23 23 # 24 24 # $Id$ 25 #26 25 # 27 26 # - fplod 20100603T134510Z aedon.locean-ipsl.upmc.fr (Darwin) … … 126 125 @echo " " 127 126 128 129 127 before : 130 128 @-mkdir -p ${DIRWWW}/html/css/ … … 227 225 $(DIRWWW)/pdf/en/sphinx/ 228 226 229 230 227 tracwiki : \ 231 228 $(DIRWIKI)/WikiStart -
trunk/docs/docs_dev/source/guides/navigation02_html.xsl
r27 r640 168 168 </xsl:if> 169 169 170 171 170 </xsl:element> 172 171 … … 247 246 </xsl:if> 248 247 249 250 248 </xsl:element> 251 249 -
trunk/docs/docs_dev/source/guides/read_me.rst
r636 r640 19 19 .. 20 20 .. - 21 22 21 23 22 .. _read_me: -
trunk/docs/docs_dev/source/guides/requirements.rst
r638 r640 86 86 .. _NCO: http://nco.sourceforge.net/ 87 87 88 89 88 Sphinx 90 89 ------ -
trunk/docs/docs_dev/source/guides/rundemotest.rst
r633 r640 113 113 To run them: 114 114 115 116 115 .. code-block:: bash 117 116 -
trunk/docs/docs_dev/source/guides/sphinx/conf.py
r634 r640 87 87 #modindex_common_prefix = [] 88 88 89 90 89 # -- Options for HTML output --------------------------------------------------- 91 90 … … 161 160 htmlhelp_basename = 'varammadoc' 162 161 163 164 162 # -- Options for LaTeX output -------------------------------------------------- 165 163 -
trunk/docs/docs_dev/source/guides/style.css
r636 r640 240 240 } 241 241 242 243 242 /* Mise en évidence : texte légÚrement plus grand */ 244 243 /* Highlights: slightly larger texts */ … … 318 317 } 319 318 320 321 319 /* Listes de renvois */ 322 320 /* Callout lists */ -
trunk/docs/docs_dev/source/guides/svnbasiccommands.xml
r635 r640 17 17 <!-- 18 18 19 20 19 /usr/home/fplod/incas/varamma/varamma_ws/adm/guides/svnbasiccommands.xml sur aedon.locean-ipsl.upmc.fr 21 20 -
trunk/docs/docs_dev/source/guides/tracwiki.xsl
r3 r640 23 23 et http://www.sagehill.net/docbookxsl/CustomDb5Xsl.html 24 24 car 1re xsl pour db5 25 25 26 --> 26 27 27 28 28 <!-- -
trunk/serieAMSU.pro
r634 r640 101 101 saveimage,'sorties/wavelet-2005-2007'+canal+'dom'+dom+'.png',/png 102 102 103 -
trunk/src/EPSATbuildfullfilename.m
r636 r640 54 54 % test if file exists and handling error 55 55 % 56 % version in fil name is not always v3.1-02c56 % version in filename is not always v3.1-02c 57 57 % 58 58 % EVOLUTIONS -
trunk/src/EPSATread.m
r635 r640 1 1 function [PluieEPSAT, Temps, Longitudes, Latitudes] = EPSATread(ncfilesset, fullfilename) 2 2 3 %EPSAT lit le fichier NetCDF de l'algorithme EPSAT du jour souhaité3 %EPSAT lit le fichier netCDF de l'algorithme EPSAT du jour souhaité 4 4 %pour en sortir les variables (temps, longitudes et latitudes) ainsi que 5 5 %les données de l'estimation des précipitations. … … 19 19 % - Données d'entrées : 20 20 % 21 % * ncfilesset (en caractÚres) : Type de fichier NetCDF que l'on cherche,21 % * ncfilesset (en caractÚres) : Type de fichier netCDF que l'on cherche, 22 22 % 'normal' ou 'extracted'. 23 23 % * fullfilename : Chemin d'accÚs au fichier pour la date désirée. … … 25 25 % - Données de sorties : 26 26 % 27 % * Temps : Variable de temps du fichier NetCDF.27 % * Temps : Variable de temps du fichier netCDF. 28 28 % * Longitudes : Variable de longitudes. 29 29 % * Latitudes : Variable de latitudes. … … 31 31 % les dimensions sont respectivement, temps, longitudes et latitudes. 32 32 % 33 % lecture du fichier NetCDF de l'algorithme EPSAT du jour souhaité pour en33 % lecture du fichier netCDF de l'algorithme EPSAT du jour souhaité pour en 34 34 % sortir les variables (temps, longitudes et latitudes) ainsi que les 35 35 % données de l'estimation des précipitations. … … 91 91 % - jaclod 2011-07-21 92 92 % 93 % * Lecture des fichiers 'extracted' sur Matlab Cratos.93 % * Lecture des fichiers 'extracted' sur Matlab cratos. 94 94 % 95 95 % - jaclod 2011-07-13 … … 116 116 % * permute array after reading 3D data if netcdf.getVar was used 117 117 % * close netcdf opened file 118 % * Plu eEPSAT int -> double118 % * PluieEPSAT int -> double 119 119 % 120 120 % - fplod 20110304T105348Z aedon.locean-ipsl.upmc.fr (Darwin) -
trunk/src/LectureLAIHDF.pro
r636 r640 105 105 ; - modifié par L.Eymard le 25/07/2010 106 106 ; 107 ; * pour l'adapter & l'environnement local (gestion des entr ees - sorties)107 ; * pour l'adapter & l'environnement local (gestion des entrées - sorties) 108 108 ; modification du formattage du fichier ASCII de sortie pour 109 ; ecriture d'une ligne pour chaque "printf" (et non de plein de109 ; écriture d'une ligne pour chaque "printf" (et non de plein de 110 110 ; lignes de 80 colonnes...) 111 111 ; … … 220 220 ; 221 221 ;*** extraction 222 ;*** attention, on compte lig/col apartir de 0222 ;*** attention, on compte lig/col à partir de 0 223 223 FOR ilig_MODIS=YDIM*v, YDIM*(v+1)-1 DO BEGIN 224 224 lat[*,ilig_MODIS-YDIM*v] = 90.D - (ilig_MODIS + 0.5D)*180./NB_LIG_MODIS -
trunk/src/LectureLAI_2000_FM.pro
r636 r640 37 37 PRO LectureLAI_2000_FM, PLOT=plot 38 38 39 40 39 FORWARD_FUNCTION LectureLAI 41 42 40 43 41 ;*** constantes … … 55 53 lonmin=-17.55 & lonmax=-15.70 56 54 57 58 55 ;*** declarations 59 56 lat=MAKE_ARRAY(XDIM,YDIM,/FLOAT,VALUE=!VALUES.F_NAN) 60 57 lon=MAKE_ARRAY(XDIM,YDIM,/FLOAT,VALUE=!VALUES.F_NAN) 61 62 58 63 59 ;*** initialisation environnement graphique … … 68 64 !P.BACKGROUND=255 69 65 70 71 66 ;*** version -> HDF4 72 67 id=HDF_OPEN(file) … … 75 70 HDF_CLOSE, id 76 71 77 78 72 ;*** lecture du LAI 79 73 LAI=LectureLAI(file) 80 74 81 82 75 ;*** extraction 83 FOR ilig_MODIS=YDIM*v, YDIM*(v+1)-1 DO BEGIN ;*** attention, on compte lig/col apartir de 076 FOR ilig_MODIS=YDIM*v, YDIM*(v+1)-1 DO BEGIN ;*** attention, on compte lig/col à partir de 0 84 77 lat[*,ilig_MODIS-YDIM*v] = 90.D - (ilig_MODIS + 0.5D)*180./NB_LIG_MODIS 85 78 FOR jcol_MODIS=XDIM*h, XDIM*(h+1)-1 DO BEGIN … … 98 91 taille=size(lai_extract) 99 92 100 101 93 ;*** affichage graphique 102 94 IF KEYWORD_SET(plot) THEN BEGIN … … 105 97 ;stop 106 98 ENDIF 107 108 109 99 110 100 ;*** boucle de lecture des fichiers LAI … … 153 143 ENDFOR 154 144 ENDIF 155 156 145 157 146 ;*** enregistrement ascii -
trunk/src/MSGread.m
r636 r640 68 68 % handling of missing values 69 69 % 70 % check unit of Temps "HOURS since 1960-01-01 00:00:00" ; in net cdf70 % check unit of Temps "HOURS since 1960-01-01 00:00:00" ; in netCDF 71 71 % 72 72 % homegenize with datajul_value, datestr_value … … 102 102 % - jaclod 2011-07-21 103 103 % 104 % * Lecture des fichiers 'extracted' sur Matlab Cratos (Attention! Comme104 % * Lecture des fichiers 'extracted' sur Matlab cratos (Attention! Comme 105 105 % pour EPSAT, ncdump nous donne temps, longitudes, latitudes comme 106 106 % dimensions pour la matrice de données qui s'avÚre être de dimensions … … 185 185 case 'extracted' 186 186 187 % Dans ce cas, longitudes et latitudes sont invers es.187 % Dans ce cas, longitudes et latitudes sont inversés. 188 188 Temps = netcdf.getVar(ncid,3); 189 189 Longitudes = netcdf.getVar(ncid,2); -
trunk/src/affiche_EOF_PC_3D.m
r635 r640 113 113 end 114 114 115 116 115 [U,S,V] = svd(mat_2D_ssnan',0); %' 117 116 clear mat_2D_ssnan … … 120 119 pc_var_exp=(diag_s.^2/sum(diag_s.^2))*100; 121 120 std_v=nstd(V,0,1); 122 123 121 124 122 for(i=1:nb_eof) … … 173 171 end 174 172 ajoute_date_fig(gcf,3); 175 176 177 173 178 174 for(i=1:3) -
trunk/src/affiche_EOF_PC_3Dsscycle.m
r635 r640 93 93 end 94 94 95 96 95 [U,S,V] = svd(mat_2D_ssnan',0); %' 97 96 %clear mat_2D_ssnan … … 101 100 std_v=nstd(V,0,1); 102 101 103 104 102 for(i=1: 1) 105 103 PC(:,i)= V(:,i)/std_v(i)*fac(i); … … 111 109 %if(nargin<5) 112 110 mat_2D_ssnan=mat_2D_ssnan-PC(:,1)*(EOF(:,1))';%' 113 114 111 115 112 [U,S,V] = svd(mat_2D_ssnan',0); %' … … 119 116 pc_var_exp=(diag_s.^2/sum(diag_s.^2))*100; 120 117 std_v=nstd(V,0,1); 121 122 118 123 119 for(i=1: 3) … … 167 163 end 168 164 ajoute_date_fig(gcf,3); 169 170 171 165 172 166 for(i=1:3) -
trunk/src/affiche_eofpc_3Drev.m
r635 r640 83 83 end 84 84 85 86 85 [U,S,V] = svd(mat_2D_ssnan',0); %' 87 86 clear mat_2D_ssnan … … 90 89 pc_var_exp=(diag_s.^2/sum(diag_s.^2))*100; 91 90 std_v=nstd(V,0,1); 92 93 91 94 92 for(i=1: 3) … … 129 127 end 130 128 131 132 129 if(anomalies) 133 130 title('std of the field') … … 137 134 ajoute_date_fig(gcf,3); 138 135 139 140 141 136 for(i=1:3) 142 137 subplot(3,2,i+1) -
trunk/src/amsu2ncdf.pro
r636 r640 111 111 ; amsu2ncdf, use_amsua, yyyy, mm, dd, lon_min, lon_max, lat_min, lat_max 112 112 ; 113 ;114 113 ; SEE ALSO 115 114 ; ======== … … 129 128 ; ==== 130 129 ; 131 ;132 130 ; EVOLUTIONS 133 131 ; ========== 134 ;135 132 ; 136 133 ;- … … 140 137 ; 141 138 @cm_project 142 143 139 144 140 testfilename='' -
trunk/src/amsu2ncdf_bis.pro
r637 r640 128 128 ; ==== 129 129 ; 130 ;131 130 ; EVOLUTIONS 132 131 ; ========== 133 132 ; 134 ;135 133 ;- 136 134 PRO amsu2ncdf, numch, yyyy, mm, dd, lon_min, lon_max, lat_min, lat_max, resol … … 139 137 ; 140 138 @cm_project 141 142 139 143 140 testfilename='' -
trunk/src/anaclimecm.pro
r637 r640 163 163 164 164 saveimage,'sorties/'+chp1+'climatology-saison.png',/png 165 166 165 167 166 ;lecture plusieurs mois pour hovmöller : champ de surface … … 250 249 saveimage,'sorties/'+chp+'smooth04-09hovxclimSahara-2meanlevels.png',/png 251 250 252 253 251 ;géopotentiel : delta entre niveau 700 (7) et 925 hPa 254 252 no=2 … … 293 291 ;saveimage,'sorties/'+chp+'smooth-2005-200704-09hovy10W5E700.png',/png 294 292 295 296 297 293 ;r (q relative) - moyenne deux niveaux adjacents 298 294 no=0 -
trunk/src/box5x5serietemp.m
r636 r640 94 94 %%%%%%%%%%%%%%%%% 95 95 96 97 96 %%%%%%%%%%%%%%%%%%%%%%%%%%%%% 98 97 % Moyennage pour chaque zone … … 118 117 119 118 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 120 121 119 122 120 %%%%%%%%%%%%%%%%%%%%% -
trunk/src/cal_hcl.m
r636 r640 152 152 [tps_tmp,lat,lon,msl_tmp]=extract_nc(sprintf('/bdd/OPERA/NETCDF/WESTAFR_025/4xdaily/AN_SF/2006/msl.2006%02d.ash.WESTAFR_025.nc',m_min+i),{'msl'},'time',NaN,tps_max_mod,'lat',latmin,latmax,'lon',lonmin,lonmax); 153 153 154 155 154 tps=[tps;tps_tmp]; 156 155 ta=[ta;ta_tmp]; -
trunk/src/choix_mois.m
r635 r640 12 12 % compris entre mois_deb et mois_fin inclus. ex: pour prendre d'avril à 13 13 % aout, on choisit mois_deb=4,mois_fin=8 , de décembre à fevrier: 14 % mois_deb=12,mois_fin=2, toute l'annee (par d efaut): mois_deb=1,mois_fin=1214 % mois_deb=12,mois_fin=2, toute l'annee (par défaut): mois_deb=1,mois_fin=12 15 15 % 16 16 % EVOLUTIONS -
trunk/src/classes_principales.tex
r636 r640 96 96 \end{figure} 97 97 98 99 98 \begin{figure} 100 99 \includegraphics[width=0.8\textwidth]{lai_2008_ferlo_vsavane_boisee.jpeg} … … 102 101 \end{figure} 103 102 104 105 103 \begin{figure} 106 104 \includegraphics[width=0.8\textwidth]{lai_2008_ferlo_vsavane_herbacee_et_arbustive.jpeg} 107 105 \caption{savane herbacée et arbustive 2008} 108 106 \end{figure} 109 110 107 111 108 Dans la classe de végétation nommée savane, les profils du LAI qui traduisent le cycle phénologique annuel des différentes classes de végétation sont trÚs variés c'est à dire la date ou bien la période des phases change d'un type de végétation à l'autre. Une étude comparative entre les classes montre qu'avec la savane herbacée arbustive la réponse de la végétation est plus rapide par rapport aux autres types de savane et que la date de début de la verdure est avancée que celle de certaine et en retard par rapport à d'autres. Par contre les savanes arbustive et boisée répond de maniÚre trÚs lente avec des phases de croissance et senescence trÚs courte c'est à dire le maximum correspondant à la maturité est atteint plus tÃŽt. Ces dates de début de la verdure et de fin de la senescence sont calculées au 10 \% de l'amplitude totale. -
trunk/src/compar_tbnoyau_disp.m
r635 r640 136 136 %- 137 137 138 139 138 usage='[Annee, Mois, Jour] = compar_tbnoyau_disp(Annee, Mois, Jour, lonmin, lonmax, latmin, latmax, Tps_SeuilInf, Tb_SeuilSup_ini, dTb_SeuilSup, Tb_LimNoyau, nbclasses)'; 140 139 if nargin~=nargin(mfilename) … … 173 172 174 173 %%%%%%%%%%%%%%%%%%%% 175 176 174 177 175 %%%%%%%%%%%%%%%%%%% -
trunk/src/compar_tbseuil_histo.m
r636 r640 176 176 177 177 %%%%%%%%%%%%%%%%%%%% 178 179 178 180 179 %%%%%%%%%%%%%%%%%%% -
trunk/src/compar_tbseuil_image.m
r636 r640 172 172 173 173 %%%%%%%%%%%%%%%%%%%%%%%% 174 175 174 176 175 %%%%%%%%%%%%%%%%%%%%%%%%%%%% -
trunk/src/correc_amsu.pro
r637 r640 68 68 s = size(SCOPE_TRACEBACK(/STRUCTURE),/DIMENSION) 69 69 routine = (SCOPE_TRACEBACK(/STRUCTURE))[s-1].Routine 70 ;71 ;72 70 ; 73 71 nparam = N_PARAMS() … … 182 180 errplot,moytb-sigtb,moytb+sigtb 183 181 184 185 182 print,moytb 186 183 -
trunk/src/correct_nadir_amsu-v2.pro
r635 r640 93 93 ; :Examples: 94 94 ; 95 ;96 95 ; canal='a5' 97 96 ; resol=1 … … 109 108 ; déplacement de l'interpolation des points dans la fauchée de 110 109 ; extract a correct_nadir 111 ; définition de la zone et ecriture du fichier "amélioré"(!)110 ; définition de la zone et écriture du fichier "amélioré"(!) 112 111 ; 113 112 ; $Id: correct_nadir_amsu.pro 436 2011-11-25 17:48:35Z lelod $ … … 185 184 tbmin=100 186 185 tbmax=350 187 188 186 189 187 ; lecture fichier land - sea (S. Masson) … … 204 202 jpie = n_elements;module lecture fichier a introduire ici 205 203 206 207 208 204 ; lecture fichier jour 209 205 ; décodage header 210 206 ; 211 207 ; à modifier en fonction du module de lecture 212 213 208 214 209 month=tab[0,*] … … 256 251 ind=where(landseamask eq 0) 257 252 258 print, '#correct_nadir_amsu# ecriture du fichier resultant ', fichier2253 print, '#correct_nadir_amsu# écriture du fichier resultant ', fichier2 259 254 260 255 for i=0L, nn-1L do begin … … 263 258 endif 264 259 endfor 265 266 267 260 268 261 ; appel à interpolswatw pour ajuster les pixels amsua sur une grille … … 299 292 endfor 300 293 301 302 294 ; écriture fichier sortie avec nbgrid valeurs dans la fauchée 303 304 305 295 306 296 printf,lun_obsa, month[i],landseamask[i],yrday[i],orbit[i],hour[i],fov[i],nosat[i],lon[i],lat[i],zen[i],$ … … 316 306 endfor 317 307 318 319 308 end -
trunk/src/correct_nadir_amsu.pro
r636 r640 80 80 ; 81 81 ; :Todo: 82 ;83 82 ; 84 83 ; vérifier la correction nadir avec comparaison de traces des orbites … … 218 217 ENDIF 219 218 220 221 219 for i=0, 8 do begin 222 220 readf, lun1, t1,t2,t3,t4,t5,$ … … 312 310 interpol_correc,cor_lb,cor_sb 313 311 endelse 314 315 312 316 313 ; fichier land - sea (S. Masson) … … 530 527 endfor 531 528 ;;correction des valeurs aberrantes 532 ;et ecriture du fichier resultant529 ;et écriture du fichier résultant 533 530 ; 534 531 ind=where(landseamask eq 0) … … 536 533 ; oplot,fov[ind],ach5_adj[ind],psym=1 537 534 538 539 print, '#correct_nadir_amsu# ecriture du fichier resultant ', fichier2 535 print, '#correct_nadir_amsu# écriture du fichier résultant ', fichier2 540 536 541 537 openw, lun_obsa, fichier2, /get_lun, ERROR=error -
trunk/src/create_gif_msg-temperature.jnl
r635 r640 61 61 LABEL 10,37,-1,0,0.12 @AR Source: EUMETSAT-CMS-IPSL 62 62 63 64 63 can reg 65 64 !! Attention cette routine redéfinit le viewport et donc tout ce qui est aprÚs ne marche pas sur -
trunk/src/cresamsu.pro
r636 r640 195 195 time1 = SYSTIME(1) 196 196 ENDIF 197 ;198 197 ; 199 198 ; Return to caller if errors … … 271 270 nbjour=jmax+1 272 271 273 274 272 ; conversions degré / radians et rayon en km 275 273 vbidon=-999. … … 295 293 msg = report(['iii : ' + routine + ' : rayon temporel en heures' + string(rtemphor)]) 296 294 297 298 295 ;____________________________________ 299 296 ; pour le fichier netCDF … … 349 346 correction_m=moytbliss-mean(moytbliss[22:24]) 350 347 endif 351 352 348 353 349 testfilename='' … … 510 506 ; fin boucle sur le temps 511 507 512 513 508 if nbfile GT 0 THEN BEGIN 514 509 good=where (kont1 ge 3,nb) -
trunk/src/cresamsu_cdf.pro
r636 r640 496 496 endfor 497 497 498 499 498 end -
trunk/src/define_area.m
r634 r640 146 146 error(usage); 147 147 end 148 %149 148 % 150 149 arg_info=whos('lat_value'); -
trunk/src/essailai.m
r635 r640 76 76 % 77 77 % $URL$ 78 %79 78 % 80 79 % - fplod 20110920T091450Z aedon.locean-ipsl.upmc.fr (Darwin) … … 254 253 255 254 % initialise les jours (les valeurs de x, les valeurs de y étant les sommes, les moyennes, 256 % les ecart typecalculées255 % les écarts-types calculées 257 256 startdate=datenum(yyyy_d,1,1); 258 257 enddate=datenum(yyyy_d,12,eomday(yyyy_d,12)); … … 301 300 clear fullfilename; 302 301 % 303 % dessin de l'évolution des écart -types des LAI302 % dessin de l'évolution des écarts-types des LAI 304 303 ifigure=ifigure+1; 305 304 figure(ifigure); -
trunk/src/extract_amsua.pro
r636 r640 519 519 + string(SYSTIME(1)-time1,format='(F12.6)')]) 520 520 ENDIF 521 522 521 523 522 ; appel à interpolswath pour ajuster les pixels amsua sur une grille -
trunk/src/extract_amsuab.pro
r636 r640 196 196 ; 197 197 ; get rid of uppercase 198 ;199 198 ; 200 199 ;- … … 249 248 STOP 250 249 ENDIF 251 252 250 253 251 ; ouverture des fichiers liste (année, mois, jour, tous satellites) pour … … 268 266 readf, lun_a, onefile 269 267 filea[index_filea] = onefile 270 ; isolate string independ ant from satellite268 ; isolate string independent from satellite 271 269 a[index_filea]=strmid(file_basename(filea[index_filea]),7) 272 270 index_filea = index_filea + 1 … … 282 280 while (not eof(lun_b)) do begin 283 281 readf,lun_b,fileb 284 ; isolate string independ ant from satellite282 ; isolate string independent from satellite 285 283 b=strmid(file_basename(fileb),7) 286 284 ; check for match between AMSU-A and AMSU-B … … 540 538 endelse 541 539 542 543 540 ; fin boucle sur les fichiers lus 544 541 labfile: -
trunk/src/extract_nc.m
r635 r640 112 112 ext_data=getnc(file,cell2str(var(k)),ind_min,ind_max,-1,-1,2,1,0); 113 113 % le dernier 0 du getnc permet de ne pas utiliser 114 % squeeze par d efaut sur la matrice de sortie114 % squeeze par défaut sur la matrice de sortie 115 115 % mais besoin d'ajustement qd le singleton est la premiÚre ou derniÚre dim 116 116 if ind_min(1)-ind_max(1)==0 | ind_min(end)-ind_max(end)==0 -
trunk/src/figure01.pro
r633 r640 27 27 ; 28 28 ; :Todo: 29 ;30 29 ; 31 30 ; :Examples: -
trunk/src/figure02.pro
r633 r640 27 27 ; 28 28 ; :Todo: 29 ;30 29 ; 31 30 ; :Examples: -
trunk/src/figure03.pro
r633 r640 27 27 ; 28 28 ; :Todo: 29 ;30 29 ; 31 30 ; :Examples: -
trunk/src/figure04.pro
r633 r640 27 27 ; 28 28 ; :Todo: 29 ;30 29 ; 31 30 ; :Examples: -
trunk/src/figure05.pro
r633 r640 28 28 ; :Todo: 29 29 ; 30 ;31 30 ; :Examples: 32 31 ; -
trunk/src/figure06.pro
r633 r640 27 27 ; 28 28 ; :Todo: 29 ;30 29 ; 31 30 ; :Examples: -
trunk/src/figure07.pro
r633 r640 25 25 ; :ref:`varamma_profile.sh` 26 26 ; 27 ;28 27 ; :Examples: 29 28 ; -
trunk/src/figure09.pro
r637 r640 31 31 ; 32 32 ; :Todo: 33 ;34 33 ; 35 34 ; :Examples: -
trunk/src/figure11.pro
r633 r640 30 30 ; 31 31 ; :Todo: 32 ;33 32 ; 34 33 ; :Examples: -
trunk/src/figure12.pro
r634 r640 28 28 ; :Todo: 29 29 ; 30 ;31 30 ; :Examples: 32 31 ; -
trunk/src/figure13.pro
r633 r640 30 30 ; :Todo: 31 31 ; 32 ;33 32 ; :Examples: 34 33 ; -
trunk/src/figure14.pro
r634 r640 24 24 ; for AMSU dataset be sure to have :file:`{c}{x}.an{yyyy}_climato.nc` 25 25 ; in the directory defined in :file:`${PROJECT_ID}/` 26 ;27 26 ; 28 27 ; :Post: -
trunk/src/forout.pro
r634 r640 36 36 ; 37 37 ; :Todo: 38 ;39 38 ; 40 39 ; :Examples: -
trunk/src/geolocation_to_string.m
r635 r640 10 10 % geolocation_to_string.m 11 11 % ======================= 12 %13 12 % 14 13 % .. function:: geolocation_to_string(lon_value, lat_value) -
trunk/src/get_msg_xrit.sh
r638 r640 223 223 # 224 224 # Autre piste de conversion /bdd/msg/software/read_xrit/readXRIT: 225 #226 225 # 227 226 # Compilation + édition de lien … … 526 525 echo "yyyymmdde=${yyyymmdde}" 1>> ${log} 527 526 echo "" 1>> ${log} 528 #529 527 # 530 528 dirref_base=/bdd/msg/images/xrit/ -
trunk/src/get_station_geolocation.m
r636 r640 139 139 end 140 140 % 141 %142 141 arg_info=whos('station_value'); 143 142 if ~strcmp(arg_info.class,'cell') … … 149 148 % 150 149 % +todo+ disp('arg station_value must be a cell array of strings'); 151 %152 150 % 153 151 nb_station=size(station_value,1); -
trunk/src/graphe_zicode.m
r639 r640 1 1 function [result]=graphe_zicode(yyyy, nbligne, zi_code, seuil, classes_select, msdnom_select) 2 3 2 4 3 %GRAPHE_ZICODE appelle les programmes read_lai_gis.m et … … 56 55 % applied to select data. 57 56 % 58 % ++ tracer de la moyenne et de l ecarttype du lai des classes de végétation59 % correspondant a chaque type de sol autour de chaque pluviomÚtrespour57 % ++ tracer de la moyenne et de l'écart-type du lai des classes de végétation 58 % correspondant à chaque type de sol autour de chaque pluviomÚtre pour 60 59 % 61 60 % ++for soil types and vegetation plot LAI data … … 236 235 for ind_classes=1:numel(list_classes_select) 237 236 disp([' ind_classes : ' num2str(ind_classes) ' ,classes_select : ' char(list_classes_select(ind_classes))]); 238 % calcule les moyennes et écart 237 % calcule les moyennes et écart-type non filtrés et filtrés 239 238 % et produit une figure 240 239 [lai_mean, lai_std]=plot_lai_gis(datejul_value, datestr_value, lai_value, classes_value, msdnom_value, char(list_classes_select(ind_classes)), char(list_msdnom_select(ind_msd)), zi_code, imageprefix, seuil); -
trunk/src/hovmuller.pro
r637 r640 91 91 ; % Type conversion error: Unable to convert given STRING to Long. 92 92 ; 93 ;94 93 ; ++ check timebox (datemin, datemax) validity 95 94 ; … … 107 106 ; 108 107 ; :Examples: 109 ;110 108 ; 111 109 ; To plot AMSU a4 hovmuller latitude plot according to :file:`${PROJECT_ID}/a4.an2006.nc` between 2006040101 and 20061031 and produce a PNG file: -
trunk/src/hovmullerlat.m
r637 r640 2 2 3 3 %+ 4 %5 4 % 6 5 % ============== -
trunk/src/hovmullerlonorlat.m
r637 r640 109 109 end 110 110 111 112 111 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 113 112 % Création des bornes sur la dimension moyenne … … 119 118 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 120 119 121 122 120 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 123 121 % Extraction des données concernées … … 147 145 148 146 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 149 150 147 151 148 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 161 158 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 162 159 163 164 160 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 165 161 % Dessin du diagramme de Hovmöller … … 176 172 nbtimestep = size(data,1); 177 173 178 % On cherche le pas d'affichage des pas de temps pour ne pas d epasser la174 % On cherche le pas d'affichage des pas de temps pour ne pas dépasser la 179 175 % limite autorisée (on prend l'arrondi supérieur avec la fonction 'ceil'). 180 176 tickstep = ceil(nbtimestep / nbtickmax); -
trunk/src/idl_amsu_netcdf.pro
r638 r640 245 245 ; si ou ne rien écrire 246 246 ; 247 ;248 247 ;--------------------------------------- 249 248 ; création du fichier netCDF de sortie -
trunk/src/int2precip.m
r636 r640 94 94 % par l'analyse des données MSG et des données EPSAT sur une zone de -20 à 95 95 % 0 en longitudes et de 0 à 20 en latitudes. Dans chaque cas, la formule 96 % moyenne est obtenue apartir des proches résultats pour un domaine donnée96 % moyenne est obtenue à partir des proches résultats pour un domaine donnée 97 97 % des régressions linéaires, cf. MSGanalysiscrossvalidation). 98 98 switch type -
trunk/src/intensityMSG_concat.m
r636 r640 20 20 % - Données d'entrées : 21 21 % 22 % * ncfilesset (en caractÚres) : Type de fichier NetCDF que l'on cherche,22 % * ncfilesset (en caractÚres) : Type de fichier netCDF que l'on cherche, 23 23 % 'normal' ou 'extracted'. 24 24 % * lonmin : Borne inférieure des longitudes du domaine étudié. … … 119 119 % - jaclod 2011-07-20 120 120 % 121 % * Passage sur Cratos, fonctionne!121 % * Passage sur cratos, fonctionne! 122 122 % * Modifications de la documentation et des commentaires. 123 123 % … … 155 155 156 156 %%%%%%%%%%%%%%%%% 157 158 157 159 158 %%%%%%%%%%%%%%%%%%%%%%%%%%% … … 243 242 disp(['iii : Analyse du fichier msg-tb108_' annee '-' mois '-' jour '_15min.nc.']); 244 243 245 % Lecture du fichier NetCDF, récupération des244 % Lecture du fichier netCDF, récupération des 246 245 % données et des variables. 247 246 [TbMSG, Temps, Longitudes, Latitudes] = MSGread(ncfilesset, fullfilename); … … 264 263 % Si la résolution demandée est différente de celle 265 264 % initiale, alors on interpole les données avec la 266 % fonction 'spatial resolution'.265 % fonction spatialresolution. 267 266 if longresol ~= size(TbMSG,2) || latresol ~= size(TbMSG,3); 268 267 [Longitudes, Latitudes, TbMSG] = spatialresolution(Longitudes, Latitudes, TbMSG, longresol, latresol); … … 307 306 disp(['iii : Analyse du fichier extracted-msg-tb108_' annee '-' mois '-' jour '_15min.nc.']); 308 307 309 % Lecture du fichier NetCDF, récupération des308 % Lecture du fichier netCDF, récupération des 310 309 % données et des variables. 311 310 [TbMSG, Temps, Longitudes, Latitudes] = MSGread(ncfilesset, fullfilename); … … 327 326 % Si la résolution demandée est différente de celle 328 327 % initiale, alors on interpole les données avec la 329 % fonction 'spatial resolution'.328 % fonction spatialresolution. 330 329 if longresol ~= size(TbMSG,2) || latresol ~= size(TbMSG,3); 331 330 [Longitudes, Latitudes, TbMSG] = spatialresolution(Longitudes, Latitudes, TbMSG, longresol, latresol); … … 404 403 mois = listdir(numdir+2).name; 405 404 406 % R epertoire du mois correspondant.405 % Répertoire du mois correspondant. 407 406 repertoireMSGmois = [repertoireMSGsaison '/' mois]; 408 407 … … 426 425 disp(['iii : Analyse du fichier ' listncfiles(numfiles).name '.']); 427 426 428 % Lecture du fichier NetCDF, récupération des données et des427 % Lecture du fichier netCDF, récupération des données et des 429 428 % variables. 430 429 [TbMSG, Temps, Longitudes, Latitudes] = MSGread(ncfilesset, fullfilename); … … 492 491 end 493 492 494 % On passe le temps en premi ere dimension pour plus de concordance avec les493 % On passe le temps en premiÚre dimension pour plus de concordance avec les 495 494 % autres variables. 496 495 ij_intensity_concat = permute(ij_intensity_concat,[3 1 2]); -
trunk/src/intensityMSG_day.m
r636 r640 27 27 % par défaut). 28 28 % * Tb_SeuilSup (Facultatif) : Seuil de température de brillance 29 % maximale délimitant un nuage convectif (en K, 235 par d efaut).29 % maximale délimitant un nuage convectif (en K, 235 par défaut). 30 30 % * Tb_LimNoyau (Facultatif) : Seuil de température de brillance 31 31 % délimitant le noyau dans notre modÚle à deux couches (en K, 190 par … … 211 211 %- 212 212 213 214 213 % Appel de la variable global 'application' (Il faut s'assurer d'avoir 215 214 % exécuté le programme 'varamma_startup' dans un premier temps). … … 252 251 idim_lon=2; 253 252 idim_lat=3; 254 255 253 256 254 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 297 295 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 298 296 299 300 297 %%%%%%%%%%%%%%%%%%%%%%%%%%% 301 298 % Lissage et interpolation … … 333 330 334 331 %%%%%%%%%%%%%%%%%%%%%%%%%%% 335 336 332 337 333 %%%%%%%%%%%%%%%%%%%%%%%%%%% … … 414 410 %%%%%%%%%%%%%%%%%%%%%%%%%%% 415 411 416 417 412 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 418 413 % ModÚle de nuages à deux couches … … 436 431 437 432 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 438 439 433 440 434 %%%%%%%%%%%%% -
trunk/src/interpol_correc.pro
r635 r640 39 39 ; } 40 40 ; 41 ;42 41 ; SEE ALSO 43 42 ; ======== … … 56 55 ; print, cor_sb 57 56 ; 58 ;59 57 ; TODO 60 58 ; ==== 61 ;62 59 ; 63 60 ; on peut changer le codage en dur des fauchées ou demi-fauchées, mais -
trunk/src/interpolswath.pro
r634 r640 37 37 ; 38 38 ; régler le pb des longitudes autour de 180 degré 39 ; l'interpolation ne les g ere pas et cree des points dans la zone AMMA...39 ; l'interpolation ne les gÚre pas et crée des points dans la zone AMMA... 40 40 ; 41 41 ; EVOLUTIONS … … 52 52 ; - fplod 20111130T120358Z cratos (Linux) 53 53 ; 54 ; * minimal re st header54 ; * minimal reStructuredText header 55 55 ; 56 56 ;- … … 115 115 grid1[i]=-grid2[ndemigrid-1-i] 116 116 endfor 117 ;alternative equivalente)117 ;alternative équivalente 118 118 ;for i=ndemigrid-1,0,-1 do grid1[i]=-(ndemigrid-i-1)*delta-fova[nnadir] 119 119 grid=[grid1,grid2] … … 145 145 ; interpolation des positions 146 146 latgrid=interpol(latlu,fova,grid) 147 ; pour la longitude, traitement sp ecifique des longitudes proches de -180147 ; pour la longitude, traitement spécifique des longitudes proches de -180 148 148 if (mean(lonlu) ge -120.) then begin 149 149 longrid=interpol(lonlu,fova,grid) -
trunk/src/map_amab.pro
r635 r640 96 96 nbpix_b = 90 ; nb pixels dans la fauchée AMSUB = nb_fov 97 97 98 99 98 ; sélection des lignes amsub qui minimisent la distance avec le centre 100 99 ; des lignes amsua … … 102 101 bsync=where(pulse_line GT 0) ; si use_amsua=0 alors pulse_line contient toutes les lignes amsub, entrees dans async 103 102 async=pulse_line[bsync] 104 105 103 106 104 fovb=indgen(nbpix_b - 1 ) + 1 … … 147 145 bmch4_1=REFORM(amb_scan.btemps[3,*,*]/100.) 148 146 bmch5_1=REFORM(amb_scan.btemps[4,*,*]/100.) 149 150 147 151 148 bmzen_1=REFORM(amb_scan.angles[0,*,*]/100.) -
trunk/src/meanandstdplot.m
r636 r640 35 35 % moyennes sur le temps de la variable de la matrice 'data'. 36 36 % * stddata : Matrice 2D (longitudes, latitudes) comportant les données 37 % d' ecart-types sur le temps de la variable de la matrice 'data'.37 % d'écarts-types sur le temps de la variable de la matrice 'data'. 38 38 % 39 39 % EXAMPLES … … 87 87 % Création de variables comportant les titres et noms des axes. 88 88 Title1 = ['Moyennes des ' nomvar ' du ' datestr(dim1(1),'dd mmm yyyy') ' au ' datestr(dim1(size(dim1,2)),'dd mmm yyyy')]; 89 Title2 = [' Ecart-types des ' nomvar ' du ' datestr(dim1(1),'dd mmm yyyy') ' au ' datestr(dim1(size(dim1,2)),'dd mmm yyyy')];89 Title2 = ['Ãcarts-types des ' nomvar ' du ' datestr(dim1(1),'dd mmm yyyy') ' au ' datestr(dim1(size(dim1,2)),'dd mmm yyyy')]; 90 90 Abs = 'Longitudes'; 91 91 Ord = 'Latitudes'; -
trunk/src/mem_to_file_amsu_t2.pro
r636 r640 166 166 ENDIF 167 167 168 print, 'ouverture pour ecriture de ', fullfilename168 print, 'ouverture pour écriture de ', fullfilename 169 169 ; 170 170 ; build header format string -
trunk/src/missingdaytonan.m
r636 r640 106 106 107 107 %%%%%%%%%%%%%%%%% 108 109 108 110 109 %%%%%%%%%%%%%%%%%%%%%% -
trunk/src/msgnuage.pro
r636 r640 16 16 ; ========== 17 17 ; 18 ; base sur :ref:`ssminuage.pro`, ecrit par L. Eymard pour l'etude des18 ; base sur :ref:`ssminuage.pro`, écrit par L. Eymard pour l'étude des 19 19 ; nuages durant FASTEX 20 20 ; … … 54 54 ; 55 55 ; * si un jour on retravaille avec SSMI 56 ; Premier nettoyage: si moins de 20 lignes dans l'image, eliminer l'image,56 ; Premier nettoyage: si moins de 20 lignes dans l'image, éliminer l'image, 57 57 ; sinon traitement 58 ;59 58 ; 60 59 ; EVOLUTIONS … … 77 76 ; detection part 78 77 ; * move time information decoding after reading or initializing fake data 79 ; because independ ant of cloud detection78 ; because independent of cloud detection 80 79 ; * 2 clouds in fake dataset 81 80 ; * loop to write each cloud informations … … 256 255 ;loadct,18 257 256 plt,nuag,window=2,/realcont,/cont_nofill,title='nuages MSG',/rempli,/nocontour,min=min(nuag),max=max(nuag) 258 ;Moyennes et ecartstypes dans les nuages257 ;Moyennes et écarts-types dans les nuages 259 258 compt=intarr(numax) 260 259 lmoy=fltarr(numax) -
trunk/src/nonprecipperiods.m
r636 r640 4 4 %pauses dans les précipitations et de tracer la carte des moyennes des 5 5 %jours de non précipitation dans une pause ainsi que la carte des 6 %écart -types.6 %écarts-types. 7 7 8 8 % … … 37 37 % * minmean (Facultatif) : Borne inférieure de l'échelle des moyennes. 38 38 % * maxmean (Facultatif) : Borne supérieure de l'échelle des moyennes. 39 % * minstd (Facultatif) : Borne inférieure de l'échelle des écart -types.40 % * maxstd (Facultatif) : Borne inférieure de l'échelle des écart -types.39 % * minstd (Facultatif) : Borne inférieure de l'échelle des écarts-types. 40 % * maxstd (Facultatif) : Borne inférieure de l'échelle des écarts-types. 41 41 % 42 42 % - Données de sorties : … … 45 45 % non précipitations et dont les dimensions sont longitudes et 46 46 % latitudes. 47 % * stdtimestepof0 : Matrice comportant les données des écart -types des47 % * stdtimestepof0 : Matrice comportant les données des écarts-types des 48 48 % jours de non précipitations et dont les dimensions sont longitudes et 49 49 % latitudes. … … 51 51 % Cette fonction permet de suivre les périodes de pauses dans les 52 52 % et de tracer la carte des moyennes des jours de non précipitation dans 53 % une pause ainsi que la carte des écart -types.53 % une pause ainsi que la carte des écarts-types. 54 54 % 55 55 % EXAMPLES … … 98 98 % 99 99 % * Ajout des données de sorties. 100 % * Ajout d'une restriction sur les échelles de moyennes et écart -types100 % * Ajout d'une restriction sur les échelles de moyennes et écarts-types 101 101 % pour une meilleure observation. 102 102 % * Ajout de la démonstration. … … 136 136 meantimestepof0 = zeros(size(ij_concat,2),size(ij_concat,3)); 137 137 138 % Initialisation de la matrice qui contiendra les ecart-types du temps en138 % Initialisation de la matrice qui contiendra les écarts-types du temps en 139 139 % jour des périodes de non précipitations. 140 140 stdtimestepof0 = zeros(size(ij_concat,2),size(ij_concat,3)); … … 214 214 215 215 %%%%%%%%%%%%%%%%%%%%%%%% 216 % Moyenne et ecart-type216 % Moyenne et écart-type 217 217 %%%%%%%%%%%%%%%%%%%%%%%% 218 218 … … 226 226 meantimestepof0(ilon,ilat) = ij_meantimestepof0; 227 227 228 % Calcul de l' ecart-type.228 % Calcul de l'écart-type. 229 229 if ~isempty(indzero) 230 230 ij_stdtimestepof0 = sqrt(sum(power(timestepof0,2))/nbperiods - power(ij_meantimestepof0,2)); … … 246 246 247 247 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 248 % Images des moyennes et ecart-types248 % Images des moyennes et écarts-types 249 249 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 250 250 … … 260 260 261 261 Title1 = ('Moyennes des jours de non précipitation'); 262 Title2 = (' Ecart-types des jours de non précipitation');262 Title2 = ('Ãcarts-types des jours de non précipitation'); 263 263 Abs = 'Longitudes'; 264 264 Ord = 'Latitudes'; … … 286 286 plot(Xgco,Ygco,'k.'); 287 287 288 % Ecart-types (la transposée est également utilise ici).288 % Ãcarts-types (la transposée est également utilise ici). 289 289 subplot(1,2,2); 290 290 imagesc(Longitudes,Latitudes,stdtimestepof0'); -
trunk/src/nuage.pro
r637 r640 115 115 ; END 116 116 ;ENDCASE 117 118 117 119 118 ;ENDIF … … 177 176 ENDFOR 178 177 179 180 178 labfin: 181 179 numax=max(nuag) -
trunk/src/obsolete/cresamsu_cdf.pro
r571 r640 101 101 ; a partir de fov=+/-11, dedoublement du pixel pixel +/- diametre tache au sol/6 102 102 ; a partir de fov=+/-14, pixel triple centre et +/- diametre tache au sol/4 103 ;104 103 ; 105 104 ; 2011/10/24 lelod … … 495 494 endfor 496 495 497 498 496 end -
trunk/src/obsolete/extract_amsuab.pro
r571 r640 187 187 ; 188 188 ; get rid of uppercase 189 ;190 189 ; 191 190 ;- … … 240 239 STOP 241 240 ENDIF 242 243 241 244 242 ; ouverture des fichiers liste (annee, mois, jour, tous satellites) pour … … 531 529 endelse 532 530 533 534 531 ; fin boucle sur les fichiers lus 535 532 labfile: -
trunk/src/obsolete/map_amab.pro
r571 r640 96 96 nbpix_b = 90 ; nb pixels dans la fauchee AMSUB = nb_fov 97 97 98 99 98 ; selection des lignes amsub qui minimisent la distance avec le centre 100 99 ; des lignes amsua … … 102 101 bsync=where(pulse_line GT 0) ; si use_amsua=0 alors pulse_line contient toutes les lignes amsub, entrees dans async 103 102 async=pulse_line[bsync] 104 105 103 106 104 fovb=indgen(nbpix_b - 1 ) + 1 … … 147 145 bmch4_1=REFORM(amb_scan.btemps[3,*,*]/100.) 148 146 bmch5_1=REFORM(amb_scan.btemps[4,*,*]/100.) 149 150 147 151 148 bmzen_1=REFORM(amb_scan.angles[0,*,*]/100.) -
trunk/src/obsolete/paper01.pro
r571 r640 9 9 ; 10 10 ; plot all figures for paper #1 11 ;12 11 ; 13 12 ; .. only:: man -
trunk/src/obsolete/ssmi_fastex.pro
r571 r640 75 75 ;imcol=[1,8,16,30,50,63,80,110,150,180,230] 76 76 loadct,couleur 77 78 77 79 78 ;**************************************************** … … 162 161 print,'3 ',lat(63,l-1),lon(63,l-1) 163 162 print,'4 ',lat(0,l-1),lon(0,l-1) 164 165 163 166 164 ;***************************************************** … … 261 259 endfor 262 260 263 264 261 ;*************************************************** 265 262 ; Visualisation des images … … 402 399 0,sx,sy), px(0),py(0) 403 400 404 405 401 if rep eq 'ps' then $ 406 402 tv,image,px(0),py(0),xsize=sx,ysize=sy … … 485 481 Endfor 486 482 487 488 483 plots,[x0,x1],[yy,yy],color=color_trait,/normal 489 484 plots,[x1,x1],[y0,y1],thick=0.5,color=color_trait,/normal -
trunk/src/paper01.pro
r633 r640 9 9 ; 10 10 ; plot all figures for paper #1 11 ;12 11 ; 13 12 ; .. only:: man … … 46 45 ; 47 46 ; } 48 ;49 47 ; 50 48 ; :Categories: -
trunk/src/peakdet.m
r422 r640 157 157 %! plot(mintab(:,1), mintab(:,2), 'g+'); 158 158 %! plot(maxtab(:,1), maxtab(:,2), 'r+'); 159 160 159 161 160 %!demo -
trunk/src/plot_amsu.pro
r636 r640 111 111 if rep eq 'oui' then goto,lab_orbite 112 112 113 114 113 ;mer=0 115 114 ;if mer ne 0 then begin -
trunk/src/plot_cartes3B42.m
r636 r640 156 156 %print -djpeg -r72 'carte3B42_2001.jpg'; 157 157 subplot1(2) 158 %title('PPT Ãcart 158 %title('PPT Ãcart-type 2000-2009','fontsize',14) 159 159 montitre=sprintf('ppt_ecartyp : %s\n',ppt_ecartyp); 160 160 montitre=sprintf('Année : %d\n',annee); … … 166 166 colorbar("SouthOutside"); 167 167 168 %title('Ãcart Type PPT en 2008','fontsize',14)168 %title('Ãcart-type PPT en 2008','fontsize',14) 169 169 % pour l'inversion de l'axe des ordonnées y (origine en bas à 170 170 % gauche, valeurs Y croissantes de bas en haut (par défaut) -
trunk/src/plot_clim3B42.m
r636 r640 75 75 axis('xy'); 76 76 axis([lonmin lonmax latmin latmax]); 77 title('PPT Ecart Type 2000-2009','fontsize',14);77 title('PPT Ãcart-type 2000-2009','fontsize',14); 78 78 xlabel('Longitude','fontsize',14); 79 79 %ylabel('Latitude','fontsize',14); -
trunk/src/plot_cumulEPSAT.m
r635 r640 2 2 3 3 %PLOT_CUMULEPSAT dessine les précipitations cumulées sur une journée d'aprÚs les données EPSAT 4 5 4 6 5 %+ … … 15 14 % =========== 16 15 % 17 % - Donn ï¿œes d'entrï¿œes :16 % - Données d'entrées : 18 17 % 19 18 % * Annee (en caractÚres) : Année (AAAA) du jour … … 62 61 % ==== 63 62 % 64 % trouver l' ï¿œquivalent de 'alphadata'63 % trouver l'équivalent de 'alphadata' 65 64 % 66 65 % strange continents contour on the plot !! … … 70 69 % cumulEPSAT_day should output metadata (date and geolocation) 71 70 % 72 % !! cumul == 0 forc ï¿œ ï¿œNaN !!! mais pourquoi ?71 % !! cumul == 0 forcé à NaN !!! mais pourquoi ? 73 72 % 74 73 % Annee, Mois Jour should be arrays or better only one datestr_value array or usage of Temps … … 93 92 % * Pour le imagesc, repassage en transposée sinon les longitudes 94 93 % (dimension 1) sont en y sur le dessin. 95 % * Fonctionne sur le Matlab de Cratos.94 % * Fonctionne sur le Matlab de cratos. 96 95 % 97 96 % - pinsard 2011-06-08T14:42:32Z loholt1.ipsl.polytechnique.fr (Linux) … … 121 120 % - jaclod 2011-03-02 122 121 % 123 % * Ajout de la doc et de la d ï¿œmonstration122 % * Ajout de la doc et de la démonstration 124 123 % 125 124 % - jaclod 2011-02-25 126 125 % 127 % * Cr ï¿œation126 % * Creation 128 127 % 129 128 %- … … 143 142 end 144 143 % 145 % Chargement des fichiers pour tracer les c ï¿œtes.144 % Chargement des fichiers pour tracer les cÃŽtes. 146 145 load('cmgco_traitdecote.mat','Xgco','Ygco'); 147 146 … … 168 167 colorbar; 169 168 170 % Ajout des c ï¿œtes.169 % Ajout des cÃŽtes. 171 170 coast = plot(Xgco,Ygco,'k.'); 172 171 … … 176 175 title(['Precipitations cumulees en mm pour le ' ([Annee '-' Mois '-' Jour])]); 177 176 178 % Ex ecution de la fonction permettant d'afficher les coordonnées en format177 % Exécution de la fonction permettant d'afficher les coordonnées en format 179 178 % °N/S et °E/W. 180 179 %status_ticklatlon=ticklatlon(); … … 234 233 %! [data, Temps, Longitudes, Latitudes] = extractedmatrixdata(data, Temps, Longitudes, Latitudes, mindim1, maxdim1, lonmin, lonmax, latmin, latmax); 235 234 %! [cumuldata] = cumulEPSAT_day(data); 236 ï¿œ%! result = plot_cumulEPSAT(cumuldata, Annee, Mois, Jour, Longitudes, Latitudes);235 %! result = plot_cumulEPSAT(cumuldata, Annee, Mois, Jour, Longitudes, Latitudes); -
trunk/src/plot_lai_gis.m
r639 r640 149 149 % 150 150 % * fix legend problem 151 % * restore plot of ecart type151 % * restore plot of écart type 152 152 % 153 153 % - lelod 20111018 … … 238 238 % - soclod 20110812 239 239 % 240 % * ajout de sorties aa_moy et aa_ecartype les moyennes et ecart types filtres%240 % * ajout de sorties aa_moy et aa_ecartype les moyennes et écart types filtres% 241 241 % - soclod/fplod 20110812T145149Z aedon.locean-ipsl.upmc.fr (Darwin) 242 242 % … … 300 300 clear arg_info 301 301 % 302 % 303 % trace des lai moyens et ecart type pour une classe de vegetation 302 % trace des lai moyens et écart-type pour une classe de végétation 304 303 % choisie avec classes_select (tri dans le tableau classes_value) et avec 305 304 % msdnom_select (tri dans le tableau msdnom_value) … … 335 334 disp(['iii : pas assez de points (' num2str(aa) ' < ' num2str(seuil) ' )']) 336 335 else 337 % calcul moyenne et ecarttype338 % filtre des valeurs moyennes et des ecarttypes336 % calcul moyenne et écart-type 337 % filtre des valeurs moyennes et des écarts-types 339 338 % index_juin=20; 340 339 for index_doy=1:nb_doy … … 367 366 malegende_p1=[malegende_p1 sprintf(' (%d points)', aa)]; 368 367 malegende_p2='spline'; 369 malegende_p3='d erivee 2de de spline';370 malegende_p4='+/- ecarttype';368 malegende_p3='dérivée 2de de spline'; 369 malegende_p4='+/- écart-type'; 371 370 % 372 371 % lissage de la moyenne par spline … … 456 455 s={s{:},get(h,'userdata')}; 457 456 end 458 h_legend=legend('{\fontsize{11}+/- ecarttype}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}')457 h_legend=legend('{\fontsize{11}+/-écart-type}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}') 459 458 % legend(hc,s) 460 % h_legend=legend('{\fontsize{11}+/- ecarttype}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}', ...459 % h_legend=legend('{\fontsize{11}+/-écart-type}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}', ... 461 460 % '{\fontsize{11}spline}','{\fontsize{11}derivee 2nde de spline}','{\fontsize{11}min}','{\fontsize{11}max%}') 462 461 %set(h_legend, 'fontsize',14) … … 655 654 %! [lai_mean, lai_std] = plot_lai_gis(datejul_value, datestr_value, lai_value, classes_value, msdnom_value, classes_select, msdnom_select, zi_code, imageprefix, seuil); 656 655 657 658 656 % x=transpose(datejul_value); 659 657 % y=transpose(lai_mean); -
trunk/src/plot_lai_parameterscycle.m
r638 r640 59 59 % automatiser sur toutes les années 60 60 % 61 %62 61 % EVOLUTIONS 63 62 % ========== … … 102 101 ind_maxtab=find(maxtab_ext(:,1)>=ddeb & maxtab_ext(:,1)<=dfin); 103 102 maxtab_fin=maxtab_ext(ind_maxtab) 104 105 103 106 104 %recherche date debut vegetation -
trunk/src/plot_lai_spline_detect_anomalies.m
r639 r640 153 153 % 154 154 % * fix legend problem 155 % * restore plot of ecarttype155 % * restore plot of écart-type 156 156 % 157 157 % - lelod 20111018 … … 242 242 % - soclod 20110812 243 243 % 244 % * ajout de sorties aa_moy et aa_ecartype les moyennes et ecart types filtres% 244 % * ajout de sorties aa_moy et aa_ecartype les moyennes et écarts-types filtrés 245 % 245 246 % - soclod/fplod 20110812T145149Z aedon.locean-ipsl.upmc.fr (Darwin) 246 247 % … … 304 305 clear arg_info 305 306 % 306 % 307 % trace des lai moyens et ecart type pour une classe de vegetation 307 % trace des lai moyens et écart-type pour une classe de végétation 308 308 % choisie avec classes_select (tri dans le tableau classes_value) et avec 309 309 % msdnom_select (tri dans le tableau msdnom_value) … … 339 339 % disp(['iii : pas assez de points (' num2str(aa) ' < ' num2str(seuil) ' )']) 340 340 % else 341 % calcul moyenne et ecarttype342 % filtre des valeurs moyennes et des ecarttypes341 % calcul moyenne et écart-type 342 % filtre des valeurs moyennes et des écarts-types 343 343 % index_juin=20; 344 344 % for index_doy=1:nb_doy … … 371 371 %malegende_p1=[malegende_p1 sprintf(' (%d points)', aa)]; 372 372 %malegende_p2='spline'; 373 %malegende_p3='d erivee 2de de spline';374 %malegende_p4='+/- ecarttype';373 %malegende_p3='dérivée 2de de spline'; 374 %malegende_p4='+/- écart-type'; 375 375 % 376 376 % lissage de la moyenne par spline … … 438 438 mintab_ext=datevec2doy(datevec(mintab)) 439 439 440 441 440 % mintab_ext1=datevec2doy(datevec(mintab)) 442 441 disp([ '1er max a : ' datestr(maxtab(6,1)) '-' num2str(maxtab(6 … … 477 476 s={s{:},get(h,'userdata')}; 478 477 end 479 h_legend=legend('{\fontsize{11}+/- ecarttype}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}')478 h_legend=legend('{\fontsize{11}+/-écart-type}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}') 480 479 % legend(hc,s) 481 % h_legend=legend('{\fontsize{11}+/- ecarttype}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}', ...480 % h_legend=legend('{\fontsize{11}+/-écart-type}','{\fontsize{11}"Savane arboree arbustive"-"FERRUGINEUX TROPICAUX"}', ... 482 481 % '{\fontsize{11}spline}','{\fontsize{11}derivee 2nde de spline}','{\fontsize{11}min}','{\fontsize{11}max%}') 483 482 %set(h_legend, 'fontsize',14) -
trunk/src/plot_zi.m
r638 r640 166 166 end 167 167 % 168 %169 168 arg_info=whos('west_value'); 170 169 if ~strcmp(arg_info.class,'double') … … 220 219 % 221 220 clear arg_info 222 %223 221 % 224 222 arg_info=whos('zi_code'); … … 231 229 % 232 230 % +todo+ disp('arg zi_code must be a cell array of strings'); 233 %234 231 % 235 232 arg_info=whos('zi_description'); -
trunk/src/plotimageMSG.m
r634 r640 2 2 3 3 %+ 4 %5 4 % 6 5 % ============== -
trunk/src/precipPlot.m
r637 r640 328 328 ylabel('precipjournaliere (mm)'); 329 329 title('précipitation journaliÚre'); 330 331 330 332 331 ifigure=1; -
trunk/src/project_init.pro
r634 r640 31 31 ; 32 32 ; @project_init 33 ;34 33 ; 35 34 ; SEE ALSO -
trunk/src/rapport_jaclod.m
r374 r640 30 30 % and inherit of visumodele1 reshape pb 31 31 % 32 % 20110801 : crash on ma lab 7.4 zeus ::32 % 20110801 : crash on matlab 7.4 zeus :: 33 33 % 34 34 % ncinquire: ncid 0: NetCDF: Not a valid ID … … 66 66 % 67 67 %- 68 %69 68 % 70 69 Annee = '2006'; -
trunk/src/readAMSU.m
r635 r640 62 62 lon=longitude(lo); 63 63 64 65 66 67 64 temperature=squeeze(nc{'moyenne_tb'}(jo,la,lo)); 68 69 65 70 66 %% mat_3D est la matrice (lat,lon,jours) que l''on utilise par la suite -
trunk/src/readAMSU_clim.m
r635 r640 67 67 lon=longitude(lo); 68 68 69 70 71 69 temperature=squeeze(nc{'moyenne_tb'}(1:length(jours),la,lo)); 72 73 70 74 71 %% mat_3D est la matrice (lat,lon,jours) que l''on utilise par la suite -
trunk/src/readERA.m
r635 r640 67 67 jours=jours(jo); 68 68 69 70 69 lom=find(longitude-360>z_deb(2) & longitude-360<z_fin(2)); 71 70 lop=find(longitude>z_deb(2) & longitude<z_fin(2)); 72 71 lon=cat(1,longitude(lom)-360,longitude(lop)); 73 74 72 75 73 %[lat,la]=sort(latitude,'ascend'); -
trunk/src/readERA_LE.m
r635 r640 7 7 % ============ 8 8 % 9 %10 9 % .. function:: readERA(file, z_deb, z_fin,nom_var,chemin) 11 10 % … … 13 12 % =========== 14 13 % 15 % routine de lecture d un fichier de données ERA40 au format net cdfsur16 % la zone d efinie par les points z_deb zfin de forme [lat lon time level]14 % routine de lecture d un fichier de données ERA40 au format netCDF sur 15 % la zone définie par les points z_deb zfin de forme [lat lon time level] 17 16 % 18 17 % EVOLUTIONS … … 64 63 jours=jours(jo); 65 64 66 67 65 lom=find(longitude-360>z_deb(2) & longitude-360<z_fin(2)); 68 66 lop=find(longitude>z_deb(2) & longitude<z_fin(2)); 69 67 lon=cat(1,longitude(lom)-360,longitude(lop)); 70 71 68 72 69 %[lat,la]=sort(latitude,'ascend'); … … 104 101 mat_3D=permute(variable,[2 3 1]); 105 102 106 107 103 else %Variable 4D dont on moyenne la dimension "niveau" (nom_var=r) ou la difference (geopt) 108 104 %extraction de la variable pour les longitudes positives et negatives -
trunk/src/read_lai_gis.m
r636 r640 28 28 % :param zi_code: zi_code ex 'barkedji' 29 29 % :type zi_code: string 30 %31 30 % 32 31 % read file :file:`${PROJECT_ID}/LAI_CSV/LAI{yyyy}F_vegfao_sol.csv` … … 430 429 end 431 430 % 432 %433 431 clear nb_line_data_max 434 432 % -
trunk/src/read_precip.m
r636 r640 55 55 % $ diff ${PROJECT_ID}/PRECIP/barkedji_pluie_2006.csv ${PROJECT_ID}/PRECIP/barkedji_pluie_2006_write_precip.csv 56 56 % 57 %58 57 % To read :file:`${PROJECT_ID}/PRECIP/pluvio_linguere_2006.csv`: 59 58 % … … 178 177 end 179 178 clear arg_info 180 %181 179 % 182 180 % build precipitation filename to be read -
trunk/src/read_rt_file.pro
r636 r640 91 91 ; TMPA-RT release 92 92 ;--------------------------------------------------------------- 93 94 93 95 94 pro read_3B4XRT_alg_ID, file, alg_ID -
trunk/src/read_station.m
r636 r640 51 51 % 52 52 % :ref:`data_precipitation_insitu` 53 %54 53 % 55 54 % :func:`simul_station` -
trunk/src/recup_camp.m
r638 r640 208 208 % ind_lat=find(15.062<=lat&lat<=15.121); 209 209 % ind_pt_fixes=find(-15.344<=lon&lon<=-15.270&15.062<=lat&lat<=15.121); 210 211 210 212 211 ind_pt_fixes=1:size(lai_value,1); -
trunk/src/scatterplot.m
r635 r640 211 211 end 212 212 213 214 213 %trace histogrammes 215 214 figure(2) -
trunk/src/script_cat_climeraint3d.sh
r635 r640 65 65 #data=${datera40} 66 66 67 68 67 # ERA INTERIM 69 68 dataerai=/bdd/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL … … 129 128 ncrcat ${chemin}/M${groupe}.2004.*.nc ${chemin}/cat${groupe}.2004.6TU.nc 130 129 131 132 130 # Moyenne des fichiers annuels 133 131 ncea ${chemin}/cat${groupe}.*.6TU.nc ${chemin}/climato.${groupe}.${an_deb}.${an_fin}.6TU.nc -
trunk/src/search_amsufiles.pro
r634 r640 100 100 return, result 101 101 ENDIF 102 ;103 102 ; 104 103 ; check required and optional parameters -
trunk/src/selectionimages.m
r635 r640 146 146 [a,indicefin]=min(abs(vecteurdatesamsu-date2)); 147 147 148 149 150 148 %lecture et enregistrement des dates et images AMSU 151 152 149 153 150 %definition du vecteur vecteurdatesamsu et de la matrice imagesamsu … … 186 183 vecteurdatesamsu=s; 187 184 188 189 190 191 185 %recherche des images MSG les plus proches temporellement des images AMSU 192 186 cd(cheminmsg) … … 205 199 end 206 200 207 208 201 %lecture des données msg, remplissage des vecteurs dates msg et heures msg 209 202 for i=1:indicefin-indicedebut+1 … … 233 226 end 234 227 235 236 228 pos=find(vecteurdatesmsg == 0); 237 229 vecteurdatesmsg(pos)=[]; … … 249 241 ylabel('TbMSG'); 250 242 251 252 243 cd .. 253 244 save donnees-msg.oct vecteurdatesmsg imagesmsg -
trunk/src/serietemp.m
r635 r640 126 126 %! clear all 127 127 %! close all 128 %! result = serietemp(rand(25,1), datenum(2011,07,25):datenum(2011,07,25)+24, 'mmm dd',7, 'valeurs aléatoires quotidiennes', 'sur 25 jours apartir du 25 juillet 2011');128 %! result = serietemp(rand(25,1), datenum(2011,07,25):datenum(2011,07,25)+24, 'mmm dd',7, 'valeurs aléatoires quotidiennes', 'sur 25 jours à partir du 25 juillet 2011'); -
trunk/src/showgrid.m
r638 r640 161 161 % octave> lon_value=load('lon.txt'); 162 162 % octave> plot(lon_value,lat_value) 163 %164 163 % 165 164 %- -
trunk/src/sigmoide3.m
r635 r640 2 2 3 3 %+ 4 %5 4 % 6 5 % =========== -
trunk/src/simul_lai_2d.m
r635 r640 240 240 clear index_datestr 241 241 % 242 %243 242 result=0; 244 243 % -
trunk/src/simul_msg_xrit_raw.py
r634 r640 41 41 42 42 Missing parameters 43 44 43 45 44 EXAMPLES -
trunk/src/simul_precip.m
r634 r640 23 23 % :type station: string 24 24 % :raise station: required 25 %26 25 % 27 26 % Charge en mémoire un tableau 2D simulé **precip_value** de précipitation dans : … … 116 115 clear arg_info 117 116 % 118 %119 117 arg_info=whos('station'); 120 118 if ~strcmp(arg_info.class,'char') -
trunk/src/simul_zi.m
r636 r640 97 97 clear arg_info 98 98 % 99 %100 99 % définition de la zone géographique 101 100 if (index_simulation == 1) -
trunk/src/simul_zi.py
r634 r640 52 52 53 53 Missing parameters 54 55 54 56 55 EXAMPLES -
trunk/src/ssmi_fastex.pro
r635 r640 75 75 ;imcol=[1,8,16,30,50,63,80,110,150,180,230] 76 76 loadct,couleur 77 78 77 79 78 ;**************************************************** … … 162 161 print,'3 ',lat(63,l-1),lon(63,l-1) 163 162 print,'4 ',lat(0,l-1),lon(0,l-1) 164 165 163 166 164 ;***************************************************** … … 261 259 ENDFOR 262 260 263 264 261 ;*************************************************** 265 262 ; Visualisation des images … … 402 399 0,sx,sy), px(0),py(0) 403 400 404 405 401 if rep eq 'ps' then $ 406 402 tv,image,px(0),py(0),xsize=sx,ysize=sy … … 485 481 ENDFOR 486 482 487 488 483 plots,[x0,x1],[yy,yy],color=color_trait,/normal 489 484 plots,[x1,x1],[y0,y1],thick=0.5,color=color_trait,/normal -
trunk/src/ssminuage.pro
r636 r640 56 56 ; 57 57 ; * si un jour on retravaille avec SSMI 58 ; Premier nettoyage: si moins de 20 lignes dans l'image, eliminer l'image,58 ; Premier nettoyage: si moins de 20 lignes dans l'image, éliminer l'image, 59 59 ; sinon traitement 60 ;61 60 ; 62 61 ; EVOLUTIONS … … 84 83 ; detection part 85 84 ; * move time information decoding after reading or initializing fake data 86 ; because independ ant of cloud detection85 ; because independent of cloud detection 87 86 ; * 2 clouds in fake dataset 88 87 ; * loop to write each cloud informations … … 264 263 ;loadct,18 265 264 plt,nuag,window=2,/realcont,/cont_nofill,title='nuages MSG',/rempli,/nocontour,min=min(nuag),max=max(nuag) 266 ;Moyennes et ecartstypes dans les nuages265 ;Moyennes et écarts-types dans les nuages 267 266 compt=intarr(numax) 268 267 lmoy=fltarr(numax) -
trunk/src/stagejasonjuillet.m
r637 r640 116 116 %%%%%%%%%% 117 117 118 119 118 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 120 119 % Moyennes et écart types des intensités sur la saison 2006 … … 132 131 133 132 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 134 135 133 136 134 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 180 178 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 181 179 182 183 180 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 184 181 % Diagramme de hovmöller des intensités et cumuls EPSAT sur 15 jours … … 199 196 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 200 197 201 202 198 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 203 199 % Distributions des intensités et cumuls EPSAT sur la saison 2006 … … 211 207 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 212 208 213 214 209 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 215 210 % Comparaison des séries temporelles sur la saison dans 5x5 degrés … … 227 222 228 223 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 229 230 224 231 225 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% … … 246 240 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 247 241 248 249 242 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 250 243 % Passage des intensités en précipitations … … 255 248 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 256 249 257 258 250 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 259 251 % Comparaison des nouvelles cartes sur la saison … … 279 271 280 272 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 281 282 273 283 274 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -
trunk/src/test_correc.pro
r633 r640 14 14 ; EXAMPLES 15 15 ; ======== 16 ;17 16 ; 18 17 ; .. code-block:: idl -
trunk/src/testinterpolswath.pro
r635 r640 53 53 for i=0,na-1 do lon(i)=0.2*(i-nnadir-1)*(i-nnadir-1)+10 54 54 55 56 55 lat=lat*deg & lon=lon*deg 57 56 … … 69 68 70 69 ;print,fova 71 72 73 74 70 75 71 dmax=max(fovtest)+max(pixa) -
trunk/src/testts.pro
r634 r640 24 24 ; :Todo: 25 25 ; 26 ;27 26 ; ++ comment ne pas faire figurer dans les labels de l'axe x l'année bidon ? 28 27 ; :Examples: 29 ;30 28 ; 31 29 ; .. code-block:: idl -
trunk/src/timeserie.pro
r634 r640 140 140 ; pltt,window=1,seriea5, 't',title='climatology time series',min=0.1,max=0.9 141 141 ; xyouts,max(time)+20,0.1,'a5',color=0,charsize=2 142 ;143 142 ; 144 143 ; :Version: -
trunk/src/utilitaires/ajoute_date_fig.m
r635 r640 48 48 end 49 49 50 51 50 text(pos_x,pos_y,date) 52 51 -
trunk/src/utilitaires/boucheNaN_3d.m
r308 r640 36 36 end 37 37 38 39 38 mat_bouche=mat; 40 39 for(i=1:ti) -
trunk/src/utilitaires/supprime_point_trous.m
r635 r640 7 7 % 8 8 % .. function:: supprime_point_trous(mat_3D) 9 %10 9 % 11 10 % DESCRIPTION -
trunk/src/varamma.kml
r308 r640 219 219 % ind_pt_fixes=find(-16.091<=lon&lon<=-16.064&15.613<=lat&lat<=15.637); 220 220 221 222 221 --> 223 222 <GroundOverlay> -
trunk/src/varamma_startup.m
r635 r640 25 25 % 26 26 % >> varamma_startup 27 %28 27 % 29 28 % If linked or copy to ./startup.m, it will be automatically -
trunk/src/visuEPSATprecipcumul.m
r636 r640 2 2 3 3 %VISUEPSATPRECIPCUMUL permet une visualisation des précipitations estimées 4 %cumulées sur une journée d'apr es EPSAT.4 %cumulées sur une journée d'aprÚs EPSAT. 5 5 6 6 % … … 81 81 % - jaclod 2011-07-28 82 82 % 83 % * Fonctionne sur le Matlab de Cratos.83 % * Fonctionne sur le Matlab de cratos. 84 84 % 85 85 % - jaclod 2011-07-24 86 86 % 87 % * Externalisation de la lecture et modification des arguments d'entr ee.88 % * Modification en cons equence de la demo.87 % * Externalisation de la lecture et modification des arguments d'entrée. 88 % * Modification en conséquence de la demo. 89 89 % 90 90 % - fplod 20110608T161752Z aedon.locean-ipsl.upmc.fr (Darwin) -
trunk/src/visuMSG.m
r635 r640 1 1 function [result] = visuMSG(TbMSG, Temps, Longitudes, Latitudes, timestep) 2 2 3 %visuMSG permet de faire une image des donn ï¿œes de tempï¿œratures de3 %visuMSG permet de faire une image des données de températures de 4 4 %brillances 5 5 … … 16 16 % =========== 17 17 % 18 % - Donn ï¿œes d'entrï¿œes :19 % 20 % * TbMSG : Données de temp erature de brillance d'apres MSG dans une18 % - Données d'entrées : 19 % 20 % * TbMSG : Données de température de brillance d'aprÚs MSG dans une 21 21 % matrice de dimensions temps, longitudes, latitudes. 22 22 % * Temps : Données de la dimension temporelle de la matrice en jour 23 % julien (premi ere dimension).24 % * Longitudes : Données de la premi ere dimension spatiale de la matrice23 % julien (premiÚre dimension). 24 % * Longitudes : Données de la premiÚre dimension spatiale de la matrice 25 25 % (seconde dimension). 26 26 % * Latitudes : Données de la seconde dimension spatiale de la matrice 27 27 % (troisiÚme dimension). 28 % * timestep : Num ï¿œro du quart d'heure correspondant (de 1 pour 00:00 ï¿œ29 % 96 pour 23:45 pour des données compl etes).28 % * timestep : Numéro du quart d'heure correspondant (de 1 pour 00:00 ï¿œ 29 % 96 pour 23:45 pour des données complÚtes). 30 30 % 31 31 % EXAMPLES … … 72 72 % - jaclod 2011-07-28 73 73 % 74 % * Fonctionne sur le Matlab de Cratos.74 % * Fonctionne sur le Matlab de cratos. 75 75 % 76 76 % - jaclod 2011-07-24 77 77 % 78 78 % * Externalisation de la lecture et modification des arguments d'entrée. 79 % * Modification en cons equence de la demo.79 % * Modification en conséquence de la demo. 80 80 % 81 81 % - pinsard 2011-06-08T11:29:08Z loholt1.ipsl.polytechnique.fr (Linux) … … 98 98 % - jaclod 2011-04-09 99 99 % 100 % * Cr ï¿œation100 % * Creation 101 101 % 102 102 %- … … 131 131 set(gca,'Ydir','normal'); 132 132 133 % Chargement des donn ï¿œes cotiï¿œres.133 % Chargement des données cÃŽtiÚres. 134 134 load('cmgco_traitdecote.mat','Xgco','Ygco'); 135 135 136 % Ajout d'une barre de couleur, born ï¿œe entre le minimum et le maximum des137 % temp ï¿œratures de brillances.136 % Ajout d'une barre de couleur, bornée entre le minimum et le maximum des 137 % températures de brillances. 138 138 colorbar 139 139 caxis([min(TbMSG(:)) max(TbMSG(:))]); … … 141 141 hold on; 142 142 143 % Normalisation de la taille de la figure pour obtenir un carr ï¿œ.143 % Normalisation de la taille de la figure pour obtenir un carré. 144 144 set(gca, 'DataAspectRatio', [1 1 1]); 145 145 … … 151 151 title(['Temperatures de brillance (K) d''apres MSG pour le ' datestr(Temps(timestep), 'dd mmm yyyy a HH:MM')]) 152 152 153 % Aper ï¿œu des donnï¿œes153 % Aperçu des données 154 154 imagesc(Longitudes,Latitudes,data); 155 155 156 % Ajout des c ï¿œtes.156 % Ajout des cÃŽtes. 157 157 plot(Xgco,Ygco,'k.'); 158 158 -
trunk/src/visuMSGsystconv.m
r636 r640 1 1 function [result] = visuMSGsystconv(TbMSG, Temps, Longitudes, Latitudes, timestep,Tb_SeuilSup) 2 2 3 %visuMSGsystconv permet de faire une image des donn ï¿œes de tempï¿œratures de3 %visuMSGsystconv permet de faire une image des données de températures de 4 4 %brillances 5 %pour les syst ï¿œmes convectifs potentiels dont la tempï¿œrature est6 %inf ï¿œrieure, en tout point, ï¿œla valeur7 %seuil d ï¿œfinie.5 %pour les systÚmes convectifs potentiels dont la température est 6 %inférieure, en tout point, à la valeur 7 %seuil définie. 8 8 9 9 % … … 21 21 % - Donnï¿œes d'entrï¿œes : 22 22 % 23 % * TbMSG : Données de température de brillance d'apr es MSG dans une23 % * TbMSG : Données de température de brillance d'aprÚs MSG dans une 24 24 % matrice de dimensions temps, longitudes, latitudes. 25 25 % * Temps : Données de la dimension temporelle de la matrice en jour … … 29 29 % * Latitudes : Données de la seconde dimension spatiale de la matrice 30 30 % (troisiÚme dimension). 31 % * timestep : Num ï¿œro du quart d'heure correspondant (de 1 pour 00:00 ï¿œ31 % * timestep : Numéro du quart d'heure correspondant (de 1 pour 00:00 à 32 32 % 96 pour 23:45 pour des données complÚtes). 33 % * Tb_SeuilSup : Seuil de temp ï¿œrature de brillance maximale dï¿œlimitant34 % un syst ï¿œme convectif.33 % * Tb_SeuilSup : Seuil de température de brillance maximale délimitant 34 % un systÚme convectif. 35 35 % 36 36 % EXAMPLES … … 66 66 % demo with simulation 67 67 % 68 % trouver l' ï¿œquivalent de 'alphadata'68 % trouver l'équivalent de 'alphadata' 69 69 % 70 70 % handling error … … 89 89 % - jaclod 2011-07-28 90 90 % 91 % * Fonctionne sur le Matlab de Cratos mais les cotes depassent.91 % * Fonctionne sur le Matlab de cratos mais les cÃŽtes dépassent. 92 92 % 93 93 % - jaclod 2011-07-24 94 94 % 95 95 % * Externalisation de la lecture et modification des arguments d'entrée. 96 % * Modification en cons equence de la demo.96 % * Modification en conséquence de la demo. 97 97 % 98 98 % - pinsard 2011-06-08T13:48:24Z loholt1.ipsl.polytechnique.fr (Linux) … … 116 116 % - jaclod 2011-04-10 117 117 % 118 % * Cr ï¿œation118 % * Création 119 119 % 120 120 %- … … 137 137 data = transpose(squeeze(TbMSG(timestep,:,:))); 138 138 139 % On ne garde que les syst ï¿œmes convectifs potentiels (dont les tempï¿œratures140 % de brillance sont inf ï¿œrieurs au seuil dï¿œfinie).139 % On ne garde que les systÚmes convectifs potentiels (dont les températures 140 % de brillance sont inférieurs au seuil définie). 141 141 data(data >= Tb_SeuilSup)=NaN; 142 142 … … 154 154 set(gca,'Ydir','normal'); 155 155 156 % Chargement des donn ï¿œes cotiï¿œres.156 % Chargement des données cÃŽtiÚres. 157 157 load('cmgco_traitdecote.mat','Xgco','Ygco'); 158 158 159 % Ajout d'une barre de couleur, born ï¿œe entre le minimum et le maximum des160 % temp ï¿œratures de brillances.159 % Ajout d'une barre de couleur, bornée entre le minimum et le maximum des 160 % températures de brillances. 161 161 colorbar 162 162 caxis([min(data(:)) max(data(:))]); … … 164 164 hold on; 165 165 166 % Normalisation de la taille de la figure pour obtenir un carr ï¿œ.166 % Normalisation de la taille de la figure pour obtenir un carré. 167 167 set(gca, 'DataAspectRatio', [1 1 1]); 168 168 … … 174 174 title(['Temperatures de brillance (K) d''apres MSG pour le ' datestr(Temps(timestep), 'dd mm yyyy a HH:MM')]); 175 175 176 % Aper ï¿œu des donnï¿œes de tempï¿œrature de brillance pour les systï¿œmes convectifs.176 % Aperçu des données de température de brillance pour les systÚmes convectifs. 177 177 my_image=imagesc(Longitudes,Latitudes,data); 178 178 switch application … … 185 185 end 186 186 187 % Ajout des c ï¿œtes.187 % Ajout des cÃŽtes. 188 188 plot(Xgco,Ygco,'k.'); 189 189 -
trunk/src/visumodele1.m
r636 r640 20 20 % - Données d'entrées : 21 21 % 22 % 23 % * TbMSG : Données de température de brillance d'apres MSG dans une 22 % * TbMSG : Données de température de brillance d'aprÚs MSG dans une 24 23 % matrice de dimensions temps, longitudes, latitudes. 25 24 % * Temps : Données de la dimension temporelle de la matrice en jour … … 71 70 % $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/varamma/svn/trunk/src/visumodele1.m $ 72 71 % 73 %74 72 % - fplod 20110830T074359Z cratos (Linux) 75 73 % … … 92 90 % 93 91 % * squeeze 'DetectionNuage' car de dimension 1*371*371 et doit etre 2D. 94 % * Fonctionne sur le Matlab de Cratos.92 % * Fonctionne sur le Matlab de cratos. 95 93 % 96 94 % - jaclod 2011-07-24 … … 191 189 ZoneNuageuse = zeros(size(TbMSG)); 192 190 193 194 191 %%%%%%%%%%%%%%%%%%%%%%%%%% 195 192 % Boucle … … 207 204 208 205 % Lissage des données pour avoir moins de détails (Une conversion 209 % en type "double" de la matrice est n ï¿œcessaire).206 % en type "double" de la matrice est nécessaire). 210 207 if ispc() 211 208 TbLiss_15min = smooth(double(TbMSG(timestep,:,:))); … … 219 216 disp('www : no smooth available using on non-pc platform'); 220 217 end 221 222 218 223 219 %%%%%%%%%%%%%%%%%%%%%%%%%% … … 287 283 end 288 284 289 290 285 %%%%%%%%%%%%%%%%%%%%%%%%%% 291 286 % Labellisation … … 312 307 end 313 308 314 315 309 %%%%%%%%%%%%%%%%%%%%%%%%%% 316 310 % Sélection des nuages … … 379 373 380 374 % On regarde les valeurs des différents nuages présents au 381 % dessus de ce point au cours de la journ ï¿œe.375 % dessus de ce point au cours de la journée. 382 376 Nuages3D_ij = Nuages3D(:,i,j); 383 377 … … 471 465 472 466 end 473 474 467 475 468 %%%%%%%%%%%%%%%%%%%%%%%%%% -
trunk/src/visumodele2.m
r637 r640 81 81 % 82 82 % * Externalisation de la lecture et modification des arguments d'entrée. 83 % * Modification en cons equence de la demo.84 % * Fonctionne sur Matlab Cratos.83 % * Modification en conséquence de la demo. 84 % * Fonctionne sur Matlab cratos. 85 85 % 86 86 % - fplod 20110704T104826Z aedon.locean-ipsl.upmc.fr (Darwin) … … 147 147 for num_nuage = 1 : numax; 148 148 149 % On cherche ï¿œquel quart d'heure le nuage apparait pour la premiÚre149 % On cherche à quel quart d'heure le nuage apparait pour la premiÚre 150 150 % fois. 151 151 timestep = 1; … … 154 154 end 155 155 156 % On vérifie que le nuage est pr ï¿œsent, depuis le temps oï¿œ il157 % apparait jusqu' a atteindre le seuil dï¿œfinit prï¿œcï¿œdemment.156 % On vérifie que le nuage est présent, depuis le temps oï¿œ il 157 % apparait jusqu'à atteindre le seuil définit précédemment. 158 158 Elimination = 0; 159 159 timestepSup = 0; 160 160 while timestepSup ~= Tps_SeuilInf-1 && Elimination == 0; 161 % Si parmi les quarts d'heure oï¿œ le nuage doit ï¿œtre prï¿œsent, il162 % n'est pas d ï¿œtectï¿œ, alors il ne satisfait pas le critï¿œre de161 % Si parmi les quarts d'heure oï¿œ le nuage doit être présent, il 162 % n'est pas détecté, alors il ne satisfait pas le critÚre de 163 163 % temps et n'est plus pris en compte. 164 164 if size(find(DetectionNuage(timestep+timestepSup,:,:) == num_nuage), 1) == 0; … … 192 192 %%%%%%%%%%%%%%%%%%%%%%%%%%% 193 193 194 195 194 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 196 % Mod ï¿œle de nuages ï¿œdeux couches195 % ModÚle de nuages à deux couches 197 196 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 198 197 199 % Cr ï¿œation d'une matrice qui contiendra les valeurs des Tb dans ce modÚle.198 % Création d'une matrice qui contiendra les valeurs des Tb dans ce modÚle. 200 199 Tb2couches = zeros(size(TbMSG)); 201 200 … … 207 206 Tb2couches(ZoneNoyau) = mean(TbMSG(ZoneNoyau)); 208 207 209 % De m ï¿œme pour la zone extï¿œrieure.208 % De même pour la zone extérieure. 210 209 ZoneExt = find(TbMSG > Tb_LimNoyau & TbMSG <= Tb_SeuilSup & DetectionNuage == num_nuage); 211 210 Tb2couches(ZoneExt) = mean(TbMSG(ZoneExt)); … … 214 213 215 214 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 216 217 215 218 216 %%%%%%%%%%%%% … … 229 227 230 228 %%%%%%%%%%%% 231 232 229 233 230 %%%%%%%%% -
trunk/src/visuscatternorm.m
r636 r640 268 268 269 269 result=0; 270 271 270 272 271 end -
trunk/src/write_lai_gis.m
r385 r640 44 44 % :rtype: integer 45 45 % 46 %47 46 % Write in a file 48 47 % file :file:`${PROJECT_ID}/LAI/LAI{yyyy}F_veg2010_sol_{zi_code}_write_lai_gis.csv` … … 89 88 % 90 89 % $URL: svn+ssh://pinsard@forge.ipsl.jussieu.fr/ipsl/forge/projets/varamma/svn/trunk/src/write_lai_gis.m $ 91 %92 90 % 93 91 % - fplod 20110920T091450Z aedon.locean-ipsl.upmc.fr (Darwin) -
trunk/src/write_msg_xrit_nc.py
r635 r640 145 145 ncdump -h /homedata/pinsard/varamma_d/MSG/2006/08/01/H-000-MSG1__-MSG1________-IR108____-CYCLE____-200608010000.write_msg_xrit_nc.nc 146 146 ncdump: /homedata/pinsard/varamma_d/MSG/2006/08/01/H-000-MSG1__-MSG1________-IR108____-CYCLE____-200608010000.write_msg_xrit_nc.nc: NetCDF: Unknown file format 147 148 147 149 148 and: … … 265 264 result = 0 266 265 267 268 266 # Run main, if called from the command line 269 267 if __name__ == '__main__': -
trunk/src/write_station.m
r624 r640 114 114 % 115 115 %- 116 %117 116 % 118 117 global application; -
trunk/src/write_zi.m
r635 r640 168 168 % 169 169 % +todo+ disp('arg zi_code must be a cell array of strings'); 170 %171 170 % 172 171 arg_info=whos('zi_description'); -
trunk/src/write_zi.py
r637 r640 69 69 >>> result = write_zi(west_value, east_value, north_value, south_value, zi_code, zi_description) 70 70 iii : opening for writing /usr/work/incas/fplod/varamma_d//maps/varamma_zi_write_zi_py.csv 71 72 71 73 72 From real file :file:`${PROJECT_ID}/maps/varamma_zi.csv`:
Note: See TracChangeset
for help on using the changeset viewer.