Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

write_histday.h

Last change on this file was 506, checked in by lmdzadmin, 21 years ago

IM: ajout lev_hist(hf,day,mth) niveaux output definis dans physiq.def et

variables ENSEMBLES + rearrangement

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

File size: 17.4 KB

Rev	Line
[416]	1	IF (ok_journe) THEN
	2	c
	3	ndex2d = 0
	4	ndex3d = 0
	5	c
	6	c Champs 2D:
	7	c
	8	zsto = dtime
	9	zout = dtime * FLOAT(ecrit_day)
	10	itau_w = itau_phy + itap
	11
	12	i = NINT(zout/zsto)
	13	CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d)
[506]	14	c
	15	IF(lev_histday.GE.1) THEN
	16	c
[416]	17	CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	18	varname = 'phis'
	19	vartitle= 'Surface geop. height'
	20	varunits= '-'
	21	c
	22	i = NINT(zout/zsto)
	23	CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d)
	24	CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	25	varname = 'aire'
	26	vartitle= 'Grid area'
	27	varunits= '-'
[506]	28	c
	29	DO i=1, klon
	30	zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
	31	ENDDO
	32	c
	33	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	34	CALL histwrite(nid_day,"contfracATM",itau_w,
	35	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	36	c
	37	CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d)
	38	CALL histwrite(nid_day,"contfracOR",itau_w,
	39	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	40	c
[416]	41	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d)
	42	CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	43	C
	44	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	45	CALL histwrite(nid_day,"t2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	46	c
[506]	47	c En attendant un eventuel debugage.
	48
	49	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	50	CALL histwrite(nid_day,"t2m_min",itau_w,zx_tmp_2d,
	51	. iim*jjmp1,ndex2d)
	52	c
	53	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d)
	54	CALL histwrite(nid_day,"t2m_max",itau_w,zx_tmp_2d,
	55	. iim*jjmp1,ndex2d)
	56	c
	57	DO i = 1, klon
	58	zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
	59	ENDDO
	60	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	61	CALL histwrite(nid_day,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	62	c
	63	DO i = 1, klon
	64	zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
	65	ENDDO
	66	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	67	CALL histwrite(nid_day,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	68	c
	69	CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_lsc,zx_tmp_2d)
	70	CALL histwrite(nid_day,"snowl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	71	c
	72	IF(1.EQ.0) THEN
	73	CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_con,zx_tmp_2d)
	74	CALL histwrite(nid_day,"snowc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	75	ENDIF
	76	c
	77	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d)
	78	CALL histwrite(nid_day,"flat",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	79	c
	80	DO i = 1, klon
	81	zx_tmp_fi2d(i) = pctsrf(i,is_sic)
	82	ENDDO
	83	c
	84	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	85	CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	86	c
[416]	87	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zq2m,zx_tmp_2d)
	88	CALL histwrite(nid_day,"q2m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	89	c
	90	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zu10m,zx_tmp_2d)
	91	CALL histwrite(nid_day,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,
	92	. ndex2d)
	93	c
	94	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zv10m,zx_tmp_2d)
	95	CALL histwrite(nid_day,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,
	96	. ndex2d)
	97	c
[506]	98	DO i=1, klon
	99	zx_tmp_fi2d(i)=SQRT(zu10m(i)zu10m(i)+zv10m(i)zv10m(i))
	100	ENDDO
	101	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	102	CALL histwrite(nid_day,"wind10m",itau_w,zx_tmp_2d,
	103	$ iim*jjmp1,ndex2d)
	104	c
	105	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	106	CALL histwrite(nid_day,"wind10max",itau_w,zx_tmp_2d,
	107	. iim*jjmp1,ndex2d)
	108	c
	109	DO i=1, klon
[416]	110	zx_tmp_fi2d(i) = paprs(i,1)
	111	ENDDO
	112	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	113	CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	114	c
[506]	115	DO i=1, klon
[416]	116	zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i))
	117	ENDDO
	118	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	119	CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	120	c
	121	CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d)
[506]	122	CALL histwrite(nid_day,"snowf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[416]	123	c
	124	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
	125	CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1,
	126	. ndex2d)
	127	c
	128	CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d)
	129	CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	130	c
	131	CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d)
	132	CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	133	c
	134	CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d)
	135	CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	136	c
	137	CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d)
	138	CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	139	c
	140	CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d)
	141	CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	142	c
