Context Navigation

source: LMDZ5/branches/testing/libf/phylmd/ini_histrac.h @ 1783

Last change on this file since 1783 was 1750, checked in by Laurent Fairhead, 12 years ago

Version testing basée sur r1745

Testing release based on r1745

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 11.0 KB

Line
1	!
2	! $Id $
3	!
4	IF (ecrit_tra>0.) THEN
5	!$OMP MASTER
6	!!! CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
7	! correction pour l heure initiale !jyg
8	! !jyg
9	CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian) !jyg
10
11	CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
12	CALL histvert(nid_tra, "presnivs", "Vertical levels", "Pa",klev, presnivs, nvert,"down")
13
14	zsto = pdtphys
15	zout = ecrit_tra
16	CALL histdef(nid_tra, "phis", "Surface geop. height", "-", &
17	iim,jj_nb,nhori, 1,1,1, -99, 32,"once", zsto,zout)
18	CALL histdef(nid_tra, "aire", "Grid area", "-", &
19	iim,jj_nb,nhori, 1,1,1, -99, 32,"once", zsto,zout)
20	CALL histdef(nid_tra, "zmasse", "column density of air in cell", &
21	"kg m-2", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
22	zsto,zout)
23	! RomP >>>
24	CALL histdef(nid_tra, "sourceBE", "source 7Be", &
25	"at/kgA/s", iim, jj_nb, nhori, klev, 1, klev, nvert, 32, "ave(X)", &
26	zsto,zout)
27	! RomP <<<
28
29	!TRACEURS
30	!----------------
31	DO it = 1,nbtr
32	iiq = niadv(it+2)
33
34	! CONCENTRATIONS
35	CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga", &
36	iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
37
38	! TD LESSIVAGE
39	IF (lessivage .AND. aerosol(it)) THEN
40	CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
41	"at/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
42	"ave(X)", zsto,zout)
43	CALL histdef(nid_tra, "d_tr_ls_"//tname(iiq), &
44	"tendance lessivage large scale"// ttext(iiq), "?",&
45	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
46	"ave(X)", zsto,zout)
47	CALL histdef(nid_tra, "d_tr_insc_"//tname(iiq), &
48	"tendance lessivage large scale"// ttext(iiq), "?",&
49	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
50	"ave(X)", zsto,zout)
51	CALL histdef(nid_tra, "d_tr_bcscav_"//tname(iiq), &
52	"tendance lessivage large scale"// ttext(iiq), "?",&
53	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
54	"ave(X)", zsto,zout)
55	CALL histdef(nid_tra, "d_tr_evls_"//tname(iiq), &
56	"tendance lessivage large scale"// ttext(iiq), "?",&
57	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
58	"ave(X)", zsto,zout)
59	! Tracer concentration in LS precipitation at surface
60	CALL histdef(nid_tra, "qpr_ls_"//tname(iiq), &
61	"concentration in LS precip"// ttext(iiq), "at/kgw", &
62	iim,jj_nb,nhori, 1,1,1, -99, 32, &
63	"ave(X)", zsto,zout)
64	END IF
65
66	! TD THERMIQUES
67	IF (iflag_thermals.gt.0) THEN
68	CALL histdef(nid_tra, "d_tr_th_"//tname(iiq), &
69	"tendance thermique"// ttext(iiq), "?", &
70	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
71	"ave(X)", zsto,zout)
72	ENDIF
73
74	! TD CONVECTION
75	IF (iflag_con.GE.2) THEN
76	CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq), &
77	"tendance convection"// ttext(iiq), "?", &
78	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
79	"ave(X)", zsto,zout)
80	ENDIF
81
82	! RomP >>>
83	IF (iflag_con.EQ.30) THEN
84	CALL histdef(nid_tra, "d_tr_cvMA_"//tname(iiq), &
85	"tendance convection"// ttext(iiq), "?",&
86	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
87	"ave(X)", zsto,zout)
88	CALL histdef(nid_tra, "d_tr_trsp_"//tname(iiq), &
89	"tendance transport "// ttext(iiq), "at/kga", &
90	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
91	"ave(X)", zsto,zout)
92	CALL histdef(nid_tra, "d_tr_sscav_"//tname(iiq), &
93	"tendance lessivage flux satures "// ttext(iiq), "at/kga", &
94	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
95	"ave(X)", zsto,zout)
96	CALL histdef(nid_tra, "d_tr_sat_"//tname(iiq), &
97	"tendance flux satures "// ttext(iiq), "at/kga", &
98	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
99	"ave(X)", zsto,zout)
100	CALL histdef(nid_tra, "d_tr_uscav_"//tname(iiq), &
101	"tendance flux insatures "// ttext(iiq), "at/kga", &
102	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
103	"ave(X)", zsto,zout)
104	CALL histdef(nid_tra, "tr_pr_"//tname(iiq), &
105	"concentration dans precip"// ttext(iiq), "at/kga", &
106	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
107	"ave(X)", zsto,zout)
108	CALL histdef(nid_tra, "tr_aa_"//tname(iiq), &
109	"concentration precip issu AA"// ttext(iiq), "at/kga", &
110	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
111	"ave(X)", zsto,zout)
112	CALL histdef(nid_tra, "tr_mel_"//tname(iiq), &
113	"concentration precip issu melange"// ttext(iiq), "at/kga", &
114	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
115	"ave(X)", zsto,zout)
116	CALL histdef(nid_tra, "tr_di_"//tname(iiq), &
117	"concentration dans descente insaturee"// ttext(iiq), "at/kga", &
118	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
119	"ave(X)", zsto,zout)
120	CALL histdef(nid_tra, "tr_trspdi_"//tname(iiq), &
121	"conc descente insaturee MA"// ttext(iiq), "at/kga", &
122	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
123	"ave(X)", zsto,zout)
124	CALL histdef(nid_tra, "zmfd1a_"//tname(iiq), &
125	"zmfd1a"// ttext(iiq), "_", &
126	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
127	"ave(X)", zsto,zout)
128	CALL histdef(nid_tra, "zmfphi2_"//tname(iiq), &
129	"zmfphi2"// ttext(iiq), "_", &
130	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
131	"ave(X)", zsto,zout)
132	CALL histdef(nid_tra, "zmfdam_"//tname(iiq), &
133	"zmfdam"// ttext(iiq), "_", &
134	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
135	"ave(X)", zsto,zout)
136	ENDIF
137	! RomP <<<
138	CALL histdef(nid_tra, "dtrdyn_"//tname(iiq), &
139	"td dyn tra"// ttext(iiq), "at/kga", &
140	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
141	"ave(X)", zsto,zout)
142	! TD decroissance radioactive
143	CALL histdef(nid_tra, "d_tr_dec_"//tname(iiq), &
144	"tendance decroi radio "// ttext(iiq), "", &
145	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
146	"ave(X)", zsto,zout)
147
148	! TD COUCHE-LIMITE
149	IF (couchelimite) THEN
150	CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq), &
151	"tendance couche limite"// ttext(iiq), "?", &
152	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
153	"ave(X)", zsto,zout)
154	! Dry deposit (1st layer and surface)
155	CALL histdef(nid_tra, "d_tr_dry_"//tname(iiq), &
156	"tendancy dry deposit"// ttext(iiq), "at/kga/step", &
157	iim,jj_nb,nhori, 1,1,1, -99, 32, &
158	"ave(X)", zsto,zout)
159	CALL histdef(nid_tra, "flux_tr_dry_"//tname(iiq), &
160	"dry deposit at surf (downward)"// ttext(iiq), "at/m2/step", &
161	iim,jj_nb,nhori, 1,1,1, -99, 32, &
162	"ave(X)", zsto,zout)
163	ENDIF
164	ENDDO
165
166	CALL histdef(nid_tra, "Mint", "Mint","", &
167	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
168	"inst(X)", zout,zout)
169	CALL histdef(nid_tra, "frac_impa", "frac_impa","", &
170	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
171	"inst(X)", zout,zout)
172	CALL histdef(nid_tra, "frac_nucl", "frac_nucl","", &
173	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
174	"inst(X)", zout,zout)
175	!---------------
176	!
177	! VENT (niveau 1)
178	CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", &
179	iim,jj_nb,nhori, 1,1,1, -99, 32, &
180	"inst(X)", zout,zout)
181	CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", &
182	iim,jj_nb,nhori, 1,1,1, -99, 32, &
183	"inst(X)", zout,zout)
184
185	! TEMPERATURE DU SOL
186	CALL histdef(nid_tra, "ftsol1", "temper sol", "-", &
187	iim,jj_nb,nhori, 1,1,1, -99, 32, &
188	"inst(X)", zout,zout)
189	CALL histdef(nid_tra, "ftsol2", "temper sol", "-", &
190	iim,jj_nb,nhori, 1,1,1, -99, 32, &
191	"inst(X)", zout,zout)
192	CALL histdef(nid_tra, "ftsol3", "temper sol", "-", &
193	iim,jj_nb,nhori, 1,1,1, -99, 32, &
194	"inst", zout,zout)
195	CALL histdef(nid_tra, "ftsol4", "temper sol", "-", &
196	iim,jj_nb,nhori, 1,1,1, -99, 32, &
197	"inst(X)", zout,zout)
198
199	! NATURE DU SOL
200	CALL histdef(nid_tra, "psrf1", "nature sol", "-", &
201	iim,jj_nb,nhori, 1,1,1, -99, 32, &
202	"inst(X)", zout,zout)
203	CALL histdef(nid_tra, "psrf2", "nature sol", "-", &
204	iim,jj_nb,nhori, 1,1,1, -99, 32, &
205	"inst(X)", zout,zout)
206	CALL histdef(nid_tra, "psrf3", "nature sol", "-", &
207	iim,jj_nb,nhori, 1,1,1, -99, 32, &
208	"inst(X)", zout,zout)
209	CALL histdef(nid_tra, "psrf4", "nature sol", "-", &
210	iim,jj_nb,nhori, 1,1,1, -99, 32, &
211	"inst(X)", zout,zout)
212	! DIVERS
213	CALL histdef(nid_tra, "pplay", "pressure","-", &
214	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
215	"inst(X)", zout,zout)
216	CALL histdef(nid_tra, "T", "temperature","K", &
217	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
218	"inst(X)", zout,zout)
219	CALL histdef(nid_tra, "mfu", "flux u mont","-", &
220	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
221	"ave(X)", zsto,zout)
222	CALL histdef(nid_tra, "mfd", "flux u decen","-", &
223	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
224	"ave(X)", zsto,zout)
225	CALL histdef(nid_tra, "en_u", "flux u mont","-", &
226	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
227	"ave(X)", zsto,zout)
228	CALL histdef(nid_tra, "en_d", "flux u mont","-", &
229	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
230	"ave(X)", zsto,zout)
231	CALL histdef(nid_tra, "de_d", "flux u mont","-", &
232	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
233	"ave(X)", zsto,zout)
234	CALL histdef(nid_tra, "de_u", "flux u decen","-", &
235	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
236	"ave(X)", zsto,zout)
237	CALL histdef(nid_tra, "coefh", "turbulent coef","-", &
238	iim,jj_nb,nhori, klev,1,klev,nvert, 32, &
239	"ave(X)", zsto,zout)
240
241	CALL histend(nid_tra)
242	!$OMP END MASTER
243	END IF ! ecrit_tra>0.
244

Note: See TracBrowser for help on using the repository browser.

Download in other formats: