c ... a quick hack to test the UW MAS Planck routines

c ... Tree of subprogram calls:
c ... MAIN
c ...    INITCF
c ...        GETCOF
c ...            BNDFIT
c ...                LINREG
c ...                WBRITE
c ...                WPLANC
c ...            CENWAV
c ...            HMWSRF
c ...                GETLUN
c ...                INTERP
c ...    MASRAD
c ...        PLANCK

      implicit none
      
      integer band, i
      real rad( 6 ), tem( 6 ), masrad
      external masrad

      integer init
      real wc( 50 ), tcs( 50 ), tci( 50 )
      common / mascof / init, wc, tcs, tci

      data tem / 200.0, 220.0, 240.0, 260.0, 280.0, 300.0 /
      
      call initcf()

      write(*,*) 'MAS IR BAND RADIANCES (WATTS/SQ.METER/SR/MICRON)'
      write(*,*) 'BAND  CENTER    SLOPE   INTRCEPT  ' //
     &  'T=200K   T=220K   T=240K   T=260K   T=280K   T=300K'

      do band = 26, 50
        do i = 1, 6
          rad( i ) = masrad( band, tem( i ) )
        end do
        write(*,'(2x,i2,2x,f9.4,8f9.6)')
     &    band, wc( band ), tcs( band ), tci( band ), rad
      end do

      end