anchor.diff:  
     integer, private, save      :: I_RHOGQstr =  7 ! tracers
     integer, private, save      :: I_RHOGQend = -1 !
  +
  +  character(len=16), private, save :: PRG_name(PRG_vmax0)
  +  data PRG_name / 'rhog', 'rhogvx', 'rhogvy', 'rhogvz', 'rhogw', 'rhoge' /
   
     integer, private, parameter :: DIAG_vmax0 = 6
   
  
positive.code:     ! <- [add]&[Mod] H.Yashiro 20110819
  
      write(ADM_LOG_FID,*)
      write(ADM_LOG_FID,*) '====== data range check : prognostic variables ======'
      do nq = 1, DIAG_vmax0
         val_max = GTL_max( DIAG_var   (:,:,:,nq),       &
                            DIAG_var_pl(:,:,:,nq),       &
  
negative.code:                                 DIAG_var(:,:,:,I_qstr:I_qend),        & !--- [INOUT]
                                  PRG_var (:,:,:,I_RHOGQstr:I_RHOGQend) ) !--- [INOUT]
      endif
  
      write(ADM_LOG_FID,*)
      write(ADM_LOG_FID,*) '====== data range check : prognostic variables ======'
      do nq = 1, DIAG_vmax0
         val_max = GTL_max( PRG_var   (:,:,:,nq),        &
                            PRG_var_pl(:,:,:,nq),        &
                            ADM_kall, ADM_kmin, ADM_kmax )
         val_min = GTL_min( PRG_var   (:,:,:,nq),        &
                            PRG_var_pl(:,:,:,nq),        &
                            ADM_kall, ADM_kmin, ADM_kmax )
  
         write(ADM_LOG_FID,'(1x,A,A16,2(A,1PE24.17))') '--- ', PRG_name(nq), ': max=', val_max, ', min=', val_min
      enddo
  
      do nq = 1, TRC_vmax
         val_max = GTL_max( PRG_var   (:,:,:,PRG_vmax0+nq), &
                            PRG_var_pl(:,:,:,PRG_vmax0+nq), &
                            ADM_kall, ADM_kmin, ADM_kmax    )
  
         if ( val_max <= 0.D0 ) then
            nonzero = .false.
         else
            nonzero = .true.
         endif
  
         val_min = GTL_min( PRG_var   (:,:,:,PRG_vmax0+nq),      &
                            PRG_var_pl(:,:,:,PRG_vmax0+nq),      &
                            ADM_kall, ADM_kmin, ADM_kmax, nonzero)
  
         write(ADM_LOG_FID,'(1x,A,A16,2(A,1PE24.17))') '--- rhog * ', TRC_name(nq),  ': max=', val_max, ', min=', val_min
      enddo
  
      return
    end subroutine restart_input
  