Context Navigation

phys_output_mod.F90 @ 1181

Last change on this file since 1181 was 1176, checked in by Laurent Fairhead, 16 years ago
Modif pour la compilation avec fcm sur Vargas SD Corrections relatives au controle de histins et de l'utilisation des options INST(X) pour l'appel au routines de IOIPSL FH Mise en place de diagnostics sur les critères CFL pour l'advection de traceurs. FH Controle dans les .def de la dépendance verticale de l'efficacite de la diffusion. Actif pour le moment uniquement avec ok_strato=y et llm=39. FH LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 60.2 KB

Rev	Line
[1176]	1	!
	2	! $Id: phys_output_mod.F90 1176 2009-06-11 08:54:10Z jghattas $
	3	!
[907]	4	! Abderrahmane 12 2007
	5	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	6	!!! Ecreture des Sorties du modele dans les fichiers Netcdf :
	7	! histmth.nc : moyennes mensuelles
	8	! histday.nc : moyennes journalieres
	9	! histhf.nc : moyennes toutes les 3 heures
	10	! histins.nc : valeurs instantanees
	11	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	12
	13	MODULE phys_output_mod
	14
	15	IMPLICIT NONE
	16
[1083]	17	private histdef2d, histdef3d, conf_physoutputs
[1065]	18
[1083]	19
[1054]	20	integer, parameter :: nfiles = 5
[929]	21	logical, dimension(nfiles), save :: clef_files
	22	integer, dimension(nfiles), save :: lev_files
	23	integer, dimension(nfiles), save :: nid_files
[1176]	24	!$OMP THREADPRIVATE(clef_files, lev_files,nid_files)
	25
[1065]	26	integer, dimension(nfiles), private, save :: nhorim, nvertm
[1176]	27	real, dimension(nfiles), private, save :: zoutm
	28	real, private, save :: zdtime
[1065]	29	CHARACTER(len=20), dimension(nfiles), private, save :: type_ecri
[1176]	30	!$OMP THREADPRIVATE(nhorim, nvertm, zoutm,zdtime,type_ecri)
[929]	31
[907]	32	! integer, save :: nid_hf3d
	33
	34	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	35	!! Definition pour chaque variable du niveau d ecriture dans chaque fichier
[1083]	36	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/ histmth, histday, histhf, histins /),'!!!!!!!!!!!!
[907]	37	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	38
[1065]	39	integer, private:: levmin(nfiles) = 1
	40	integer, private:: levmax(nfiles)
	41
[1083]	42	TYPE ctrl_out
	43	integer,dimension(5) :: flag
	44	character(len=20) :: name
	45	END TYPE ctrl_out
	46
	47
[907]	48	!!! 1D
[1083]	49	type(ctrl_out) :: o_phis = ctrl_out((/ 1, 1, 10, 1, 1 /), 'phis')
	50	type(ctrl_out) :: o_aire = ctrl_out((/ 1, 1, 10, 1, 1 /),'aire')
	51	type(ctrl_out) :: o_contfracATM = ctrl_out((/ 10, 1, 1, 10, 10 /),'contfracATM')
	52	type(ctrl_out) :: o_contfracOR = ctrl_out((/ 10, 1, 1, 10, 10 /),'contfracOR')
	53	type(ctrl_out) :: o_aireTER = ctrl_out((/ 10, 10, 1, 10, 10 /),'aireTER')
[907]	54
	55	!!! 2D
[1083]	56	type(ctrl_out) :: o_flat = ctrl_out((/ 10, 1, 10, 10, 1 /),'flat')
	57	type(ctrl_out) :: o_slp = ctrl_out((/ 1, 1, 1, 10, 1 /),'slp')
	58	type(ctrl_out) :: o_tsol = ctrl_out((/ 1, 1, 1, 1, 1 /),'tsol')
	59	type(ctrl_out) :: o_t2m = ctrl_out((/ 1, 1, 1, 1, 1 /),'t2m')
	60	type(ctrl_out) :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_min')
	61	type(ctrl_out) :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10 /),'t2m_max')
	62	type(ctrl_out),dimension(4) :: o_t2m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_ter'), &
	63	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_lic'), &
	64	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_oce'), &
	65	ctrl_out((/ 10, 4, 10, 10, 10 /),'t2m_sic') /)
[907]	66
[1083]	67	type(ctrl_out) :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10 /),'wind10m')
	68	type(ctrl_out) :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10 /),'wind10max')
	69	type(ctrl_out) :: o_sicf = ctrl_out((/ 1, 1, 10, 10, 10 /),'sicf')
	70	type(ctrl_out) :: o_q2m = ctrl_out((/ 1, 1, 1, 1, 1 /),'q2m')
	71	type(ctrl_out) :: o_u10m = ctrl_out((/ 1, 1, 1, 1, 1 /),'u10m')
	72	type(ctrl_out) :: o_v10m = ctrl_out((/ 1, 1, 1, 1, 1 /),'v10m')
	73	type(ctrl_out) :: o_psol = ctrl_out((/ 1, 1, 1, 1, 1 /),'psol')
	74	type(ctrl_out) :: o_qsurf = ctrl_out((/ 1, 10, 10, 10, 10 /),'qsurf')
[907]	75
[1083]	76	type(ctrl_out),dimension(4) :: o_u10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_ter'), &
	77	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_lic'), &
	78	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_oce'), &
	79	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10m_sic') /)
[907]	80
[1083]	81	type(ctrl_out),dimension(4) :: o_v10m_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_ter'), &
	82	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_lic'), &
	83	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_oce'), &
	84	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10m_sic') /)
[907]	85
[1083]	86	type(ctrl_out) :: o_qsol = ctrl_out((/ 1, 10, 10, 1, 1 /),'qsol')
[907]	87
[1083]	88	type(ctrl_out) :: o_ndayrain = ctrl_out((/ 1, 10, 10, 10, 10 /),'ndayrain')
	89	type(ctrl_out) :: o_precip = ctrl_out((/ 1, 1, 1, 1, 1 /),'precip')
	90	type(ctrl_out) :: o_plul = ctrl_out((/ 1, 1, 1, 1, 10 /),'plul')
	91
	92	type(ctrl_out) :: o_pluc = ctrl_out((/ 1, 1, 1, 1, 10 /),'pluc')
	93	type(ctrl_out) :: o_snow = ctrl_out((/ 1, 1, 10, 1, 10 /),'snow')
	94	type(ctrl_out) :: o_evap = ctrl_out((/ 1, 1, 10, 1, 10 /),'evap')
	95	type(ctrl_out) :: o_tops = ctrl_out((/ 1, 1, 10, 10, 10 /),'tops')
	96	type(ctrl_out) :: o_tops0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'tops0')
	97	type(ctrl_out) :: o_topl = ctrl_out((/ 1, 1, 10, 1, 10 /),'topl')
	98	type(ctrl_out) :: o_topl0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'topl0')
	99	type(ctrl_out) :: o_SWupTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupTOA')
	100	type(ctrl_out) :: o_SWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupTOAclr')
	101	type(ctrl_out) :: o_SWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnTOA')
	102	type(ctrl_out) :: o_SWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnTOAclr')
	103	type(ctrl_out) :: o_SWup200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'SWup200')
	104	type(ctrl_out) :: o_SWup200clr = ctrl_out((/ 10, 1, 10, 10, 10 /),'SWup200clr')
	105	type(ctrl_out) :: o_SWdn200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'SWdn200')
	106	type(ctrl_out) :: o_SWdn200clr = ctrl_out((/ 10, 1, 10, 10, 10 /),'SWdn200clr')
	107
[907]	108	! arajouter
[1083]	109	! type(ctrl_out) :: o_LWupTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupTOA')
	110	! type(ctrl_out) :: o_LWupTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupTOAclr')
	111	! type(ctrl_out) :: o_LWdnTOA = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnTOA')
	112	! type(ctrl_out) :: o_LWdnTOAclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnTOAclr')
[907]	113
[1083]	114	type(ctrl_out) :: o_LWup200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWup200')
	115	type(ctrl_out) :: o_LWup200clr = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWup200clr')
	116	type(ctrl_out) :: o_LWdn200 = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWdn200')
	117	type(ctrl_out) :: o_LWdn200clr = ctrl_out((/ 1, 10, 10, 10, 10 /),'LWdn200clr')
	118	type(ctrl_out) :: o_sols = ctrl_out((/ 1, 1, 10, 1, 10 /),'sols')
	119	type(ctrl_out) :: o_sols0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'sols0')
	120	type(ctrl_out) :: o_soll = ctrl_out((/ 1, 1, 10, 1, 10 /),'soll')
	121	type(ctrl_out) :: o_soll0 = ctrl_out((/ 1, 5, 10, 10, 10 /),'soll0')
	122	type(ctrl_out) :: o_radsol = ctrl_out((/ 1, 1, 10, 10, 10 /),'radsol')
	123	type(ctrl_out) :: o_SWupSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupSFC')
	124	type(ctrl_out) :: o_SWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWupSFCclr')
	125	type(ctrl_out) :: o_SWdnSFC = ctrl_out((/ 1, 1, 10, 10, 10 /),'SWdnSFC')
	126	type(ctrl_out) :: o_SWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'SWdnSFCclr')
	127	type(ctrl_out) :: o_LWupSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupSFC')
	128	type(ctrl_out) :: o_LWupSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWupSFCclr')
	129	type(ctrl_out) :: o_LWdnSFC = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnSFC')
	130	type(ctrl_out) :: o_LWdnSFCclr = ctrl_out((/ 1, 4, 10, 10, 10 /),'LWdnSFCclr')
	131	type(ctrl_out) :: o_bils = ctrl_out((/ 1, 2, 10, 1, 10 /),'bils')
	132	type(ctrl_out) :: o_sens = ctrl_out((/ 1, 1, 10, 1, 1 /),'sens')
	133	type(ctrl_out) :: o_fder = ctrl_out((/ 1, 2, 10, 1, 10 /),'fder')
	134	type(ctrl_out) :: o_ffonte = ctrl_out((/ 1, 10, 10, 10, 10 /),'ffonte')
	135	type(ctrl_out) :: o_fqcalving = ctrl_out((/ 1, 10, 10, 10, 10 /),'fqcalving')
	136	type(ctrl_out) :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10 /),'fqfonte')
