Changeset 41 for trunk/src/mode_sahelien/time_serie_histog_phase_num.m
- Timestamp:
- 01/28/09 11:09:56 (15 years ago)
- File:
-
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/time_serie_histog_phase_num.m
r38 r41 1 % 2 3 clear; 4 status=load('phase1.txt'); 5 phase1=status; 1 function time_serie_histog_phase_num(phase_num) 2 % time_serie_histog_phase_num ++explication between 1 and 8 3 4 %+ 5 % module 6 % ====== 7 % 8 % ++ 9 % 10 % DESCRIPTION 11 % =========== 12 % 13 % From occu.txt and phase``phase_num``.txt, time_serie_histog_phase_num plots 14 % for each WT ++. 15 % 16 % EXAMPLES 17 % ======== 18 % 19 % Following line read phase1.txt and occu.txt and plots for each WT ++ 20 % :: 21 % 22 % >> time_serie_histog_phase_num(1); 23 % 24 % SEE ALSO 25 % ======== 26 % 27 % 28 % composite_olr_phase1.m_ 29 % 30 % .. _composite_olr_phase1.m : composite_olr_phase1.m.html 31 % 32 % time_serie_histog_phase.m_ 33 % 34 % .. _time_serie_histog_phase.m : time_serie_histog_phase.m.html 35 % 36 % TODO 37 % ==== 38 % 39 % improve description 40 % 41 % improve file pb 42 % 43 % use return function 44 % 45 % EVOLUTIONS 46 % ========== 47 % 48 % $Id$ 49 % 50 % - fplod 2009-01-28T09:01:45Z aedon.locean-ipsl.upmc.fr (Darwin) 51 % 52 % * created from time_serie_histog_phase1.m to replace it as well as 53 % time_serie_histog_phase[2-8].m 54 % 55 %- 56 57 global IRCAAM_ID; 58 59 phase_num_min=1; 60 phase_num_max=8; 61 62 if nargin==1 63 % phase_num must be greater or equal to 1 and lower or equal to 8 64 if ((phase_num < phase_num_min) || (phase_num > phase_num_max)) 65 disp(['phase_num_min = ', int2str(phase_num_min)]); 66 disp(['phase_num_max = ', int2str(phase_num_max)]); 67 error('eee : phase_num must be between phase_num_min and phase_num_max'); 68 end 69 else 70 usage='usage : time_serie_histog_phase_num(phase_num)'; 71 disp(usage); 72 error('eee : wrong arguments numbers') 73 end 74 75 fullfilename=[IRCAAM_ID,'phase',int2str(phase_num),'.txt']; 76 status=load(fullfilename); 77 clear fullfilename; 78 79 phase=status; 6 80 clear status; 7 phase 1=reshape(phase1,122,28,3);81 phase=reshape(phase,122,28,3); 8 82 9 83 %%%%%%%%%%%% FICHIER TYPE DE TEMPS 10 84 % on ouvre le fichier 11 status=load('occu.txt'); 85 fullfilename=[IRCAAM_ID,'occu.txt']; 86 status=load(fullfilename); 87 clear fullfilename; 12 88 occu=status; 13 89 clear status; … … 32 108 % condition 1 on est a un maximum deltas apres sont negetifs et 33 109 % delta avant postifs 34 if phase 1(b,a,1)==1;110 if phase(b,a,1)==1; 35 111 compteurmax=compteurmax+1; 36 112 type_temps_compomax(compteurmax,1,:)=occu_type_tps(b-10,a,:); … … 95 171 end; 96 172 173 offset=mod(phase_num-1,4)*4; 97 174 for a=1:4; 98 subplot(4,4,a );175 subplot(4,4,a+offset); 99 176 diff=kara_diff(:,a)*100; 100 177 jour=[-10:10]'; … … 116 193 if a==1; 117 194 title('Atl low','Fontsize',8,'Fontname','Arial'); 118 text(-20,0,'Phase1'); 195 ylabel=['Phase ' int2str(phase_num)]; 196 text(-20,0,ylabel); 197 clear ylabel; 119 198 elseif a==2; 120 199 title('Atl Ridge','Fontsize',8,'Fontname','Arial'); … … 125 204 end; 126 205 end; 206 clear offset;
Note: See TracChangeset
for help on using the changeset viewer.