API tools faq
paste
Advertisement
Guest User

Untitled

a guest
Jun 26th, 2017
50
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 23.78 KB | None | 0 0
raw download clone embed print report
  1. %-- Unknown date --%
  2. eeze(EOL60(5,:,:));squeeze(EOL60(6,:,:));
  3. squeeze(EOL60(7,:,:));squeeze(EOL60(8,:,:));squeeze(EOL60(9,:,:));squeeze(EOL60(10,:,:));
  4. squeeze(EOL60(11,:,:));squeeze(EOL60(12,:,:));squeeze(EOL60(13,:,:));squeeze(EOL60(14,:,:));
  5. squeeze(EOL60(15,:,:));squeeze(EOL60(16,:,:));squeeze(EOL60(17,:,:));squeeze(EOL60(18,:,:));
  6. squeeze(EOL60(19,:,:));squeeze(EOL60(20,:,:));squeeze(EOL60(21,:,:));squeeze(EOL60(22,:,:));
  7. squeeze(EOL60(23,:,:));squeeze(EOL60(24,:,:));squeeze(EOL60(25,:,:));squeeze(EOL60(26,:,:));
  8. squeeze(EOL60(27,:,:));squeeze(EOL60(28,:,:));squeeze(EOL60(29,:,:))];
  9. combarr40=[squeeze(EOL40(1,:,:));squeeze(EOL40(2,:,:));
  10. squeeze(EOL40(3,:,:));squeeze(EOL40(4,:,:));squeeze(EOL40(5,:,:));squeeze(EOL40(6,:,:));
  11. squeeze(EOL40(7,:,:));squeeze(EOL40(8,:,:));squeeze(EOL40(9,:,:));squeeze(EOL40(10,:,:));
  12. squeeze(EOL40(11,:,:));squeeze(EOL40(12,:,:));squeeze(EOL40(13,:,:));squeeze(EOL40(14,:,:));
  13. squeeze(EOL40(15,:,:));squeeze(EOL40(16,:,:));squeeze(EOL40(17,:,:));squeeze(EOL40(18,:,:));
  14. squeeze(EOL40(19,:,:));squeeze(EOL40(20,:,:));squeeze(EOL40(21,:,:));squeeze(EOL40(22,:,:));
  15. squeeze(EOL40(23,:,:));squeeze(EOL40(24,:,:));squeeze(EOL40(25,:,:));squeeze(EOL40(26,:,:));
  16. squeeze(EOL40(27,:,:));squeeze(EOL40(28,:,:));squeeze(EOL40(29,:,:))];
  17. cd ..
  18. % figure(1)
  19. MODELarr80=combarr80(7:649,:); %80m
  20. MODELarr60=combarr60(7:649,:); %60m
  21. MODELarr40=combarr40(7:649,:); %40m
  22. [mobs,mest,sdobs,sdest,R2,imae,rmsee,E,Elog,E1,d,d1]=goodfit2(RAWarr80(:,7),MODELarr80(:,10)*0.514444445);
  23. rmsetext=num2str(rmsee);iagreetxt=num2str(d);
  24. plot(MODELarr80(:,10)*0.514444445,'*b');hold on;plot(RAWarr80(:,7),'*r');
  25. legend('Model Simulation (m/s)','Anemometer (m/s)');
  26. xlabel('Elapsed time in hours starting at 00Z 9 Feb 2007','fontsize',16);
  27. ylabel('Wind speed (m/s)','fontsize',16);
  28. title('00Z Initialized [NARR] - Wind speed comparison at 80m height [2km RES]','fontsize',16);
  29. textul({['RMSE = ' rmsetext] ['IoA = ' iagreetxt]});
  30. screen_size = get(0, 'ScreenSize');
  31. set(gcf,'paperpositionmode','auto')
  32. f1 = figure(1);
  33. set(f1, 'Position', [0 0 screen_size(3) screen_size(4) ] );
  34. saveas(gcf,['00Z Initialized [NARR] - Wind speed comparison at 80m height [2km RES]'],'png');
  35. [mobs,mest,sdobs,sdest,R2,imae,rmsee,E,Elog,E1,d,d1]=goodfit2(RAWarr60(:,7),MODELarr60(:,10)*0.514444445);
  36. rmsetext=num2str(rmsee);iagreetxt=num2str(d);
  37. plot(MODELarr60(:,10)*0.514444445,'*b');hold on;plot(RAWarr60(:,7),'*r');
  38. legend('Model Simulation (m/s)','Anemometer (m/s)');
  39. xlabel('Elapsed time in hours starting at 00Z 9 Feb 2007','fontsize',16);
  40. ylabel('Wind speed (m/s)','fontsize',16);
  41. title('00Z Initialized [NARR] - Wind speed comparison at 60m height [2km RES]','fontsize',16);
  42. textul({['RMSE = ' rmsetext] ['IoA = ' iagreetxt]});
  43. screen_size = get(0, 'ScreenSize');
  44. set(gcf,'paperpositionmode','auto')
  45. f1 = figure(1);
  46. set(f1, 'Position', [0 0 screen_size(3) screen_size(4) ] );
  47. saveas(gcf,['00Z Initialized [NARR] - Wind speed comparison at 60m height [2km RES]'],'png');
  48. [mobs,mest,sdobs,sdest,R2,imae,rmsee,E,Elog,E1,d,d1]=goodfit2(RAWarr40(:,7),MODELarr40(:,10)*0.514444445);
  49. rmsetext=num2str(rmsee);iagreetxt=num2str(d);
  50. plot(MODELarr40(:,10)*0.514444445,'*b');hold on;plot(RAWarr40(:,7),'*r');
  51. legend('Model Simulation (m/s)','Anemometer (m/s)');
  52. xlabel('Elapsed time in hours starting at 00Z 9 Feb 2007','fontsize',16);
  53. ylabel('Wind speed (m/s)','fontsize',16);
  54. title('00Z Initialized [NARR] - Wind speed comparison at 40m height [2km RES]','fontsize',16);
  55. textul({['RMSE = ' rmsetext] ['IoA = ' iagreetxt]});
  56. screen_size = get(0, 'ScreenSize');
  57. set(gcf,'paperpositionmode','auto')
  58. f1 = figure(1);
  59. set(f1, 'Position', [0 0 screen_size(3) screen_size(4) ] );
  60. saveas(gcf,['00Z Initialized [NARR] - Wind speed comparison at 40m height [2km RES]'],'png');
  61. save temp
  62. clear all
  63. %-- 11/4/10 3:54 PM --%
  64. ncload UWND.CDF
  65. ncload VWND.cdf
  66. [lon,lat]=meshgrid[X,Y];
  67. [lon,lat]=meshgrid(X,Y);
  68. lon
  69. size(u)
  70. size(v)
  71. v=reshape(v,2760/30,30,17,29,67)
  72. v=reshape(v,(2760/30),30,17,29,67)
  73. 2760/30
  74. v=reshape(v,(2760/30),30,17,29,57)
  75. v=reshape(v,(2760/30),30,17,29,57);
  76. size(v)
  77. vclim=squeeze(nanmean(vclim,2));
  78. vclim=squeeze(nanmean(v,2));
  79. size(vclim)
  80. vevent=v(29:40,25,:,:,:);
  81. vevent=squeeze(vevent);
  82. size(vevent)
  83. for i=1:12 vanom(i,:,:,:)=vevent(i,:,:,:)-vclim(28+i,:,:,:); endl
  84. for i=1:12 vanom(i,:,:,:)=vevent(i,:,:,:)-vclim(28+i,:,:,:); end;
  85. size(vanom)
  86. hist(vanom)
  87. hist(vanom(:))
  88. hist(vanom(:),30)
  89. a = peaks;
  90. [trash c]=contourf(a);
  91. set(get(c,'Children'),'LineStyle','none')
  92. a = peaks;
  93. [trash c]=contourf(a);
  94. A=peaks;
  95. minA=min(A(:));
  96. level=linspace(minA,0,10);
  97. contourf(A,level);
  98. A=peaks;
  99. [C,h] = contourf(A,20);
  100. p=get(h,'Children');
  101. c=get(p,'Cdata');
  102. set(p([c{:}]<=0),'LineStyle','-.')
  103. A=peaks;
  104. [C,h] = contourf(A,20);
  105. p=get(h,'Children');
  106. c=get(p,'Cdata');
  107. set(p([c{:}]<=0),'LineStyle','.')
  108. set(p([c{:}]<=0),'LineStyle','-.')
  109. A=peaks;
  110. [C,h] = contourf(A,20);
  111. p=get(h,'Children');
  112. c=get(p,'Cdata');
  113. set(p([c{:}]<=0),'LineStyle','--')
  114. size(vanom)
  115. Pre=300;
  116. abs(min(min(min(min(vanom(:,find(Pr==Pre),:,:))))))
  117. abs(min(min(min(min(vanom(:,find(P==Pre)
  118. abs(min(min(min(min(vanom(:,find(P==Pre),:,:))))))
  119. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  120. Pr=300;
  121. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  122. for i=1:12
  123. subplot(4,3,i);Pr=300;
  124. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  125. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  126. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  127. else
  128. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  129. end
  130. makemap(lon,lat,vanom(i,find(P==Pr),:,:),[-lim lim]);
  131. end
  132. for i=1:12
  133. subplot(4,3,i);Pr=300;
  134. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  135. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  136. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  137. else
  138. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  139. end
  140. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  141. end
  142. for i=1:12
  143. subplot(4,3,i);Pr=300;
  144. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  145. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  146. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  147. else
  148. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  149. end
  150. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  151. subplotcolorbar;
  152. end
  153. for i=1:12
  154. subplot(4,3,i);Pr=300;
  155. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  156. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  157. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  158. else
  159. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  160. end
  161. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  162. end
  163. subplotcolorbar;
  164. for i=1:12
  165. subplot(4,3,i);Pr=300;
  166. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  167. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  168. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  169. else
  170. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  171. end
  172. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  173. end
  174. subplotcolorbar('V-wind Anom at 300hPa');
  175. clf
  176. for i=1:12
  177. subplot(4,3,i);Pr=300;
  178. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  179. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  180. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  181. else
  182. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  183. end
  184. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  185. end
  186. subindex;subplotcolorbar('V-wind Anom at 300hPa');
  187. edit subindex
  188. for i=1:12
  189. subplot(4,3,i);Pr=300;
  190. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  191. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  192. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  193. else
  194. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  195. end
  196. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  197. end
  198. subplotspace('horizontal',45);subplotspace('vertical',15);
  199. subindex;subplotcolorbar('V-wind Anom (m/s)');
  200. clf
  201. for i=1:12
  202. subplot(4,3,i);Pr=300;
  203. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  204. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  205. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  206. else
  207. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  208. end
  209. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  210. end
  211. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  212. subindex;subplotcolorbar('V-wind Anom (m/s)');
  213. for i=1:12
  214. subplot(4,3,i);Pr=300;
  215. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  216. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  217. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  218. else
  219. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  220. end
  221. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  222. end
  223. mtit([num2str(Pr) 'hPa V-Wind Anomalies');
  224. subindex;subplotcolorbar('V-wind Anom (m/s)');
  225. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  226. clf
  227. for i=1:12
  228. subplot(4,3,i);Pr=300;
  229. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  230. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  231. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  232. else
  233. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  234. end
  235. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  236. end
  237. mtit([num2str(Pr) 'hPa V-Wind Anomalies');
  238. subindex;subplotcolorbar('V-wind Anom (m/s)');
  239. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  240. for i=1:12
  241. subplot(4,3,i);Pr=300;
  242. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  243. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  244. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  245. else
  246. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  247. end
  248. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  249. end
  250. subindex;subplotcolorbar('V-wind Anom (m/s)');
  251. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  252. mtit([num2str(Pr) 'hPa V-Wind Anomalies');
  253. edit mtit
  254. mtit([num2str(Pr) 'hPa V-Wind Anomalies']);
  255. mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);
  256. for i=1:12
  257. subplot(4,3,i);Pr=300;
  258. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  259. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  260. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  261. else
  262. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  263. end
  264. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  265. end
  266. subindex;subplotcolorbar('V-wind Anom (m/s)');
  267. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  268. mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);
  269. for i=1:12
  270. subplot(4,3,i);Pr=300;
  271. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  272. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  273. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  274. else
  275. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  276. end
  277. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  278. end
  279. subindex;
  280. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  281. mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);subplotcolorbar('V-wind Anom (m/s)');
  282. figure(2)
  283. for i=1:12
  284. subplot(4,3,i);Pr=250;
  285. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  286. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  287. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  288. else
  289. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  290. end
  291. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  292. end
  293. subindex;
  294. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  295. mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);subplotcolorbar('V-wind Anom (m/s)');
  296. figure(3)
  297. for i=1:12
  298. subplot(4,3,i);Pr=500;
  299. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  300. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  301. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  302. else
  303. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  304. end
  305. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  306. end
  307. subindex;
  308. subplotspace('horizontal',-45);subplotspace('vertical',-15);
  309. mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);subplotcolorbar('V-wind Anom (m/s)');
  310. cd 2km_00Z_NARR
  311. x=['20070209_80m_2km_00Z_narr.xls';'20070210_80m_2km_00Z_narr.xls';'20070211_80m_2km_00Z_narr.xls';'20070212_80m_2km_00Z_narr.xls';'20070213_80m_2km_00Z_narr.xls';
  312. '20070214_80m_2km_00Z_narr.xls';'20070215_80m_2km_00Z_narr.xls';'20070216_80m_2km_00Z_narr.xls';'20070217_80m_2km_00Z_narr.xls';
  313. '20070218_80m_2km_00Z_narr.xls';'20070219_80m_2km_00Z_narr.xls';'20070220_80m_2km_00Z_narr.xls';'20070221_80m_2km_00Z_narr.xls';
  314. '20070222_80m_2km_00Z_narr.xls';'20070223_80m_2km_00Z_narr.xls';'20070224_80m_2km_00Z_narr.xls';'20070225_80m_2km_00Z_narr.xls';
  315. '20070226_80m_2km_00Z_narr.xls';'20070227_80m_2km_00Z_narr.xls';'20070228_80m_2km_00Z_narr.xls';'20070301_80m_2km_00Z_narr.xls';
  316. '20070302_80m_2km_00Z_narr.xls';'20070303_80m_2km_00Z_narr.xls';'20070304_80m_2km_00Z_narr.xls';'20070305_80m_2km_00Z_narr.xls';
  317. '20070306_80m_2km_00Z_narr.xls';'20070307_80m_2km_00Z_narr.xls';'20070308_80m_2km_00Z_narr.xls';'20070309_80m_2km_00Z_narr.xls';];
  318. y=['20070209_60m_2km_00Z_narr.xls';'20070210_60m_2km_00Z_narr.xls';'20070211_60m_2km_00Z_narr.xls';'20070212_60m_2km_00Z_narr.xls';'20070213_60m_2km_00Z_narr.xls';
  319. '20070214_60m_2km_00Z_narr.xls';'20070215_60m_2km_00Z_narr.xls';'20070216_60m_2km_00Z_narr.xls';'20070217_60m_2km_00Z_narr.xls';
  320. '20070218_60m_2km_00Z_narr.xls';'20070219_60m_2km_00Z_narr.xls';'20070220_60m_2km_00Z_narr.xls';'20070221_60m_2km_00Z_narr.xls';
  321. '20070222_60m_2km_00Z_narr.xls';'20070223_60m_2km_00Z_narr.xls';'20070224_60m_2km_00Z_narr.xls';'20070225_60m_2km_00Z_narr.xls';
  322. '20070226_60m_2km_00Z_narr.xls';'20070227_60m_2km_00Z_narr.xls';'20070228_60m_2km_00Z_narr.xls';'20070301_60m_2km_00Z_narr.xls';
  323. '20070302_60m_2km_00Z_narr.xls';'20070303_60m_2km_00Z_narr.xls';'20070304_60m_2km_00Z_narr.xls';'20070305_60m_2km_00Z_narr.xls';
  324. '20070306_60m_2km_00Z_narr.xls';'20070307_60m_2km_00Z_narr.xls';'20070308_60m_2km_00Z_narr.xls';'20070309_60m_2km_00Z_narr.xls';];
  325. z=['20070209_60m_2km_00Z_narr.xls';'20070210_40m_2km_00Z_narr.xls';'20070211_40m_2km_00Z_narr.xls';'20070212_40m_2km_00Z_narr.xls';'20070213_40m_2km_00Z_narr.xls';
  326. '20070214_40m_2km_00Z_narr.xls';'20070215_40m_2km_00Z_narr.xls';'20070216_40m_2km_00Z_narr.xls';'20070217_40m_2km_00Z_narr.xls';
  327. '20070218_40m_2km_00Z_narr.xls';'20070219_40m_2km_00Z_narr.xls';'20070220_40m_2km_00Z_narr.xls';'20070221_40m_2km_00Z_narr.xls';
  328. '20070222_40m_2km_00Z_narr.xls';'20070223_40m_2km_00Z_narr.xls';'20070224_40m_2km_00Z_narr.xls';'20070225_40m_2km_00Z_narr.xls';
  329. '20070226_40m_2km_00Z_narr.xls';'20070227_40m_2km_00Z_narr.xls';'20070228_40m_2km_00Z_narr.xls';'20070301_40m_2km_00Z_narr.xls';
  330. '20070302_40m_2km_00Z_narr.xls';'20070303_40m_2km_00Z_narr.xls';'20070304_40m_2km_00Z_narr.xls';'20070305_40m_2km_00Z_narr.xls';
  331. '20070306_40m_2km_00Z_narr.xls';'20070307_40m_2km_00Z_narr.xls';'20070308_40m_2km_00Z_narr.xls';'20070309_40m_2km_00Z_narr.xls';];
  332. for i=1:29
  333. M80(i,:,:) = dlmread(x(i,:),'',1,0);
  334. M60(i,:,:) = dlmread(y(i,:),'',1,0);
  335. M40(i,:,:) = dlmread(z(i,:),'',1,0);
  336. end
  337. EOL80=M80(:,1:2:48,:);
  338. EOL60=M60(:,1:2:48,:);
  339. EOL40=M40(:,1:2:48,:);
  340. combarr80=[squeeze(EOL80(1,:,:));squeeze(EOL80(2,:,:));
  341. squeeze(EOL80(3,:,:));squeeze(EOL80(4,:,:));squeeze(EOL80(5,:,:));squeeze(EOL80(6,:,:));
  342. squeeze(EOL80(7,:,:));squeeze(EOL80(8,:,:));squeeze(EOL80(9,:,:));squeeze(EOL80(10,:,:));
  343. squeeze(EOL80(11,:,:));squeeze(EOL80(12,:,:));squeeze(EOL80(13,:,:));squeeze(EOL80(14,:,:));
  344. squeeze(EOL80(15,:,:));squeeze(EOL80(16,:,:));squeeze(EOL80(17,:,:));squeeze(EOL80(18,:,:));
  345. squeeze(EOL80(19,:,:));squeeze(EOL80(20,:,:));squeeze(EOL80(21,:,:));squeeze(EOL80(22,:,:));
  346. squeeze(EOL80(23,:,:));squeeze(EOL80(24,:,:));squeeze(EOL80(25,:,:));squeeze(EOL80(26,:,:));
  347. squeeze(EOL80(27,:,:));squeeze(EOL80(28,:,:));squeeze(EOL80(29,:,:))];
  348. combarr60=[squeeze(EOL60(1,:,:));squeeze(EOL60(2,:,:));
  349. squeeze(EOL60(3,:,:));squeeze(EOL60(4,:,:));squeeze(EOL60(5,:,:));squeeze(EOL60(6,:,:));
  350. squeeze(EOL60(7,:,:));squeeze(EOL60(8,:,:));squeeze(EOL60(9,:,:));squeeze(EOL60(10,:,:));
  351. squeeze(EOL60(11,:,:));squeeze(EOL60(12,:,:));squeeze(EOL60(13,:,:));squeeze(EOL60(14,:,:));
  352. squeeze(EOL60(15,:,:));squeeze(EOL60(16,:,:));squeeze(EOL60(17,:,:));squeeze(EOL60(18,:,:));
  353. squeeze(EOL60(19,:,:));squeeze(EOL60(20,:,:));squeeze(EOL60(21,:,:));squeeze(EOL60(22,:,:));
  354. squeeze(EOL60(23,:,:));squeeze(EOL60(24,:,:));squeeze(EOL60(25,:,:));squeeze(EOL60(26,:,:));
  355. squeeze(EOL60(27,:,:));squeeze(EOL60(28,:,:));squeeze(EOL60(29,:,:))];
  356. combarr40=[squeeze(EOL40(1,:,:));squeeze(EOL40(2,:,:));
  357. squeeze(EOL40(3,:,:));squeeze(EOL40(4,:,:));squeeze(EOL40(5,:,:));squeeze(EOL40(6,:,:));
  358. squeeze(EOL40(7,:,:));squeeze(EOL40(8,:,:));squeeze(EOL40(9,:,:));squeeze(EOL40(10,:,:));
  359. squeeze(EOL40(11,:,:));squeeze(EOL40(12,:,:));squeeze(EOL40(13,:,:));squeeze(EOL40(14,:,:));
  360. squeeze(EOL40(15,:,:));squeeze(EOL40(16,:,:));squeeze(EOL40(17,:,:));squeeze(EOL40(18,:,:));
  361. squeeze(EOL40(19,:,:));squeeze(EOL40(20,:,:));squeeze(EOL40(21,:,:));squeeze(EOL40(22,:,:));
  362. squeeze(EOL40(23,:,:));squeeze(EOL40(24,:,:));squeeze(EOL40(25,:,:));squeeze(EOL40(26,:,:));
  363. squeeze(EOL40(27,:,:));squeeze(EOL40(28,:,:));squeeze(EOL40(29,:,:))];
  364. cd ..
  365. % figure(1)
  366. MODELarr80=combarr80(7:649,:); %80m
  367. MODELarr60=combarr60(7:649,:); %60m
  368. MODELarr40=combarr40(7:649,:); %40m
  369. for i=1:12
  370. % subplot(4,3,i);
  371. Pr=500;
  372. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  373. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  374. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  375. else
  376. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  377. end
  378. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  379. colorbar;
  380. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');
  381. end
  382. for i=1:12
  383. % subplot(4,3,i);
  384. Pr=500;
  385. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  386. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  387. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  388. else
  389. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  390. end
  391. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  392. colorbar;
  393. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');
  394. end
  395. for i=1:12
  396. % subplot(4,3,i);
  397. Pr=500;
  398. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  399. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  400. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  401. else
  402. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  403. end
  404. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  405. colorbar;
  406. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');pause(.5);clf;
  407. end
  408. for i=1:12
  409. % subplot(4,3,i);
  410. Pr=300;
  411. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  412. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  413. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  414. else
  415. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  416. end
  417. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  418. colorbar;
  419. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');pause(.5);clf;
  420. end
  421. for i=1:12
  422. % subplot(4,3,i);
  423. Pr=250;
  424. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  425. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  426. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  427. else
  428. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  429. end
  430. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  431. colorbar;
  432. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');pause(.5);clf;
  433. end
  434. % subindex;
  435. % subplotspace('horizontal',-45);subplotspace('vertical',-15);
  436. % mtit([num2str(Pr) 'hPa V-Wind Anomalies'],'fontsize',20);
  437. % subplotcolorbar('V-wind Anom (m/s)');
  438. for i=1
  439. % subplot(4,3,i);
  440. Pr=500;
  441. if abs(min(min(min(min(vanom(:,find(P==Pr),:,:)))))) >...,
  442. abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))))
  443. lim=abs(min(min(min(min(vanom(:,find(P==Pr),:,:))))));
  444. else
  445. lim=abs(max(max(max(max(vanom(:,find(P==Pr),:,:))))));
  446. end
  447. makemap(lon,lat,squeeze(vanom(i,find(P==Pr),:,:)),[-lim lim]);
  448. colorbar;
  449. saveas(gcf,[num2str(i-1+29) num2str(Pr)],'png');pause(.5);clf;
  450. end
  451. save temp X Y P lat lon v vanom vclim vevent
  452. %-- 11/11/10 12:00 PM --%
  453. 360/12
  454. precip=precip(1:360,:,:);
  455. precip=nmean(precip(:,fin),2);
  456. plot(ao,precip)
  457. plot(ao,precip,'r*')
  458. set(gca,'fontsize',16)
  459. hist(precip)
  460. clear all
  461. precip=precip(1:360,:,:);
  462. for i=1:360 panom(i,:,:)=precip(i,:,:)-precipmean(:,:); end
  463. for i=1:360 panom(i,:,:)=squeeze(precip(i,:,:))-precipmean(:,:); end
  464. hist(panom)
  465. hist(panom(:))
  466. panomdom=panom(:,fin);
  467. panomdom=squeeze(nmean(panomdom,2));
  468. hist(panomdom(:))
  469. plot(ao,panomdom,'r.')
  470. set(gca,'fontsize',16)
  471. plot(ao,panomdom,'r*')
  472. plot(ao,panomdom,'b*')
  473. plot(ao,panomdom,'o*')
  474. plot(ao,panomdom,'bo')
  475. xlabel('Monthly AO Index')
  476. ylabel('Monthly Precipitation Anomalies for AZ/NM Domain')
  477. set(gca,'fontsize',18)
  478. set(gca,'fontsize',12)
  479. set(gca,'fontsize',16)
  480. xlabel('Monthly AO Index')
  481. hist(panomdom(:))
  482. plot(ao,panomdom,'bo')
  483. xlabel('Monthly AO Index')
  484. ylabel('Monthly Precipitation Anomalies for AZ/NM Domain')
  485. set(gca,'fontsize',16)
  486. precipdom=precip(360,fin);
  487. precipdom=precip(:,fin);
  488. precipdom=squeeze(nmean(precipdom(:),2))
  489. precipdom=precip(:,fin);
  490. precipdom=squeeze(nmean(precipdom,2));
  491. twit('Yea thats right I just updated my status via MATLAB... not nerdy enough','andyknows','najnaj12'
  492. twit('Yea thats right I just updated my status via MATLAB... not nerdy enough','andyknows','najnaj12');
  493. twit('Yea thats right I just updated my status via MATLAB... not nerdy enough',andyknows,najnaj12);
  494. figure(2);
  495. plot(ao,precipdom,'bo')
  496. set(gca,'fontsize',16)
  497. xlabel('Monthly AO Index')
  498. ylabel('Monthly Precipitation for AZ/NM Domain')
  499. twitter('Yea thats right I just updated my status via MATLAB... not nerdy enough');
  500. ncdump('precip.cdf')
  501. ncdump('ALLprecip.cdf')
  502. x
  503. y
  504. lon
  505. ncdump('temp.cdf')
  506. save temp
  507. clear all
  508. ncload temp.cdf
  509. X
  510. Y
  511. load temp2 fin
  512. load temp fin
  513. [lon,lat]=meshgrid(X,Y);
  514. lon
  515. BX=[-103 -114.5 -114.5 -103 -103];BY=[32.5 32.5 38 38 32.5];
  516. clear fin
  517. in=inpolygon(lon,lat,BX,BY);
  518. fin=find(in==1);
  519. tempdom=tempdom(:,fin);
  520. tempdom=temp(:,fin);
  521. tempavg=nmean(tempdom,2);
  522. figure(3);
  523. plot(ao,tempavg,'bo')
  524. xlabel('Monthly AO Index')
  525. ylabel('Monthly Temperature for AZ/NM Domain')
  526. plot(ao,tempavg-273,'bo')
  527. plot(ao(1:359),tempavg(2:360)-273,'bo')
  528. corrcoef(ao(1:359),tempavg(2:360)-273)
  529. corrcoef(ao(:),tempavg(:)-273)
  530. corrcoef(ao(1:358),tempavg(3:360)-273)
  531. corrcoef(ao(1:357),tempavg(4:360)-273)
  532. corrcoef(ao(1:356),tempavg(5:360)-273)
  533. corrcoef(ao(1:355),tempavg(6:360)-273)
  534. [c,p]=corrcoef(ao(1:355),tempavg(6:360)-273)
  535. [c,p]=corrcoef(ao(1:354),tempavg(7:360)-273)
  536. [c,p]=corrcoef(ao(1:353),tempavg(8:360)-273)
  537. [c,p]=corrcoef(ao(:),tempavg(:)-273)
  538. [c,p]=corrcoef(ao(1:359),tempavg(2:360)-273)
  539. figure(3)
  540. plot(ao(1:354),tempavg(7:360)-273,'bo')
  541. save tempTEMPERATURE
  542. clear all
  543. figure(4)
  544. bigarr=[ao precipdom];
  545. p=anova1(bigarr)
  546. p=anova2(bigarr,2)
  547. p=anova1(bigarr)
  548. [c,p]=corrcoef(ao(1:359),precipdom(2:360))
  549. [c,p]=corrcoef(ao(1:358),precipdom(3:360))
  550. [c,p]=corrcoef(ao(1:357),precipdom(4:360))
  551. edit xcorr.m
  552. clear c
  553. c=xcorr(ao,precipdom);
  554. c
  555. plot(c)
  556. c=xcorr(ao(:),precipdom(:),[0 8]);
  557. c=xcorr(ao(:),precipdom(:),[0:8]);
  558. c=xcorr(ao,precipdom,[0:8]);
  559. c=xcorr(ao,precipdom,1);
  560. size(c)
  561. c=xcorr(ao,precipdom,5);
  562. size(c)
  563. c=xcorr(ao,precipdom,7);
  564. figure(6)
  565. plot(c)
  566. c=xcorr(ao,precipdom,20);
  567. figure(6)
  568. plot(c)
  569. c=xcorr(ao,precipdom,60);
  570. plot(c)
  571. c=xcorr(ao,precipdom,60,'coeff');
  572. plot(c)
  573. clf
  574. figure(6)
  575. clf
  576. c=xcorr(ao,precipdom,60,'coeff');
  577. plot(c)
  578. test2=[1 2 3 4 5 3 9 8 7 6 20 9 8 6 89];test3=[1 2 3 4 5 6 7 8 9 2 3 4 6 1 123 0];
  579. test2=[1 2 3 4 5 3 9 8 7 6 20 9 8 6 89 2];test3=[1 2 3 4 5 6 7 8 9 2 3 4 6 1 123 0];
  580. size(c)
  581. c=xcorr(test2,test3,5,'coeff');
  582. figure(6)
  583. plot(c)
  584. c=xcorr(test2,test3,5);
  585. figure(6)
  586. c=xcorr(test2,test3,5);
  587. plot(c)
  588. c=xcorr(test2,test3,5,'coeff');
  589. plot(c)
  590. c=xcorr(ao,precipdom,60,'coeff');
  591. subplot(2,1,1);plot(c)
  592. c=xcorr(ao,precipdom,10,'coeff');
  593. c=xcorr(ao,precipdom,12,'coeff');
  594. figure(6)
  595. c=xcorr(ao,precipdom,10,'coeff');
  596. plot(c)
  597. subplot(2,1,2);
  598. c=xcorr(precipdom,ao,10,'coeff');
  599. plot(c)
  600. c=xcorr(ao,precipdom,3,'coeff');
  601. subplot(2,1,1);plot(c)
  602. c=xcorr(ao,precipdom,6,'coeff');
  603. subplot(2,1,1);plot(c)
  604. clf
  605. c=xcorr(ao,precipdom,7,'coeff');
  606. figure(6)
  607. plot(c)
  608. clf
  609. figure(6)
  610. crosscorr(ao,precipdom,10)
  611. grid off
  612. corrcoef(ao(1:359),precipdom(2:360))
  613. corrcoef(ao(1:358),precipdom(3:360))
  614. save temp
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!