/[lmdze]/trunk/dyn3d/Dissipation/gradiv2.f
ViewVC logotype

Contents of /trunk/dyn3d/Dissipation/gradiv2.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 54 - (show annotations)
Tue Dec 6 15:07:04 2011 UTC (12 years, 5 months ago) by guez
Original Path: trunk/libf/dyn3d/gradiv2.f90
File size: 1699 byte(s) 
Removed Numerical Recipes procedure "ran1". Replaced calls to "ran1"
in "inidissip" by calls to intrinsic procedures.

Split file "interface_surf.f90" into a file with a module containing
only variables, "interface_surf", and single-procedure files. Gathered
files into directory "Interface_surf".

Added argument "cdivu" to "gradiv" and "gradiv2", "cdivh" to
"divgrad2" and "divgrad", and "crot" to "nxgraro2" and
"nxgrarot". "dissip" now uses variables "cdivu", "cdivh" and "crot"
from module "inidissip_m", so it can pass them to "gradiv2",
etc. Thanks to this modification, we avoid a circular dependency
betwwen "inidissip.f90" and "gradiv2.f90", etc. The value -1. used by
"gradiv2", for instance, during computation of eigenvalues is not the
value "cdivu" computed by "inidissip".

Extracted procedure "start_inter_3d" from module "startdyn", to its
own module.

In "inidissip", unrolled loop on "ii". I find it clearer now.

Moved variables "matriceun", "matriceus", "matricevn", "matricevs",
"matrinvn" and "matrinvs" from module "parafilt" to module
"inifilr_m". Moved variables "jfiltnu", "jfiltnv", "jfiltsu",
"jfiltsv" from module "coefils" to module "inifilr_m".
1 module gradiv2_m
2
3 IMPLICIT NONE
4
5 contains
6
7 SUBROUTINE gradiv2(klevel, xcov, ycov, ld, gdx, gdy, cdivu)
8
9 ! From LMDZ4/libf/dyn3d/gradiv2.F, version 1.1.1.1 2004/05/19 12:53:07
10
11 ! P. Le Van
12 ! calcul de grad div du vecteur v
13 ! xcov et ycov etant les composantes covariantes de v
14 ! xcont, ycont et ld sont des arguments d'entree pour le sous-programme
15 ! gdx et gdy sont des arguments de sortie pour le sous-programme
16
17 use dimens_m
18 use paramet_m
19 use comgeom
20 use filtreg_m, only: filtreg
21
22 ! variables en arguments
23
24 INTEGER klevel
25 REAL xcov( ip1jmp1,klevel), ycov( ip1jm,klevel)
26 integer, intent(in):: ld
27 REAL gdx( ip1jmp1,klevel), gdy( ip1jm,klevel)
28 real, intent(in):: cdivu
29
30 ! variables locales
31
32 REAL div(ip1jmp1,llm)
33 REAL nugrads
34 INTEGER l,ij,iter
35
36 !--------------------------------------------------------------
37
38 CALL SCOPY( ip1jmp1 * klevel, xcov, 1, gdx, 1)
39 CALL SCOPY( ip1jm * klevel, ycov, 1, gdy, 1)
40
41 CALL divergf( klevel, gdx, gdy, div)
42
43 IF( ld.GT.1) THEN
44 CALL laplacien ( klevel, div, div)
45
46 ! Iteration de l'operateur laplacien_gam
47 DO iter = 1, ld -2
48 CALL laplacien_gam ( klevel,cuvscvgam1,cvuscugam1,unsair_gam1, &
49 unsapolnga1, unsapolsga1, div, div)
50 ENDDO
51 ENDIF
52
53 CALL filtreg( div, jjp1, klevel, 2, 1, .TRUE., 1)
54 CALL grad ( klevel, div, gdx, gdy)
55 nugrads = (-1.)**ld * cdivu
56
57 DO l = 1, klevel
58 DO ij = 1, ip1jmp1
59 gdx( ij,l) = gdx( ij,l) * nugrads
60 ENDDO
61 DO ij = 1, ip1jm
62 gdy( ij,l) = gdy( ij,l) * nugrads
63 ENDDO
64 ENDDO
65
66 END SUBROUTINE gradiv2
67
68 end module gradiv2_m
  ViewVC Help
Powered by ViewVC 1.1.21