source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.F @ 5099

Last change on this file since 5099 was 5099, checked in by abarral, 8 weeks ago

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line
  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.1 KB
Line 
1
2! $Header$
3
4      SUBROUTINE nxgradst (klevel,rot, x, y )
5c
6      IMPLICIT NONE
7c     Auteur :  P. Le Van
8c
9c   ********************************************************************
10c      calcul du gradient tourne de pi/2 du rotationnel du vect.v
11c   ********************************************************************
12c       rot          est un argument  d'entree pour le s-prog
13c       x  et y    sont des arguments de sortie pour le s-prog
14c
15      INCLUDE "dimensions.h"
16      INCLUDE "paramet.h"
17      INCLUDE "comgeom.h"
18
19      INTEGER klevel
20      REAL rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
21      INTEGER l,ij
22c
23      DO l = 1,klevel
24c
25      DO ij = 2, ip1jm
26      y(ij,l)=( rot(ij,l) - rot(ij-1,l))
27      END DO
28c
29c    ..... correction pour  y ( 1,j,l )  ......
30c
31c    ....    y(1,j,l)= y(iip1,j,l) ....
32
33      DO ij = 1, ip1jm, iip1
34      y( ij,l ) = y( ij +iim,l )
35      END DO
36c
37      DO ij = iip2,ip1jm
38      x(ij,l)= rot(ij,l)-rot(ij-iip1,l)
39      END DO
40      DO ij = 1,iip1
41      x(    ij    ,l ) = 0.
42      x( ij +ip1jm,l ) = 0.
43      END DO
44c
45      END DO
46      RETURN
47      END
Note: See TracBrowser for help on using the repository browser.