Context Navigation

source: LMDZ4/trunk/libf/phylmd/write_histrac.h @ 717

Last change on this file since 717 was 684, checked in by lmdzadmin, 19 years ago
Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def + Lecture freq. ecriture en nombre de jours IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 21.9 KB

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	ndex = 0
	5	ndex2d = 0
	6	ndex3d = 0
	7	c
	8	itau_w = itau_phy + nstep
	9
	10	CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
	11	CALL histwrite(nid_tra,"phis",itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	12	C
	13	CALL gr_fi_ecrit(1,klon,iim,jjm+1,airephy,zx_tmp_2d)
	14	CALL histwrite(nid_tra,"aire",itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	15	#ifdef INCA
	16	CALL gr_fi_ecrit(1, klon,iim,jjm+1, ps,zx_tmp_2d)
	17	CALL histwrite(nid_tra,"ps",itau_w,zx_tmp_2d,
	18	. iim*(jjm+1),ndex2d)
	19
	20	CALL gr_fi_ecrit(1, klon,iim,jjm+1, ptrop,zx_tmp_2d)
	21	CALL histwrite(nid_tra,"ptrop",itau_w,zx_tmp_2d,
	22	. iim*(jjm+1),ndex2d)
	23
	24	C 3d FIELDS
	25
	26	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,t_seri, zx_tmp_3d)
	27	CALL histwrite(nid_tra,"temp",itau_w,zx_tmp_3d,
	28	. iim(jjm+1)klev,ndex3d)
	29
	30	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,u, zx_tmp_3d)
	31	CALL histwrite(nid_tra,"u",itau_w,zx_tmp_3d,
	32	. iim(jjm+1)klev,ndex3d)
	33
	34	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,v, zx_tmp_3d)
	35	CALL histwrite(nid_tra,"v",itau_w,zx_tmp_3d,
	36	. iim(jjm+1)klev,ndex3d)
	37
	38	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,sh, zx_tmp_3d)
	39	CALL histwrite(nid_tra,"h2o",itau_w,zx_tmp_3d,
	40	. iim(jjm+1)klev,ndex3d)
	41
	42	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pdel, zx_tmp_3d)
	43	CALL histwrite(nid_tra,"pdel",itau_w,zx_tmp_3d,
	44	. iim(jjm+1)klev,ndex3d)
	45
	46	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pplay, zx_tmp_3d)
	47	CALL histwrite(nid_tra,"pmid",itau_w,zx_tmp_3d,
	48	. iim(jjm+1)klev,ndex3d)
	49
[616]	50	! Ajout Anne
[635]	51	#ifdef INCA_AER
[624]	52	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,airm, zx_tmp_3d)
	53	CALL histwrite(nid_tra,"airm",itau_w,zx_tmp_3d,
	54	. iim(jjm+1)klev,ndex)
[635]	55	#endif
[646]	56
[616]	57	! Fin ajout Anne
	58
	59	#ifdef INCA
[524]	60	#ifdef INCAINFO
[673]	61	#ifdef INCA_NMHC
[524]	62	DO it=1, phtcnt
	63	WRITE(str2,'(i2.2)') it
	64	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it),
	65	. zx_tmp_3d)
	66	CALL histwrite(nid_tra,"j"//str2,itau_w,zx_tmp_3d,
	67	. iim(jjm+1)klev,ndex3d)
	68	ENDDO
	69
	70	DO it=1, hetcnt
	71	WRITE(str2,'(i2.2)') it
	72	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it),
	73	. zx_tmp_3d)
	74	CALL histwrite(nid_tra,"w"//str2,itau_w,zx_tmp_3d,
	75	. iim(jjm+1)klev,ndex3d)
	76	ENDDO
	77
	78	DO it=1, extcnt
	79	WRITE(str2,'(i2.2)') it
	80	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it),
	81	. zx_tmp_3d)
	82	CALL histwrite(nid_tra,"ext"//str2,itau_w,zx_tmp_3d,
	83	. iim(jjm+1)klev,ndex3d)
	84	ENDDO
[673]	85	#endif
	86	#ifdef INCA_CH4
	87	DO it=1, phtcnt
	88	WRITE(str2,'(i2.2)') it
	89	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it),
	90	. zx_tmp_3d)
	91	CALL histwrite(nid_tra,"j"//str2,itau_w,zx_tmp_3d,
	92	. iim(jjm+1)klev,ndex3d)
	93	ENDDO
