/[lmdze]/trunk/libf/dyn3d/Dissipation/divgrad.f
ViewVC logotype

Contents of /trunk/libf/dyn3d/Dissipation/divgrad.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 57 - (show annotations)
Mon Jan 30 12:54:02 2012 UTC (12 years, 5 months ago) by guez
File size: 1425 byte(s) 
Write used namelists to file "" instead of standard output.

Avoid aliasing in "inidissip" in calls to "divgrad2", "divgrad",
"gradiv2", "gradiv", "nxgraro2" and "nxgrarot". Add a degenerate
dimension to arrays so they have rank 3, like the dummy arguments in
"divgrad2", "divgrad", "gradiv2", "gradiv", "nxgraro2" and "nxgrarot".

Extract the initialization part from "bilan_dyn" and make a separate
procedure, "init_dynzon", from it.

Move variables from modules "iniprint" and "logic" to module
"conf_gcm_m".

Promote internal procedures of "fxy" to private procedures of module
"fxy_m".

Extracted documentation from "inigeom". Removed useless "save"
attributes. Removed useless intermediate variables. Extracted
processing of poles from loop on latitudes. Write coordinates to file
"longitude_latitude.txt" instead of standard output.

Do not use ozone tracer for radiative transfer.
1 !
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/divgrad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3 !
4 SUBROUTINE divgrad (klevel,h, lh, divgra, cdivh )
5 use dimens_m
6 use paramet_m
7 use conf_gcm_m
8 use comgeom
9 use filtreg_m, only: filtreg
10 IMPLICIT NONE
11 c
12 c=======================================================================
13 c
14 c Auteur : P. Le Van
15 c ----------
16 c
17 c lh
18 c calcul de (div( grad )) de h .....
19 c h et lh sont des arguments d'entree pour le s-prog
20 c divgra est un argument de sortie pour le s-prog
21 c
22 c=======================================================================
23 c
24 c declarations:
25 c -------------
26 c
27 c
28 INTEGER klevel
29 REAL, intent(in):: h( ip1jmp1,klevel )
30 real, intent(out):: divgra( ip1jmp1,klevel )
31 real, intent(in):: cdivh
32 c
33 REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
34
35 INTEGER l,ij,iter
36 integer, intent(in):: lh
37 c
38 c
39 c
40 divgra = h
41 c
42 DO 10 iter = 1,lh
43
44 CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1 )
45
46 CALL grad (klevel,divgra, ghx , ghy )
47 CALL diverg (klevel, ghx , ghy , divgra )
48
49 CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)
50
51 DO 5 l = 1,klevel
52 DO 4 ij = 1, ip1jmp1
53 divgra( ij,l ) = - cdivh * divgra( ij,l )
54 4 CONTINUE
55 5 CONTINUE
56 c
57 10 CONTINUE
58 RETURN
59 END
  ViewVC Help
Powered by ViewVC 1.1.21