anchor.diff:        alpha = 1.0D0/3.0D0
         beta = 1.0D0/15.0D0
   
  -      if(testcase == 8)then
  +      if(testcase == 5)then
           lat2 = -25.d0*deg2rad
         else
           lat2 = 0.5D0*piby2
  
positive.code:     !======================================================================================
      ! Galewski - Jet in Southern Hemisphere
      !======================================================================================
      case(7,8)
        !Parameters
        !mu1    = 0.05_r8
        mu1    = 0.00002_r8
  
negative.code:       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
          qv_error%f = qv_exact%f - qv%f
  