PRO calmodel

RsperDN=[130.,65.,75.,65.,101.,115.,105.,90.,72.,55.]
temperature=2225.
attenuation=4.6E-11
rdcwlrpb,cwl,rpb
planck,10.*cwl,[temperature,attenuation],bb
trans=[0.844,0.88,0.61,0.795,0.61,0.79,0.768,0.796,0.70,0.86]
bias=[62.,62.,62.,62.,62.,135.,135.,135.,135.,135.]
rate=[0.463,0.463,0.463,0.463,0.463,0.409,0.409,0.409,0.409,0.409]
maxsignal=2L^14-1
maxexp=(REPLICATE(maxsignal,10)-bias)/(rate+1000.*(bb*trans*10.*rpb)/RsperDN)
PRINT,maxexp,FORMAT='(5F7.2)'
STOP
RETURN
END