source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrarot.f90 @ 5221

! $Header$

SUBROUTINE nxgrarot(klevel, xcov, ycov, lr, grx, gry)
  !   ***********************************************************

  !    Auteur :  P.Le Van

  !                             lr
  !  calcul de  ( nXgrad (rot) )   du vect. v  ....

  !   xcov et ycov  etant les compos. covariantes de  v
  !   ***********************************************************
  ! xcov , ycov et lr  sont des arguments  d'entree pour le s-prog
  !  grx   et  gry     sont des arguments de sortie pour le s-prog


  USE lmdz_filtreg, ONLY: filtreg
  USE lmdz_ssum_scopy, ONLY: scopy
  USE lmdz_comdissipn, ONLY: tetaudiv, tetaurot, tetah, cdivu, crot, cdivh

USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
  USE lmdz_paramet
  IMPLICIT NONE

  !



  INTEGER :: klevel
  REAL :: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)
  REAL :: grx(ip1jmp1, klevel), gry(ip1jm, klevel)

  REAL :: rot(ip1jm, llm)

  INTEGER :: l, ij, iter, lr



  CALL SCOPY (ip1jmp1 * klevel, xcov, 1, grx, 1)
  CALL SCOPY (ip1jm * klevel, ycov, 1, gry, 1)

  DO iter = 1, lr
    CALL  rotat (klevel, grx, gry, rot)
    CALL filtreg(rot, jjm, klevel, 2, 1, .FALSE., 2)
    CALL nxgrad (klevel, rot, grx, gry)

    DO l = 1, klevel
      DO ij = 1, ip1jm
        gry(ij, l) = - gry(ij, l) * crot
      END DO
      DO ij = 1, ip1jmp1
        grx(ij, l) = - grx(ij, l) * crot
      END DO
    END DO

  END DO
  RETURN
END SUBROUTINE nxgrarot