[907]	137
[1083]	138	type(ctrl_out),dimension(4) :: o_taux_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_ter'), &
	139	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_lic'), &
	140	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_oce'), &
	141	ctrl_out((/ 1, 4, 10, 1, 10 /),'taux_sic') /)
[907]	142
[1083]	143	type(ctrl_out),dimension(4) :: o_tauy_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_ter'), &
	144	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_lic'), &
	145	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_oce'), &
	146	ctrl_out((/ 1, 4, 10, 1, 10 /),'tauy_sic') /)
[907]	147
[1083]	148
	149	type(ctrl_out),dimension(4) :: o_pourc_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_ter'), &
	150	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_lic'), &
	151	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_oce'), &
	152	ctrl_out((/ 1, 4, 10, 1, 10 /),'pourc_sic') /)
	153
	154	type(ctrl_out),dimension(4) :: o_fract_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_ter'), &
	155	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_lic'), &
	156	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_oce'), &
	157	ctrl_out((/ 1, 4, 10, 1, 10 /),'fract_sic') /)
	158
	159	type(ctrl_out),dimension(4) :: o_tsol_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_ter'), &
	160	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_lic'), &
	161	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_oce'), &
	162	ctrl_out((/ 1, 4, 10, 1, 10 /),'tsol_sic') /)
	163
	164	type(ctrl_out),dimension(4) :: o_sens_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_ter'), &
	165	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_lic'), &
	166	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_oce'), &
	167	ctrl_out((/ 1, 4, 10, 1, 10 /),'sens_sic') /)
	168
	169	type(ctrl_out),dimension(4) :: o_lat_srf = (/ ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_ter'), &
	170	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_lic'), &
	171	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_oce'), &
	172	ctrl_out((/ 1, 4, 10, 1, 10 /),'lat_sic') /)
	173
	174	type(ctrl_out),dimension(4) :: o_flw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_ter'), &
	175	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_lic'), &
	176	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_oce'), &
	177	ctrl_out((/ 1, 10, 10, 10, 10 /),'flw_sic') /)
	178
	179	type(ctrl_out),dimension(4) :: o_fsw_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_ter'), &
	180	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_lic'), &
	181	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_oce'), &
	182	ctrl_out((/ 1, 10, 10, 10, 10 /),'fsw_sic') /)
	183
	184	type(ctrl_out),dimension(4) :: o_wbils_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_ter'), &
	185	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_lic'), &
	186	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_oce'), &
	187	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbils_sic') /)
	188
	189	type(ctrl_out),dimension(4) :: o_wbilo_srf = (/ ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_ter'), &
	190	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_lic'), &
	191	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_oce'), &
	192	ctrl_out((/ 1, 10, 10, 10, 10 /),'wbilo_sic') /)
	193
	194
	195	type(ctrl_out) :: o_cdrm = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrm')
	196	type(ctrl_out) :: o_cdrh = ctrl_out((/ 1, 10, 10, 1, 10 /),'cdrh')
	197	type(ctrl_out) :: o_cldl = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldl')
	198	type(ctrl_out) :: o_cldm = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldm')
	199	type(ctrl_out) :: o_cldh = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldh')
	200	type(ctrl_out) :: o_cldt = ctrl_out((/ 1, 1, 2, 10, 10 /),'cldt')
	201	type(ctrl_out) :: o_cldq = ctrl_out((/ 1, 1, 10, 10, 10 /),'cldq')
	202	type(ctrl_out) :: o_lwp = ctrl_out((/ 1, 5, 10, 10, 10 /),'lwp')
	203	type(ctrl_out) :: o_iwp = ctrl_out((/ 1, 5, 10, 10, 10 /),'iwp')
	204	type(ctrl_out) :: o_ue = ctrl_out((/ 1, 10, 10, 10, 10 /),'ue')
	205	type(ctrl_out) :: o_ve = ctrl_out((/ 1, 10, 10, 10, 10 /),'ve')
	206	type(ctrl_out) :: o_uq = ctrl_out((/ 1, 10, 10, 10, 10 /),'uq')
	207	type(ctrl_out) :: o_vq = ctrl_out((/ 1, 10, 10, 10, 10 /),'vq')
[907]	208
[1083]	209	type(ctrl_out) :: o_cape = ctrl_out((/ 1, 10, 10, 10, 10 /),'cape')
	210	type(ctrl_out) :: o_pbase = ctrl_out((/ 1, 10, 10, 10, 10 /),'pbase')
	211	type(ctrl_out) :: o_ptop = ctrl_out((/ 1, 4, 10, 10, 10 /),'ptop')
	212	type(ctrl_out) :: o_fbase = ctrl_out((/ 1, 10, 10, 10, 10 /),'fbase')
	213	type(ctrl_out) :: o_prw = ctrl_out((/ 1, 1, 10, 10, 10 /),'prw')
[907]	214
[1083]	215	type(ctrl_out) :: o_s_pblh = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_pblh')
	216	type(ctrl_out) :: o_s_pblt = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_pblt')
	217	type(ctrl_out) :: o_s_lcl = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_lcl')
	218	type(ctrl_out) :: o_s_capCL = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_capCL')
	219	type(ctrl_out) :: o_s_oliqCL = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_oliqCL')
	220	type(ctrl_out) :: o_s_cteiCL = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_cteiCL')
	221	type(ctrl_out) :: o_s_therm = ctrl_out((/ 1, 10, 10, 1, 1 /),'s_therm')
	222	type(ctrl_out) :: o_s_trmb1 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb1')
	223	type(ctrl_out) :: o_s_trmb2 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb2')
	224	type(ctrl_out) :: o_s_trmb3 = ctrl_out((/ 1, 10, 10, 1, 10 /),'s_trmb3')
[907]	225
[1083]	226	type(ctrl_out) :: o_slab_bils = ctrl_out((/ 1, 1, 10, 10, 10 /),'slab_bils_oce')
[907]	227
[1083]	228	type(ctrl_out) :: o_ale_bl = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale_bl')
	229	type(ctrl_out) :: o_alp_bl = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp_bl')
	230	type(ctrl_out) :: o_ale_wk = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale_wk')
	231	type(ctrl_out) :: o_alp_wk = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp_wk')
[945]	232
[1083]	233	type(ctrl_out) :: o_ale = ctrl_out((/ 1, 1, 1, 1, 10 /),'ale')
	234	type(ctrl_out) :: o_alp = ctrl_out((/ 1, 1, 1, 1, 10 /),'alp')
	235	type(ctrl_out) :: o_cin = ctrl_out((/ 1, 1, 1, 1, 10 /),'cin')
	236	type(ctrl_out) :: o_wape = ctrl_out((/ 1, 1, 1, 1, 10 /),'wape')
[973]	237
	238
[907]	239	! Champs interpolles sur des niveaux de pression ??? a faire correctement
	240	! if=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
	241	! if=2 on ecrit w et ph 500 seulement au niv 1
	242	! et u v sur 850 700 500 200
	243	! if=3 on ecrit ph a 500 seulement au niv 1
	244	! on ecrit u v t q a 850 700 500 200 au niv 3
	245	! on ecrit ph a 500 au niv 3
	246
[1083]	247
	248	type(ctrl_out),dimension(4) :: o_uSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'u850'), &
	249	ctrl_out((/ 1, 1, 3, 10, 10 /),'u700'), &
	250	ctrl_out((/ 1, 1, 3, 10, 10 /),'u500'), &
	251	ctrl_out((/ 1, 1, 3, 10, 10 /),'u200') /)
[907]	252
[1083]	253	type(ctrl_out),dimension(4) :: o_vSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'v850'), &
	254	ctrl_out((/ 1, 1, 3, 10, 10 /),'v700'), &
	255	ctrl_out((/ 1, 1, 3, 10, 10 /),'v500'), &
	256	ctrl_out((/ 1, 1, 3, 10, 10 /),'v200') /)
[907]	257
[1083]	258	type(ctrl_out),dimension(4) :: o_wSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'w850'), &
	259	ctrl_out((/ 1, 1, 3, 10, 10 /),'w700'), &
	260	ctrl_out((/ 1, 1, 3, 10, 10 /),'w500'), &
	261	ctrl_out((/ 1, 1, 3, 10, 10 /),'w200') /)
[907]	262
[1083]	263	type(ctrl_out),dimension(4) :: o_tSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'t850'), &
	264	ctrl_out((/ 1, 1, 3, 10, 10 /),'t700'), &
	265	ctrl_out((/ 1, 1, 3, 10, 10 /),'t500'), &
	266	ctrl_out((/ 1, 1, 3, 10, 10 /),'t200') /)
[907]	267
[1083]	268	type(ctrl_out),dimension(4) :: o_qSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'q850'), &
	269	ctrl_out((/ 1, 1, 3, 10, 10 /),'q700'), &
	270	ctrl_out((/ 1, 1, 3, 10, 10 /),'q500'), &
	271	ctrl_out((/ 1, 1, 3, 10, 10 /),'q200') /)
[907]	272
[1083]	273	type(ctrl_out),dimension(4) :: o_phiSTDlevs = (/ ctrl_out((/ 1, 1, 3, 10, 10 /),'phi850'), &
	274	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi700'), &
	275	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi500'), &
	276	ctrl_out((/ 1, 1, 3, 10, 10 /),'phi200') /)
	277
	278
	279	type(ctrl_out) :: o_t_oce_sic = ctrl_out((/ 1, 10, 10, 10, 10 /),'t_oce_sic')
	280
	281	type(ctrl_out) :: o_weakinv = ctrl_out((/ 10, 1, 10, 10, 10 /),'weakinv')
	282	type(ctrl_out) :: o_dthmin = ctrl_out((/ 10, 1, 10, 10, 10 /),'dthmin')
	283	type(ctrl_out),dimension(4) :: o_u10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_ter'), &
