source: LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/inter_barxy.F @ 1074

Last change on this file since 1074 was 774, checked in by Laurent Fairhead, 17 years ago

Suite du merge entre la version et la HEAD: quelques modifications de
Yann sur le

LF

  • 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!
4       SUBROUTINE inter_barxy ( interfd,jnterfd,dlonid,dlatid ,
5     ,        champ,imod,jmod,rlonimod,rlatimod, jsort,champint )
6
7c    Auteur :   P. Le Van
8c
9       INTEGER interfd,jnterfd,imod,jmod
10       REAL champ(interfd,jnterfd +1 ),dlonid(interfd),dlatid(jnterfd),
11     ,      champint(imod,jsort)
12       REAL rlonimod(imod),rlatimod(jmod)
13
14#include "dimensions.h"
15#include "paramet.h"
16#include "comgeom2.h"
17
18       REAL champx(imod),champy(jnterfd +1,imod),chpn(imod),chps(imod)
19       REAL chhpn,chhps
20       REAL fmody(jjp1)
21c
22
23       DO j = 1, jnterfd + 1
24        CALL inter_barx( interfd, dlonid, champ( 1,j ),
25     ,                       imod, rlonimod , champx )
26         DO i = 1,imod
27           champy(j,i) = champx(i)
28         ENDDO
29       ENDDO
30
31       DO i = 1, imod
32        CALL inter_bary( jjm,jnterfd,dlatid,champy(1,i),
33     ,                     jmod ,rlatimod,  fmody     )
34          DO j = 1, jsort
35           champint(i,j) = fmody(j)
36          ENDDO
37       ENDDO
38
39       IF( jsort.EQ.jjp1)  THEN
40
41c   ....  Valeurs uniques  aux  poles ....
42c
43         DO i =  1,imod
44          chpn(i)  = aire( i,  1   ) * champint( i, 1   )
45          chps(i)  = aire( i, jjp1 ) * champint( i,jjp1 )
46         ENDDO
47          chhpn  = SSUM(imod,chpn,1)/apoln
48          chhps  = SSUM(imod,chps,1)/apols
49
50         DO i = 1, imod
51          champint( i,  1  ) = chhpn
52          champint( i, jjp1) = chhps
53         ENDDO
54c
55       ENDIF
56
57       RETURN
58       END
59
Note: See TracBrowser for help on using the repository browser.