source: LMDZ5/trunk/libf/dyn3dmem/laplacien_rot_loc.F

      SUBROUTINE laplacien_rot_loc ( klevel, rotin, rotout,ghx,ghy )
c
c    P. Le Van
c
c   ************************************************************
c    ...  calcul de  ( rotat x nxgrad )  du rotationnel rotin  .
c   ************************************************************
c
c     klevel et rotin  sont des arguments  d'entree pour le s-prog
c      rotout           est  un argument  de sortie pour le s-prog
c
      USE parallel_lmdz
      USE mod_filtreg_p
      IMPLICIT NONE
c
#include "dimensions.h"
#include "paramet.h"
#include "comgeom.h"

c 
c   ..........    variables  en  arguments     .............
c
      INTEGER klevel
      REAL rotin(ijb_v:ije_v,klevel ), rotout( ijb_v:ije_v,klevel )
c
c   ..........    variables   locales       ................
c
      REAL ghy(ijb_v:ije_v,klevel), ghx(ijb_u:ije_u,klevel)
c   ........................................................
c
c
      INTEGER :: ijb,ije,jjb,jje
      
      jjb=jj_begin-1
      jje=jj_end+1
      
      if (pole_nord) jjb=jj_begin
      if (pole_sud) jje=jj_end-1
      
      CALL  filtreg_p ( rotin ,jjb_v,jje_v,jjb,jje,jjm,
     &                  klevel,2, 1, .FALSE., 1)

      CALL   nxgrad_loc ( klevel, rotin,   ghx ,  ghy            )
      CALL   rotatf_loc  ( klevel, ghx  ,   ghy , rotout         )
c
      RETURN
      END