Context Navigation

phystokenc.F @ 1085

Last change on this file since 1085 was 1067, checked in by Laurent Fairhead, 16 years ago

Modifications lie au premier niveau du modele pour la diffusion turbulent

du vent.

Preparation pour un couplage des courrant oceaniques.

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

File size: 13.3 KB

Rev	Line
[524]	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,
[541]	7	I pfm_therm,pentr_therm,
[1067]	8	I cdragh, pcoefh,yu1,yv1,ftsol,pctsrf,
[541]	9	I frac_impa,frac_nucl,
[524]	10	I pphis,paire,dtime,itap)
	11	USE ioipsl
[766]	12	USE dimphy
	13	USE iophy
[524]	14	IMPLICIT none
	15
	16	c======================================================================
	17	c Auteur(s) FH
	18	c Objet: Moniteur general des tendances traceurs
	19	c
	20
	21	c======================================================================
	22	#include "dimensions.h"
[766]	23	cym#include "dimphy.h"
[524]	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
[541]	41	integer physid, itap
	42	save physid
[766]	43	c$OMP THREADPRIVATE(physid)
[541]	44	integer ndex2d(iim(jjm+1)),ndex3d(iim(jjm+1)*klev)
[524]	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
[766]	55	real pt(klon,klev)
	56	REAL,allocatable,save :: t(:,:)
	57	c$OMP THREADPRIVATE(t)
[524]	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
[1067]	65	REAL cdragh(klon) ! cdrag
	66	REAL pcoefh(klon,klev) ! coeff melange CL
	67	REAL pcoefh_buf(klon,klev) ! coeff melange CL + cdrag
[524]	68	REAL yv1(klon)
	69	REAL yu1(klon),pphis(klon),paire(klon)
[541]	70
	71	c Les Thermiques : (Abderr 25 11 02)
	72	c ---------------
	73	REAL pfm_therm(klon,klev+1)
[766]	74	real fm_therm1(klon,klev)
[541]	75	REAL pentr_therm(klon,klev)
[766]	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)
[524]	81	c
	82	c Lessivage:
	83	c ----------
	84	c
[541]	85	REAL frac_impa(klon,klev)
	86	REAL frac_nucl(klon,klev)
[524]	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
[766]	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
[524]	103
[766]	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)
[524]	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
[541]	117	logical ok_sync
	118
[766]	119	save dtcum
[524]	120	save iadvtr,irec
[766]	121	c$OMP THREADPRIVATE(dtcum,iadvtr,irec)
[524]	122	data iadvtr,irec/0,1/
[766]	123	logical,save :: first=.true.
	124	c$OMP THREADPRIVATE(first)
[524]	125	c
	126	c Couche limite:
	127	c======================================================================
	128
[1067]	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
[541]	133	ok_sync = .true.
	134	print*,'Dans phystokenc.F'
[524]	135	print*,'iadvtr= ',iadvtr
	136	print*,'istphy= ',istphy
	137	print*,'istdyn= ',istdyn
	138
[766]	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
[524]	159	IF (iadvtr.eq.0) THEN
	160
	161	CALL initphysto('phystoke',
	162	. rlon,rlat,dtime, dtimeistphy,dtimeistphy,nqmx,physid)
	163
	164	write(,) 'apres initphysto ds phystokenc'
	165
	166
	167	ENDIF
	168	c
[541]	169	ndex2d = 0
	170	ndex3d = 0
	171	i=itap
[766]	172	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
	173	CALL histwrite_phy(physid,"phis",i,pphis)
[541]	174	c
	175	i=itap
[766]	176	cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
	177	CALL histwrite_phy(physid,"aire",i,paire)
[541]	178
[524]	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.
[541]	193	t(i,k)=0.
	194	fm_therm(i,k)=0.
	195	entr_therm(i,k)=0.
[524]	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
[1067]	220	coefh(i,k)=coefh(i,k)+pcoefh_buf(i,k)*pdtphys
[541]	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
[524]	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
[541]	238
	239	IF(mod(iadvtr,istphy).eq.0) THEN
[524]	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
[541]	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
[524]	255	enddo
	256	enddo
	257	do i=1,klon
	258	pyv1(i)=pyv1(i)/dtcum
	259	pyu1(i)=pyu1(i)/dtcum
	260	end do
[541]	261	do k=1,nbsrf
[524]	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
[541]	269	ppsrf(i,k)=ppsrf(i,k)/dtcum
[524]	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
[541]	276	enddo
[524]	277	c
	278	c ecriture des champs
	279	c
	280	irec=irec+1
	281
	282	ccccc
[766]	283	cym CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)
	284	CALL histwrite_phy(physid,"t",itap,t)
[524]	285
[766]	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)
[541]	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
[766]	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)
[541]	310	c
[766]	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)
[524]	313	cccc
[766]	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)
[524]	316
[766]	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)
[541]	319
[766]	320	cym CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d)
	321	CALL histwrite_phy(physid,"pyu1",itap,pyu1)
[541]	322
[766]	323	cym CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d)
	324	CALL histwrite_phy(physid,"pyv1",itap,pyv1)
[541]	325
[766]	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)
[524]	334
[766]	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)
[524]	343
[766]	344	c$OMP MASTER
[541]	345	if (ok_sync) call histsync(physid)
[766]	346	c$OMP END MASTER
[541]	347	c if (ok_sync) call histsync
	348
[524]	349	c
[541]	350	cAA Test sur la valeur des coefficients de lessivage
[524]	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
[541]	372	ENDIF
[524]	373
[541]	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
[524]	401
[541]	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
[1067]	413	coefh(i,k)=coefh(i,k)+pcoefh_buf(i,k)*pdtphys
[541]	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
[524]	433	RETURN
	434	END

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

Download in other formats:

Original Format