function Out= Tape5Station2(I,Temp,esp,Angle) %I代表读入的4种大气廓线,Temp代表地温(K)esp代表1-地表发射率,Angle代表180℃-扫描角 % 生成tape5文件 cd F:\Modran4 fid=fopen('tape5','w'); for i=1:4 N=I{i}; [row,col]=size(N); fprintf(fid,'TM 7 2 1 0 6 6 6 6 6 6 1 1 0%8.3f%7.2f\n',Temp,esp); % fprintf(fid,'TF 8T 0 365.00000 0 0 f t t t 0.000\n'); fprintf(fid,'TF 8T 5 365.000 1 1 f t t f\n'); fprintf(fid,'DATA/BMP99_01.BIN\n'); fprintf(fid,'F:/Modran4/CBTIRband_FB.dat\n'); fprintf(fid,' 1 1 1 10 0 0 0 0.000 0.000 0.000'); fprintf(fid,'%10.5f\n',N(1,1)); % fprintf(fid,' 4 2 0 0 0 0 0.000 0.000 0.000 0.000 0.000\n'); fprintf(fid,'%5d 0 0*I%2d\n',row,i); fprintf(fid,'%10.3f%10.3E%10.2E%10.3E%10.3E%10.3EABH66666\n',N'); fprintf(fid,' 100.000%10.3f%10.3f 0.000 0.000 0.000 0 0.000\n',N(1,1),Angle); fprintf(fid,' 700.000 1250.000 1.000 10.000RM W2aa\n'); % fprintf(fid,' 700.000 1250.000 1 2rw w'); if i<4 fprintf(fid,' 1\n'); else fprintf(fid,' 0\n'); end end fclose(fid); cd D:\Soft\Science\Matlab7\work\CE312