/[lmdze]/trunk/dyn3d/tourabs.f
ViewVC logotype

Annotation of /trunk/dyn3d/tourabs.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 27 - (hide annotations)
Thu Mar 25 14:29:07 2010 UTC (14 years, 2 months ago) by guez
Original Path: trunk/libf/dyn3d/tourabs.f
File size: 2490 byte(s) 
"dyn3d" and "filtrez" do not contain any included file so make rules
have been updated.

"comdissip.f90" was useless, removed it.

"dynredem0" wrote undefined value in "controle(31)", that was
overwritten by "dynredem1". Now "dynredem0" just writes 0 to
"controle(31)".

Removed arguments of "inidissip". "inidissip" now accesses the
variables by use association.

In program "etat0_lim", "itaufin" is not defined so "dynredem1" wrote
undefined value to "controle(31)". Added argument "itau" of
"dynredem1" to correct that.

"itaufin" does not need to be a module variable (of "temps"), made it
a local variable of "leapfrog".

Removed calls to "diagedyn" from "leapfrog".
1 guez 3 SUBROUTINE tourabs ( ntetaSTD,vcov, ucov, vorabs )
2    
3     c=======================================================================
4     c
5     c Modif: I. Musat (28/10/04)
6     c -------
7     c adaptation du code tourpot.F pour le calcul de la vorticite absolue
8     c cf. P. Le Van
9     c
10     c Objet:
11     c ------
12     c
13     c *******************************************************************
14     c ............. calcul de la vorticite absolue .................
15     c *******************************************************************
16     c
17     c ntetaSTD, vcov,ucov sont des argum. d'entree pour le s-pg .
18     c vorabs est un argum.de sortie pour le s-pg .
19     c
20     c=======================================================================
21    
22     use dimens_m
23     use paramet_m
24     use comconst
25     use logic
26     use comgeom
27 guez 27 use filtreg_m, only: filtreg
28    
29 guez 3 IMPLICIT NONE
30     c
31     INTEGER ntetaSTD
32     REAL vcov( ip1jm,ntetaSTD ), ucov( ip1jmp1,ntetaSTD )
33     REAL vorabs( ip1jm,ntetaSTD )
34     c
35     c variables locales
36     INTEGER l, ij, i, j
37     REAL rot( ip1jm,ntetaSTD )
38    
39    
40    
41     c ... vorabs = ( Filtre( d(vcov)/dx - d(ucov)/dy ) + fext ) ..
42    
43    
44    
45     c ........ Calcul du rotationnel du vent V puis filtrage ........
46    
47     DO 5 l = 1,ntetaSTD
48    
49     DO 2 i = 1, iip1
50     DO 2 j = 1, jjm
51     c
52     ij=i+(j-1)*iip1
53     IF(ij.LE.ip1jm - 1) THEN
54     c
55     IF(cv(ij).EQ.0..OR.cv(ij+1).EQ.0..OR.
56     $ cu(ij).EQ.0..OR.cu(ij+iip1).EQ.0.) THEN
57     rot( ij,l ) = 0.
58     continue
59     ELSE
60     rot( ij,l ) = (vcov(ij+1,l)/cv(ij+1)-vcov(ij,l)/cv(ij))/
61     $ (2.*pi*RAD*cos(rlatv(j)))*float(iim)
62     $ +(ucov(ij+iip1,l)/cu(ij+iip1)-ucov(ij,l)/cu(ij))/
63     $ (pi*RAD)*(float(jjm)-1.)
64     c
65     ENDIF
66     ENDIF !(ij.LE.ip1jm - 1) THEN
67     2 CONTINUE
68    
69     c .... correction pour rot( iip1,j,l ) .....
70     c .... rot(iip1,j,l) = rot(1,j,l) .....
71    
72     CDIR$ IVDEP
73    
74     DO 3 ij = iip1, ip1jm, iip1
75     rot( ij,l ) = rot( ij -iim, l )
76     3 CONTINUE
77    
78     5 CONTINUE
79    
80    
81     CALL filtreg( rot, jjm, ntetaSTD, 2, 1, .FALSE., 1 )
82    
83    
84     DO 10 l = 1, ntetaSTD
85    
86     DO 6 ij = 1, ip1jm - 1
87     vorabs( ij,l ) = ( rot(ij,l) + fext(ij)*unsairez(ij) )
88     6 CONTINUE
89    
90     c ..... correction pour vorabs( iip1,j,l) .....
91     c .... vorabs(iip1,j,l)= vorabs(1,j,l) ....
92     CDIR$ IVDEP
93     DO 8 ij = iip1, ip1jm, iip1
94     vorabs( ij,l ) = vorabs( ij -iim,l )
95     8 CONTINUE
96    
97     10 CONTINUE
98    
99     RETURN
100     END
  ViewVC Help
Powered by ViewVC 1.1.21