[1123]	284	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_lic'), &
	285	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_oce'), &
	286	ctrl_out((/ 10, 4, 10, 10, 10 /),'u10_sic') /)
[1083]	287
	288	type(ctrl_out),dimension(4) :: o_v10_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_ter'), &
[1123]	289	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_lic'), &
	290	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_oce'), &
	291	ctrl_out((/ 10, 4, 10, 10, 10 /),'v10_sic') /)
	292
[1083]	293	type(ctrl_out) :: o_cldtau = ctrl_out((/ 10, 5, 10, 10, 10 /),'cldtau')
	294	type(ctrl_out) :: o_cldemi = ctrl_out((/ 10, 5, 10, 10, 10 /),'cldemi')
	295	type(ctrl_out) :: o_rh2m = ctrl_out((/ 10, 5, 10, 10, 10 /),'rh2m')
	296	type(ctrl_out) :: o_qsat2m = ctrl_out((/ 10, 5, 10, 10, 10 /),'qsat2m')
	297	type(ctrl_out) :: o_tpot = ctrl_out((/ 10, 5, 10, 10, 10 /),'tpot')
	298	type(ctrl_out) :: o_tpote = ctrl_out((/ 10, 5, 10, 10, 10 /),'tpote')
	299	type(ctrl_out) :: o_tke = ctrl_out((/ 4, 10, 10, 10, 10 /),'tke ')
	300	type(ctrl_out) :: o_tke_max = ctrl_out((/ 4, 10, 10, 10, 10 /),'tke_max')
[1123]	301
[1083]	302	type(ctrl_out),dimension(4) :: o_tke_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_ter'), &
[1123]	303	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_lic'), &
	304	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_oce'), &
	305	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_sic') /)
[1083]	306
[1123]	307	type(ctrl_out),dimension(4) :: o_tke_max_srf = (/ ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_ter'), &
	308	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_lic'), &
	309	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_oce'), &
	310	ctrl_out((/ 10, 4, 10, 10, 10 /),'tke_max_sic') /)
[1083]	311
	312	type(ctrl_out) :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10 /),'kz')
	313	type(ctrl_out) :: o_kz_max = ctrl_out((/ 4, 10, 10, 10, 10 /),'kz_max')
	314	type(ctrl_out) :: o_SWnetOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'SWnetOR')
	315	type(ctrl_out) :: o_SWdownOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'SWdownOR')
	316	type(ctrl_out) :: o_LWdownOR = ctrl_out((/ 10, 10, 2, 10, 10 /),'LWdownOR')
	317
	318	type(ctrl_out) :: o_snowl = ctrl_out((/ 10, 1, 10, 10, 10 /),'snowl')
	319	type(ctrl_out) :: o_cape_max = ctrl_out((/ 10, 1, 10, 10, 10 /),'cape_max')
	320	type(ctrl_out) :: o_solldown = ctrl_out((/ 10, 1, 10, 1, 10 /),'solldown')
	321
	322	type(ctrl_out) :: o_dtsvdfo = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfo')
	323	type(ctrl_out) :: o_dtsvdft = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdft')
	324	type(ctrl_out) :: o_dtsvdfg = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfg')
	325	type(ctrl_out) :: o_dtsvdfi = ctrl_out((/ 10, 10, 10, 1, 10 /),'dtsvdfi')
	326	type(ctrl_out) :: o_rugs = ctrl_out((/ 10, 10, 10, 1, 1 /),'rugs')
	327
[1095]	328	type(ctrl_out) :: o_topswad = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswad')
	329	type(ctrl_out) :: o_topswai = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswai')
	330	type(ctrl_out) :: o_solswad = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswad')
	331	type(ctrl_out) :: o_solswai = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswai')
[907]	332	!!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1083]	333	type(ctrl_out) :: o_lwcon = ctrl_out((/ 2, 5, 10, 10, 1 /),'lwcon')
	334	type(ctrl_out) :: o_iwcon = ctrl_out((/ 2, 5, 10, 10, 10 /),'iwcon')
	335	type(ctrl_out) :: o_temp = ctrl_out((/ 2, 3, 4, 1, 1 /),'temp')
	336	type(ctrl_out) :: o_theta = ctrl_out((/ 2, 3, 4, 1, 1 /),'theta')
	337	type(ctrl_out) :: o_ovap = ctrl_out((/ 2, 3, 4, 1, 1 /),'ovap')
[1123]	338	type(ctrl_out) :: o_ovapinit = ctrl_out((/ 2, 3, 4, 1, 1 /),'ovapinit')
[1083]	339	type(ctrl_out) :: o_wvapp = ctrl_out((/ 2, 10, 10, 10, 10 /),'wvapp')
	340	type(ctrl_out) :: o_geop = ctrl_out((/ 2, 3, 10, 1, 1 /),'geop')
	341	type(ctrl_out) :: o_vitu = ctrl_out((/ 2, 3, 4, 1, 1 /),'vitu')
	342	type(ctrl_out) :: o_vitv = ctrl_out((/ 2, 3, 4, 1, 1 /),'vitv')
	343	type(ctrl_out) :: o_vitw = ctrl_out((/ 2, 3, 10, 10, 1 /),'vitw')
	344	type(ctrl_out) :: o_pres = ctrl_out((/ 2, 3, 10, 1, 1 /),'pres')
	345	type(ctrl_out) :: o_rneb = ctrl_out((/ 2, 5, 10, 10, 1 /),'rneb')
	346	type(ctrl_out) :: o_rnebcon = ctrl_out((/ 2, 5, 10, 10, 1 /),'rnebcon')
	347	type(ctrl_out) :: o_rhum = ctrl_out((/ 2, 10, 10, 10, 10 /),'rhum')
	348	type(ctrl_out) :: o_ozone = ctrl_out((/ 2, 10, 10, 10, 10 /),'ozone')
	349	type(ctrl_out) :: o_upwd = ctrl_out((/ 2, 10, 10, 10, 10 /),'upwd')
	350	type(ctrl_out) :: o_dtphy = ctrl_out((/ 2, 10, 10, 10, 1 /),'dtphy')
	351	type(ctrl_out) :: o_dqphy = ctrl_out((/ 2, 10, 10, 10, 1 /),'dqphy')
	352	type(ctrl_out) :: o_pr_con_l = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_con_l')
	353	type(ctrl_out) :: o_pr_con_i = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_con_i')
	354	type(ctrl_out) :: o_pr_lsc_l = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_l')
	355	type(ctrl_out) :: o_pr_lsc_i = ctrl_out((/ 2, 10, 10, 10, 10 /),'pr_lsc_i')
[929]	356	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[907]	357
[1083]	358	type(ctrl_out),dimension(4) :: o_albe_srf = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_ter'), &
[1123]	359	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_lic'), &
	360	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_oce'), &
	361	ctrl_out((/ 3, 4, 10, 1, 10 /),'albe_sic') /)
[907]	362
[1083]	363	type(ctrl_out),dimension(4) :: o_ages_srf = (/ ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_ter'), &
[1123]	364	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_lic'), &
	365	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_oce'), &
	366	ctrl_out((/ 3, 10, 10, 10, 10 /),'ages_sic') /)
[907]	367
[1083]	368	type(ctrl_out),dimension(4) :: o_rugs_srf = (/ ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_ter'), &
[1123]	369	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_lic'), &
	370	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_oce'), &
	371	ctrl_out((/ 3, 4, 10, 1, 10 /),'rugs_sic') /)
