anchor.diff:              return
             end subroutine grtrans_main
   
  +!          subroutine init_pgrtrans_data(npix,)
  +!            integer, intent(in) :: npix,nmdots,nmu,nfreqs
  +!            allocate(ivals(n)); allocate(ab(n)); allocate(freqs(nfreqs))
  +!            allocate(mdots(nmdots)); allocate(mu0(nmu))
  +!          end subroutine init_pgrtrans_data
  +
             subroutine del_pgrtrans_data()
               deallocate(ivals); deallocate(ab); deallocate(freqs)
  +            deallocate(mdots); deallocate(mu0)
             end subroutine del_pgrtrans_data
             
           end module pgrtrans
  
positive.code: 
         real(8), allocatable, dimension(:,:,:) :: ivals!,ab
         real(8), allocatable, dimension(:,:) :: ab
         real(8), allocatable, dimension(:) :: freqs
  
         contains
  
  
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
  