anchor.diff:        Train=sumf_areas(qr)
         Tcloud=sumf_areas(qc)
         Tvapour=sumf_areas(qv)
  -      
  +
         if(testcase==2)then
           h_error%f = h_exact%f - h%f
           u_error%f = u_exact%f - u%f
  
positive.code:     hqr_old=hqr
      htracer_old=htracer
   
      !Time loop
      do k=1, ntime
        !Calculate u and h for time:
  
negative.code:     
      ! High order advection scheme variable
      character (len=6):: advmtd
  
  contains
  
  