[907]	372
[1083]	373	type(ctrl_out) :: o_albs = ctrl_out((/ 3, 10, 10, 1, 10 /),'albs')
	374	type(ctrl_out) :: o_albslw = ctrl_out((/ 3, 10, 10, 1, 10 /),'albslw')
	375
	376	type(ctrl_out) :: o_clwcon = ctrl_out((/ 4, 10, 10, 10, 10 /),'clwcon')
	377	type(ctrl_out) :: o_Ma = ctrl_out((/ 4, 10, 10, 10, 10 /),'Ma')
	378	type(ctrl_out) :: o_dnwd = ctrl_out((/ 4, 10, 10, 10, 10 /),'dnwd')
	379	type(ctrl_out) :: o_dnwd0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dnwd0')
	380	type(ctrl_out) :: o_dtdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtdyn')
	381	type(ctrl_out) :: o_dqdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dqdyn')
	382	type(ctrl_out) :: o_dudyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dudyn') !AXC
	383	type(ctrl_out) :: o_dvdyn = ctrl_out((/ 4, 10, 10, 10, 1 /),'dvdyn') !AXC
	384	type(ctrl_out) :: o_dtcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtcon')
	385	type(ctrl_out) :: o_ducon = ctrl_out((/ 4, 10, 10, 10, 10 /),'ducon')
	386	type(ctrl_out) :: o_dqcon = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqcon')
	387	type(ctrl_out) :: o_dtwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dtwak')
	388	type(ctrl_out) :: o_dqwak = ctrl_out((/ 4, 5, 10, 10, 10 /),'dqwak')
	389	type(ctrl_out) :: o_wake_h = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_h')
	390	type(ctrl_out) :: o_wake_s = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_s')
	391	type(ctrl_out) :: o_wake_deltat = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_deltat')
	392	type(ctrl_out) :: o_wake_deltaq = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_deltaq')
	393	type(ctrl_out) :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10 /),'wake_omg')
	394	type(ctrl_out) :: o_Vprecip = ctrl_out((/ 10, 10, 10, 10, 10 /),'Vprecip')
	395	type(ctrl_out) :: o_ftd = ctrl_out((/ 4, 5, 10, 10, 10 /),'ftd')
	396	type(ctrl_out) :: o_fqd = ctrl_out((/ 4, 5, 10, 10, 10 /),'fqd')
	397	type(ctrl_out) :: o_dtlsc = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlsc')
	398	type(ctrl_out) :: o_dtlschr = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlschr')
	399	type(ctrl_out) :: o_dqlsc = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqlsc')
	400	type(ctrl_out) :: o_dtvdf = ctrl_out((/ 4, 10, 10, 1, 10 /),'dtvdf')
	401	type(ctrl_out) :: o_dqvdf = ctrl_out((/ 4, 10, 10, 1, 10 /),'dqvdf')
	402	type(ctrl_out) :: o_dteva = ctrl_out((/ 4, 10, 10, 10, 10 /),'dteva')
	403	type(ctrl_out) :: o_dqeva = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqeva')
	404	type(ctrl_out) :: o_ptconv = ctrl_out((/ 4, 10, 10, 10, 10 /),'ptconv')
	405	type(ctrl_out) :: o_ratqs = ctrl_out((/ 4, 10, 10, 10, 10 /),'ratqs')
	406	type(ctrl_out) :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtthe')
	407	type(ctrl_out) :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'f_th')
	408	type(ctrl_out) :: o_e_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'e_th')
	409	type(ctrl_out) :: o_w_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'w_th')
	410	type(ctrl_out) :: o_lambda_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'lambda_th')
	411	type(ctrl_out) :: o_q_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'q_th')
	412	type(ctrl_out) :: o_a_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'a_th')
	413	type(ctrl_out) :: o_d_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'d_th')
	414	type(ctrl_out) :: o_f0_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'f0_th')
	415	type(ctrl_out) :: o_zmax_th = ctrl_out((/ 4, 10, 10, 10, 10 /),'zmax_th')
	416	type(ctrl_out) :: o_dqthe = ctrl_out((/ 4, 10, 10, 10, 1 /),'dqthe')
	417	type(ctrl_out) :: o_dtajs = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtajs')
	418	type(ctrl_out) :: o_dqajs = ctrl_out((/ 4, 10, 10, 10, 10 /),'dqajs')
	419	type(ctrl_out) :: o_dtswr = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtswr')
	420	type(ctrl_out) :: o_dtsw0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtsw0')
	421	type(ctrl_out) :: o_dtlwr = ctrl_out((/ 4, 10, 10, 10, 1 /),'dtlwr')
	422	type(ctrl_out) :: o_dtlw0 = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtlw0')
	423	type(ctrl_out) :: o_dtec = ctrl_out((/ 4, 10, 10, 10, 10 /),'dtec')
	424	type(ctrl_out) :: o_duvdf = ctrl_out((/ 4, 10, 10, 10, 10 /),'duvdf')
	425	type(ctrl_out) :: o_dvvdf = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvvdf')
	426	type(ctrl_out) :: o_duoro = ctrl_out((/ 4, 10, 10, 10, 10 /),'duoro')
	427	type(ctrl_out) :: o_dvoro = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvoro')
	428	type(ctrl_out) :: o_dulif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dulif')
	429	type(ctrl_out) :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10 /),'dvlif')
	430
	431	! Attention a refaire correctement
	432	type(ctrl_out),dimension(2) :: o_trac = (/ ctrl_out((/ 4, 10, 10, 10, 10 /),'trac01'), &
	433	ctrl_out((/ 4, 10, 10, 10, 10 /),'trac02') /)
[907]	434	CONTAINS
	435
	436	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	437	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
	438	!! histbeg, histvert et histdef
	439	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[929]	440
[1114]	441	SUBROUTINE phys_output_open(jjmp1,nlevSTD,clevSTD,nbteta, &
[929]	442	ctetaSTD,dtime, presnivs, ok_veget, &
[1123]	443	type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
[1095]	444	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie)
[907]	445
[929]	446	USE iophy
	447	USE dimphy
[1114]	448	USE infotrac
[929]	449	USE ioipsl
	450	USE mod_phys_lmdz_para
[907]	451
[929]	452	IMPLICIT NONE
	453	include "dimensions.h"
	454	include "temps.h"
	455	include "indicesol.h"
	456	include "clesphys.h"
[964]	457	include "thermcell.h"
[907]	458
[1114]	459	integer :: jjmp1
[929]	460	integer :: nbteta, nlevSTD, radpas
	461	logical :: ok_mensuel, ok_journe, ok_hf, ok_instan
[1095]	462	logical :: ok_LES,ok_ade,ok_aie
[929]	463	real :: dtime
	464	integer :: idayref
	465	real :: zjulian
	466	real, dimension(klev) :: presnivs
	467	character(len=4), dimension(nlevSTD) :: clevSTD
[1083]	468	integer :: nsrf, k, iq, iiq, iff, i, j, ilev
[929]	469	logical :: ok_veget
	470	integer :: iflag_pbl
[1145]	471	CHARACTER(len=4) :: bb2
[929]	472	CHARACTER(len=2) :: bb3
[1123]	473	character(len=6) :: type_ocean
[929]	474	CHARACTER(len=3) :: ctetaSTD(nbteta)
	475	real, dimension(nfiles) :: ecrit_files
	476	CHARACTER(len=20), dimension(nfiles) :: name_files
	477	INTEGER, dimension(iim*jjmp1) :: ndex2d
	478	INTEGER, dimension(iimjjmp1klev) :: ndex3d
	479	integer :: imin_ins, imax_ins
	480	integer :: jmin_ins, jmax_ins
[1176]	481	CHARACTER(len=20), dimension(nfiles) :: type_ecri_files
[907]	482
[1065]	483	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	484	! entre [lonmin_reg,lonmax_reg] et [latmin_reg,latmax_reg]
[1054]	485
[1083]	486	logical, dimension(nfiles), save :: ok_reglim = (/ .false., .false., .false., .false., .true. /)
[1065]	487	real, dimension(nfiles), save :: lonmin_reg = (/ 0., -45., 0., 0., -162. /)
	488	real, dimension(nfiles), save :: lonmax_reg = (/ 90., 45., 90., 90., -144. /)
	489	real, dimension(nfiles), save :: latmin_reg = (/ 0., -45., 0., 0., 7. /)
	490	real, dimension(nfiles), save :: latmax_reg = (/ 90., 90., 90., 90., 21. /)
	491
[1168]	492	levmax = (/ klev, klev, klev, klev, klev /)
[1065]	493
[929]	494	name_files(1) = 'histmth'
	495	name_files(2) = 'histday'
	496	name_files(3) = 'histhf'
	497	name_files(4) = 'histins'
[1054]	498	name_files(5) = 'histLES'
[907]	499
[1176]	500	type_ecri_files(1) = 'ave(X)'
	501	type_ecri_files(2) = 'ave(X)'
	502	type_ecri_files(3) = 'inst(X)'
	503	type_ecri_files(4) = 'inst(X)'
	504	type_ecri_files(5) = 'ave(X)'
[907]	505
[929]	506	clef_files(1) = ok_mensuel
	507	clef_files(2) = ok_journe
	508	clef_files(3) = ok_hf
	509	clef_files(4) = ok_instan
[1054]	510	clef_files(5) = ok_LES
[907]	511
[929]	512	lev_files(1) = lev_histmth
	513	lev_files(2) = lev_histday
	514	lev_files(3) = lev_histhf
	515	lev_files(4) = 1
[1054]	516	lev_files(5) = 1
[907]	517
[929]	518	ecrit_files(1) = ecrit_mth
	519	ecrit_files(2) = ecrit_day
	520	ecrit_files(3) = ecrit_hf
	521	ecrit_files(4) = ecrit_ins
[1054]	522	ecrit_files(5) = ecrit_LES
[907]	523
[1168]	524	!! Lectures des parametres de sorties dans physiq.def
	525	call getin('phys_out_regfkey',ok_reglim)
	526	call getin('phys_out_lonmin',lonmin_reg)
	527	call getin('phys_out_lonmax',lonmax_reg)
	528	call getin('phys_out_latmin',latmin_reg)
	529	call getin('phys_out_latmax',latmax_reg)
	530	call getin('phys_out_levmin',levmin)
	531	call getin('phys_out_levmax',levmax)
	532	call getin('phys_out_filenames',name_files)
[1176]	533	call getin('phys_out_filetypes',type_ecri_files)
[1168]	534	call getin('phys_out_filekeys',clef_files)
	535	call getin('phys_out_filelevels',lev_files)
	536	call getin('phys_out_filetimesteps',ecrit_files)
	537
[1176]	538	type_ecri(:) = type_ecri_files(:)
	539
[907]	540	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	541	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
	542	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
	543	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	544
[1176]	545	zdtime = dtime ! Frequence ou l on moyenne
	546
[929]	547	DO iff=1,nfiles
	548
	549	IF (clef_files(iff)) THEN
[907]	550
[929]	551	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit
[907]	552
[929]	553	idayref = day_ref
	554	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
[907]	555
	556	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
	557	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[929]	558	if (ok_reglim(iff)) then
[907]	559
[929]	560	imin_ins=1
	561	imax_ins=iim
	562	jmin_ins=1
	563	jmax_ins=jjmp1
[907]	564
[1062]	565	! correction abderr do i=1,iim-1
	566	do i=1,iim
[929]	567	print*,'io_lon(i)=',io_lon(i)
	568	if (io_lon(i).le.lonmin_reg(iff)) imin_ins=i
	569	if (io_lon(i).le.lonmax_reg(iff)) imax_ins=i
	570	enddo
[907]	571
[929]	572	do j=1,jjmp1-1
	573	print*,'io_lat(j)=',io_lat(j)
[932]	574	if (io_lat(j).ge.latmin_reg(iff)) jmax_ins=j+1
[929]	575	if (io_lat(j).ge.latmax_reg(iff)) jmin_ins=j
	576	enddo
[907]	577
	578	print*,'On stoke le fichier hist sur, ', &
