anchor.diff:                        / ( CNST_RAIR * ( t2m(:,K0,l) + 0.5D0 * 0.005D0 * ( GRD_zs(:,K0,l,GRD_ZSFC) + 2.D0 ) ) ) )
   
             call history_in( 'sl_slp', slp(:,:,l) )
  +       enddo
  +    endif
  +
  +    if (out_850hPa) then   ! [add] 20130705 R.Yoshida
  +       do l = 1, ADM_lall
  +          call sv_plev_uvwt( ADM_gall,              & ! [IN]
  +                             pre    (:,:,l),        & ! [IN]
  +                             u      (:,:,l),        & ! [IN]
  +                             v      (:,:,l),        & ! [IN]
  +                             w      (:,:,l),        & ! [IN]
  +                             tem    (:,:,l),        & ! [IN]
  +                             85000.d0,              & ! [IN]
  +                             u_850  (:,K0,l),       & ! [OUT]
  +                             v_850  (:,K0,l),       & ! [OUT]
  +                             w_850  (:,K0,l),       & ! [OUT]
  +                             t_850  (:,K0,l)        ) ! [OUT]
  +
  +          call history_in( 'sl_u850', u_850(:,:,l) )
  +          call history_in( 'sl_v850', v_850(:,:,l) )
  +          call history_in( 'sl_w850', w_850(:,:,l) )
  +          call history_in( 'sl_t850', t_850(:,:,l) )
          enddo
       endif
   
  
positive.code:     return
    end subroutine sv_pre_sfc
  
  end module mod_history_vars
  !-------------------------------------------------------------------------------
  
negative.code:     real(8) :: vz10m(ADM_gall,ADM_KNONE,ADM_lall)
      real(8) :: t2m(ADM_gall,ADM_KNONE,ADM_lall)
      real(8) :: q2m(ADM_gall,ADM_KNONE,ADM_lall)
      real(8) :: u_850(ADM_gall,ADM_KNONE,ADM_lall)   ! [add] 20130705 R.Yoshida
      real(8) :: v_850(ADM_gall,ADM_KNONE,ADM_lall)
      real(8) :: w_850(ADM_gall,ADM_KNONE,ADM_lall)
      real(8) :: t_850(ADM_gall,ADM_KNONE,ADM_lall)
  
      real(8) :: albedo_sfc(ADM_gall,ADM_KNONE,ADM_lall,1:NRDIR,1:NRBND)
      real(8) :: slp(ADM_gall,ADM_KNONE,ADM_lall)
  