/[lmdze]/trunk/Sources/dyn3d/Dissipation/divgrad2.f
ViewVC logotype

Contents of /trunk/Sources/dyn3d/Dissipation/divgrad2.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/divgrad2.f
File size: 2138 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 !
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad2.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3 !
4 SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra, cdivh )
5 c
6 c P. Le Van
7 c
8 c ***************************************************************
9 c
10 c ..... calcul de (div( grad )) de ( pext * h ) .....
11 c ****************************************************************
12 c h ,klevel,lh et pext sont des arguments d'entree pour le s-prg
13 c divgra est un argument de sortie pour le s-prg
14 c
15 use dimens_m
16 use paramet_m
17 use comgeom
18 IMPLICIT NONE
19 c
20
21 c ....... variables en arguments .......
22 c
23 INTEGER klevel
24 REAL h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel )
25 REAL divgra( ip1jmp1,klevel)
26 real, intent(in):: cdivh
27 c
28 c ....... variables locales ..........
29 c
30 REAL signe, nudivgrs, sqrtps( ip1jmp1,llm )
31 INTEGER l,ij,iter
32 integer, intent(in):: lh
33 c ...................................................................
34
35 c
36 signe = (-1.)**lh
37 nudivgrs = signe * cdivh
38
39 CALL SCOPY ( ip1jmp1 * klevel, h, 1, divgra, 1 )
40
41 c
42 CALL laplacien( klevel, divgra, divgra )
43
44 DO l = 1, klevel
45 DO ij = 1, ip1jmp1
46 sqrtps( ij,l ) = SQRT( deltapres(ij,l) )
47 ENDDO
48 ENDDO
49 c
50 DO l = 1, klevel
51 DO ij = 1, ip1jmp1
52 divgra(ij,l) = divgra(ij,l) * sqrtps(ij,l)
53 ENDDO
54 ENDDO
55
56 c ........ Iteration de l'operateur laplacien_gam ........
57 c
58 DO iter = 1, lh - 2
59 CALL laplacien_gam ( klevel,cuvscvgam2,cvuscugam2,unsair_gam2,
60 * unsapolnga2, unsapolsga2, divgra, divgra )
61 ENDDO
62 c
63 c ...............................................................
64
65 DO l = 1, klevel
66 DO ij = 1, ip1jmp1
67 divgra(ij,l) = divgra(ij,l) * sqrtps(ij,l)
68 ENDDO
69 ENDDO
70 c
71 CALL laplacien ( klevel, divgra, divgra )
72 c
73 DO l = 1,klevel
74 DO ij = 1,ip1jmp1
75 divgra(ij,l) = nudivgrs * divgra(ij,l) / deltapres(ij,l)
76 ENDDO
77 ENDDO
78
79 RETURN
80 END
  ViewVC Help
Powered by ViewVC 1.1.21