anchor.diff:    ep2baro(:,:) = 0.
     ep3baro(:,:) = 0.
     ep4(:,:) = 0.
  -  dc = dz/prefac
  +  dc = dz/6745.348
   
  -  do k = 2,kmax-1
  +  do k = 2,96
       zk = dz*float(k-1)
       do i = 1,imax
         do j = 6,nd-1            ! 5N and higher latitude
  
positive.code: 
          ! low-level meridional eddy heat flux
          if(k.eq.2) then     ! (26) of SI-HN17
            ep41 = 2.*om*sin0*cos0*dz/prefac       ! prefactor
            ep42 = exp(-dz/hh)*vv(i,j+90,2)*(pt(i,j+90,2)-tref(j-5,2))
            ep42 = ep42/(tg(3)-tg(1))
            ep43 = vv(i,j+90,1)*(pt(i,j+90,1)-tref(j-5,1))
            ep43 = 0.5*ep43/(tg(2)-tg(1))
  
negative.code:         ep1(i,j) = ep1(i,j)+0.5*vv(i,j+90,k)**2    !F3a+b
          ep11 = 0.5*(pt(i,j+90,k)-tref(j-5,k))**2         !F3c
          zz = dz*float(k-1)
          ep11 = ep11*(r/h)*exp(-rkappa*zz/h)
          ep11 = ep11*2.*dz/(tg(k+1)-tg(k-1))
          ep1(i,j) = ep1(i,j)-ep11                   !F3
          phip = dp*float(j)
  