LMDZ
dteta1.F
Go to the documentation of this file.
1 !
2 ! $Header$
3 !
4  SUBROUTINE dteta1 ( teta, pbaru, pbarv, dteta)
5  IMPLICIT NONE
6 
7 c=======================================================================
8 c
9 c Auteur: P. Le Van
10 c -------
11 c Modif F.Forget 03/94 (on retire q et dq pour construire dteta1)
12 c
13 c ********************************************************************
14 c ... calcul du terme de convergence horizontale du flux d'enthalpie
15 c potentielle ......
16 c ********************************************************************
17 c .. teta,pbaru et pbarv sont des arguments d'entree pour le s-pg ....
18 c dteta sont des arguments de sortie pour le s-pg ....
19 c
20 c=======================================================================
21 
22 
23 #include "dimensions.h"
24 #include "paramet.h"
25 #include "logic.h"
26 
27  REAL teta( ip1jmp1,llm ),pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm)
28  REAL dteta( ip1jmp1,llm )
29  INTEGER l,ij
30 
31  REAL hbyv( ip1jm,llm ), hbxu( ip1jmp1,llm )
32 
33 c
34 
35  DO 5 l = 1,llm
36 
37  DO 1 ij = iip2, ip1jm - 1
38  hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( teta(ij,l) + teta(ij+1,l) )
39  1 CONTINUE
40 
41 c .... correction pour hbxu(iip1,j,l) .....
42 c .... hbxu(iip1,j,l)= hbxu(1,j,l) ....
43 
44 CDIR$ IVDEP
45  DO 2 ij = iip1+ iip1, ip1jm, iip1
46  hbxu( ij, l ) = hbxu( ij - iim, l )
47  2 CONTINUE
48 
49 
50  DO 3 ij = 1,ip1jm
51  hbyv(ij,l)= pbarv(ij, l)* 0.5 * ( teta(ij, l)+ teta(ij +iip1,l) )
52  3 CONTINUE
53 
54  5 CONTINUE
55 
56 
57  CALL convflu ( hbxu, hbyv, llm, dteta )
58 
59 
60 c stockage dans dh de la convergence horizont. filtree' du flux
61 c .... ...........
62 c d'enthalpie potentielle .
63 
64  CALL filtreg( dteta, jjp1, llm, 2, 2, .true., 1)
65 
66 c
67  RETURN
68  END
iip2
!$Header iip2
Definition: paramet.h:14
ip1jmp1
!$Header llmm1 INTEGER ip1jmp1
Definition: paramet.h:14
llm
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm
Definition: 1D_interp_cases.h:103
convflu
subroutine convflu(xflu, yflu, nbniv, convfl)
Definition: convflu.F:5
ip1jm
!$Header llmm1 INTEGER ip1jm
Definition: paramet.h:14
jjp1
!$Header jjp1
Definition: paramet.h:14
true
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
Definition: calcul_STDlev.h:139
dteta1
subroutine dteta1(teta, pbaru, pbarv, dteta)
Definition: dteta1.F:5
iim
c c zjulian c cym CALL iim cym klev iim
Definition: ini_bilKP_ave.h:24
filtreg
subroutine filtreg(champ, nlat, nbniv, ifiltre, iaire, griscal, iter)
Definition: filtreg.F:6