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

Annotation of /trunk/dyn3d/Dissipation/dissip.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 57 - (hide annotations)
Mon Jan 30 12:54:02 2012 UTC (12 years, 4 months ago) by guez
Original Path: trunk/libf/dyn3d/Dissipation/dissip.f90
File size: 3008 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 guez 47 module dissip_m
2 guez 3
3 guez 47 IMPLICIT NONE
4 guez 3
5 guez 47 contains
6 guez 3
7 guez 47 SUBROUTINE dissip(vcov, ucov, teta, p, dv, du, dh)
8 guez 3
9 guez 47 ! From dyn3d/dissip.F, version 1.1.1.1 2004/05/19 12:53:05
10     ! Author: P. Le Van
11 guez 56 ! Objet : calcul de la dissipation horizontale
12     ! Avec opĂ©rateurs star : gradiv2, divgrad2, nxgraro2
13 guez 3
14 guez 55 USE dimens_m, ONLY: iim, jjm, llm
15     USE comdissnew, ONLY: lstardis, nitergdiv, nitergrot, niterh
16     USE inidissip_m, ONLY: dtdiss, tetah, tetaudiv, tetaurot, cdivu, crot, cdivh
17 guez 54 use gradiv2_m, only: gradiv2
18 guez 56 use nr_util, only: assert
19 guez 3
20 guez 55 REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)
21     REAL, intent(in):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)
22 guez 56 REAL, intent(in):: teta(:, :, :) ! (iim + 1, jjm + 1, llm)
23     REAL, INTENT(IN):: p(:, :, :) ! (iim + 1, jjm + 1, llm + 1)
24 guez 55 REAL, intent(out):: dv(:, :, :) ! (iim + 1, jjm, llm)
25     REAL, intent(out):: du(:, :, :) ! (iim + 1, jjm + 1, llm)
26     REAL, intent(out):: dh(:, :, :) ! (iim + 1, jjm + 1, llm)
27 guez 3
28 guez 47 ! Local:
29 guez 55 REAL gdx(iim + 1, jjm + 1, llm), gdy(iim + 1, jjm, llm)
30     REAL grx(iim + 1, jjm + 1, llm), gry(iim + 1, jjm, llm)
31 guez 56 REAL tedt(llm)
32     REAL deltapres(iim + 1, jjm + 1, llm)
33 guez 55 INTEGER l
34 guez 3
35 guez 47 !-----------------------------------------------------------------------
36 guez 3
37 guez 56 call assert((/size(vcov, 1), size(ucov, 1), size(teta, 1), size(p, 1), &
38     size(dv, 1), size(du, 1), size(dh, 1)/) == iim + 1, "dissip iim")
39     call assert((/size(vcov, 2), size(ucov, 2) - 1, size(teta, 2) - 1, &
40     size(p, 2) - 1, size(dv, 2), size(du, 2) - 1, size(dh, 2) - 1/) &
41     == jjm, "dissip jjm")
42     call assert((/size(vcov, 3), size(ucov, 3), size(teta, 3), size(p, 3) - 1, &
43     size(dv, 3), size(du, 3), size(dh, 3)/) == llm, "dissip llm")
44 guez 3
45 guez 56 du(:, 1, :) = 0.
46     du(:, jjm + 1, :) = 0.
47 guez 3
48 guez 47 ! Calcul de la partie grad (div) :
49 guez 3
50 guez 47 IF (lstardis) THEN
51 guez 57 CALL gradiv2(ucov, vcov, nitergdiv, gdx, gdy, cdivu)
52 guez 47 ELSE
53 guez 54 CALL gradiv(llm, ucov, vcov, nitergdiv, gdx, gdy, cdivu)
54 guez 47 END IF
55 guez 3
56 guez 56 tedt = tetaudiv * dtdiss
57 guez 55 forall (l = 1: llm)
58 guez 56 du(:, 2: jjm, l) = - tedt(l) * gdx(:, 2: jjm, l)
59     dv(:, :, l) = - tedt(l) * gdy(:, :, l)
60 guez 55 END forall
61 guez 3
62 guez 56 ! Calcul de la partie n X grad (rot) :
63 guez 3
64 guez 47 IF (lstardis) THEN
65 guez 54 CALL nxgraro2(llm, ucov, vcov, nitergrot, grx, gry, crot)
66 guez 47 ELSE
67 guez 54 CALL nxgrarot(llm, ucov, vcov, nitergrot, grx, gry, crot)
68 guez 47 END IF
69 guez 3
70 guez 56 tedt = tetaurot * dtdiss
71 guez 55 forall (l = 1: llm)
72 guez 56 du(:, 2: jjm, l) = du(:, 2: jjm, l) - tedt(l) * grx(:, 2: jjm, l)
73     dv(:, :, l) = dv(:, :, l) - tedt(l) * gry(:, :, l)
74 guez 55 END forall
75 guez 3
76 guez 47 ! calcul de la partie div (grad) :
77 guez 3
78 guez 47 IF (lstardis) THEN
79 guez 56 forall (l = 1: llm) &
80     deltapres(:, :, l) = max(0., p(:, :, l) - p(:, :, l + 1))
81 guez 54 CALL divgrad2(llm, teta, deltapres, niterh, gdx, cdivh)
82 guez 47 ELSE
83 guez 54 CALL divgrad(llm, teta, niterh, gdx, cdivh)
84 guez 47 END IF
85 guez 3
86 guez 56 forall (l = 1: llm) dh(:, :, l) = - tetah(l) * dtdiss * gdx(:, :, l)
87 guez 47
88     END SUBROUTINE dissip
89    
90     end module dissip_m
  ViewVC Help
Powered by ViewVC 1.1.21