Data Format
character(len=30) :: institution ! your institution
integer :: wri_unit ! output file number
integer :: i ! loop index
real :: reso ! spectral resolution:
! 1.0 cm**-1 or 0.01 cm**-1
integer :: atmosphere ! 1 = midlat day
! 2 = midlat night
! 3 = polar winter
! 4 = polar summer
character(len=4) :: nlte_flag ! ' LTE' or 'NLTE'
real :: mwl ! left boundary of microwindow
integer :: n_pts ! number of radiances
real :: ztan ! tangent altitude in km
real :: wdel ! wavenumber increment (0.01 or 1.00 cm**-1)
real, dimension (:) :: spec ! spectral radiance nW/(cm**2 sr cm**-1)
write(wri_unit,'(a2,a30)') '# ', institution
write(wri_unit,'(a2,i2,2x,a4,f5.0,f7.2,f7.2)')&
'# ', atmosphere, nlte_flag, z_tan, mwl,reso
do i = 1, n_pts ! on output wavenumber points
write(wri_unit,'(f18.10,1pe16.7)') mwl + (i - 1) * wdel, spec(i)
end do