Context Navigation

phystokenc.F @ 5209

Last change on this file since 5209 was 1114, checked in by jghattas, 16 years ago

Creation du module infotrac:

contient les variables de advtrac.h
contient la subroutine iniadvtrac renommer en infotrac_init
le nombre des traceurs est lu dans tracer.def en dynamique (ou par default ou recu par INCA)
ce module est utilise dans la dynamique et la physique
contient aussi la variable nbtr qui avant etait stockee dans dimphy

Le fichier advtrac.h n'existe plus.
La compilation ne prend plus en compte le nombre de traceur.

/JG

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

File size: 13.3 KB

Line
1	!
2	c
3	c
4	SUBROUTINE phystokenc (
5	I nlon,nlev,pdtphys,rlon,rlat,
6	I pt,pmfu, pmfd, pen_u, pde_u, pen_d, pde_d,
7	I pfm_therm,pentr_therm,
8	I cdragh, pcoefh,yu1,yv1,ftsol,pctsrf,
9	I frac_impa,frac_nucl,
10	I pphis,paire,dtime,itap)
11	USE ioipsl
12	USE dimphy
13	USE infotrac, ONLY : nqtot
14	USE iophy
15	IMPLICIT none
16
17	c======================================================================
18	c Auteur(s) FH
19	c Objet: Moniteur general des tendances traceurs
20	c
21
22	c======================================================================
23	#include "dimensions.h"
24	#include "tracstoke.h"
25	#include "indicesol.h"
26	#include "control.h"
27	c======================================================================
28
29	c Arguments:
30	c
31	c EN ENTREE:
32	c ==========
33	c
34	c divers:
35	c -------
36	c
37	integer nlon ! nombre de points horizontaux
38	integer nlev ! nombre de couches verticales
39	real pdtphys ! pas d'integration pour la physique (seconde)
40	c
41	integer physid, itap
42	save physid
43	c$OMP THREADPRIVATE(physid)
44	integer ndex2d(iim(jjm+1)),ndex3d(iim(jjm+1)*klev)
45
46	c convection:
47	c -----------
48	c
49	REAL pmfu(klon,klev) ! flux de masse dans le panache montant
50	REAL pmfd(klon,klev) ! flux de masse dans le panache descendant
51	REAL pen_u(klon,klev) ! flux entraine dans le panache montant
52	REAL pde_u(klon,klev) ! flux detraine dans le panache montant
53	REAL pen_d(klon,klev) ! flux entraine dans le panache descendant
54	REAL pde_d(klon,klev) ! flux detraine dans le panache descendant
55	real pt(klon,klev)
56	REAL,allocatable,save :: t(:,:)
57	c$OMP THREADPRIVATE(t)
58	c
59	REAL rlon(klon), rlat(klon), dtime
60	REAL zx_tmp_3d(iim,jjm+1,klev),zx_tmp_2d(iim,jjm+1)
61
62	c Couche limite:
63	c --------------
64	c
65	REAL cdragh(klon) ! cdrag
66	REAL pcoefh(klon,klev) ! coeff melange CL
67	REAL pcoefh_buf(klon,klev) ! coeff melange CL + cdrag
68	REAL yv1(klon)
69	REAL yu1(klon),pphis(klon),paire(klon)
70
71	c Les Thermiques : (Abderr 25 11 02)
72	c ---------------
73	REAL pfm_therm(klon,klev+1)
74	real fm_therm1(klon,klev)
75	REAL pentr_therm(klon,klev)
76
77	REAL,allocatable,save :: entr_therm(:,:)
78	REAL,allocatable,save :: fm_therm(:,:)
79	c$OMP THREADPRIVATE(entr_therm)
80	c$OMP THREADPRIVATE(fm_therm)
81	c
82	c Lessivage:
83	c ----------
84	c
85	REAL frac_impa(klon,klev)
86	REAL frac_nucl(klon,klev)
87	c
88	c Arguments necessaires pour les sources et puits de traceur
89	C
90	real ftsol(klon,nbsrf) ! Temperature du sol (surf)(Kelvin)
91	real pctsrf(klon,nbsrf) ! Pourcentage de sol f(nature du sol)
92	c======================================================================
93	c
94	INTEGER i, k
95	c
96	REAL,allocatable,save :: mfu(:,:) ! flux de masse dans le panache montant
97	REAL,allocatable,save :: mfd(:,:) ! flux de masse dans le panache descendant
98	REAL,allocatable,save :: en_u(:,:) ! flux entraine dans le panache montant
99	REAL,allocatable,save :: de_u(:,:) ! flux detraine dans le panache montant
100	REAL,allocatable,save :: en_d(:,:) ! flux entraine dans le panache descendant
101	REAL,allocatable,save :: de_d(:,:) ! flux detraine dans le panache descendant
102	REAL,allocatable,save :: coefh(:,:) ! flux detraine dans le panache descendant
103
104	REAL,allocatable,save :: pyu1(:)
105	REAL,allocatable,save :: pyv1(:)
106	REAL,allocatable,save :: pftsol(:,:)
107	REAL,allocatable,save :: ppsrf(:,:)
108	c$OMP THREADPRIVATE(mfu,mfd,en_u,de_u,en_d,de_d,coefh)
109	c$OMP THREADPRIVATE(pyu1,pyv1,pftsol,ppsrf)
110	real pftsol1(klon),pftsol2(klon),pftsol3(klon),pftsol4(klon)
111	real ppsrf1(klon),ppsrf2(klon),ppsrf3(klon),ppsrf4(klon)
112
113	REAL dtcum
114
115	integer iadvtr,irec
116	real zmin,zmax
117	logical ok_sync
118
119	save dtcum
120	save iadvtr,irec
121	c$OMP THREADPRIVATE(dtcum,iadvtr,irec)
122	data iadvtr,irec/0,1/
123	logical,save :: first=.true.
124	c$OMP THREADPRIVATE(first)
125	c
126	c Couche limite:
127	c======================================================================
128
129	c Dans le meme vecteur on recombine le drag et les coeff d'echange
130	pcoefh_buf(:,1) = cdragh(:)
131	pcoefh_buf(:,2:klev) = pcoefh(:,2:klev)
132
133	ok_sync = .true.
134	print*,'Dans phystokenc.F'
135	print*,'iadvtr= ',iadvtr
136	print*,'istphy= ',istphy
137	print*,'istdyn= ',istdyn
138
139	if (first) then
140
141	allocate( t(klon,klev))
142	allocate( mfu(klon,klev))
143	allocate( mfd(klon,klev))
144	allocate( en_u(klon,klev))
145	allocate( de_u(klon,klev))
146	allocate( en_d(klon,klev))
147	allocate( de_d(klon,klev))
148	allocate( coefh(klon,klev))
149	allocate( entr_therm(klon,klev))
150	allocate( fm_therm(klon,klev))
151	allocate( pyu1(klon))
152	allocate( pyv1(klon))
153	allocate( pftsol(klon,nbsrf))
154	allocate( ppsrf(klon,nbsrf))
155
156	first=.false.
157	endif
158
159	IF (iadvtr.eq.0) THEN
160
161	CALL initphysto('phystoke',
162	. rlon,rlat,dtime, dtimeistphy,dtimeistphy,nqtot,physid)
163
164	write(,) 'apres initphysto ds phystokenc'
165
166
167	ENDIF
168	c
169	ndex2d = 0
170	ndex3d = 0
171	i=itap
172	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
173	CALL histwrite_phy(physid,"phis",i,pphis)
174	c
175	i=itap
176	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
177	CALL histwrite_phy(physid,"aire",i,paire)
178
179	iadvtr=iadvtr+1
180	c
181	if (mod(iadvtr,istphy).eq.1.or.istphy.eq.1) then
182	print*,'reinitialisation des champs cumules
183	s a iadvtr=',iadvtr
184	do k=1,klev
185	do i=1,klon
186	mfu(i,k)=0.
187	mfd(i,k)=0.
188	en_u(i,k)=0.
189	de_u(i,k)=0.
190	en_d(i,k)=0.
191	de_d(i,k)=0.
192	coefh(i,k)=0.
193	t(i,k)=0.
194	fm_therm(i,k)=0.
195	entr_therm(i,k)=0.
196	enddo
197	enddo
198	do i=1,klon
199	pyv1(i)=0.
200	pyu1(i)=0.
201	end do
202	do k=1,nbsrf
203	do i=1,klon
204	pftsol(i,k)=0.
205	ppsrf(i,k)=0.
206	enddo
207	enddo
208
209	dtcum=0.
210	endif
211
212	do k=1,klev
213	do i=1,klon
214	mfu(i,k)=mfu(i,k)+pmfu(i,k)*pdtphys
215	mfd(i,k)=mfd(i,k)+pmfd(i,k)*pdtphys
216	en_u(i,k)=en_u(i,k)+pen_u(i,k)*pdtphys
217	de_u(i,k)=de_u(i,k)+pde_u(i,k)*pdtphys
218	en_d(i,k)=en_d(i,k)+pen_d(i,k)*pdtphys
219	de_d(i,k)=de_d(i,k)+pde_d(i,k)*pdtphys
220	coefh(i,k)=coefh(i,k)+pcoefh_buf(i,k)*pdtphys
221	t(i,k)=t(i,k)+pt(i,k)*pdtphys
222	fm_therm(i,k)=fm_therm(i,k)+pfm_therm(i,k)*pdtphys
223	entr_therm(i,k)=entr_therm(i,k)+pentr_therm(i,k)*pdtphys
224	enddo
225	enddo
226	do i=1,klon
227	pyv1(i)=pyv1(i)+yv1(i)*pdtphys
228	pyu1(i)=pyu1(i)+yu1(i)*pdtphys
229	end do
230	do k=1,nbsrf
231	do i=1,klon
232	pftsol(i,k)=pftsol(i,k)+ftsol(i,k)*pdtphys
233	ppsrf(i,k)=ppsrf(i,k)+pctsrf(i,k)*pdtphys
234	enddo
235	enddo
236
237	dtcum=dtcum+pdtphys
238
239	IF(mod(iadvtr,istphy).eq.0) THEN
240	c
241	c normalisation par le temps cumule
242	do k=1,klev
243	do i=1,klon
244	mfu(i,k)=mfu(i,k)/dtcum
245	mfd(i,k)=mfd(i,k)/dtcum
246	en_u(i,k)=en_u(i,k)/dtcum
247	de_u(i,k)=de_u(i,k)/dtcum
248	en_d(i,k)=en_d(i,k)/dtcum
249	de_d(i,k)=de_d(i,k)/dtcum
250	coefh(i,k)=coefh(i,k)/dtcum
251	c Unitel a enlever
252	t(i,k)=t(i,k)/dtcum
253	fm_therm(i,k)=fm_therm(i,k)/dtcum
254	entr_therm(i,k)=entr_therm(i,k)/dtcum
255	enddo
256	enddo
257	do i=1,klon
258	pyv1(i)=pyv1(i)/dtcum
259	pyu1(i)=pyu1(i)/dtcum
260	end do
261	do k=1,nbsrf
262	do i=1,klon
263	pftsol(i,k)=pftsol(i,k)/dtcum
264	pftsol1(i) = pftsol(i,1)
265	pftsol2(i) = pftsol(i,2)
266	pftsol3(i) = pftsol(i,3)
267	pftsol4(i) = pftsol(i,4)
268
269	ppsrf(i,k)=ppsrf(i,k)/dtcum
270	ppsrf1(i) = ppsrf(i,1)
271	ppsrf2(i) = ppsrf(i,2)
272	ppsrf3(i) = ppsrf(i,3)
273	ppsrf4(i) = ppsrf(i,4)
274
275	enddo
276	enddo
277	c
278	c ecriture des champs
279	c
280	irec=irec+1
281
282	ccccc
283	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)
284	CALL histwrite_phy(physid,"t",itap,t)
285
286	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d)
287	CALL histwrite_phy(physid,"mfu",itap,mfu)
288	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d)
289	CALL histwrite_phy(physid,"mfd",itap,mfd)
290	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d)
291	CALL histwrite_phy(physid,"en_u",itap,en_u)
292	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d)
293	CALL histwrite_phy(physid,"de_u",itap,de_u)
294	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d)
295	CALL histwrite_phy(physid,"en_d",itap,en_d)
296	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d)
297	CALL histwrite_phy(physid,"de_d",itap,de_d)
298	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d)
299	CALL histwrite_phy(physid,"coefh",itap,coefh)
300
301	c ajou...
302	do k=1,klev
303	do i=1,klon
304	fm_therm1(i,k)=fm_therm(i,k)
305	enddo
306	enddo
307
308	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d)
309	CALL histwrite_phy(physid,"fm_th",itap,fm_therm1)
310	c
311	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d)
312	CALL histwrite_phy(physid,"en_th",itap,entr_therm)
313	cccc
314	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d)
315	CALL histwrite_phy(physid,"frac_impa",itap,frac_impa)
316
317	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d)
318	CALL histwrite_phy(physid,"frac_nucl",itap,frac_nucl)
319
320	cym CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d)
321	CALL histwrite_phy(physid,"pyu1",itap,pyu1)
322
323	cym CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d)
324	CALL histwrite_phy(physid,"pyv1",itap,pyv1)
325
326	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d)
327	CALL histwrite_phy(physid,"ftsol1",itap,pftsol1)
328	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d)
329	CALL histwrite_phy(physid,"ftsol2",itap,pftsol2)
330	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d)
331	CALL histwrite_phy(physid,"ftsol3",itap,pftsol3)
332	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d)
333	CALL histwrite_phy(physid,"ftsol4",itap,pftsol4)
334
335	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d)
336	CALL histwrite_phy(physid,"psrf1",itap,ppsrf1)
337	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d)
338	CALL histwrite_phy(physid,"psrf2",itap,ppsrf2)
339	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d)
340	CALL histwrite_phy(physid,"psrf3",itap,ppsrf3)
341	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d)
342	CALL histwrite_phy(physid,"psrf4",itap,ppsrf4)
343
344	c$OMP MASTER
345	if (ok_sync) call histsync(physid)
346	c$OMP END MASTER
347	c if (ok_sync) call histsync
348
349	c
350	cAA Test sur la valeur des coefficients de lessivage
351	c
352	zmin=1e33
353	zmax=-1e33
354	do k=1,klev
355	do i=1,klon
356	zmax=max(zmax,frac_nucl(i,k))
357	zmin=min(zmin,frac_nucl(i,k))
358	enddo
359	enddo
360	Print*,'------ coefs de lessivage (min et max) --------'
361	Print*,'facteur de nucleation ',zmin,zmax
362	zmin=1e33
363	zmax=-1e33
364	do k=1,klev
365	do i=1,klon
366	zmax=max(zmax,frac_impa(i,k))
367	zmin=min(zmin,frac_impa(i,k))
368	enddo
369	enddo
370	Print*,'facteur d impaction ',zmin,zmax
371
372	ENDIF
373
374	c reinitialisation des champs cumules
375	go to 768
376	if (mod(iadvtr,istphy).eq.1) then
377	do k=1,klev
378	do i=1,klon
379	mfu(i,k)=0.
380	mfd(i,k)=0.
381	en_u(i,k)=0.
382	de_u(i,k)=0.
383	en_d(i,k)=0.
384	de_d(i,k)=0.
385	coefh(i,k)=0.
386	t(i,k)=0.
387	fm_therm(i,k)=0.
388	entr_therm(i,k)=0.
389	enddo
390	enddo
391	do i=1,klon
392	pyv1(i)=0.
393	pyu1(i)=0.
394	end do
395	do k=1,nbsrf
396	do i=1,klon
397	pftsol(i,k)=0.
398	ppsrf(i,k)=0.
399	enddo
400	enddo
401
402	dtcum=0.
403	endif
404
405	do k=1,klev
406	do i=1,klon
407	mfu(i,k)=mfu(i,k)+pmfu(i,k)*pdtphys
408	mfd(i,k)=mfd(i,k)+pmfd(i,k)*pdtphys
409	en_u(i,k)=en_u(i,k)+pen_u(i,k)*pdtphys
410	de_u(i,k)=de_u(i,k)+pde_u(i,k)*pdtphys
411	en_d(i,k)=en_d(i,k)+pen_d(i,k)*pdtphys
412	de_d(i,k)=de_d(i,k)+pde_d(i,k)*pdtphys
413	coefh(i,k)=coefh(i,k)+pcoefh_buf(i,k)*pdtphys
414	t(i,k)=t(i,k)+pt(i,k)*pdtphys
415	fm_therm(i,k)=fm_therm(i,k)+pfm_therm(i,k)*pdtphys
416	entr_therm(i,k)=entr_therm(i,k)+pentr_therm(i,k)*pdtphys
417	enddo
418	enddo
419	do i=1,klon
420	pyv1(i)=pyv1(i)+yv1(i)*pdtphys
421	pyu1(i)=pyu1(i)+yu1(i)*pdtphys
422	end do
423	do k=1,nbsrf
424	do i=1,klon
425	pftsol(i,k)=pftsol(i,k)+ftsol(i,k)*pdtphys
426	ppsrf(i,k)=ppsrf(i,k)+pctsrf(i,k)*pdtphys
427	enddo
428	enddo
429
430	dtcum=dtcum+pdtphys
431	768 continue
432
433	RETURN
434	END

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

Download in other formats:

Original Format