source: trunk/LMDZ.TITAN/Tools/moyzon2.F @ 826

Last change on this file since 826 was 816, checked in by slebonnois, 12 years ago

SL: tools for postprocessing (Veznus and Titan); see DOC/documentation/vt-tools.pdf

File size: 1.2 KB
RevLine 
[816]1      SUBROUTINE moyzon2(iim,jjp1,nl,indefini,x,xbar)
2c=======================================================================
3c
4c
5c   Subject:
6c   ------
7c   Calcul de la moyenne zonale de la variable au point scalaire x
8c
9c MEME CHOSE QUE moyzon, SAUF LONGUEUR DIM LON DE X
10c
11c=======================================================================
12      IMPLICIT NONE
13c-----------------------------------------------------------------------
14c   Declararations:
15c   ---------------
16
17c   Arguments:
18c   ----------
19
20      INTEGER iim,jjp1,nl
21      real indefini
22      REAL x(iim,jjp1,nl)
23      REAL xbar(jjp1,nl)
24
25c   Local:
26c   ------
27
28      iNTEGER i,j,l , n
29
30c------------------------------------------------------------------------
31      do l=1,nl
32        do j=1,jjp1
33           xbar(j,l)=0.
34           n = 0
35           do i=1,iim
36              if (x(i,j,l).lt.indefini) then
37                 xbar(j,l) = xbar(j,l) + x(i,j,l)
38                 n = n+1
39              end if
40           end do
41           if (n.ne.0) then
42              xbar(j,l) = xbar(j,l)/float(n)
43           else
44              xbar(j,l) = indefini
45           end if 
46        end do
47      end do
48
49      RETURN
50      END
Note: See TracBrowser for help on using the repository browser.