anchor.diff:         call hs_init ( ADM_gall, ADM_kall, ADM_lall, DIAG_var(:,:,:,:) )
       elseif( trim(init_type) == "Jablonowski" ) then
          call jbw_init( ADM_gall, ADM_kall, ADM_lall, DIAG_var(:,:,:,:) )
  -    elseif( trim(init_type) == "Tracer" ) then ! tentative test
  -       call jbw_init   ( ADM_gall, ADM_kall, ADM_lall, DIAG_var(:,:,:,:) )
  +    elseif( trim(init_type) == "Traceradvection" ) then ! tentative test
          call tracer_init( ADM_gall, ADM_kall, ADM_lall, DIAG_var(:,:,:,:) )
       else
          write(ADM_LOG_FID,*) 'xxx Invalid input_io_mode. STOP.'
  
positive.code:     K0 = ADM_KNONE
      e = exp(1.D0)
  
      !DIAG_var(:,:,:,1) = 1000.D2 ! pressure
      !DIAG_var(:,:,:,2) =  300.D0 ! temperature
      DIAG_var(:,:,:,3) =    0.D0 ! v(x)
      DIAG_var(:,:,:,4) =    0.D0 ! v(y)
      DIAG_var(:,:,:,5) =    0.D0 ! v(z)
      DIAG_var(:,:,:,6) =    0.D0 ! w
  
      do l = 1, lall
      do n = 1, ijdim
  
negative.code:     integer, intent(in)    :: ijdim
      integer, intent(in)    :: kdim
      integer, intent(in)    :: lall
      real(8), intent(inout) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(8) :: lon_center   = 120.D0 ! [degree]
      real(8) :: lat_center   =  30.D0 ! [degree]
  