clear?all;?close?all;?clc;
tic
for?yearnum=1979:2017
????cd（[‘F:\pythondownload\separate\‘num2str（yearnum）]）???????????????????%打開年份文件夾
????for?monthnum=6:12
????cd（[‘F:\pythondownload\separate\‘num2str（yearnum）]）;?
???????if?monthnum<10??????????????????????????????????????????????????????
??????????filename{monthnum}=[num2str（yearnum）‘0‘num2str（monthnum）];???
???????else
??????????filename{monthnum}=[num2str（yearnum）num2str（monthnum）];
???????end
???????filen1=strcat（[filename{monthnum}‘.nc‘]）;??????????????????????????%讀入文件名
???????ncid?=?netcdf.open（filen1‘NC_NOWRITE‘）;????????????????????????????%打開文件
???????lat=double（ncread（filen1‘latitude‘））;??????????????????????????????%讀入‘2009年1月有效波高.nc‘文件中的緯度，變量名字不能改
???????lon=double（ncread（filen1‘longitude‘））;?????????????????????????????%讀入經(jīng)度
???????[lon1lat1]=meshgrid（lonlat）;
???????%?????[lat1lon1]=meshgrid（latlon）;
???????U10=ncread（filen1‘u10‘）;???????????????????????????????????????????%讀入有效波高
???????V10=ncread（filen1‘v10‘）;???????????????????????????????????????????%讀入有效波高
???????
???????lonswn?=?[118:0.125:124];?%29?118?0?64?48?0.125?0.125
???????latswn?=?[29:0.125:37];
???????[lon2lat2]?=?meshgrid（lonswnlatswn）;
???????
???????TimeData=ncread（[filename{monthnum}‘.nc‘]‘time‘）;
???????disp（strcat（‘Done?reading‘filen1））;????????????????????????????????%顯示功能顯示所需要的數(shù)據(jù)變量等已經(jīng)輸入完畢
???????%????A=u10（5050:）
???????%????B=shiftdim（A2）
???????T=double（TimeData）;
???????L=（T-T（1））/24;
???????R=datenum（yearnummonthnum1000）;
???????format?long;
???????%????BEGIN=（datenum（20111213000）-R）*24;
???????%????END=（datenum（201211201200）-R）*24;
???????K=R+L;
???????DateString?=?datestr（K‘yyyymmddHH‘）;
???????if?mod（yearnum4）==0
???????????monthFeb=29;
???????else
???????????monthFeb=28;
???????end
???????????if?monthnum==1?||?monthnum==3?||?monthnum==5?||?monthnum==7?||?monthnum==8?||?monthnum==10?||?monthnum==12
???????????????daynum=31;
???????????elseif?monthnum==2
???????????????daynum=monthFeb;
???????????else
???????????????daynum=30;
???????????end
???????for?i=1:（daynum*4）
???????????aa（i:）=strcat（DateString（i:）‘.txt‘）;
???????end
???????%????aaa=DateString（1:）
???????%????DateString?=?datestr（989550/24+R‘yyyymmddHH‘）
???????
???????cd（[‘F:\pythondownload\windtxt\‘num2str（yearnum）‘\‘filename{monthnum}]）;
???????for?i=1:（daynum*4）
???????????w=[];
???????????name=strcat（DateString（i:）‘.txt‘）;
???????????curfold?=pwd;
???????????fid=fopen（strcat（curfold‘\‘name）‘wt+‘）;
???????????u=double（U10（::i）‘）;
???????????v=double（V10（::i）‘）;
???????????
???????????uq?=?griddata（lon1lat1ulon2lat2）;?uq（isnan（uq））=0;
???????????vq?=?griddata（lon1lat1vlon2lat2）;?vq（isnan（vq））=0;
???????????w?=?[uq;vq];
???????????dlmwrite（strcat（curfold‘\‘name）w‘delimiter‘‘\t‘）;
???????????fclose（fid）;
???????????x=w‘;
???????????%bb（i1:2）=x;
???????????
???????end
???????toc
???????clearvars?-exce

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????3425??2018-08-08?18:58??duqu_ecmwf.m

?????文件???????2069??2018-07-25?18:14??duqu_ncep.m

?????文件???????1700??2018-09-01?17:30??usemakedir.m

?????文件???????1365??2018-08-15?17:27??windin.m

-----------?---------??----------?-----??----

?????????????????8559????????????????????4