12345678910111213141516171819202122232425262728293031 |
- 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
|