anchor.diff:      lond = 0._r8*deg2rad
   
       select case(testcase)
  -    !-------------------------------------------------------------------------------  
  -    ! Steady state 
  -    !-------------------------------------------------------------------------------  
  -
  +    !======================================================================================
  +    ! Steady state - from Zerroukat and Allen JCP 2015
  +    !======================================================================================
       case(2)
         !Field at cell's center
         u0    = 20._r8
  
positive.code: 
      !Set a standart name for files
      write(atmp,'(i8)') int(testcase)
      swmname="swm_physics_tc"//trim(adjustl(trim(atmp)))
  
      !Print information
      print*, "Test Case     : ", testcase
  
negative.code:     !Variables at Voronoi centers
      do i=1, mesh%nv
          vectmp = mesh%v(i)%p
          lon = mesh%v(i)%lon
          lat = mesh%v(i)%lat
          sinn = sin(lat)
  