[929]	579	imin_ins,imax_ins,jmin_ins,jmax_ins
[907]	580	print*,'On stoke le fichier instantanne sur, ', &
[929]	581	io_lon(imin_ins),io_lon(imax_ins), &
	582	io_lat(jmin_ins),io_lat(jmax_ins)
[907]	583
[929]	584	CALL histbeg(name_files(iff),iim,io_lon,jjmp1,io_lat, &
	585	imin_ins,imax_ins-imin_ins+1, &
	586	jmin_ins,jmax_ins-jmin_ins+1, &
[1065]	587	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
[907]	588	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	589	else
[1065]	590	CALL histbeg_phy(name_files(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
[907]	591	endif
	592
[1065]	593	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "mb", &
	594	levmax(iff) - levmin(iff) + 1, &
	595	presnivs(levmin(iff):levmax(iff))/100., nvertm(iff))
[907]	596
	597	!!!!!!!!!!!!! Traitement des champs 3D pour histhf !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	598	!!!!!!!!!!!!!!! A Revoir plus tard !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	599	! IF (iff.eq.3.and.lev_files(iff).ge.4) THEN
	600	! CALL histbeg_phy("histhf3d",itau_phy, &
	601	! & zjulian, dtime, &
[929]	602	! & nhorim, nid_hf3d)
[907]	603
[1083]	604	! CALL histvert(nid_hf3d, "presnivs", &
	605	! & "Vertical levels", "mb", &
[929]	606	! & klev, presnivs/100., nvertm)
[907]	607	! ENDIF
	608
	609	!!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1083]	610	CALL histdef2d(iff,o_phis%flag,o_phis%name,"Surface geop.height", "m2/s2")
[1036]	611	type_ecri(1) = 'once'
	612	type_ecri(2) = 'once'
	613	type_ecri(3) = 'once'
	614	type_ecri(4) = 'once'
[1123]	615	type_ecri(5) = 'once'
[1083]	616	CALL histdef2d(iff,o_aire%flag,o_aire%name,"Grid area", "-")
	617	CALL histdef2d(iff,o_contfracATM%flag,o_contfracATM%name,"% sfce ter+lic", "-")
[1176]	618	type_ecri(:) = type_ecri_files(:)
[1036]	619
	620	!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1083]	621	CALL histdef2d(iff,o_contfracOR%flag,o_contfracOR%name,"% sfce terre OR", "-" )
	622	CALL histdef2d(iff,o_aireTER%flag,o_aireTER%name,"Grid area CONT", "-" )
	623	CALL histdef2d(iff,o_flat%flag,o_flat%name, "Latent heat flux", "W/m2")
	624	CALL histdef2d(iff,o_slp%flag,o_slp%name, "Sea Level Pressure", "Pa" )
	625	CALL histdef2d(iff,o_tsol%flag,o_tsol%name, "Surface Temperature", "K")
	626	CALL histdef2d(iff,o_t2m%flag,o_t2m%name, "Temperature 2m", "K" )
[1123]	627	type_ecri(1) = 't_min(X)'
	628	type_ecri(2) = 't_min(X)'
	629	type_ecri(3) = 't_min(X)'
	630	type_ecri(4) = 't_min(X)'
	631	type_ecri(5) = 't_min(X)'
[1083]	632	CALL histdef2d(iff,o_t2m_min%flag,o_t2m_min%name, "Temp 2m min", "K" )
[1123]	633	type_ecri(1) = 't_max(X)'
	634	type_ecri(2) = 't_max(X)'
	635	type_ecri(3) = 't_max(X)'
	636	type_ecri(4) = 't_max(X)'
	637	type_ecri(5) = 't_max(X)'
[1083]	638	CALL histdef2d(iff,o_t2m_max%flag,o_t2m_max%name, "Temp 2m max", "K" )
[1176]	639	type_ecri(:) = type_ecri_files(:)
[1083]	640	CALL histdef2d(iff,o_wind10m%flag,o_wind10m%name, "10-m wind speed", "m/s")
	641	CALL histdef2d(iff,o_wind10max%flag,o_wind10max%name, "10m wind speed max", "m/s")
	642	CALL histdef2d(iff,o_sicf%flag,o_sicf%name, "Sea-ice fraction", "-" )
	643	CALL histdef2d(iff,o_q2m%flag,o_q2m%name, "Specific humidity 2m", "kg/kg")
	644	CALL histdef2d(iff,o_u10m%flag,o_u10m%name, "Vent zonal 10m", "m/s" )
	645	CALL histdef2d(iff,o_v10m%flag,o_v10m%name, "Vent meridien 10m", "m/s")
	646	CALL histdef2d(iff,o_psol%flag,o_psol%name, "Surface Pressure", "Pa" )
	647	CALL histdef2d(iff,o_qsurf%flag,o_qsurf%name, "Surface Air humidity", "kg/kg")
[907]	648
[929]	649	if (.not. ok_veget) then
[1083]	650	CALL histdef2d(iff,o_qsol%flag,o_qsol%name, "Soil watter content", "mm" )
[929]	651	endif
[907]	652
[1083]	653	CALL histdef2d(iff,o_ndayrain%flag,o_ndayrain%name, "Number of dayrain(liq+sol)", "-")
	654	CALL histdef2d(iff,o_precip%flag,o_precip%name, "Precip Totale liq+sol", "kg/(s*m2)" )
	655	CALL histdef2d(iff,o_plul%flag,o_plul%name, "Large-scale Precip.", "kg/(s*m2)")
	656	CALL histdef2d(iff,o_pluc%flag,o_pluc%name, "Convective Precip.", "kg/(s*m2)")
	657	CALL histdef2d(iff,o_snow%flag,o_snow%name, "Snow fall", "kg/(s*m2)" )
	658	CALL histdef2d(iff,o_evap%flag,o_evap%name, "Evaporat", "kg/(s*m2)" )
	659	CALL histdef2d(iff,o_tops%flag,o_tops%name, "Solar rad. at TOA", "W/m2")
	660	CALL histdef2d(iff,o_tops0%flag,o_tops0%name, "CS Solar rad. at TOA", "W/m2")
	661	CALL histdef2d(iff,o_topl%flag,o_topl%name, "IR rad. at TOA", "W/m2" )
	662	CALL histdef2d(iff,o_topl0%flag,o_topl0%name, "IR rad. at TOA", "W/m2")
	663	CALL histdef2d(iff,o_SWupTOA%flag,o_SWupTOA%name, "SWup at TOA", "W/m2")
	664	CALL histdef2d(iff,o_SWupTOAclr%flag,o_SWupTOAclr%name, "SWup clear sky at TOA", "W/m2")
	665	CALL histdef2d(iff,o_SWdnTOA%flag,o_SWdnTOA%name, "SWdn at TOA", "W/m2" )
	666	CALL histdef2d(iff,o_SWdnTOAclr%flag,o_SWdnTOAclr%name, "SWdn clear sky at TOA", "W/m2")
	667	CALL histdef2d(iff,o_SWup200%flag,o_SWup200%name, "SWup at 200mb", "W/m2" )
	668	CALL histdef2d(iff,o_SWup200clr%flag,o_SWup200clr%name, "SWup clear sky at 200mb", "W/m2")
	669	CALL histdef2d(iff,o_SWdn200%flag,o_SWdn200%name, "SWdn at 200mb", "W/m2" )
	670	CALL histdef2d(iff,o_SWdn200clr%flag,o_SWdn200clr%name, "SWdn clear sky at 200mb", "W/m2")
	671	CALL histdef2d(iff,o_LWup200%flag,o_LWup200%name, "LWup at 200mb", "W/m2")
	672	CALL histdef2d(iff,o_LWup200clr%flag,o_LWup200clr%name, "LWup clear sky at 200mb", "W/m2")
	673	CALL histdef2d(iff,o_LWdn200%flag,o_LWdn200%name, "LWdn at 200mb", "W/m2")
	674	CALL histdef2d(iff,o_LWdn200clr%flag,o_LWdn200clr%name, "LWdn clear sky at 200mb", "W/m2")
	675	CALL histdef2d(iff,o_sols%flag,o_sols%name, "Solar rad. at surf.", "W/m2")
	676	CALL histdef2d(iff,o_sols0%flag,o_sols0%name, "Solar rad. at surf.", "W/m2")
	677	CALL histdef2d(iff,o_soll%flag,o_soll%name, "IR rad. at surface", "W/m2")
	678	CALL histdef2d(iff,o_radsol%flag,o_radsol%name, "Rayonnement au sol", "W/m2")
	679	CALL histdef2d(iff,o_soll0%flag,o_soll0%name, "IR rad. at surface", "W/m2")
	680	CALL histdef2d(iff,o_SWupSFC%flag,o_SWupSFC%name, "SWup at surface", "W/m2")
	681	CALL histdef2d(iff,o_SWupSFCclr%flag,o_SWupSFCclr%name, "SWup clear sky at surface", "W/m2")
	682	CALL histdef2d(iff,o_SWdnSFC%flag,o_SWdnSFC%name, "SWdn at surface", "W/m2")
	683	CALL histdef2d(iff,o_SWdnSFCclr%flag,o_SWdnSFCclr%name, "SWdn clear sky at surface", "W/m2")
	684	CALL histdef2d(iff,o_LWupSFC%flag,o_LWupSFC%name, "Upwd. IR rad. at surface", "W/m2")
	685	CALL histdef2d(iff,o_LWdnSFC%flag,o_LWdnSFC%name, "Down. IR rad. at surface", "W/m2")
	686	CALL histdef2d(iff,o_LWupSFCclr%flag,o_LWupSFCclr%name, "CS Upwd. IR rad. at surface", "W/m2")
	687	CALL histdef2d(iff,o_LWdnSFCclr%flag,o_LWdnSFCclr%name, "Down. CS IR rad. at surface", "W/m2")
	688	CALL histdef2d(iff,o_bils%flag,o_bils%name, "Surf. total heat flux", "W/m2")
	689	CALL histdef2d(iff,o_sens%flag,o_sens%name, "Sensible heat flux", "W/m2")
	690	CALL histdef2d(iff,o_fder%flag,o_fder%name, "Heat flux derivation", "W/m2")
	691	CALL histdef2d(iff,o_ffonte%flag,o_ffonte%name, "Thermal flux for snow melting", "W/m2")
	692	CALL histdef2d(iff,o_fqcalving%flag,o_fqcalving%name, "Ice Calving", "kg/m2/s")
	693	CALL histdef2d(iff,o_fqfonte%flag,o_fqfonte%name, "Land ice melt", "kg/m2/s")
