anchor.diff:  
   
       !Standard parameters file
  -    filename=trim(pardir)//"swm_physics.par"
  +    filename=trim(pardir)//"moist_swm.par"
       print*,"Shallow Water Model parameters (file): ", trim(filename)
       print*
       call getunit(fileunit)
  
positive.code:         h%f(i) = (phi0 -w*sinn**2)*gravi
  
          r  = dsqrt((lon-lon0)**2+(lat-lat0)**2)
          !r1 = dsqrt((lon-lon1)**2+(lat-lat0)**2)
  
          if(r<rmax)then
            bt%f(i)=2000._r8*(1._r8-r/rmax)
          !else if(r1<rmax)then
          !  bt%f(i)=2000._r8*(1._r8-r1/rmax)
          else
            bt%f(i)=0.
          endif
  
negative.code:       Qc_exact = Qc
        Qr_exact = Qr
  
      !======================================================================================
      ! Flow over a mountain - from Zerroukat and Allen JCP 2015
      !======================================================================================
      case(4)
  
      !Parameters
      u0     = 20._r8
  