[524]	94
[673]	95	DO it=1, hetcnt
	96	WRITE(str2,'(i2.2)') it
	97	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it),
	98	. zx_tmp_3d)
	99	CALL histwrite(nid_tra,"w"//str2,itau_w,zx_tmp_3d,
	100	. iim(jjm+1)klev,ndex3d)
	101	ENDDO
	102
	103	DO it=1, extcnt
	104	WRITE(str2,'(i2.2)') it
	105	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it),
	106	. zx_tmp_3d)
	107	CALL histwrite(nid_tra,"ext"//str2,itau_w,zx_tmp_3d,
	108	. iim(jjm+1)klev,ndex3d)
	109	ENDDO
	110	#endif
	111
[524]	112	DO it=1, nfs
	113	WRITE(str2,'(i2.2)') it
	114	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,invariants(1,1,it),
	115	. zx_tmp_3d)
	116	CALL histwrite(nid_tra,"INV"//str2,itau_w,zx_tmp_3d,
	117	. iim(jjm+1)klev,ndex3d)
	118	ENDDO
	119	#else
[673]	120	#ifdef INCA_NMHC
[524]	121	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
	122	. zx_tmp_3d)
	123	CALL histwrite(nid_tra,"jO3",itau_w,zx_tmp_3d,
	124	. iim(jjm+1)klev,ndex3d)
	125
	126	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4),
	127	. zx_tmp_3d)
	128	CALL histwrite(nid_tra,"jNO2",itau_w,zx_tmp_3d,
	129	. iim(jjm+1)klev,ndex3d)
	130
	131	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13),
	132	. zx_tmp_3d)
	133	CALL histwrite(nid_tra,"jH2O2",itau_w,zx_tmp_3d,
	134	. iim(jjm+1)klev,ndex3d)
	135
	136	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1),
	137	. zx_tmp_3d)
	138	CALL histwrite(nid_tra,"wHNO3",itau_w,zx_tmp_3d,
	139	. iim(jjm+1)klev,ndex3d)
	140
	141	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1),
	142	. zx_tmp_3d)
	143	CALL histwrite(nid_tra,"kN2O5",itau_w,zx_tmp_3d,
	144	. iim(jjm+1)klev,ndex3d)
	145
	146	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1),
	147	. zx_tmp_3d)
	148	CALL histwrite(nid_tra,"LghtNO",itau_w,zx_tmp_3d,
	149	. iim(jjm+1)klev,ndex3d)
	150	#endif
[673]	151	#ifdef INCA_CH4
	152	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
	153	. zx_tmp_3d)
	154	CALL histwrite(nid_tra,"jO3",itau_w,zx_tmp_3d,
	155	. iim(jjm+1)klev,ndex3d)
[524]	156
[673]	157	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4),
	158	. zx_tmp_3d)
	159	CALL histwrite(nid_tra,"jNO2",itau_w,zx_tmp_3d,
[524]	160	. iim(jjm+1)klev,ndex3d)
[673]	161
	162	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13),
	163	. zx_tmp_3d)
	164	CALL histwrite(nid_tra,"jH2O2",itau_w,zx_tmp_3d,
	165	. iim(jjm+1)klev,ndex3d)
	166
	167	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1),
	168	. zx_tmp_3d)
	169	CALL histwrite(nid_tra,"wHNO3",itau_w,zx_tmp_3d,
	170	. iim(jjm+1)klev,ndex3d)
	171
	172	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1),
	173	. zx_tmp_3d)
	174	CALL histwrite(nid_tra,"kN2O5",itau_w,zx_tmp_3d,
	175	. iim(jjm+1)klev,ndex3d)
	176	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1),
	177	. zx_tmp_3d)
	178	CALL histwrite(nid_tra,"LghtNO",itau_w,zx_tmp_3d,
	179	. iim(jjm+1)klev,ndex3d)
[524]	180	#endif
[673]	181	#endif
	182	! DO it=1, grpcnt
	183	!
	184	! CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
	185	! zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
	186	! CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
	187	! . iim(jjm+1)klev,ndex3d)
	188	! ENDDO
	189	#endif
