anchor.diff:      real (r8):: utmp
       real (r8):: vtmp
       real (r8):: u0
  -    
  +
       u0 = 2._r8*pi*erad/(12._r8*day2sec)
       call cart2sph(p(1), p(2), p(3), lon, lat)
       utmp = u0*cos(lat)
  
positive.code:         l1 = long - 0.93d0*pi
  
          clat = COS(l2)
          if(testcase ==8)then
             l1 = l1 + 120.d0*deg2rad
          end if
          e1 = EXP(-(l1/alpha)**2)
  
negative.code:         print '(a22, 3e16.8)',' rain = ',minval(qr%f),maxval(qr%f), Train
        end if
  
        print '(a22, 2e16.8)',' rain = ',minval(htracer%f),maxval(htracer%f)
        if (rel_error_h>1000.d0) stop
        !print*,'CFL = ',cfl
        !Plot fields
        call plotfields_mswm(k, time)
  