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., 1)

END SUBROUTINE dteta1
1	guez	44	SUBROUTINE dteta1(teta, pbaru, pbarv, dteta)
2	guez	27
3	guez	44	! 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	guez	3
6	guez	44	! Calcul du terme de convergence horizontale du flux d'enthalpie
7			! potentielle.
8	guez	3
9	guez	44	! dteta est un argument de sortie pour le s-pg
10	guez	3
11	guez	44	use dimens_m
12			use paramet_m
13	guez	57	use conf_gcm_m
14	guez	44	use filtreg_m, only: filtreg
15	guez	3
16	guez	44	IMPLICIT NONE
17	guez	3
18	guez	44	REAL, intent(in):: teta(ip1jmp1, llm), pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
19			REAL dteta(ip1jmp1, llm)
20			INTEGER l, ij
21	guez	3
22	guez	44	REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm)
23	guez	3
24	guez	44	!----------------------------------------------------------------
25	guez	3
26	guez	44	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	guez	3
31	guez	44	DO ij = iip1+ iip1, ip1jm, iip1
32			hbxu(ij, l) = hbxu(ij - iim, l)
33			end DO
34	guez	3
35	guez	44	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	guez	3
40	guez	44	CALL convflu(hbxu, hbyv, llm, dteta)
41	guez	3
42	guez	44	! stockage dans dh de la convergence horizont. filtree' du flux
43			! d'enthalpie potentielle
44			CALL filtreg(dteta, jjp1, llm, 2, 2, .true., 1)
45	guez	3
46	guez	44	END SUBROUTINE dteta1