[524]	190
	191	#ifdef INCA_AER
	192
	193	it = id_CIDUSTM
	194	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_st(1,1,it),
	195	. zx_tmp_3d)
[616]	196	CALL histwrite(nid_tra2,"scavcoef_st",itau_w,zx_tmp_3d,
[524]	197	. iim(jjm+1)klev,ndex3d)
	198	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_cv(1,1,it),
	199	. zx_tmp_3d)
[616]	200	CALL histwrite(nid_tra2,"scavcoef_cv",itau_w,zx_tmp_3d,
[524]	201	. iim(jjm+1)klev,ndex3d)
	202
	203	CALL gr_fi_ecrit(1, klon,iim,jjm+1,angst(1),zx_tmp_2d)
	204	CALL histwrite(nid_tra2,"AngstroemComp",itau_w,zx_tmp_2d,
	205	. iim*(jjm+1),ndex2d)
	206
[646]	207	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,totaerh2o(1,1), zx_tmp_3d)
	208	CALL histwrite(nid_tra2,"TOTAERH2O",itau_w,zx_tmp_3d,
	209	. iim(jjm+1)klev,ndex3d)
	210
[524]	211	#endif
	212	#endif
	213
	214	DO it=1,nqmax
	215	C champs 2D
	216	#ifdef INCA
[679]	217	IF ( prt_flag_ts(it) == 0 ) CYCLE
[673]	218
[524]	219	CALL gr_fi_ecrit(1, klon,iim,jjm+1, eflux(1,it),zx_tmp_2d)
	220	CALL histwrite(nid_tra,"Emi_"//solsym(it),itau_w,zx_tmp_2d,
	221	. iim*(jjm+1),ndex2d)
	222
	223	CALL gr_fi_ecrit(1, klon,iim,jjm+1, dvel(1,it),zx_tmp_2d)
	224	CALL histwrite(nid_tra,"Dep_"//solsym(it),itau_w,zx_tmp_2d,
	225	. iim*(jjm+1),ndex2d)
	226	#ifdef INCA_AER
