source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.f90 @ 5185

Last change on this file since 5185 was 5159, checked in by abarral, 5 months ago

Put dimensions.h and paramet.h into modules
  • 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.0 KB
RevLine 
[5099]1
[524]2! $Header$
[5099]3
[5106]4SUBROUTINE nxgradst(klevel,rot, x, y )
[5136]5
6  USE lmdz_comgeom
7
[5159]8  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
9  USE lmdz_paramet
[5105]10  IMPLICIT NONE
11  ! Auteur :  P. Le Van
[5159]12
[5105]13  !   ********************************************************************
14  !  calcul du gradient tourne de pi/2 du rotationnel du vect.v
15  !   ********************************************************************
16  !   rot          est un argument  d'entree pour le s-prog
17  !   x  et y    sont des arguments de sortie pour le s-prog
18  !
[524]19
[5159]20
21
[5105]22  INTEGER :: klevel
23  REAL :: rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
24  INTEGER :: l,ij
[5159]25
[5105]26  DO l = 1,klevel
[5159]27
[5105]28  DO ij = 2, ip1jm
29  y(ij,l)=( rot(ij,l) - rot(ij-1,l))
30  END DO
[5159]31
[5105]32  !    ..... correction pour  y ( 1,j,l )  ......
[5159]33
[5105]34  !    ....    y(1,j,l)= y(iip1,j,l) ....
[524]35
[5105]36  DO ij = 1, ip1jm, iip1
37  y( ij,l ) = y( ij +iim,l )
38  END DO
[5159]39
[5105]40  DO ij = iip2,ip1jm
41  x(ij,l)= rot(ij,l)-rot(ij-iip1,l)
42  END DO
43  DO ij = 1,iip1
44  x(    ij    ,l ) = 0.
45  x( ij +ip1jm,l ) = 0.
46  END DO
[5159]47
[5105]48  END DO
49  RETURN
50END SUBROUTINE nxgradst
Note: See TracBrowser for help on using the repository browser.