Changeset 16 for trunk/src/mode_sahelien/read_pluie_ird.m
- Timestamp:
- 01/06/09 10:53:18 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/mode_sahelien/read_pluie_ird.m
r10 r16 1 % 2 1 3 %%%%%%%%%%%%%%% FICHIER pluie IRD 2 4 % ouverture / lecture du fichier … … 19 21 for c=1:365; 20 22 if v(a,b,c)>99; 21 v(a,b,c)=NaN; 23 v(a,b,c)=NaN; 22 24 elseif v(a,b,c)<0; 23 25 v(a,b,c)=NaN; … … 36 38 for c=1:366; 37 39 if v(a,b,c)>99; 38 v(a,b,c)=NaN; 40 v(a,b,c)=NaN; 39 41 elseif v(a,b,c)<0; 40 42 v(a,b,c)=NaN; … … 53 55 for c=1:365; 54 56 if v(a,b,c)>99; 55 v(a,b,c)=NaN; 57 v(a,b,c)=NaN; 56 58 elseif v(a,b,c)<0; 57 59 v(a,b,c)=NaN; … … 70 72 for c=1:365; 71 73 if v(a,b,c)>99; 72 v(a,b,c)=NaN; 74 v(a,b,c)=NaN; 73 75 elseif v(a,b,c)<0; 74 76 v(a,b,c)=NaN; … … 87 89 for c=1:365; 88 90 if v(a,b,c)>99; 89 v(a,b,c)=NaN; 91 v(a,b,c)=NaN; 90 92 elseif v(a,b,c)<0; 91 93 v(a,b,c)=NaN; … … 106 108 for c=1:366; 107 109 if v(a,b,c)==999.99; 108 v(a,b,c)=NaN; 110 v(a,b,c)=NaN; 109 111 elseif v(a,b,c)<0; 110 112 v(a,b,c)=NaN; … … 123 125 for c=1:366; 124 126 if v(a,b,c)==999.99; 125 v(a,b,c)=NaN; 127 v(a,b,c)=NaN; 126 128 elseif v(a,b,c)<0; 127 129 v(a,b,c)=NaN; … … 140 142 for c=1:365; 141 143 if v(a,b,c)==999.99; 142 v(a,b,c)=NaN; 144 v(a,b,c)=NaN; 143 145 elseif v(a,b,c)<0; 144 146 v(a,b,c)=NaN; … … 157 159 for c=1:365; 158 160 if v(a,b,c)==999.99; 159 v(a,b,c)=NaN; 161 v(a,b,c)=NaN; 160 162 elseif v(a,b,c)<0; 161 163 v(a,b,c)=NaN; … … 174 176 for c=1:366; 175 177 if v(a,b,c)==999.99; 176 v(a,b,c)=NaN; 178 v(a,b,c)=NaN; 177 179 elseif v(a,b,c)<0; 178 180 v(a,b,c)=NaN; … … 191 193 for c=1:365; 192 194 if v(a,b,c)==999.99; 193 v(a,b,c)=NaN; 195 v(a,b,c)=NaN; 194 196 elseif v(a,b,c)<0; 195 197 v(a,b,c)=NaN; … … 208 210 for c=1:365; 209 211 if v(a,b,c)==999.99; 210 v(a,b,c)=NaN; 212 v(a,b,c)=NaN; 211 213 elseif v(a,b,c)<0; 212 214 v(a,b,c)=NaN; … … 221 223 cd('..'); 222 224 % essai de visualisation d'une carte de pluie annuelle 223 orient('landscape'); 224 lat=[0:23]; lon=[-20:26]'; 225 orient('landscape'); 226 lat=[0:23]; 227 lon=[-20:26]'; 225 228 for a=1:12; 226 229 subplot(3,4,a);
Note: See TracChangeset
for help on using the changeset viewer.