[646]	227	call diag(airephy,tr_seri,.false.)
[524]	228
[646]	229	CALL gr_fi_ecrit(1, klon,iim,jjm+1,dflux(1,it),zx_tmp_2d)
	230	CALL histwrite(nid_tra2,"Dry_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	231
[524]	232	IF ((it .ge. trmx) .and. (it .le. trnx)) then
	233	CALL gr_fi_ecrit(1, klon,iim,jjm+1,sflux(1,it),zx_tmp_2d)
[646]	234	CALL histwrite(nid_tra2,"Sed_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	235
[524]	236	CALL gr_fi_ecrit(1, klon,iim,jjm+1,wflux(1,it),zx_tmp_2d)
[646]	237	CALL histwrite(nid_tra2,"Wet_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
[524]	238	CALL gr_fi_ecrit(1, klon,iim,jjm+1,wsflux(1,it),zx_tmp_2d)
[646]	239	CALL histwrite(nid_tra2,"WetST_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
[524]	240	CALL gr_fi_ecrit(1, klon,iim,jjm+1,wcflux(1,it),zx_tmp_2d)
[646]	241	CALL histwrite(nid_tra2,"WetCV_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	242
[524]	243	CALL gr_fi_ecrit(klev, klon,iim,jjm+1,eflux_alt(1,1,it),zx_tmp_3d)
[646]	244	CALL histwrite(nid_tra2,"Emi_alt_"//solsym(it),itau_w,zx_tmp_3d,iim(jjm+1)klev,ndex3d)
	245
[524]	246	CALL gr_fi_ecrit(1, klon,iim,jjm+1,aload(1,it),zx_tmp_2d)
[646]	247	CALL histwrite(nid_tra2,"Load_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	248	CALL histwrite(nid_tra3,"Inst_Load_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
[524]	249
	250	CALL gr_fi_ecrit(1, klon,iim,jjm+1,sconc(1,it),zx_tmp_2d)
[646]	251	CALL histwrite(nid_tra2,"SConc_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
[524]	252
	253	do la=1,las
	254	CALL gr_fi_ecrit(1, klon,iim,jjm+1,tausum(1,la,it),zx_tmp_2d)
[646]	255	CALL histwrite(nid_tra2,"OD"//cla(la)//"_"//solsym(it),itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)
	256	CALL gr_fi_ecrit(klev, klon,iim,jjm+1,tau(1,1,la,it),zx_tmp_3d)
	257	CALL histwrite(nid_tra2,"OD3D"//cla(la)//"_"//solsym(it),itau_w,zx_tmp_3d,iim(jjm+1)klev,ndex3d)
[524]	258	enddo
	259
	260	CALL gr_fi_ecrit(klev, klon,iim,jjm+1,md(1,1,it),zx_tmp_3d)
[646]	261	CALL histwrite(nid_tra2,"MD_"//solsym(it),itau_w,zx_tmp_3d,iim(jjm+1)klev,ndex3d)
[524]	262
[646]	263	CALL gr_fi_ecrit(klev, klon,iim,jjm+1,mdw(1,1,it),zx_tmp_3d)
	264	CALL histwrite(nid_tra2,"MDW_"//solsym(it),itau_w,zx_tmp_3d,iim(jjm+1)klev,ndex3d)
	265
	266	CALL gr_fi_ecrit(klev, klon,iim,jjm+1,cload(1,1,it),zx_tmp_3d)
	267	CALL histwrite(nid_tra2,"CLOAD_"//solsym(it),itau_w,zx_tmp_3d,iim(jjm+1)klev,ndex3d)
	268
	269
[524]	270	endif
[646]	271
	272
[524]	273	#endif
	274	C champs 3D
	275	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
	276
	277	!Prefer vmr to mmr for transported species
	278	if( adv_mass(it) /= 0. ) then
	279	#ifdef INCA_AER
	280	if (it .lt. trmx) then
	281	#endif
	282	zx_tmp_3d = zx_tmp_3d * dry_mass / adv_mass(it)
	283	#ifdef INCA_AER
	284	endif
	285	#endif
	286	else
[673]	287	#ifdef INCA_CH4
[524]	288	if ( solsym(it) == 'OX' ) then
	289	zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
	290	end if
	291	#endif
[673]	292	#ifdef INCA_NMHC
	293	if ( solsym(it) == 'OX' ) then
	294	zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
[524]	295	end if
[673]	296	#endif
	297	end if
[524]	298
	299	CALL histwrite(nid_tra,solsym(it),itau_w,zx_tmp_3d,
	300	. iim(jjm+1)klev,ndex3d)
	301	#else
	302
	303	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
	304	CALL histwrite(nid_tra,tnom(it+2),itau_w,zx_tmp_3d,
	305	. iim(jjm+1)klev,ndex3d)
	306	if (lessivage) THEN
	307	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,flestottr(1,1,it),zx_tmp_3d)
	308	CALL histwrite(nid_tra,"fl"//tnom(it+2),itau_w,zx_tmp_3d,
	309	. iim(jjm+1)klev,ndex3d)
	310	endif
[541]	311
	312	c----Olivia
	313	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_th(1,1,it),zx_tmp_3d)
	314	CALL histwrite(nid_tra,"d_tr_th_"//tnom(it+2),itau_w,zx_tmp_3d,
	315	. iim(jjm+1)klev,ndex3d)
[619]	316	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cv(1,1,it),zx_tmp_3d)
[541]	317	CALL histwrite(nid_tra,"d_tr_cv_"//tnom(it+2),itau_w,zx_tmp_3d,
	318	. iim(jjm+1)klev,ndex3d)
[619]	319	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,d_tr_cl(1,1,it),zx_tmp_3d)
[541]	320	CALL histwrite(nid_tra,"d_tr_cl_"//tnom(it+2),itau_w,zx_tmp_3d,
	321	. iim(jjm+1)klev,ndex3d)
	322	c---fin Olivia
	323
[524]	324	#endif
	325	ENDDO
	326
	327	#ifdef INCA
[616]	328	!#ifdef INCA_CH4
	329	CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_col(1), zx_tmp_2d)
[524]	330	CALL histwrite(nid_tra,"O3_column",itau_w,zx_tmp_2d,
	331	. iim*(jjm+1),ndex2d)
	332
[616]	333	CALL gr_fi_ecrit(1, klon,iim,jjm+1, co_col(1), zx_tmp_2d)
[524]	334	CALL histwrite(nid_tra,"CO_column",itau_w,zx_tmp_2d,
	335	. iim*(jjm+1),ndex2d)
	336
[616]	337	CALL gr_fi_ecrit(1, klon,iim,jjm+1, ch4_col(1), zx_tmp_2d)
[524]	338	CALL histwrite(nid_tra,"CH4_column",itau_w,zx_tmp_2d,
	339	. iim*(jjm+1),ndex2d)
	340
[616]	341	CALL gr_fi_ecrit(1, klon,iim,jjm+1, no2_col(1), zx_tmp_2d)
[524]	342	CALL histwrite(nid_tra,"NO2_column",itau_w,zx_tmp_2d,
	343	. iim*(jjm+1),ndex2d)
	344
	345	CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_st_flx(1), zx_tmp_2d)
	346	CALL histwrite(nid_tra,"O3_ste",itau_w,zx_tmp_2d,
	347	. iim*(jjm+1),ndex2d)
	348
	349	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,o3_prod(1,1),
	350	. zx_tmp_3d)
	351	CALL histwrite(nid_tra,"O3_prod",itau_w,zx_tmp_3d,
	352	. iim(jjm+1)klev,ndex3d)
	353
	354	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,o3_loss(1,1),
	355	. zx_tmp_3d)
	356	CALL histwrite(nid_tra,"O3_loss",itau_w,zx_tmp_3d,
	357	. iim(jjm+1)klev,ndex3d)
[624]	358	! Ajout Anne
	359	#ifdef INCA_AER
	360	! for sulfur cycle
[524]	361
[624]	362	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsoh(1,1),
	363	. zx_tmp_3d)
	364	CALL histwrite(nid_tra,"SO2_p_dmsoh",itau_w,zx_tmp_3d,
	365	. iim(jjm+1)klev,ndex)
	366
	367	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsno3(1,1),
	368	. zx_tmp_3d)
	369	CALL histwrite(nid_tra,"SO2_p_dmsno3",itau_w,zx_tmp_3d,
	370	. iim(jjm+1)klev,ndex)
	371
	372	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_h2soh(1,1),
	373	. zx_tmp_3d)
	374	CALL histwrite(nid_tra,"SO2_p_h2soh",itau_w,zx_tmp_3d,
	375	. iim(jjm+1)klev,ndex)
	376
	377	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsooh(1,1),
	378	. zx_tmp_3d)
	379	CALL histwrite(nid_tra,"SO2_p_dmsooh",itau_w,zx_tmp_3d,
	380	. iim(jjm+1)klev,ndex)
	381
	382	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,DMSO_p_dmsoh(1,1),
	383	. zx_tmp_3d)
	384	CALL histwrite(nid_tra,"DMSO_p_dmsoh",itau_w,zx_tmp_3d,
	385	. iim(jjm+1)klev,ndex)
	386
	387	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASMSAM_p_dmsooh(1,1),
	388	. zx_tmp_3d)
	389	CALL histwrite(nid_tra,"ASMSAM_p_dmsooh",itau_w,zx_tmp_3d,
	390	. iim(jjm+1)klev,ndex)
	391
	392	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2oh(1,1),
	393	. zx_tmp_3d)
	394	CALL histwrite(nid_tra,"ASSO4M_p_so2oh",itau_w,zx_tmp_3d,
	395	. iim(jjm+1)klev,ndex)
	396	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2h2o2(1,1),
	397	. zx_tmp_3d)
	398	CALL histwrite(nid_tra,"ASSO4M_p_so2h2o2",itau_w,zx_tmp_3d,
	399	. iim(jjm+1)klev,ndex)
	400
	401	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2o3(1,1),
	402	. zx_tmp_3d)
	403	CALL histwrite(nid_tra,"ASSO4M_p_so2o3",itau_w,zx_tmp_3d,
	404	. iim(jjm+1)klev,ndex)
	405
	406	c closing the sulfur budget
	407	CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_so2(1,1), zx_tmp_3d)
	408	CALL histwrite(nid_tra,"Wet3D_SO2",itau_w,zx_tmp_3d,
	409	. iim(jjm+1)klev,ndex)
	410	CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_dms(1,1), zx_tmp_3d)
	411	CALL histwrite(nid_tra,"Wet3D_DMS",itau_w,zx_tmp_3d,
	412	. iim(jjm+1)klev,ndex)
	413	CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_hno3(1,1), zx_tmp_3d)
	414	CALL histwrite(nid_tra,"Wet3D_HNO3",itau_w,zx_tmp_3d,
	415	. iim(jjm+1)klev,ndex)
	416	CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_h2o2(1,1), zx_tmp_3d)
	417	CALL histwrite(nid_tra,"Wet3D_H2O2",itau_w,zx_tmp_3d,
	418	. iim(jjm+1)klev,ndex)
	419
	420
	421	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,PH_HIST(1,1),
	422	. zx_tmp_3d)
	423	CALL histwrite(nid_tra,"PH_HIST",itau_w,zx_tmp_3d,
	424	. iim(jjm+1)klev,ndex)
	425	#endif
	426
