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

Contents of /trunk/dyn3d/laplacien.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 57 - (show annotations)
Mon Jan 30 12:54:02 2012 UTC (12 years, 4 months ago) by guez
Original Path: trunk/libf/dyn3d/laplacien.f
File size: 1247 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/laplacien.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3 !
4 SUBROUTINE laplacien ( klevel, teta, divgra )
5 c
6 c P. Le Van
7 c
8 c ************************************************************
9 c .... calcul de (div( grad )) de teta .....
10 c ************************************************************
11 c klevel et teta sont des arguments d'entree pour le s-prog
12 c divgra est un argument de sortie pour le s-prog
13 c
14 use dimens_m
15 use paramet_m
16 use comgeom
17 use filtreg_m, only: filtreg
18 use divergf_m, only: divergf
19
20 IMPLICIT NONE
21 c
22
23 c
24 c ......... variables en arguments ..............
25 c
26 INTEGER, intent(in):: klevel
27 REAL teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
28 c
29 c ............ variables locales ..............
30 c
31 REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
32 c .......................................................
33
34
35 c
36 CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 )
37
38 CALL filtreg( divgra, jjp1, klevel, 2, 1, .TRUE., 1 )
39 CALL grad ( klevel,divgra, ghx , ghy )
40 CALL divergf ( klevel, ghx , ghy , divgra )
41
42 RETURN
43 END
  ViewVC Help
Powered by ViewVC 1.1.21