anchor.diff:      !Check for blow ups
       integer(i4)::blowup=0
   
  -    !Error variables - for tc2 and tc3
  +    !Error variables - for tc2
       real(r8):: rel_error_h
       real(r8):: rel_error_u
       real(r8):: rel_error_Qv
  
positive.code:         !lon = lon+1.d0*pi + 0.52777d0*pi !-95graus
          lon = lon+pi+pi + 30.d0*deg2rad !+ 0.52777d0*pi !-95graus
          theta%f(i)=F_quad(theta_sp,(1._r8-mu1)*theta_eq,theta_np,lat) + mu1*theta_eq*dcos(lat)*dsin(lon)
  
          !Vapour
          Qv%f(i) = mu2*qsat(theta%f(i),h%f(i),bt%f(i),1._r8)
          
  
negative.code:       call ode_rk4_moist_swm(time, h_old, u_old, htheta_old, hQv_old, hQc_old, hQr_old, &
                           h, u, htheta, hQv, hQc, hQr, dt)
  
        !Apply the monotonic filter for tracers
        call monotonic_filter(hQv)             
        call monotonic_filter(hQr)
        call monotonic_filter(hQc)
  
        !call mass_fixer(hQv,hQr,hQc, iniwater)
        !call scalar_elem_divide(htheta, h, theta)
  