trunk/dyn3d/dteta1.f

SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)

  ! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1 2004/05/19 12:53:06
  ! Auteurs : P. Le Van, F. Forget

  ! Calcul du terme de convergence horizontale du flux d'enthalpie
  ! potentielle.

  ! dteta est un argument de sortie pour le s-pg

  use dimens_m
  use paramet_m
  use conf_gcm_m
  use filtreg_m, only: filtreg

  IMPLICIT NONE

  REAL, intent(in):: teta(ip1jmp1, llm), pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
  REAL dteta(ip1jmp1, llm)
  INTEGER l, ij

  REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)

  !----------------------------------------------------------------

  DO l = 1, llm
     DO ij = iip2, ip1jm - 1
        hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
     end DO

     DO ij = iip1+ iip1, ip1jm, iip1
        hbxu(ij, l) = hbxu(ij - iim, l)
     end DO

     DO ij = 1, ip1jm
        hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
     end DO
  end DO

  CALL convflu(hbxu, hbyv, llm, dteta)

  ! stockage dans dh de la convergence horizont. filtree' du flux
  ! d'enthalpie potentielle
  CALL filtreg(dteta, jjp1, llm, 2, 2, .true.)

END SUBROUTINE dteta1
1	SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)
2
3	! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1 2004/05/19 12:53:06
4	! Auteurs : P. Le Van, F. Forget
5
6	! Calcul du terme de convergence horizontale du flux d'enthalpie
7	! potentielle.
8
9	! dteta est un argument de sortie pour le s-pg
10
11	use dimens_m
12	use paramet_m
13	use conf_gcm_m
14	use filtreg_m, only: filtreg
15
16	IMPLICIT NONE
17
18	REAL, intent(in):: teta(ip1jmp1, llm), pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
19	REAL dteta(ip1jmp1, llm)
20	INTEGER l, ij
21
22	REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)
23
24	!----------------------------------------------------------------
25
26	DO l = 1, llm
27	DO ij = iip2, ip1jm - 1
28	hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l))
29	end DO
30
31	DO ij = iip1+ iip1, ip1jm, iip1
32	hbxu(ij, l) = hbxu(ij - iim, l)
33	end DO
34
35	DO ij = 1, ip1jm
36	hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l))
37	end DO
38	end DO
39
40	CALL convflu(hbxu, hbyv, llm, dteta)
41
42	! stockage dans dh de la convergence horizont. filtree' du flux
43	! d'enthalpie potentielle
44	CALL filtreg(dteta, jjp1, llm, 2, 2, .true.)
45
46	END SUBROUTINE dteta1