;
; simrad calls tomrad to simulate radiances for a finite
; bandwidth using solar flux weighting across a triangular
; slit function with fwhm in angstrom to the nearest half A,
; and band center wavelength in angstrom.
;
cartrn=''
cnvrt=0.017453
;
;  define circumstances of the simulated measurement
;
sza      = [7.8]
satza    = [0.0]           ;leave spacecraft zenith angle 0.0 for nadir view
phi      = [0.0]           ;leave relative azimuth angle 0.0 for nadir view
refl      = [15.0]         ;%
;
trprs= 1.0
toz = 225.0
getprf21,'l225',prfoz,prftemp
;
;  sbuv/2 wavelengths
;
fwhm=11.32
bcw=[2556.5,2736.1,2831.0,2877.0,2923.0,2976.0,3020.0,3059.0, $
     3126.0,3176.0,3313.0,3399.0,3430.0]
;
radtot=fltarr(13)
nvalue=fltarr(13)
;
;  calculate simulated radiances for each wavelength of interest
;
  for i=0,12 do begin
    simrad,prfoz,prftemp,trprs,sza,satza,phi,refl,bcw(i),fwhm,rad
    radtot(i)=rad
    nvalue(i)=-100.0*alog10(rad)
  endfor
;
  print,bcw
  print,radtot
  print,nvalue
;
;  write results to a simulated pmf record
;
openw,3,'simpmf.dat',/f77_unformatted
irec=0
;
;  make some stuff up
;
yyyy=2002
doy=365
cldprs=0.4                   ;atm
lat=45.0
lon=0.0
;
;  write output simulated pmf data record
;
buf=fltarr(207)
buf(1)=irec+2
buf(3)=yyyy*1000.0+doy
for i=10,13 do buf(i)=nvalue(12)
for i=8,11 do buf(25-i)=nvalue(i)
buf(21)=cldprs
buf(43)=trprs
buf(47)=lat
buf(48)=lon
buf(49)=sza
for i=50,57 do buf(i)=nvalue(12)
for i=58,65 do buf(i)=nvalue(i-58)
buf(205)=sza*cnvrt*10000.0
buf(206)=sza*cnvrt*10000.0
writeu,3,buf
;
close,3
;
end