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

Annotation of /trunk/Sources/dyn3d/Dissipation/divgrad2.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 56 - (hide annotations)
Tue Jan 10 19:02:02 2012 UTC (12 years, 4 months ago) by guez
Original Path: trunk/libf/dyn3d/Dissipation/divgrad2.f
File size: 2127 byte(s) 
Imported "writehist.f" from LMDZ.

Moved module variable "histaveid" from "com_io_dyn" to "initdynav_m".

In "inithist", access directly module variables from "com_io_dyn"
instead of going through the arguments. Copying from LMDZ, write "u"
and scalar variables to separate files. Create a new variable for the
new file in "com_io_dyn". Copying from LMDZ, change the vertical axes
of the three files.

Removed some useless initializations in "dissip".

In "bilan_dyn", removed useless variable "time". Avoiding the
approximate test on "dt_cum" being a multiple of "dt_app", just
compute "ncum" from known usage of "bilan_dyn" and compute "dt_cum"
from "ncum". Change "periodav" from real to integer in
"conf_gcm_m". Since "day_step" is required to be a multiple of
"iperiod", so is "ncum".
1 guez 3 !
2     ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad2.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3     !
4 guez 54 SUBROUTINE divgrad2 ( klevel, h, deltapres, lh, divgra, cdivh )
5 guez 3 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 guez 55 REAL, intent(in):: h( ip1jmp1,klevel ), deltapres( ip1jmp1,klevel)
25 guez 56 REAL, intent(out):: divgra( ip1jmp1,klevel)
26 guez 54 real, intent(in):: cdivh
27 guez 3 c
28     c ....... variables locales ..........
29     c
30     REAL signe, nudivgrs, sqrtps( ip1jmp1,llm )
31 guez 26 INTEGER l,ij,iter
32     integer, intent(in):: lh
33 guez 3 c ...................................................................
34    
35     c
36     signe = (-1.)**lh
37     nudivgrs = signe * cdivh
38 guez 55 divgra = h
39 guez 3
40     c
41     CALL laplacien( klevel, divgra, divgra )
42    
43     DO l = 1, klevel
44     DO ij = 1, ip1jmp1
45     sqrtps( ij,l ) = SQRT( deltapres(ij,l) )
46     ENDDO
47     ENDDO
48     c
49     DO l = 1, klevel
50     DO ij = 1, ip1jmp1
51     divgra(ij,l) = divgra(ij,l) * sqrtps(ij,l)
52     ENDDO
53     ENDDO
54    
55     c ........ Iteration de l'operateur laplacien_gam ........
56     c
57     DO iter = 1, lh - 2
58     CALL laplacien_gam ( klevel,cuvscvgam2,cvuscugam2,unsair_gam2,
59     * unsapolnga2, unsapolsga2, divgra, divgra )
60     ENDDO
61     c
62     c ...............................................................
63    
64     DO l = 1, klevel
65     DO ij = 1, ip1jmp1
66     divgra(ij,l) = divgra(ij,l) * sqrtps(ij,l)
67     ENDDO
68     ENDDO
69     c
70     CALL laplacien ( klevel, divgra, divgra )
71     c
72     DO l = 1,klevel
73     DO ij = 1,ip1jmp1
74     divgra(ij,l) = nudivgrs * divgra(ij,l) / deltapres(ij,l)
75     ENDDO
76     ENDDO
77    
78     RETURN
79     END
  ViewVC Help
Powered by ViewVC 1.1.21