[506]	143	CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d)
	144	CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	145	c
	146	CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d)
	147	CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	148	c
	149	CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d)
	150	CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	151	c
	152	CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d)
	153	CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	154	c
	155	CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d)
	156	CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	157	c
	158	CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d)
	159	CALL histwrite(nid_day,"prw",itau_w,zx_tmp_2d,
	160	. iim*jjmp1,ndex2d)
	161
	162	c Ecriture de champs dynamiques sur des niveaux de pression
	163	DO k=1, nlevENS
	164	c
	165	CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d)
	166	CALL histwrite(nid_day,"u"//clev(k),itau_w,zx_tmp_2d,
	167	$ iim*jjmp1,ndex2d)
	168	c
	169	CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d)
	170	CALL histwrite(nid_day,"v"//clev(k),itau_w,zx_tmp_2d,
	171	$ iim*jjmp1,ndex2d)
	172	c
	173	ENDDO !nlevENS
	174	c
	175	DO k=1, nlevENS
	176	IF(clev(k).EQ."500") THEN
	177	c
	178	CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlev(:,indENS(k)),zx_tmp_2d)
	179	CALL histwrite(nid_day,"w500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	180	c
	181	c DO i=1, klon
	182	c print*,' ind, phi 500hPa',i,indENS(k),philev(:,indENS(k))
	183	c print*,' ind, phi 500hPa',i,indENS(k)
	184	c ENDDO
	185	c
	186	cBAD CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,indENS(k)),zx_tmp_2d)
	187	CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d)
	188	CALL histwrite(nid_day,"phi500",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	189	c
	190	ENDIF !clev(k).EQ.500
	191	ENDDO !k=1, nlevENS
	192	c
	193	CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d)
	194	CALL histwrite(nid_day,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	195
	196	CALL gr_fi_ecrit(1,klon,iim,jjmp1,cape,zx_tmp_2d)
	197	CALL histwrite(nid_day,"cape_max",itau_w,zx_tmp_2d,
	198	. iim*jjmp1,ndex2d)
	199	c
[416]	200	CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d)
	201	CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1,
	202	. ndex2d)
	203	c
[506]	204	DO i=1, klon
	205	zx_tmp_fi2d(i)=-1*sens(i)
	206	ENDDO
	207	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	208	CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	209	c
	210	zx_tmp_fi2d(1 : klon) = swdn(1 : klon, 1)
	211
	212	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	213	CALL histwrite(nid_day, "SWdnSFC",itau_w,zx_tmp_2d,
	214	. iim*jjmp1,ndex2d)
	215	c
	216	ENDIF !lev_histday.GE.1
	217	c
	218	IF(lev_histday.GE.2) THEN
	219	c
[416]	220	CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d)
	221	CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	222	c
	223	CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d)
	224	CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
	225	c
[506]	226	ENDIF !lev_histday.GE.2
[416]	227	c
[506]	228	IF(lev_histday.GE.3) THEN
	229	c=================================================================
	230	c ECRITURE DES CHAMPS 3D
	231	c=================================================================
	232
[433]	233	c
[506]	234	c Champs 3D:
[433]	235	c
[506]	236	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d)
	237	CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d,
	238	. iimjjmp1klev,ndex3d)
[433]	239	c
[506]	240	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d)
	241	CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d,
	242	. iimjjmp1klev,ndex3d)
[433]	243	c
[506]	244	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d)
	245	CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d,
	246	. iimjjmp1klev,ndex3d)
[433]	247	c
[506]	248	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d)
	249	CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d,
	250	. iimjjmp1klev,ndex3d)
[433]	251	c
[506]	252	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d)
	253	CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d,
	254	. iimjjmp1klev,ndex3d)
	255	c
	256	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d)
	257	CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d,
	258	. iimjjmp1klev,ndex3d)
	259	c
	260	CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d)
	261	CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d,
	262	. iimjjmp1klev,ndex3d)
	263	c
	264	ENDIF !lev_histday.GE.3
	265	c=================================================================
	266	c FIN ECRITURE DES CHAMPS 3D
	267	c=================================================================
	268	IF(lev_histday.ge.4) THEN
	269	c=================================================================
	270	c
	271	c ECRITURE DES CHAMPS SUR LES SOUS SURFACES
	272	c
	273	c=================================================================
	274	zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 )
[433]	275	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	276	CALL histwrite(nid_day, "SWupTOA",itau_w,zx_tmp_2d,
	277	. iim*jjmp1,ndex2d)
	278	c
[467]	279	zx_tmp_fi2d(1 : klon) = swup( 1 : klon, 1 )
[433]	280	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	281	CALL histwrite(nid_day, "SWupSFC",itau_w,zx_tmp_2d,
	282	. iim*jjmp1,ndex2d)
	283	c
[506]	284	zx_tmp_fi2d(1 : klon) = swdn( 1 : klon, klevp1 )
[433]	285	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	286	CALL histwrite(nid_day, "SWdnTOA",itau_w,zx_tmp_2d,
	287	. iim*jjmp1,ndex2d)
	288	c
[506]	289	zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, klevp1 )
[433]	290	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
[467]	291	CALL histwrite(nid_day, "SWupTOAclr",itau_w,zx_tmp_2d,
	292	. iim*jjmp1,ndex2d)
[506]	293	c
[467]	294	zx_tmp_fi2d(1 : klon) = swup0( 1 : klon, 1 )
	295	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	296	CALL histwrite(nid_day, "SWupSFCclr",itau_w,zx_tmp_2d,
	297	. iim*jjmp1,ndex2d)
