source: LMDZ5/trunk/libf/dyn3dmem/nxgrad.F @ 1934

Last change on this file since 1934 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.2 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE nxgrad (klevel, rot, x, y )
5c
6c     P. Le Van
7c
8c   ********************************************************************
9c      calcul du gradient tourne de pi/2 du rotationnel du vect.v
10c   ********************************************************************
11c       rot          est un argument  d'entree pour le s-prog
12c       x  et y    sont des arguments de sortie pour le s-prog
13c
14      IMPLICIT NONE
15c
16#include "dimensions.h"
17#include "paramet.h"
18#include "comgeom.h"
19      INTEGER klevel
20      REAL rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
21      INTEGER   l,ij
22c
23c
24      DO 10 l = 1,klevel
25c
26      DO 1  ij = 2, ip1jm
27      y( ij,l ) = (  rot( ij,l ) - rot( ij-1,l )  ) * cvsurcuv( ij )
28   1  CONTINUE
29c
30c    ..... correction pour  y ( 1,j,l )  ......
31c
32c    ....    y(1,j,l)= y(iip1,j,l) ....
33CDIR$ IVDEP
34      DO 2  ij = 1, ip1jm, iip1
35      y( ij,l ) = y( ij +iim,l )
36   2  CONTINUE
37c
38      DO 4  ij = iip2,ip1jm
39      x( ij,l ) = (  rot( ij,l ) - rot( ij -iip1,l )  ) * cusurcvu( ij )
40   4  CONTINUE
41      DO 6 ij = 1,iip1
42      x(    ij    ,l ) = 0.
43      x( ij +ip1jm,l ) = 0.
44   6  CONTINUE
45c
46  10  CONTINUE
47      RETURN
48      END
Note: See TracBrowser for help on using the repository browser.