[673]	427	#ifdef INCA_CH4
	428	DO it=1, grpcnt
	429
	430	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
	431	zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
	432	CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
	433	. iim(jjm+1)klev,ndex3d)
	434	ENDDO
	435
	436	#endif
	437
	438
[624]	439	#ifdef INCA_NMHC
	440	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1),
	441	. zx_tmp_3d)
	442	CALL histwrite(nid_tra,"CO2_basprod",itau_w,zx_tmp_3d,
	443	. iim(jjm+1)klev,ndex)
	444
	445	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_nmhcprod(1,1),
	446	. zx_tmp_3d)
	447	CALL histwrite(nid_tra,"CO2_nmhcprod",itau_w,zx_tmp_3d,
	448	. iim(jjm+1)klev,ndex)
	449
	450	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_radicalprod(1,1),
	451	. zx_tmp_3d)
	452	CALL histwrite(nid_tra,"CO2_radicalprod",itau_w,zx_tmp_3d,
	453	. iim(jjm+1)klev,ndex)
	454
	455	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_prod(1,1),
	456	. zx_tmp_3d)
	457	CALL histwrite(nid_tra,"HNO3_prod",itau_w,zx_tmp_3d,
	458	. iim(jjm+1)klev,ndex)
	459
	460	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_loss(1,1),
	461	. zx_tmp_3d)
	462	CALL histwrite(nid_tra,"HNO3_loss",itau_w,zx_tmp_3d,
	463	. iim(jjm+1)klev,ndex)
	464
	465	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_prod(1,1),
	466	. zx_tmp_3d)
	467	CALL histwrite(nid_tra,"CO_prod",itau_w,zx_tmp_3d,
	468	. iim(jjm+1)klev,ndex)
	469
	470	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_loss(1,1),
	471	. zx_tmp_3d)
	472	CALL histwrite(nid_tra,"CO_loss",itau_w,zx_tmp_3d,
	473	. iim(jjm+1)klev,ndex)
	474
