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:                sparams(iii)%p2=p2
                 call assign_source_params_type(sparams(iii),stype)
              enddo
              NCAMS=nparams*nfreq*nmu*nt
              allocate(c(NCAMS))
              if(outfile.ne."") write(6,*) 'outfile grtrans: ',outfile, NCAMS
              do m=1,NCAMS
  
negative.code: 
         real(8), allocatable, dimension(:,:,:) :: ivals!,ab
         real(8), allocatable, dimension(:,:) :: ab
         real(8), allocatable, dimension(:) :: freqs!,mdots,mu0
  
         contains
  
  