anchor.diff:          order=4
   
       case default
  -        print*, "Invalid advection method. Please select a proper method."
  +        print*, "Invalid advection method order. Please select a proper order."
           stop
       end select
   
  -    swmname=trim(swmname)//"_advorder"//trim(advmtd)
  +    if (method == 'O') then
  +        swmname=trim(swmname)//"_advmethod"//trim(method)//"_advorder"//trim(advmtd)
  +    else
  +        order = 3
  +        swmname=trim(swmname)//"_advmethod"//trim(method)
  +    endif
       print*, "SWM Name for Plots: ", trim(swmname)
       return
     end subroutine swm_phys_pars
  
positive.code:       if(testcase==2 .or. testcase==3)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
          Qc_error%f = Qc_exact%f - Qc%f
          Qr_error%f = Qr_exact%f - Qr%f
          theta_error%f = theta_exact%f - theta%f
          rel_error_h = maxval(abs(h_error%f))/maxval(abs(h_exact%f))
          rel_error_u = maxval(abs(u_error%f))/maxval(abs(u_exact%f))
          rel_error_Qv = maxval(abs(Qv_error%f))/maxval(abs(Qv_exact%f))
          rel_error_Qc = maxval(abs(Qc_error%f))
          rel_error_Qr = maxval(abs(Qr_error%f))
          rel_error_theta = maxval(abs(theta_error%f))/maxval(abs(theta_exact%f))
          !print*, k, ntime
          print*, "Time (dys) :",   k*dt*sec2day, " of ", ntime*dt*sec2day
          print*, "Step = ", k, " of ", ntime
          print '(a33, 3e16.8)','linf errors of (h, u, theta) = ',rel_error_h,rel_error_u,rel_error_theta
          print '(a33, 3e16.8)','linf errors of  (qv, qc, qr) = ',rel_error_Qv,rel_error_Qc,rel_error_Qr
          !print '(a22, 2e16.8)',' height = ',minval(h%f),maxval(h%f)
          !print '(a22, 2e16.8)',' velocity = ',minval(u%f),maxval(u%f)
          !print '(a22, 2e16.8)',' temperature = ',minval(theta%f),maxval(theta%f)
  
negative.code:     !Error variables - for tc2
      real(r8):: rel_error_h
      real(r8):: rel_error_u
      real(r8):: rel_error_qv
      real(r8):: rel_error_qc
      real(r8):: rel_error_qr
      real(r8):: rel_error_theta
  
      !Total mass of cloud, rain and vapour
  