doxy/html/laplacien__rot__p_8_f_source.html

      SUBROUTINE laplacien_rot_p ( 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

      IMPLICIT NONE

c

#include "dimensions.h"

#include "paramet.h"

#include "comgeom.h"


c

c   ..........    variables  en  arguments     .............

c

      INTEGER klevel

      REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )

c

c   ..........    variables   locales       ................

c

      REAL ghy(ip1jm,klevel), ghx(ip1jmp1,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,jje,jjm, klevel,2, 1, .false., 1)


      CALL   nxgrad_p( klevel, rotin,   ghx ,  ghy               )

      CALL   rotatf_p( klevel, ghx  ,   ghy , rotout             )

c

      RETURN

      END