anchor.diff:  
         !Calculate erngies
         call calc_energies(Penergy, Kenergy, Tenergy, Availenergy)
  -      print*,iniwater,Twater,(Twater-iniwater)/iniwater
  -
  +      print '(a33, 2e16.8)','Change in mass of h*(total water):', (Twater-iniwater)/iniwater
  +      print*,''
         !update fields
         u_old=u
         h_old=h
  
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:       Qv%f = q0*Qv%f
        hQv%f = q0*hQv%f
  
        !print*,q0,maxval(qv%f),minval(qv%f)
        !stop 
  
      case default
  