anchor.diff:  
          real(8), allocatable, dimension(:,:,:) :: ivals!,ab
          real(8), allocatable, dimension(:,:) :: ab
  -       real(8), allocatable, dimension(:) :: freqs
  +       real(8), allocatable, dimension(:) :: freqs!,mdots,mu0
   
          contains
   
  
positive.code:             return
            end subroutine grtrans_main
  
            subroutine del_pgrtrans_data()
              deallocate(ivals); deallocate(ab); deallocate(freqs)
            end subroutine del_pgrtrans_data
            
          end module pgrtrans
  
negative.code:        !          write(6,*) 'read inputs nup: ',spin
              ! Compute frequencies w/ log spacing between and including fmin, fmax:
              allocate(freqs(nfreq)); allocate(mu0(nmu)); allocate(mdots(nmdot))
              NCAMS=nmdot*nfreq*nmu*nt
  !            call init_pgrtrans_data(nro*nphi,NCAMS,nfreq,nmdot,nmu)
              if(nfreq==1) then
                 freqs=(/fmin/)
              else
  