Tape5Station2.asv 1.3 KB

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