anchor.diff:  !     lower - o2, p = 706.2720 mbar, t = 278.94 k
   !     upper - o2, p = 4.758820 mbarm t = 250.85 k
   
  -! Compute the optical depth by interpolating in ln(pressure) and 
  +! Compute the optical depth by interpolating in ln(pressure) and
   ! temperature.  Below laytrop, the water vapor self-continuum and
   ! foreign continuum is interpolated (in temperature) separately.
   
  
positive.code:    REAL, DIMENSION( kts:kte+2 ),INTENT(IN   )    ::      Plev
  
  ! LOCAL VAR
    
     INTEGER :: k
  
  !                                                                                
  !  COMPUTE OZONE MIXING RATIO DISTRIBUTION                                       
  !                                                                                
     DO K=kts,kte+1
        O3PROF(K)=0.                                                       
     ENDDO
                                                                                   
     CALL O3DATA(O3PROF, Plev, kts, kte)
  
     END SUBROUTINE INIRAD
                                                                                   
  !-------------------------------------------------------------------------
     SUBROUTINE O3DATA (O3PROF, Plev, kts, kte)
  !-------------------------------------------------------------------------
  
negative.code:          enddo
  
  !  For band 16, if radiative transfer will be performed on just
  !  this band, use integrated Planck values up to 3250 cm-1.
  !  If radiative transfer will be performed across all 16 bands,
  !  then include in the integrated Planck values for this band
  !  contributions from 2600 cm-1 to infinity.
  