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

Contents of /trunk/dyn3d/addfi.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 91 - (show annotations)
Wed Mar 26 17:18:58 2014 UTC (10 years, 1 month ago) by guez
File size: 2616 byte(s) 
Removed unused variables lock_startdate and time_stamp of module
calendar.

Noticed that physiq does not change the surface pressure. So removed
arguments ps and dpfi of subroutine addfi. dpfi was always 0. The
computation of ps in addfi included some averaging at the poles. In
principle, this does not change ps but in practice it does because of
finite numerical precision. So the results of the simulation are
changed. Removed arguments ps and dpfi of calfis. Removed argument
d_ps of physiq.

du at the poles is not computed by dudv1, so declare only the
corresponding latitudes in dudv1. caldyn passes only a section of the
array dudyn as argument.

Removed variable niadv of module iniadvtrac_m.

Declared arguments of exner_hyb as assumed-shape arrays and made all
other horizontal sizes in exner_hyb dynamic. This allows the external
program test_disvert to use exner_hyb at a single horizontal position.
1 module addfi_m
2
3 IMPLICIT NONE
4
5 contains
6
7 SUBROUTINE addfi(ucov, vcov, teta, q, dufi, dvfi, dtetafi, dqfi)
8
9 ! From dyn3d/addfi.F, v 1.1.1.1 2004/05/19 12:53:06
10
11 ! Addition of the physical tendencies
12
13 USE comconst, ONLY: dtphys
14 USE comgeom, ONLY: aire, apoln, apols
15 USE dimens_m, ONLY: iim, jjm, llm, nqmx
16
17 ! First and second components of the covariant velocity:
18 REAL, intent(inout):: ucov((iim + 1) * (jjm + 1), llm)
19 REAL, intent(inout):: vcov((iim + 1) * jjm, llm)
20
21 REAL, intent(inout):: teta((iim + 1) * (jjm + 1), llm)
22 ! potential temperature
23
24 real, intent(inout):: q((iim + 1) * (jjm + 1), llm, nqmx)
25
26 ! Tendencies:
27 REAL, intent(in):: dufi((iim + 1) * (jjm + 1), llm)
28 REAL, intent(in):: dvfi((iim + 1) * jjm, llm)
29 real, intent(in):: dtetafi((iim + 1) * (jjm + 1), llm)
30 REAL, intent(in):: dqfi((iim + 1) * (jjm + 1), llm, nqmx)
31
32 ! Local variables :
33 REAL xpn(iim), xps(iim), tpn, tps
34 INTEGER j, k, iq, ij
35 REAL, PARAMETER:: qtestw = 1e-15, qtestt = 1e-40
36
37 !-----------------------------------------------------------------------
38
39 teta = teta + dtetafi * dtphys
40
41 DO k = 1, llm
42 DO ij = 1, iim
43 xpn(ij) = aire(ij) * teta(ij , k)
44 xps(ij) = aire(ij+(iim + 1) * jjm) * teta(ij+(iim + 1) * jjm, k)
45 ENDDO
46 tpn = SUM(xpn)/ apoln
47 tps = SUM(xps)/ apols
48
49 DO ij = 1, iim + 1
50 teta(ij , k) = tpn
51 teta(ij+(iim + 1) * jjm, k) = tps
52 ENDDO
53 ENDDO
54
55 DO k = 1, llm
56 DO j = iim + 2, (iim + 1) * jjm
57 ucov(j, k)= ucov(j, k) + dufi(j, k) * dtphys
58 ENDDO
59 ENDDO
60
61 vcov = vcov + dvfi * dtphys
62
63 DO iq = 1, 2
64 DO k = 1, llm
65 DO j = 1, (iim + 1) * (jjm + 1)
66 q(j, k, iq)= q(j, k, iq) + dqfi(j, k, iq) * dtphys
67 q(j, k, iq)= MAX(q(j, k, iq), qtestw)
68 ENDDO
69 ENDDO
70 ENDDO
71
72 DO iq = 3, nqmx
73 DO k = 1, llm
74 DO j = 1, (iim + 1) * (jjm + 1)
75 q(j, k, iq)= q(j, k, iq) + dqfi(j, k, iq) * dtphys
76 q(j, k, iq)= MAX(q(j, k, iq), qtestt)
77 ENDDO
78 ENDDO
79 ENDDO
80
81 DO iq = 1, nqmx
82 DO k = 1, llm
83 DO ij = 1, iim
84 xpn(ij) = aire(ij) * q(ij , k, iq)
85 xps(ij) = aire(ij+(iim + 1) * jjm) * q(ij+(iim + 1) * jjm, k, iq)
86 ENDDO
87 tpn = SUM(xpn)/apoln
88 tps = SUM(xps)/apols
89
90 DO ij = 1, iim + 1
91 q(ij , k, iq) = tpn
92 q(ij+(iim + 1) * jjm, k, iq) = tps
93 ENDDO
94 ENDDO
95 ENDDO
96
97 END SUBROUTINE addfi
98
99 end module addfi_m
  ViewVC Help
Powered by ViewVC 1.1.21