source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.f90 @ 5105

Last change on this file since 5105 was 5105, checked in by abarral, 3 months ago

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines
  • 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.6 KB
Line 
1
2! $Header$
3
4SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
5  !
6  !  P. Le Van
7  !
8  !
9  !    *******************************************************************
10  !  ... calcule la (convergence horiz. * aire locale)du flux ayant pour
11  !  composantes xflu et yflu ,variables extensives .  ......
12  !    *******************************************************************
13  !  xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
14  !  convfl                est  un argument de sortie pour le s-pg .
15  !
16  ! njxflu  est le nombre de lignes de latitude de xflu,
17  ! ( = jjm ou jjp1 )
18  ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
19  !
20  IMPLICIT NONE
21  !
22  INCLUDE "dimensions.h"
23  INCLUDE "paramet.h"
24  REAL :: xflu,yflu,convfl,convpn,convps
25  INTEGER :: l,ij,nbniv
26  DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) , &
27        convfl( ip1jmp1,nbniv )
28  !
29  REAL :: SSUM
30  !
31  !
32  INCLUDE "comgeom.h"
33  !
34  DO l = 1,nbniv
35  !
36  DO ij = iip2, ip1jm - 1
37  convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   + &
38        yflu(ij +1,l ) - yflu( ij -iim,l )
39  END DO
40  !
41  !
42
43  ! ....  correction pour  convfl( 1,j,l)  ......
44  ! ....   convfl(1,j,l)= convfl(iip1,j,l) ...
45  !
46  !DIR$ IVDEP
47  DO ij = iip2,ip1jm,iip1
48  convfl( ij,l ) = convfl( ij + iim,l )
49  END DO
50  !
51  ! ......  calcul aux poles  .......
52  !
53  convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
54  convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
55  DO ij = 1,iip1
56  convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
57  convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
58  END DO
59  !
60  END DO
61  RETURN
62END SUBROUTINE convflu
Note: See TracBrowser for help on using the repository browser.