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

Annotation of /trunk/dyn3d/dteta1.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 104 - (hide annotations)
Thu Sep 4 10:05:52 2014 UTC (9 years, 8 months ago) by guez
File size: 1333 byte(s) 
Removed procedure sortvarc0. Called sortvarc with an additional
argument resetvarc instead. (Following LMDZ.) Moved current time
computations and some printing statements from sortvarc to
caldyn. Could then remove arguments itau and time_0 of sortvarc, and
could remove "use dynetat0". Better to keep "dynetat0.f" as a gcm-only
file.

Moved some variables from module ener to module sortvarc.

Split file "mathelp.f" into single-procedure files.

Removed unused argument nadv of adaptdt. Removed dimension arguments
of bernoui.

Removed unused argument nisurf of interfoce_lim. Changed the size of
argument lmt_sst of interfoce_lim from klon to knon. Removed case when
newlmt is false.

dynredem1 is called only once in each run, either ce0l or gcm. So
variable nb in call to nf95_put_var was always 1. Removed variable nb.

Removed dimension arguments of calcul_fluxs. Removed unused arguments
precip_rain, precip_snow, snow of calcul_fluxs. Changed the size of
all the arrays in calcul_fluxs from klon to knon.

Removed dimension arguments of fonte_neige. Changed the size of all
the arrays in fonte_neige from klon to knon.

Changed the size of arguments tsurf and tsurf_new of interfsurf_hq
from klon to knon. Changed the size of argument ptsrf of soil from
klon to knon.
1 guez 88 module dteta1_m
2 guez 27
3 guez 88 IMPLICIT NONE
4 guez 3
5 guez 88 contains
6 guez 3
7 guez 88 SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)
8 guez 3
9 guez 88 ! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1, 2004/05/19 12:53:06
10     ! Authors: P. Le Van, F. Forget
11 guez 3
12 guez 88 ! Calcul du terme de convergence horizontale du flux d'enthalpie
13     ! potentielle.
14 guez 3
15 guez 104 USE dimens_m, ONLY: iim, jjm, llm
16 guez 88 USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1, jjp1
17     USE filtreg_m, ONLY: filtreg
18 guez 3
19 guez 88 REAL, intent(in):: teta(ip1jmp1, llm)
20     REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
21 guez 104 REAL, intent(out):: dteta(iim + 1, jjm + 1, llm)
22 guez 3
23 guez 88 ! Local:
24     INTEGER l, ij
25     REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)
26 guez 3
27 guez 88 !----------------------------------------------------------------
28 guez 3
29 guez 88 DO l = 1, llm
30     DO ij = iip2, ip1jm - 1
31     hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
32     end DO
33 guez 3
34 guez 88 DO ij = iip1+ iip1, ip1jm, iip1
35     hbxu(ij, l) = hbxu(ij - iim, l)
36     end DO
37 guez 3
38 guez 88 DO ij = 1, ip1jm
39     hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
40     end DO
41     end DO
42 guez 3
43 guez 88 CALL convflu(hbxu, hbyv, llm, dteta)
44 guez 3
45 guez 88 ! stockage dans dh de la convergence horizontale filtrée du flux
46     ! d'enthalpie potentielle
47     CALL filtreg(dteta, jjp1, llm, 2, 2, .true.)
48    
49     END SUBROUTINE dteta1
50    
51     end module dteta1_m
  ViewVC Help
Powered by ViewVC 1.1.21