c-----------------------------------------------------------------------
c        Included code to initialize non-scalar field vars.
c        Simply calls 'calc_al'---note that RNPL wants 'nm1'
c        level initialized.
c-----------------------------------------------------------------------

         logical     ltrace
         parameter ( ltrace = .true. )

         if( ltrace ) then
            write(0,*) 'init_al: In routine with g1_nr = ', g1_nr
            write(0,*) 'init_al: t = ', t
            write(0,*) 'init_al: hcdlnatol = ', hcdlnatol
            write(0,*) 'init_al: hcmxiter  = ', hcmxiter
            write(0,*) 'init_al: flatspace  = ', flatspace
            write(0,*)
         end if

         call calc_al(
     &                pp_nm1,   pp_n,     pp_np1,
     &                pi_nm1,   pi_n,     pi_np1,
     &                a_nm1,    a_n,      a_np1, 
     &                l_nm1,    l_n,      l_np1, 
     &                tmr,
     &                f0,      f2,   g0,      
     &                g1_Nr,   r,    t,
     &                flatspace,
     &                hcdlnatol,     hcmxiter,         hcrc,
     &                tmrbh
     &               )