[506]	298	c
	299	zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, klevp1 )
[467]	300	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	301	CALL histwrite(nid_day, "SWdnTOAclr",itau_w,zx_tmp_2d,
	302	. iim*jjmp1,ndex2d)
[506]	303	c
[467]	304	zx_tmp_fi2d(1 : klon) = swdn0( 1 : klon, 1 )
	305	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	306	CALL histwrite(nid_day, "SWdnSFCclr",itau_w,zx_tmp_2d,
	307	. iim*jjmp1,ndex2d)
[433]	308	c
	309	zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter)
	310	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d)
	311	CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[506]	312	c
[433]	313	zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic)
	314	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	315	CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[506]	316	c
[433]	317	zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce)
	318	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	319	CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[506]	320	c
[433]	321	zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic)
	322	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
	323	CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)
[506]	324	c
[433]	325	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	326	CALL histwrite(nid_day,"t2mter",itau_w,zx_tmp_2d,iim*jjmp1,
	327	. ndex2d)
	328	c
	329	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
	330	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	331	CALL histwrite(nid_day,"t2mter_min",itau_w,zx_tmp_2d,
	332	. iim*jjmp1,ndex2d)
	333	c
	334	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_ter)
	335	CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
	336	CALL histwrite(nid_day,"t2mter_max",itau_w,zx_tmp_2d,
	337	. iim*jjmp1,ndex2d)
	338	c
	339	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_lic)
	340	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	341	CALL histwrite(nid_day,"t2mlic",itau_w,zx_tmp_2d,iim*jjmp1,
	342	. ndex2d)
	343	c
	344	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_oce)
	345	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	346	CALL histwrite(nid_day,"t2moce",itau_w,zx_tmp_2d,iim*jjmp1,
	347	. ndex2d)
	348	c
	349	zx_tmp_fi2d(1 : klon) = t2m(1 : klon, is_sic)
	350	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	351	CALL histwrite(nid_day,"t2msic",itau_w,zx_tmp_2d,iim*jjmp1,
	352	. ndex2d)
	353	c
	354	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_ter)
	355	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	356	CALL histwrite(nid_day,"u10mter",itau_w,zx_tmp_2d,iim*jjmp1,
	357	. ndex2d)
	358	c
	359	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_ter)
	360	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	361	CALL histwrite(nid_day,"v10mter",itau_w,zx_tmp_2d,iim*jjmp1,
	362	. ndex2d)
	363	c
	364	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_lic)
	365	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	366	CALL histwrite(nid_day,"u10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
	367	. ndex2d)
	368	c
	369	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_lic)
	370	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	371	CALL histwrite(nid_day,"v10mlic",itau_w,zx_tmp_2d,iim*jjmp1,
	372	. ndex2d)
	373	c
	374	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_oce)
	375	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	376	CALL histwrite(nid_day,"u10moce",itau_w,zx_tmp_2d,iim*jjmp1,
	377	. ndex2d)
	378	c
	379	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_oce)
	380	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	381	CALL histwrite(nid_day,"v10moce",itau_w,zx_tmp_2d,iim*jjmp1,
	382	. ndex2d)
	383	c
	384	zx_tmp_fi2d(1 : klon) = u10m(1 : klon, is_sic)
	385	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	386	CALL histwrite(nid_day,"u10msic",itau_w,zx_tmp_2d,iim*jjmp1,
	387	. ndex2d)
	388	C
	389	zx_tmp_fi2d(1 : klon) = v10m(1 : klon, is_sic)
	390	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
	391	CALL histwrite(nid_day,"v10msic",itau_w,zx_tmp_2d,iim*jjmp1,
	392	. ndex2d)
	393	C
[416]	394	DO nsrf = 1, nbsrf
[506]	395	C
[416]	396	zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
	397	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	398	CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w,
	399	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	400	C
	401	zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
	402	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	403	CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w,
	404	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	405	C
	406	zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
	407	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	408	CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w,
	409	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	410	C
	411	zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
	412	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	413	CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w,
	414	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	415	C
	416	zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
	417	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	418	CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w,
	419	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	420	C
	421	zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
	422	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	423	CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w,
	424	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	425	C
	426	zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf)
	427	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	428	CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w,
	429	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	430	C
	431	zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf)
	432	CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d)
	433	CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w,
	434	$ zx_tmp_2d,iim*jjmp1,ndex2d)
	435	C
	436	END DO
[433]	437	c=================================================================
	438	c FIN ECRITURE DES CHAMPS SUR LES SOUS SURFACES
	439	c=================================================================
[506]	440	ENDIF !lev_histday.GE.4
[416]	441	c
[433]	442	c=================================================================
	443	c=================================================================
	444	c=================================================================
[506]	445	c
[416]	446	if (ok_sync) then
	447	call histsync(nid_day)
	448	endif
[433]	449
[416]	450	ENDIF

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

Context Navigation

source: LMDZ.3.3/branches/rel-LF/libf/phylmd/write_histday.h

Download in other formats: