source: LMDZ6/branches/Amaury_posttrusting/libf/dyn3d/dudv1.f90 @ 5333

Last change on this file since 5333 was 5312, checked in by abarral, 3 weeks ago

.f90 <-> .F90

  • 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.4 KB
Line 
1!
2! $Header$
3!
4SUBROUTINE dudv1 ( vorpot, pbaru, pbarv, du, dv )
5  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
6USE paramet_mod_h
7IMPLICIT NONE
8  !
9  !-----------------------------------------------------------------------
10  !
11  !   Auteur:   P. Le Van
12  !   -------
13  !
14  !   Objet:
15  !   ------
16  !   calcul du terme de  rotation
17  !   ce terme est ajoute a  d(ucov)/dt et a d(vcov)/dt  ..
18  !   vorpot, pbaru et pbarv sont des arguments d'entree  pour le s-pg ..
19  !   du  et dv              sont des arguments de sortie pour le s-pg ..
20  !
21  !-----------------------------------------------------------------------
22
23
24
25
26  REAL :: vorpot( ip1jm,llm ) ,pbaru( ip1jmp1,llm ) , &
27        pbarv( ip1jm,llm ) ,du( ip1jmp1,llm ) ,dv( ip1jm,llm )
28  INTEGER :: l,ij
29  !
30  !
31  DO l = 1,llm
32  !
33  DO  ij = iip2, ip1jm - 1
34  du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) * &
35        (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) + &
36        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
37  END DO
38  !
39  DO ij = 1, ip1jm - 1
40  dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) * &
41        (   pbaru(ij, l)  +  pbaru(ij+1   , l) + &
42        pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
43  END DO
44  !
45  !    .... correction  pour  dv( 1,j,l )  .....
46  !    ....   dv(1,j,l)= dv(iip1,j,l) ....
47  !
48  !DIR$ IVDEP
49  DO ij = 1, ip1jm, iip1
50  dv( ij,l ) = dv( ij + iim, l )
51  END DO
52  !
53  END DO
54  RETURN
55END SUBROUTINE dudv1
Note: See TracBrowser for help on using the repository browser.