6S模型IDL调用

最新推荐文章于 2023-05-17 11:31:43 发布

ρ∧ρ

最新推荐文章于 2023-05-17 11:31:43 发布

阅读量938

点赞数 7

文章标签：开发语言

本文链接：https://blog.csdn.net/weixin_71224572/article/details/130477686

版权

pro s6
n=0
m=0
y=0
for j=1,3 do begin
if j eq 1 then begin
  aod_toa=fltarr(13)
  aod_toa_1=fltarr(13)
  for i=0.0,3.0,0.25 do begin
    aod_toa[n]=output(i,0.1,0.0)
    aod_toa_1[n]=i
;    print,aod_toa_1[n],aod_toa[n]
    n=n+1
  endfor
p=plot(aod_toa_1,aod_toa,xtitle='AOD',ytitle='TOA',title='大气厚度对表观反射率的影响')
endif
  if j eq 2 then begin
    arc_toa=fltarr(15)
    arc_toa_1=fltarr(15)
    for i=0,70,5 do begin
      arc_toa[m]=output(1,0.1,i)
      arc_toa_1[m]=i
      m=m+1
    endfor
    q=plot(arc_toa_1,arc_toa,xtitle='ARC',ytitle='TOA',title='天顶角对表观反射率的影响')
  endif
  if j eq 3 then begin
    ref_toa=fltarr(11)
    ref_toa_1=fltarr(11)
    for i=0.0,0.5,0.05 do begin
      ref_toa[y]=output(1,i,45)
      ref_toa_1[y]=i
      y=y+1
    endfor
    o=plot(ref_toa_1,ref_toa,xtitle='REF',ytitle='TOA',title='地表反射率对表观反射率的影响')
  endif
endfor
end

;x=[1,2,3,4]
;y=[2,3,4,5]
;plot,x,y,xtitle='x',ytitle='y',title='xy'
;plot(x,y,xtitle='x',ytitle='y',title='xy')


;为参数定义初值
;AOD=0.25
;REF=0.1
;ARC=0.0
function output,AOD,REF,ARC
  outfilename='E:\yaoganhomework\ch02_6s_test\ch02_6s_test\input.txt'
  openw,var_lun,outfilename,/get_lun
  printf,var_lun,7,format='(i-2)'
  printf,var_lun,ARC,24,2.72,116.381,39.977,format='(i-3.2,i-3,f4.2,f8.3,f7.3)'
  printf,var_lun,8,format='(i1)'
  printf,var_lun,1.366,0.29,format='(f5.3,f5.2)'
  printf,var_lun,1,format='(i1)'
  printf,var_lun,0,format='(i1)'
  printf,var_lun,AOD,format='(f5.3)'
  printf,var_lun,-0.092,format='(f6.3)'
  printf,var_lun,-1000,format='(i5)'
  printf,var_lun,27,format='(i2)'
  printf,var_lun,0,format='(i1)'
  printf,var_lun,0,format='(i1)'
  printf,var_lun,0,format='(i1)'
  printf,var_lun,REF,format='(f3.1)'
  printf,var_lun,-1,format='(i2)'
  close,var_lun
  free_lun,var_lun
  SPAWN, '"E:\yaoganhomework\ch02_6s_test\ch02_6s_test\6SV-run.bat"' ,/hide
  filepath='E:\yaoganhomework\ch02_6s_test\ch02_6s_test\output.txt'
  openr,lun,filepath,/get_lun
  data_line=''
  flag=0
  i=0.0
  while flag eq 0 do begin
    readf,lun,data_line
    if strpos(data_line,'*       apparent reflectance') ne -1 then begin
      tdata=strsplit(data_line,'  ',/extract)
      i=float(tdata[3])
      flag=1
    endif
  endwhile
  free_lun, lun
  file_delete,'E:\yaoganhomework\ch02_6s_test\ch02_6s_test\output.txt'
  return,i
end