anchor.diff:  !
   
   USE sdata, ONLY: mtem, cden, ftem, tin, xyz, cflow, nyy, nzz, cf, ent, heatf, &
  -                 ystag, tfm, nt, rpos, rdel, rf, rg, rc, pi, zdel, dia, ystag
  +                 ystag, tfm, nt, rpos, rdel, rf, rg, rc, pi, zdel, dia, ystag, &
  +                 farea
   
   IMPLICIT NONE
   
  
positive.code: 
              mdens = cden(xyz(i,j,k)) * 1000.                                    ! Coolant density (kg/m3)
              cpline = heatf(xyz(i,j,k)) * pi * dia  &
                    + cf * xpline(xyz(i,j,k)) * 100.                              ! Coolant Linear power densisty (W/m)
              vol   = farea * zdel(k) * 0.01
              IF (k == 1) THEN                                                    ! Calculate coolant enthalpy
                  eps = 0.5 * mdens * vol
  
negative.code:     DO j = 1, nyy
          DO i = ystag(j)%smin, ystag(j)%smax
  
              cpline = cpline = heatf(xyz(i,j,k)) * pi * dia  &
                     + cf * xpline(xyz(i,j,k)) * 100. * farea / (pi *rf**2)       ! Coolant Linear power densisty (W/m)                             ! Coolant Linear power densisty (W/m)
  
              IF (k == 1) THEN                                                    ! Calculate coolant enthalpy and
                  ent(xyz(i,j,k)) = enti + cpline * zdel(k) * 0.01 / cflow        ! corresponding temp and density
  