source: LMDZ5/trunk/libf/dyn3d/traceurpole.F @ 1941

Last change on this file since 1941 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
Line 
1!
2! $Id: traceurpole.F 1907 2013-11-26 13:10:46Z jghattas $
3!
4          subroutine traceurpole(q,masse)
5
6      USE control_mod
7
8          implicit none
9     
10#include "dimensions.h"
11c#include "paramr2.h"
12#include "paramet.h"
13#include "comconst.h"
14#include "comdissip.h"
15#include "comvert.h"
16#include "comgeom2.h"
17#include "logic.h"
18#include "temps.h"
19#include "ener.h"
20#include "description.h"
21
22
23c   Arguments
24       integer iq
25       real masse(iip1,jjp1,llm)
26       real q(iip1,jjp1,llm)
27       
28
29c   Locals
30      integer i,j,l
31      real sommemassen(llm)
32      real sommemqn(llm)
33      real sommemasses(llm)
34      real sommemqs(llm)
35      real qpolen(llm),qpoles(llm)
36
37   
38c On impose une seule valeur au pôle Sud j=jjm+1=jjp1       
39      sommemasses=0
40      sommemqs=0
41          do l=1,llm
42             do i=1,iip1         
43                 sommemasses(l)=sommemasses(l)+masse(i,jjp1,l)
44                 sommemqs(l)=sommemqs(l)+masse(i,jjp1,l)*q(i,jjp1,l)
45             enddo         
46          qpoles(l)=sommemqs(l)/sommemasses(l)
47          enddo
48
49c On impose une seule valeur du traceur au pôle Nord j=1
50      sommemassen=0
51      sommemqn=0 
52         do l=1,llm
53           do i=1,iip1             
54               sommemassen(l)=sommemassen(l)+masse(i,1,l)
55               sommemqn(l)=sommemqn(l)+masse(i,1,l)*q(i,1,l)
56           enddo
57           qpolen(l)=sommemqn(l)/sommemassen(l)
58         enddo
59   
60c On force le traceur à prendre cette valeur aux pôles
61        do l=1,llm
62            do i=1,iip1
63               q(i,1,l)=qpolen(l)
64               q(i,jjp1,l)=qpoles(l)
65             enddo
66        enddo
67
68     
69      return
70      end           
Note: See TracBrowser for help on using the repository browser.