source: LMDZ4/branches/LMDZ4_V2_patch/libf/dyn3d/traceurpole.F @ 1798

Last change on this file since 1798 was 524, checked in by lmdzadmin, 20 years ago

Initial revision

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