[907]	694
[929]	695	DO nsrf = 1, nbsrf
[1083]	696	CALL histdef2d(iff,o_pourc_srf(nsrf)%flag,o_pourc_srf(nsrf)%name,"% "//clnsurf(nsrf),"%")
	697	CALL histdef2d(iff,o_fract_srf(nsrf)%flag,o_fract_srf(nsrf)%name,"Fraction "//clnsurf(nsrf),"1")
	698	CALL histdef2d(iff,o_taux_srf(nsrf)%flag,o_taux_srf(nsrf)%name,"Zonal wind stress"//clnsurf(nsrf),"Pa")
	699	CALL histdef2d(iff,o_tauy_srf(nsrf)%flag,o_tauy_srf(nsrf)%name,"Meridional wind stress "//clnsurf(nsrf),"Pa")
	700	CALL histdef2d(iff,o_tsol_srf(nsrf)%flag,o_tsol_srf(nsrf)%name,"Temperature "//clnsurf(nsrf),"K")
	701	CALL histdef2d(iff,o_u10m_srf(nsrf)%flag,o_u10m_srf(nsrf)%name,"Vent Zonal 10m "//clnsurf(nsrf),"m/s")
	702	CALL histdef2d(iff,o_v10m_srf(nsrf)%flag,o_v10m_srf(nsrf)%name,"Vent meredien 10m "//clnsurf(nsrf),"m/s")
	703	CALL histdef2d(iff,o_t2m_srf(nsrf)%flag,o_t2m_srf(nsrf)%name,"Temp 2m "//clnsurf(nsrf),"K")
	704	CALL histdef2d(iff,o_sens_srf(nsrf)%flag,o_sens_srf(nsrf)%name,"Sensible heat flux "//clnsurf(nsrf),"W/m2")
	705	CALL histdef2d(iff,o_lat_srf(nsrf)%flag,o_lat_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
	706	CALL histdef2d(iff,o_flw_srf(nsrf)%flag,o_flw_srf(nsrf)%name,"LW "//clnsurf(nsrf),"W/m2")
	707	CALL histdef2d(iff,o_fsw_srf(nsrf)%flag,o_fsw_srf(nsrf)%name,"SW "//clnsurf(nsrf),"W/m2")
	708	CALL histdef2d(iff,o_wbils_srf(nsrf)%flag,o_wbils_srf(nsrf)%name,"Bilan sol "//clnsurf(nsrf),"W/m2" )
	709	CALL histdef2d(iff,o_wbilo_srf(nsrf)%flag,o_wbilo_srf(nsrf)%name,"Bilan eau "//clnsurf(nsrf),"kg/(m2*s)")
[929]	710	if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then
[1083]	711	CALL histdef2d(iff,o_tke_srf(nsrf)%flag,o_tke_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
[1123]	712	type_ecri(1) = 't_max(X)'
	713	type_ecri(2) = 't_max(X)'
	714	type_ecri(3) = 't_max(X)'
	715	type_ecri(4) = 't_max(X)'
	716	type_ecri(5) = 't_max(X)'
[1083]	717	CALL histdef2d(iff,o_tke_max_srf(nsrf)%flag,o_tke_max_srf(nsrf)%name,"Max Turb. Kinetic Energy "//clnsurf(nsrf),"-")
[1176]	718	type_ecri(:) = type_ecri_files(:)
[929]	719	endif
[1083]	720	CALL histdef2d(iff,o_albe_srf(nsrf)%flag,o_albe_srf(nsrf)%name,"Albedo surf. "//clnsurf(nsrf),"-")
	721	CALL histdef2d(iff,o_rugs_srf(nsrf)%flag,o_rugs_srf(nsrf)%name,"Latent heat flux "//clnsurf(nsrf),"W/m2")
	722	CALL histdef2d(iff,o_ages_srf(nsrf)%flag,o_ages_srf(nsrf)%name,"Snow age", "day")
[929]	723	END DO
[907]	724
[1095]	725	IF (ok_ade) THEN
	726	CALL histdef2d(iff,o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2")
	727	CALL histdef2d(iff,o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2")
	728	ENDIF
	729
	730	IF (ok_aie) THEN
	731	CALL histdef2d(iff,o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2")
	732	CALL histdef2d(iff,o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2")
	733	ENDIF
	734
	735
[1083]	736	CALL histdef2d(iff,o_albs%flag,o_albs%name, "Surface albedo", "-")
	737	CALL histdef2d(iff,o_albslw%flag,o_albslw%name, "Surface albedo LW", "-")
	738	CALL histdef2d(iff,o_cdrm%flag,o_cdrm%name, "Momentum drag coef.", "-")
	739	CALL histdef2d(iff,o_cdrh%flag,o_cdrh%name, "Heat drag coef.", "-" )
	740	CALL histdef2d(iff,o_cldl%flag,o_cldl%name, "Low-level cloudiness", "-")
	741	CALL histdef2d(iff,o_cldm%flag,o_cldm%name, "Mid-level cloudiness", "-")
	742	CALL histdef2d(iff,o_cldh%flag,o_cldh%name, "High-level cloudiness", "-")
	743	CALL histdef2d(iff,o_cldt%flag,o_cldt%name, "Total cloudiness", "%")
	744	CALL histdef2d(iff,o_cldq%flag,o_cldq%name, "Cloud liquid water path", "kg/m2")
	745	CALL histdef2d(iff,o_lwp%flag,o_lwp%name, "Cloud water path", "kg/m2")
	746	CALL histdef2d(iff,o_iwp%flag,o_iwp%name, "Cloud ice water path", "kg/m2" )
	747	CALL histdef2d(iff,o_ue%flag,o_ue%name, "Zonal energy transport", "-")
	748	CALL histdef2d(iff,o_ve%flag,o_ve%name, "Merid energy transport", "-")
	749	CALL histdef2d(iff,o_uq%flag,o_uq%name, "Zonal humidity transport", "-")
	750	CALL histdef2d(iff,o_vq%flag,o_vq%name, "Merid humidity transport", "-")
[907]	751
[929]	752	IF(iflag_con.GE.3) THEN ! sb
[1083]	753	CALL histdef2d(iff,o_cape%flag,o_cape%name, "Conv avlbl pot ener", "J/kg")
	754	CALL histdef2d(iff,o_pbase%flag,o_pbase%name, "Cld base pressure", "mb")
	755	CALL histdef2d(iff,o_ptop%flag,o_ptop%name, "Cld top pressure", "mb")
	756	CALL histdef2d(iff,o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s")
	757	CALL histdef2d(iff,o_prw%flag,o_prw%name, "Precipitable water", "kg/m2")
[1123]	758	type_ecri(1) = 't_max(X)'
	759	type_ecri(2) = 't_max(X)'
	760	type_ecri(3) = 't_max(X)'
	761	type_ecri(4) = 't_max(X)'
	762	type_ecri(5) = 't_max(X)'
	763	CALL histdef2d(iff,o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg")
[1176]	764	type_ecri(:) = type_ecri_files(:)
[1123]	765	CALL histdef3d(iff,o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s")
	766	CALL histdef3d(iff,o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s")
	767	CALL histdef3d(iff,o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s")
	768	CALL histdef3d(iff,o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s")
[929]	769	ENDIF !iflag_con .GE. 3
[907]	770
[1083]	771	CALL histdef2d(iff,o_s_pblh%flag,o_s_pblh%name, "Boundary Layer Height", "m")
	772	CALL histdef2d(iff,o_s_pblt%flag,o_s_pblt%name, "t at Boundary Layer Height", "K")
	773	CALL histdef2d(iff,o_s_lcl%flag,o_s_lcl%name, "Condensation level", "m")
	774	CALL histdef2d(iff,o_s_capCL%flag,o_s_capCL%name, "Conv avlbl pot enerfor ABL", "J/m2" )
	775	CALL histdef2d(iff,o_s_oliqCL%flag,o_s_oliqCL%name, "Liq Water in BL", "kg/m2")
	776	CALL histdef2d(iff,o_s_cteiCL%flag,o_s_cteiCL%name, "Instability criteria(ABL)", "K")
	777	CALL histdef2d(iff,o_s_therm%flag,o_s_therm%name, "Exces du thermique", "K")
	778	CALL histdef2d(iff,o_s_trmb1%flag,o_s_trmb1%name, "deep_cape(HBTM2)", "J/m2")
	779	CALL histdef2d(iff,o_s_trmb2%flag,o_s_trmb2%name, "inhibition (HBTM2)", "J/m2")
	780	CALL histdef2d(iff,o_s_trmb3%flag,o_s_trmb3%name, "Point Omega (HBTM2)", "m")
[907]	781
	782	! Champs interpolles sur des niveaux de pression
	783	! iif=1 on ecrit u v w phi sur 850 700 500 200 au niv 1
	784	! iif=2 on ecrit w et ph 500 seulement au niv 1
	785	! et u v sur 850 700 500 200
	786	! iif=3 on ecrit ph a 500 seulement au niv 1
	787	! on ecrit u v t q a 850 700 500 200 au niv 3
	788
[1055]	789	type_ecri(1) = 'inst(X)'
	790	type_ecri(2) = 'inst(X)'
	791	type_ecri(3) = 'inst(X)'
	792	type_ecri(4) = 'inst(X)'
[1123]	793	type_ecri(5) = 'inst(X)'
[1083]	794
	795	! Attention a reverifier
	796
	797	ilev=0
[929]	798	DO k=1, nlevSTD
	799	IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)
[1083]	800	! IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)
[929]	801	IF(bb2.EQ."850".OR.bb2.EQ."700".OR.bb2.EQ."500".OR.bb2.EQ."200")THEN
[1083]	802	ilev=ilev+1
	803	print*,'ilev bb2 flag name ',ilev,bb2,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name
	804	CALL histdef2d(iff,o_uSTDlevs(ilev)%flag,o_uSTDlevs(ilev)%name,"Zonal wind "//bb2//"mb", "m/s")
	805	CALL histdef2d(iff,o_vSTDlevs(ilev)%flag,o_vSTDlevs(ilev)%name,"Meridional wind "//bb2//"mb", "m/s")
	806	CALL histdef2d(iff,o_wSTDlevs(ilev)%flag,o_wSTDlevs(ilev)%name,"Vertical wind "//bb2//"mb", "Pa/s")
	807	CALL histdef2d(iff,o_phiSTDlevs(ilev)%flag,o_phiSTDlevs(ilev)%name,"Geopotential "//bb2//"mb", "m")
	808	CALL histdef2d(iff,o_qSTDlevs(ilev)%flag,o_qSTDlevs(ilev)%name,"Specific humidity "//bb2//"mb", "kg/kg" )
	809	CALL histdef2d(iff,o_tSTDlevs(ilev)%flag,o_tSTDlevs(ilev)%name,"Temperature "//bb2//"mb", "K")
	810	ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR."500".OR.bb2.EQ."200")
[907]	811	ENDDO
[1176]	812	type_ecri(:) = type_ecri_files(:)
[907]	813
[1083]	814	CALL histdef2d(iff,o_t_oce_sic%flag,o_t_oce_sic%name, "Temp mixte oce-sic", "K")
[907]	815
[1123]	816	IF (type_ocean=='slab') &
[1083]	817	CALL histdef2d(iff,o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2")
[907]	818
[1123]	819	! Couplage conv-CL
	820	IF (iflag_con.GE.3) THEN
	821	IF (iflag_coupl.EQ.1) THEN
	822	CALL histdef2d(iff,o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2")
	823	CALL histdef2d(iff,o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2")
	824	ENDIF
	825	ENDIF !(iflag_con.GE.3)
[945]	826
[1083]	827	CALL histdef2d(iff,o_weakinv%flag,o_weakinv%name, "Weak inversion", "-")
	828	CALL histdef2d(iff,o_dthmin%flag,o_dthmin%name, "dTheta mini", "K/m")
	829	CALL histdef2d(iff,o_rh2m%flag,o_rh2m%name, "Relative humidity at 2m", "%" )
	830	CALL histdef2d(iff,o_qsat2m%flag,o_qsat2m%name, "Saturant humidity at 2m", "%")
	831	CALL histdef2d(iff,o_tpot%flag,o_tpot%name, "Surface air potential temperature", "K")
	832	CALL histdef2d(iff,o_tpote%flag,o_tpote%name, "Surface air equivalent potential temperature", "K")
	833	CALL histdef2d(iff,o_SWnetOR%flag,o_SWnetOR%name, "Sfce net SW radiation OR", "W/m2")
	834	CALL histdef2d(iff,o_SWdownOR%flag,o_SWdownOR%name, "Sfce incident SW radiation OR", "W/m2")
	835	CALL histdef2d(iff,o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2")
	836	CALL histdef2d(iff,o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)")
[1123]	837
[1083]	838	CALL histdef2d(iff,o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2")
	839	CALL histdef2d(iff,o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s")
	840	CALL histdef2d(iff,o_dtsvdft%flag,o_dtsvdft%name, "Boundary-layer dTs(t)", "K/s")
	841	CALL histdef2d(iff,o_dtsvdfg%flag,o_dtsvdfg%name, "Boundary-layer dTs(g)", "K/s")
	842	CALL histdef2d(iff,o_dtsvdfi%flag,o_dtsvdfi%name, "Boundary-layer dTs(g)", "K/s")
	843	CALL histdef2d(iff,o_rugs%flag,o_rugs%name, "rugosity", "-" )
[907]	844
	845	! Champs 3D:
[1083]	846	CALL histdef3d(iff,o_lwcon%flag,o_lwcon%name, "Cloud liquid water content", "kg/kg")
	847	CALL histdef3d(iff,o_iwcon%flag,o_iwcon%name, "Cloud ice water content", "kg/kg")
	848	CALL histdef3d(iff,o_temp%flag,o_temp%name, "Air temperature", "K" )
	849	CALL histdef3d(iff,o_theta%flag,o_theta%name, "Potential air temperature", "K" )
[1123]	850	CALL histdef3d(iff,o_ovap%flag,o_ovap%name, "Specific humidity + dqphy", "kg/kg" )
	851	CALL histdef3d(iff,o_ovapinit%flag,o_ovapinit%name, "Specific humidity", "kg/kg" )
[1083]	852	CALL histdef3d(iff,o_geop%flag,o_geop%name, "Geopotential height", "m2/s2")
	853	CALL histdef3d(iff,o_vitu%flag,o_vitu%name, "Zonal wind", "m/s" )
	854	CALL histdef3d(iff,o_vitv%flag,o_vitv%name, "Meridional wind", "m/s" )
	855	CALL histdef3d(iff,o_vitw%flag,o_vitw%name, "Vertical wind", "Pa/s" )
	856	CALL histdef3d(iff,o_pres%flag,o_pres%name, "Air pressure", "Pa" )
	857	CALL histdef3d(iff,o_rneb%flag,o_rneb%name, "Cloud fraction", "-")
	858	CALL histdef3d(iff,o_rnebcon%flag,o_rnebcon%name, "Convective Cloud Fraction", "-")
	859	CALL histdef3d(iff,o_rhum%flag,o_rhum%name, "Relative humidity", "-")
	860	CALL histdef3d(iff,o_ozone%flag,o_ozone%name, "Ozone concentration", "ppmv")
	861	CALL histdef3d(iff,o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s")
	862	CALL histdef3d(iff,o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s")
	863	CALL histdef3d(iff,o_cldtau%flag,o_cldtau%name, "Cloud optical thickness", "1")
	864	CALL histdef3d(iff,o_cldemi%flag,o_cldemi%name, "Cloud optical emissivity", "1")
[973]	865	!IM: bug ?? dimensionnement variables (klon,klev+1) pmflxr, pmflxs, prfl, psfl
[1083]	866	! CALL histdef3d(iff,o_pr_con_l%flag,o_pmflxr%name, "Convective precipitation lic", " ")
	867	! CALL histdef3d(iff,o_pr_con_i%flag,o_pmflxs%name, "Convective precipitation ice", " ")
	868	! CALL histdef3d(iff,o_pr_lsc_l%flag,o_prfl%name, "Large scale precipitation lic", " ")
	869	! CALL histdef3d(iff,o_pr_lsc_i%flag,o_psfl%name, "Large scale precipitation ice", " ")
[907]	870
[929]	871	!FH Sorties pour la couche limite
	872	if (iflag_pbl>1) then
[1083]	873	CALL histdef3d(iff,o_tke%flag,o_tke%name, "TKE", "m2/s2")
[1123]	874	type_ecri(1) = 't_max(X)'
	875	type_ecri(2) = 't_max(X)'
	876	type_ecri(3) = 't_max(X)'
	877	type_ecri(4) = 't_max(X)'
	878	type_ecri(5) = 't_max(X)'
[1083]	879	CALL histdef3d(iff,o_tke_max%flag,o_tke_max%name, "TKE max", "m2/s2")
[1176]	880	type_ecri(:) = type_ecri_files(:)
[929]	881	endif
[907]	882
[1083]	883	CALL histdef3d(iff,o_kz%flag,o_kz%name, "Kz melange", "m2/s")
[1123]	884	type_ecri(1) = 't_max(X)'
	885	type_ecri(2) = 't_max(X)'
	886	type_ecri(3) = 't_max(X)'
	887	type_ecri(4) = 't_max(X)'
	888	type_ecri(5) = 't_max(X)'
[1083]	889	CALL histdef3d(iff,o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" )
[1176]	890	type_ecri(:) = type_ecri_files(:)
[1123]	891	CALL histdef3d(iff,o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg")
[1083]	892	CALL histdef3d(iff,o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s")
	893	CALL histdef3d(iff,o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s")
	894	CALL histdef3d(iff,o_dudyn%flag,o_dudyn%name, "Dynamics dU", "m/s2")
	895	CALL histdef3d(iff,o_dvdyn%flag,o_dvdyn%name, "Dynamics dV", "m/s2")
	896	CALL histdef3d(iff,o_dtcon%flag,o_dtcon%name, "Convection dT", "K/s")
	897	CALL histdef3d(iff,o_ducon%flag,o_ducon%name, "Convection du", "m/s2")
	898	CALL histdef3d(iff,o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s")
[1123]	899
	900	! Wakes
	901	IF(iflag_con.EQ.3) THEN
	902	IF (iflag_wake == 1) THEN
	903	CALL histdef2d(iff,o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2")
	904	CALL histdef2d(iff,o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2")
	905	CALL histdef2d(iff,o_ale%flag,o_ale%name, "ALE", "m2/s2")
	906	CALL histdef2d(iff,o_alp%flag,o_alp%name, "ALP", "W/m2")
	907	CALL histdef2d(iff,o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2")
	908	CALL histdef2d(iff,o_wape%flag,o_WAPE%name, "WAPE", "m2/s2")
	909	CALL histdef2d(iff,o_wake_h%flag,o_wake_h%name, "wake_h", "-")
	910	CALL histdef2d(iff,o_wake_s%flag,o_wake_s%name, "wake_s", "-")
	911	CALL histdef3d(iff,o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s")
	912	CALL histdef3d(iff,o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s")
	913	CALL histdef3d(iff,o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ")
	914	CALL histdef3d(iff,o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ")
	915	CALL histdef3d(iff,o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ")
	916	ENDIF
	917	CALL histdef3d(iff,o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-")
	918	CALL histdef3d(iff,o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-")
	919	CALL histdef3d(iff,o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-")
	920	ENDIF !(iflag_con.EQ.3)
	921
[1083]	922	CALL histdef3d(iff,o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s")
	923	CALL histdef3d(iff,o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s")
	924	CALL histdef3d(iff,o_dqlsc%flag,o_dqlsc%name, "Condensation dQ", "(kg/kg)/s")
	925	CALL histdef3d(iff,o_dtvdf%flag,o_dtvdf%name, "Boundary-layer dT", "K/s")
	926	CALL histdef3d(iff,o_dqvdf%flag,o_dqvdf%name, "Boundary-layer dQ", "(kg/kg)/s")
	927	CALL histdef3d(iff,o_dteva%flag,o_dteva%name, "Reevaporation dT", "K/s")
	928	CALL histdef3d(iff,o_dqeva%flag,o_dqeva%name, "Reevaporation dQ", "(kg/kg)/s")
	929	CALL histdef3d(iff,o_ptconv%flag,o_ptconv%name, "POINTS CONVECTIFS", " ")
	930	CALL histdef3d(iff,o_ratqs%flag,o_ratqs%name, "RATQS", " ")
	931	CALL histdef3d(iff,o_dtthe%flag,o_dtthe%name, "Dry adjust. dT", "K/s")
[1036]	932
	933	if(iflag_thermals.gt.1) THEN
[1083]	934	CALL histdef3d(iff,o_f_th%flag,o_f_th%name, "Thermal plume mass flux", "K/s")
	935	CALL histdef3d(iff,o_e_th%flag,o_e_th%name,"Thermal plume entrainment","K/s")
	936	CALL histdef3d(iff,o_w_th%flag,o_w_th%name,"Thermal plume vertical velocity","m/s")
	937	CALL histdef3d(iff,o_lambda_th%flag,o_lambda_th%name,"Thermal plume vertical velocity","m/s")
	938	CALL histdef3d(iff,o_q_th%flag,o_q_th%name, "Thermal plume total humidity", "kg/kg")
	939	CALL histdef3d(iff,o_a_th%flag,o_a_th%name, "Thermal plume fraction", "")
	940	CALL histdef3d(iff,o_d_th%flag,o_d_th%name, "Thermal plume detrainment", "K/s")
[1036]	941	endif !iflag_thermals.gt.1
[1083]	942	CALL histdef2d(iff,o_f0_th%flag,o_f0_th%name, "Thermal closure mass flux", "K/s")
	943	CALL histdef2d(iff,o_zmax_th%flag,o_zmax_th%name, "Thermal plume height", "K/s")
	944	CALL histdef3d(iff,o_dqthe%flag,o_dqthe%name, "Dry adjust. dQ", "(kg/kg)/s")
	945	CALL histdef3d(iff,o_dtajs%flag,o_dtajs%name, "Dry adjust. dT", "K/s")
	946	CALL histdef3d(iff,o_dqajs%flag,o_dqajs%name, "Dry adjust. dQ", "(kg/kg)/s")
	947	CALL histdef3d(iff,o_dtswr%flag,o_dtswr%name, "SW radiation dT", "K/s")
	948	CALL histdef3d(iff,o_dtsw0%flag,o_dtsw0%name, "CS SW radiation dT", "K/s")
	949	CALL histdef3d(iff,o_dtlwr%flag,o_dtlwr%name, "LW radiation dT", "K/s")
	950	CALL histdef3d(iff,o_dtlw0%flag,o_dtlw0%name, "CS LW radiation dT", "K/s")
	951	CALL histdef3d(iff,o_dtec%flag,o_dtec%name, "Cinetic dissip dT", "K/s")
	952	CALL histdef3d(iff,o_duvdf%flag,o_duvdf%name, "Boundary-layer dU", "m/s2")
	953	CALL histdef3d(iff,o_dvvdf%flag,o_dvvdf%name, "Boundary-layer dV", "m/s2")
[907]	954
[929]	955	IF (ok_orodr) THEN
[1083]	956	CALL histdef3d(iff,o_duoro%flag,o_duoro%name, "Orography dU", "m/s2")
	957	CALL histdef3d(iff,o_dvoro%flag,o_dvoro%name, "Orography dV", "m/s2")
[929]	958	ENDIF
[907]	959
[929]	960	IF (ok_orolf) THEN
[1083]	961	CALL histdef3d(iff,o_dulif%flag,o_dulif%name, "Orography dU", "m/s2")
	962	CALL histdef3d(iff,o_dvlif%flag,o_dvlif%name, "Orography dV", "m/s2")
[929]	963	ENDIF
[907]	964
[1114]	965	if (nqtot>=3) THEN
	966	!Attention DO iq=3,nqtot
[1083]	967	DO iq=3,4
[929]	968	iiq=niadv(iq)
[1083]	969	! CALL histdef3d (iff, o_trac%flag,'o_'//tnom(iq)%name,ttext(iiq), "-" )
	970	CALL histdef3d (iff, o_trac(iq-2)%flag,o_trac(iq-2)%name,ttext(iiq), "-" )
[929]	971	ENDDO
[907]	972	endif
	973
[929]	974	CALL histend(nid_files(iff))
[907]	975
[929]	976	ndex2d = 0
	977	ndex3d = 0
[907]	978
[929]	979	ENDIF ! clef_files
[907]	980
[929]	981	ENDDO !
	982	end subroutine phys_output_open
[907]	983
[929]	984	SUBROUTINE histdef2d (iff,flag_var,nomvar,titrevar,unitvar)
	985
	986	use ioipsl
	987	USE dimphy
	988	USE mod_phys_lmdz_para
[907]	989
[929]	990	IMPLICIT NONE
	991
	992	include "dimensions.h"
	993	include "temps.h"
	994	include "indicesol.h"
	995	include "clesphys.h"
[907]	996
[929]	997	integer :: iff
	998	integer, dimension(nfiles) :: flag_var
[1083]	999	character(len=20) :: nomvar
[929]	1000	character(len=*) :: titrevar
	1001	character(len=*) :: unitvar
[1083]	1002
[1176]	1003	real zstophym
	1004
	1005	if (type_ecri(iff)=='inst(X)') then
	1006	zstophym=zoutm(iff)
	1007	else
	1008	zstophym=zdtime
	1009	endif
	1010
[1123]	1011	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
	1012	call conf_physoutputs(nomvar,flag_var)
[929]	1013
	1014	if ( flag_var(iff)<=lev_files(iff) ) then
	1015	call histdef (nid_files(iff),nomvar,titrevar,unitvar, &
	1016	iim,jj_nb,nhorim(iff), 1,1,1, -99, 32, &
[1176]	1017	type_ecri(iff), zstophym,zoutm(iff))
[929]	1018	endif
	1019	end subroutine histdef2d
[907]	1020
[929]	1021	SUBROUTINE histdef3d (iff,flag_var,nomvar,titrevar,unitvar)
[907]	1022
[929]	1023	use ioipsl
	1024	USE dimphy
	1025	USE mod_phys_lmdz_para
[907]	1026
[929]	1027	IMPLICIT NONE
[907]	1028
[929]	1029	include "dimensions.h"
	1030	include "temps.h"
	1031	include "indicesol.h"
	1032	include "clesphys.h"
[907]	1033
[929]	1034	integer :: iff
	1035	integer, dimension(nfiles) :: flag_var
[1083]	1036	character(len=20) :: nomvar
[929]	1037	character(len=*) :: titrevar
	1038	character(len=*) :: unitvar
[907]	1039
[1176]	1040	real zstophym
	1041
[1123]	1042	! Appel a la lecture des noms et niveau d'ecriture des variables dans output.def
	1043	call conf_physoutputs(nomvar,flag_var)
[1083]	1044
[1176]	1045	if (type_ecri(iff)=='inst(X)') then
	1046	zstophym=zoutm(iff)
	1047	else
	1048	zstophym=zdtime
	1049	endif
	1050
[929]	1051	if ( flag_var(iff)<=lev_files(iff) ) then
[1065]	1052	call histdef (nid_files(iff), nomvar, titrevar, unitvar, &
	1053	iim, jj_nb, nhorim(iff), klev, levmin(iff), &
	1054	levmax(iff)-levmin(iff)+1, nvertm(iff), 32, type_ecri(iff), &
[1176]	1055	zstophym, zoutm(iff))
[929]	1056	endif
	1057	end subroutine histdef3d
[907]	1058
[1083]	1059	SUBROUTINE conf_physoutputs(nam_var,flag_var)
[1123]	1060	!!! Lecture des noms et niveau de sortie des variables dans output.def
	1061	! en utilisant les routines getin de IOIPSL
[1083]	1062	use ioipsl
	1063
	1064	IMPLICIT NONE
	1065
[1134]	1066	include 'iniprint.h'
	1067
[1083]	1068	character(len=20) :: nam_var
	1069	integer, dimension(nfiles) :: flag_var
	1070	integer, dimension(nfiles),save :: flag_var_omp
	1071	character(len=20),save :: nam_var_omp
	1072
	1073	flag_var_omp = flag_var
	1074	nam_var_omp = nam_var
[1134]	1075	IF(prt_level>10) WRITE(lunout,*)'Avant getin: nam_var flag_var ',nam_var,flag_var(:)
[1083]	1076	call getin('flag_'//nam_var,flag_var_omp)
	1077	flag_var = flag_var_omp
	1078	call getin('name_'//nam_var,nam_var_omp)
	1079	nam_var=nam_var_omp
	1080
[1134]	1081	IF(prt_level>10) WRITE(lunout,*)'Apres getin: nam_var flag_var ',nam_var,flag_var(:)
[1083]	1082
	1083	END SUBROUTINE conf_physoutputs
	1084
[907]	1085	END MODULE phys_output_mod
	1086

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

Download in other formats:

Original Format