write(unitid,'(f12.4)') prof_info(nf_time)   ! +/- around analysis time (00Z, 06Z, 12Z, or 18Z)
   write(unitid,'(2f12.4)') prof_info(nf_lat), prof_info(nf_lon)  !degree
   write(unitid,'(f12.4)') prof_info(nf_saza)  !degree
   write(unitid,'(f12.4)') prof_info(nf_soza)  !degree
   write(unitid,'(f12.4)') prof_2d(nf_ps)*0.01 !hPa
   write(unitid,'(f12.4)') prof_2d(nf_ts)      !K

   do i= 1, nlayers     !nlayers=72
      write(unitid,'(i4,f12.2,f12.2,f12.6)') i, Atmosphere(1)%Pressure(i), &
                    Atmosphere(1)%Temperature(i), Atmosphere(1)%Absorber(i,1)    !hPa, K, g/kg
   enddo

  !n_Channels = 2960 (see gxs_channel_wavenum_info.txt)
  write(unitid,'(3000f12.4)') (btchan(i),i=1,ChannelInfo(1)%n_Channels)    !K
  write(unitid,'(3000f12.4)') (anoise(i),i=1,ChannelInfo(1)%n_Channels)    !K
  write(unitid,'(3000f8.3)') (btemissivity(i),i=1,ChannelInfo(1)%n_Channels)