[673]	475	DO it=1, grpcnt
	476
	477	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
	478	zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
	479	CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
	480	. iim(jjm+1)klev,ndex3d)
	481	ENDDO
	482
	483
[624]	484	#endif
	485
	486	! Fin ajout Anne
	487
[524]	488	! ... Special section for daytime averaging
	489	! CALL gr_fi_ecrit(klev,klon,iim,jjm+1,day_cnt(1,1),
	490	! . zx_tmp_3d)
	491	! CALL histwrite(nid_tra,"day_cnt",itau_w,zx_tmp_3d,
	492	! . iim(jjm+1)klev,ndex3d)
	493	! CALL gr_fi_ecrit(klev,klon,iim,jjm+1,no_daytime(1,1),
	494	! . zx_tmp_3d)
	495	! CALL histwrite(nid_tra,"NO_day",itau_w,zx_tmp_3d,
	496	! . iim(jjm+1)klev,ndex3d)
	497
[616]	498	!#endif
[524]	499	#else
	500
	501	C abder
	502	CALL gr_fi_ecrit(1,klon,iim,jjm+1,yu1,zx_tmp_2d)
	503	CALL histwrite(nid_tra,"pyu1",itau_w,zx_tmp_2d,
	504	s iim*(jjm+1),ndex2d)
	505
	506	CALL gr_fi_ecrit(1,klon,iim,jjm+1,yv1,zx_tmp_2d)
	507	CALL histwrite(nid_tra,"pyv1",itau_w,zx_tmp_2d,
	508	s iim*(jjm+1),ndex2d)
	509
	510	CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol1,zx_tmp_2d)
	511	CALL histwrite(nid_tra,"ftsol1",itau_w,zx_tmp_2d,
	512	s iim*(jjm+1),ndex2d)
	513
	514	CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol2,zx_tmp_2d)
	515	CALL histwrite(nid_tra,"ftsol2",itau_w,zx_tmp_2d,
	516	s iim*(jjm+1),ndex2d)
	517
	518	CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol3,zx_tmp_2d)
	519	CALL histwrite(nid_tra,"ftsol3",itau_w,zx_tmp_2d,
	520	s iim*(jjm+1),ndex2d)
	521
	522	CALL gr_fi_ecrit(1,klon,iim,jjm+1,pftsol4,zx_tmp_2d)
	523	CALL histwrite(nid_tra,"ftsol4",itau_w,zx_tmp_2d,
	524	s iim*(jjm+1),ndex2d)
	525
	526	CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf1,zx_tmp_2d)
	527	CALL histwrite(nid_tra,"psrf1",itau_w,zx_tmp_2d,
	528	s iim*(jjm+1),ndex2d)
	529
	530	CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf2,zx_tmp_2d)
	531	CALL histwrite(nid_tra,"psrf2",itau_w,zx_tmp_2d,
	532	s iim*(jjm+1),ndex2d)
	533
	534	CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf3,zx_tmp_2d)
	535	CALL histwrite(nid_tra,"psrf3",itau_w,zx_tmp_2d,
	536	s iim*(jjm+1),ndex2d)
	537
	538	CALL gr_fi_ecrit(1,klon,iim,jjm+1,ppsrf4,zx_tmp_2d)
	539	CALL histwrite(nid_tra,"psrf4",itau_w,zx_tmp_2d,
	540	s iim*(jjm+1),ndex2d)
	541	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pplay,zx_tmp_3d)
	542	CALL histwrite(nid_tra,"pplay",itau_w,zx_tmp_3d,
	543	. iim(jjm+1)klev,ndex3d)
	544
	545	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,t_seri,zx_tmp_3d)
	546	CALL histwrite(nid_tra,"t",itau_w,zx_tmp_3d,
	547	. iim(jjm+1)klev,ndex3d)
	548	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfu,zx_tmp_3d)
	549	CALL histwrite(nid_tra,"mfu",itau_w,zx_tmp_3d,
	550	. iim(jjm+1)klev,ndex3d)
	551	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pmfd,zx_tmp_3d)
	552	CALL histwrite(nid_tra,"mfd",itau_w,zx_tmp_3d,
	553	. iim(jjm+1)klev,ndex3d)
	554	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_u,zx_tmp_3d)
	555	CALL histwrite(nid_tra,"en_u",itau_w,zx_tmp_3d,
	556	. iim(jjm+1)klev,ndex3d)
	557	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pen_d,zx_tmp_3d)
	558	CALL histwrite(nid_tra,"en_d",itau_w,zx_tmp_3d,
	559	. iim(jjm+1)klev,ndex3d)
	560	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_d,zx_tmp_3d)
	561	CALL histwrite(nid_tra,"de_d",itau_w,zx_tmp_3d,
	562	. iim(jjm+1)klev,ndex3d)
	563	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pde_u,zx_tmp_3d)
	564	CALL histwrite(nid_tra,"de_u",itau_w,zx_tmp_3d,
	565	. iim(jjm+1)klev,ndex3d)
	566	CALL gr_fi_ecrit(klev,klon,iim,jjm+1,coefh,zx_tmp_3d)
	567	CALL histwrite(nid_tra,"coefh",itau_w,zx_tmp_3d,
	568	. iim(jjm+1)klev,ndex3d)
	569
	570
	571	c abder
	572	#endif
	573
	574	if (ok_sync) then
	575	call histsync(nid_tra)
	576	#ifdef INCA_AER
	577	call histsync(nid_tra2)
	578	call histsync(nid_tra3)
	579	#endif
	580	endif
	581
	582

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

Download in other formats: