anchor.diff:            DIAG_var(n,k,l,3) = vx_local(k)
             DIAG_var(n,k,l,4) = vy_local(k)
             DIAG_var(n,k,l,5) = vz_local(k)
  -          DIAG_var(n,k,l,6) = 0.D0
          enddo
   
       enddo
  
positive.code:        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(:,:,:,:) )
         call tracer_init( ADM_gall, ADM_kall, ADM_lall, DIAG_var(:,:,:,:) )
      else
         write(ADM_LOG_FID,*) 'xxx Invalid input_io_mode. STOP.'
  
negative.code: 
    !-----------------------------------------------------------------------------
    subroutine tracer_init( &
       ijdim,   &
       kdim,    &
       lall,    &
       DIAG_var )
      use mod_misc, only: &
         MISC_get_latlon, &
         MISC_get_distance
      use mod_adm, only: &
         ADM_KNONE
      use mod_cnst, only: &
         CNST_PI, &
         CNST_ERADIUS
      use mod_grd, only: &
         GRD_vz,         &
         GRD_x,          &
         GRD_x_pl,       &
         GRD_XDIR,       &
         GRD_YDIR,       &
         GRD_ZDIR,       &
         GRD_Z,          &
         GRD_ZH
      use mod_gmtr, only: &
         GMTR_lon, &
         GMTR_lat
  