Context Navigation

phys_output_mod.F90 @ 3603

Last change on this file since 3603 was 3003, checked in by Laurent Fairhead, 7 years ago
Modifications to the code and xml files to output Ap and B, the coefficients of the hybrid coordinates as requested by the CMIP6 DataRequest? LF (with guidance from A. Caubel and S. Senesi)
Property copyright set to Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 25.0 KB

Line
1	! $Id: phys_output_mod.F90 3003 2017-10-03 14:45:26Z musat $
2	!
3
4	MODULE phys_output_mod
5	USE indice_sol_mod
6	USE phys_output_var_mod
7	USE phys_output_write_mod, ONLY : phys_output_write
8	REAL, DIMENSION(nfiles),SAVE :: ecrit_files
9
10	! Abderrahmane 12 2007
11	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
12	!!! Ecreture des Sorties du modele dans les fichiers Netcdf :
13	! histmth.nc : moyennes mensuelles
14	! histday.nc : moyennes journalieres
15	! histhf.nc : moyennes toutes les 3 heures
16	! histins.nc : valeurs instantanees
17	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
18
19	CONTAINS
20
21	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
22	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
23	!! histbeg, histvert et histdef
24	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
25
26	SUBROUTINE phys_output_open(rlon,rlat,pim,tabij,ipt,jpt,plon,plat, &
27	jjmp1,nlevSTD,clevSTD,rlevSTD, dtime, ok_veget, &
28	type_ocean, iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
29	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
30	phys_out_filestations, &
31	new_aod, aerosol_couple, flag_aerosol_strat, &
32	pdtphys, paprs, pphis, pplay, lmax_th, ptconv, ptconvth, ivap, &
33	d_u, d_t, qx, d_qx, zmasse, ok_sync)
34
35	USE iophy
36	USE dimphy
37	USE infotrac_phy, ONLY: nqtot, nqo, niadv, tname, ttext
38	USE ioipsl
39	USE phys_cal_mod, only : hour, calend
40	USE mod_phys_lmdz_para
41	!Martin
42	USE surface_data, ONLY : ok_snow
43	USE phys_output_ctrlout_mod
44	USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat
45	USE print_control_mod, ONLY: prt_level,lunout
46	USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt
47	USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
48	#ifdef CPP_XIOS
49	! ug Pour les sorties XIOS
50	USE wxios
51	#endif
52
53	IMPLICIT NONE
54	include "clesphys.h"
55	include "thermcell.h"
56	include "YOMCST.h"
57
58	! ug Nouveaux arguments n\'ecessaires au histwrite_mod:
59	INTEGER, INTENT(IN) :: ivap
60	INTEGER, DIMENSION(klon), INTENT(IN) :: lmax_th
61	LOGICAL, INTENT(IN) :: ok_sync
62	LOGICAL, DIMENSION(klon, klev), INTENT(IN) :: ptconv, ptconvth
63	REAL, INTENT(IN) :: pdtphys
64	REAL, DIMENSION(klon), INTENT(IN) :: pphis
65	REAL, DIMENSION(klon, klev), INTENT(IN) :: pplay, d_u, d_t
66	REAL, DIMENSION(klon, klev+1), INTENT(IN) :: paprs
67	REAL, DIMENSION(klon,klev,nqtot), INTENT(IN):: qx, d_qx
68	REAL, DIMENSION(klon, klev), INTENT(IN) :: zmasse
69
70
71	REAL,DIMENSION(klon),INTENT(IN) :: rlon
72	REAL,DIMENSION(klon),INTENT(IN) :: rlat
73	INTEGER, INTENT(IN) :: pim
74	INTEGER, DIMENSION(pim) :: tabij
75	INTEGER,DIMENSION(pim), INTENT(IN) :: ipt, jpt
76	REAL,DIMENSION(pim), INTENT(IN) :: plat, plon
77	REAL,DIMENSION(pim,2) :: plat_bounds, plon_bounds
78
79	INTEGER :: jjmp1
80	INTEGER :: nlevSTD, radpas
81	LOGICAL :: ok_mensuel, ok_journe, ok_hf, ok_instan
82	LOGICAL :: ok_LES,ok_ade,ok_aie
83	INTEGER :: flag_aerosol_strat
84	LOGICAL :: new_aod, aerosol_couple
85	INTEGER, INTENT(IN):: read_climoz ! read ozone climatology
86	! Allowed values are 0, 1 and 2
87	! 0: do not read an ozone climatology
88	! 1: read a single ozone climatology that will be used day and night
89	! 2: read two ozone climatologies, the average day and night
90	! climatology and the daylight climatology
91
92	REAL :: dtime
93	INTEGER :: idayref
94	REAL :: zjulian_start, zjulian
95	CHARACTER(LEN=4), DIMENSION(nlevSTD) :: clevSTD
96	REAL, DIMENSION(nlevSTD) :: rlevSTD
97	INTEGER :: nsrf, k, iq, iiq, iff, i, j, ilev
98	INTEGER :: naero
99	LOGICAL :: ok_veget
100	INTEGER :: iflag_pbl
101	INTEGER :: iflag_pbl_split
102	CHARACTER(LEN=4) :: bb2
103	CHARACTER(LEN=2) :: bb3
104	CHARACTER(LEN=6) :: type_ocean
105	INTEGER, DIMENSION(nbp_lon*jjmp1) :: ndex2d
106	INTEGER, DIMENSION(nbp_lonjjmp1klev) :: ndex3d
107	INTEGER :: imin_ins, imax_ins
108	INTEGER :: jmin_ins, jmax_ins
109	INTEGER, DIMENSION(nfiles) :: phys_out_levmin, phys_out_levmax
110	INTEGER, DIMENSION(nfiles) :: phys_out_filelevels
111	CHARACTER(LEN=20), DIMENSION(nfiles) :: chtimestep = (/ 'Default', 'Default', 'Default', 'Default', 'Default', &
112	'Default', 'Default', 'Default', 'Default', 'Default' /)
113	LOGICAL, DIMENSION(nfiles) :: phys_out_filekeys
114	LOGICAL, DIMENSION(nfiles) :: phys_out_filestations
115
116	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
117	! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]
118
119	LOGICAL, DIMENSION(nfiles), SAVE :: phys_out_regfkey = (/ .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., &
120	.FALSE., .FALSE., .FALSE., .FALSE., .FALSE. /)
121	REAL, DIMENSION(nfiles), SAVE :: phys_out_lonmin = (/ -180., -180., -180., -180., -180., &
122	-180., -180., -180., -180., -180. /)
123	REAL, DIMENSION(nfiles), SAVE :: phys_out_lonmax = (/ 180., 180., 180., 180., 180., &
124	180., 180., 180., 180., 180. /)
125	REAL, DIMENSION(nfiles), SAVE :: phys_out_latmin = (/ -90., -90., -90., -90., -90., &
126	-90., -90., -90., -90., -90. /)
127	REAL, DIMENSION(nfiles), SAVE :: phys_out_latmax = (/ 90., 90., 90., 90., 90., &
128	90., 90., 90., 90., 90. /)
129	REAL, DIMENSION(klev,2) :: Ahyb_bounds, Bhyb_bounds
130	REAL, DIMENSION(klev) :: lev_index
131
132	#ifdef CPP_XIOS
133	! ug Variables utilis\'ees pour r\'ecup\'erer le calendrier pour xios
134	INTEGER :: x_an, x_mois, x_jour
135	REAL :: x_heure
136	INTEGER :: ini_an, ini_mois, ini_jour
137	REAL :: ini_heure
138	#endif
139
140	WRITE(lunout,*) 'Debut phys_output_mod.F90'
141	! Initialisations (Valeurs par defaut
142
143	DO ilev=1,klev
144	Ahyb_bounds(ilev,1) = ap(ilev)
145	Ahyb_bounds(ilev,2) = ap(ilev+1)
146	Bhyb_bounds(ilev,1) = bp(ilev)
147	Bhyb_bounds(ilev,2) = bp(ilev+1)
148	lev_index(ilev) = REAL(ilev)
149	END DO
150
151	IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
152	IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot))
153	ALLOCATE(o_dtr_the(nqtot),o_dtr_con(nqtot),o_dtr_lessi_impa(nqtot))
154	ALLOCATE(o_dtr_lessi_nucl(nqtot),o_dtr_insc(nqtot),o_dtr_bcscav(nqtot))
155	ALLOCATE(o_dtr_evapls(nqtot),o_dtr_ls(nqtot),o_dtr_trsp(nqtot))
156	ALLOCATE(o_dtr_sscav(nqtot),o_dtr_sat(nqtot),o_dtr_uscav(nqtot))
157	ALLOCATE(o_dtr_dry(nqtot),o_dtr_vdf(nqtot))
158
159	levmax = (/ klev, klev, klev, klev, klev, klev, nlevSTD, nlevSTD, nlevSTD, klev /)
160
161	phys_out_filenames(1) = 'histmth'
162	phys_out_filenames(2) = 'histday'
163	phys_out_filenames(3) = 'histhf6h'
164	phys_out_filenames(4) = 'histhf3h'
165	phys_out_filenames(5) = 'histhf3hm'
166	phys_out_filenames(6) = 'histstn'
167	phys_out_filenames(7) = 'histmthNMC'
168	phys_out_filenames(8) = 'histdayNMC'
169	phys_out_filenames(9) = 'histhfNMC'
170	phys_out_filenames(10)= 'histstrataer'
171
172	type_ecri(1) = 'ave(X)'
173	type_ecri(2) = 'ave(X)'
174	type_ecri(3) = 'inst(X)'
175	type_ecri(4) = 'inst(X)'
176	type_ecri(5) = 'ave(X)'
177	type_ecri(6) = 'inst(X)'
178	type_ecri(7) = 'inst(X)'
179	type_ecri(8) = 'inst(X)'
180	type_ecri(9) = 'inst(X)'
181	type_ecri(10)= 'ave(X)'
182
183	clef_files(1) = ok_mensuel
184	clef_files(2) = ok_journe
185	clef_files(3) = ok_hf
186	clef_files(4) = ok_instan
187	clef_files(5) = ok_LES
188	clef_files(6) = ok_instan
189	clef_files(7) = ok_histNMC(1)
190	clef_files(8) = ok_histNMC(2)
191	clef_files(9) = ok_histNMC(3)
192	#ifdef CPP_StratAer
193	clef_files(10)= .TRUE.
194	#else
195	clef_files(10)= .FALSE.
196	#endif
197
198	!sortir des fichiers "stations" si clef_stations(:)=.TRUE.
199	clef_stations(1) = .FALSE.
200	clef_stations(2) = .FALSE.
201	clef_stations(3) = .FALSE.
202	clef_stations(4) = .FALSE.
203	clef_stations(5) = .FALSE.
204	clef_stations(6) = .FALSE.
205	clef_stations(7) = .FALSE.
206	clef_stations(8) = .FALSE.
207	clef_stations(9) = .FALSE.
208	clef_stations(10)= .FALSE.
209
210	lev_files(1) = lev_histmth
211	lev_files(2) = lev_histday
212	lev_files(3) = lev_histhf
213	lev_files(4) = lev_histins
214	lev_files(5) = lev_histLES
215	lev_files(6) = lev_histins
216	lev_files(7) = levout_histNMC(1)
217	lev_files(8) = levout_histNMC(2)
218	lev_files(9) = levout_histNMC(3)
219	lev_files(10)= 5
220
221	ecrit_files(1) = ecrit_mth
222	ecrit_files(2) = ecrit_day
223	ecrit_files(3) = ecrit_hf
224	ecrit_files(4) = ecrit_ins
225	ecrit_files(5) = ecrit_LES
226	ecrit_files(6) = ecrit_ins
227	ecrit_files(7) = freq_outNMC(1)
228	ecrit_files(8) = freq_outNMC(2)
229	ecrit_files(9) = freq_outNMC(3)
230	ecrit_files(10)= ecrit_mth
231
232	!! Lectures des parametres de sorties dans physiq.def
233
234	CALL getin('phys_out_regfkey',phys_out_regfkey)
235	CALL getin('phys_out_lonmin',phys_out_lonmin)
236	CALL getin('phys_out_lonmax',phys_out_lonmax)
237	CALL getin('phys_out_latmin',phys_out_latmin)
238	CALL getin('phys_out_latmax',phys_out_latmax)
239	phys_out_levmin(:)=levmin(:)
240	CALL getin('phys_out_levmin',levmin)
241	phys_out_levmax(:)=levmax(:)
242	CALL getin('phys_out_levmax',levmax)
243	CALL getin('phys_out_filenames',phys_out_filenames)
244	phys_out_filekeys(:)=clef_files(:)
245	CALL getin('phys_out_filekeys',clef_files)
246	phys_out_filestations(:)=clef_stations(:)
247	CALL getin('phys_out_filestations',clef_stations)
248	phys_out_filelevels(:)=lev_files(:)
249	CALL getin('phys_out_filelevels',lev_files)
250	CALL getin('phys_out_filetimesteps',chtimestep)
251	phys_out_filetypes(:)=type_ecri(:)
252	CALL getin('phys_out_filetypes',type_ecri)
253
254	type_ecri_files(:)=type_ecri(:)
255
256	! if (ok_all_xml) phys_out_filelevels = 999
257
258	WRITE(lunout,*)'phys_out_lonmin=',phys_out_lonmin
259	WRITE(lunout,*)'phys_out_lonmax=',phys_out_lonmax
260	WRITE(lunout,*)'phys_out_latmin=',phys_out_latmin
261	WRITE(lunout,*)'phys_out_latmax=',phys_out_latmax
262	WRITE(lunout,*)'phys_out_filenames=',phys_out_filenames
263	WRITE(lunout,*)'phys_out_filetypes=',type_ecri
264	WRITE(lunout,*)'phys_out_filekeys=',clef_files
265	WRITE(lunout,*)'phys_out_filestations=',clef_stations
266	WRITE(lunout,*)'phys_out_filelevels=',lev_files
267	WRITE(lunout,*)'phys_out_regfkey=',phys_out_regfkey
268
269	! A noter pour
270	! l heure initiale - dans les fichiers histoire hist* - on met comme
271	! heure de debut soit la vraie heure (pour le 1D) soit 0h (pour le 3D)
272	! afin d avoir une seule sortie mensuelle par mois lorsque l on tourne
273	! par annee (IM).
274	!
275	idayref = day_ref
276	IF (klon_glo==1) THEN
277	! current_time (used to compute hour) is updated at the begining of
278	! the physics; to set the correct outputs "initial time" we thus
279	! have to use (hour-dtphys).
280	CALL ymds2ju(annee_ref, 1, idayref, hour-pdtphys, zjulian)
281	print *,'phys_output_mod: annee,iday,hour,zjulian=',annee_ref,idayref, hour, zjulian
282	ELSE
283	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
284	CALL ymds2ju(annee_ref, 1, day_ini, start_time*rday, zjulian_start)
285	ENDIF
286
287	#ifdef CPP_XIOS
288	! ug R\'eglage du calendrier xios
289	!Temps julian => an, mois, jour, heure
290	CALL ju2ymds(zjulian, x_an, x_mois, x_jour, x_heure)
291	CALL ju2ymds(zjulian_start, ini_an, ini_mois, ini_jour, ini_heure)
292	CALL wxios_set_cal(dtime, calend, x_an, x_mois, x_jour, x_heure, ini_an, &
293	ini_mois, ini_jour, ini_heure )
294	#endif
295
296	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
297	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
298	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
299	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
300
301	zdtime_moy = dtime ! Frequence ou l on moyenne
302
303
304	ecrit_files(7) = ecrit_files(1)
305	ecrit_files(8) = ecrit_files(2)
306	ecrit_files(9) = ecrit_files(3)
307
308	DO iff=1,nfiles
309
310	! Calculate ecrit_files for all files
311	IF ( chtimestep(iff).eq.'Default' ) THEN
312	! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf
313	! ...)*86400.
314	ecrit_files(iff)=ecrit_files(iff)*86400.
315	ELSE IF (chtimestep(iff).eq.'-1') THEN
316	PRINT*,'ecrit_files(',iff,') < 0 so IOIPSL work on different'
317	PRINT*,'months length'
318	ecrit_files(iff)=-1.
319	ELSE
320	CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
321	ENDIF
322
323	WRITE(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
324	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
325
326
327	#ifdef CPP_XIOS
328	!!! Ouverture de chaque fichier XIOS !!!!!!!!!!!
329	IF (.not. ok_all_xml) THEN
330	IF (prt_level >= 10) THEN
331	print*,'phys_output_open: call wxios_add_file with phys_out_filenames(iff)=',trim(phys_out_filenames(iff))
332	ENDIF
333	CALL wxios_add_file(phys_out_filenames(iff),chtimestep(iff),lev_files(iff))
334	ENDIF
335
336	!!! Declaration des axes verticaux de chaque fichier:
337	IF (prt_level >= 10) THEN
338	print*,'phys_output_open: Declare vertical axes for each file'
339	ENDIF
340
341	IF (iff.LE.6.OR.iff.EQ.10) THEN
342	CALL wxios_add_vaxis("presnivs", &
343	levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
344	CALL wxios_add_vaxis("Ahyb", &
345	levmax(iff) - levmin(iff) + 1, aps(levmin(iff):levmax(iff)), positif='down', &
346	bnds=Ahyb_bounds(levmin(iff):levmax(iff),:))
347	CALL wxios_add_vaxis("Bhyb", &
348	levmax(iff) - levmin(iff) + 1, bps(levmin(iff):levmax(iff)), positif='down', &
349	bnds=Bhyb_bounds(levmin(iff):levmax(iff),:))
350	CALL wxios_add_vaxis("klev", levmax(iff) - levmin(iff) + 1, &
351	lev_index(levmin(iff):levmax(iff)))
352	CALL wxios_add_vaxis("bnds", 2, (/1.,2./))
353
354	CALL wxios_add_vaxis("Alt", &
355	levmax(iff) - levmin(iff) + 1, pseudoalt)
356	ELSE
357	! NMC files
358	CALL wxios_add_vaxis("plev", &
359	levmax(iff) - levmin(iff) + 1, rlevSTD(levmin(iff):levmax(iff)))
360	ENDIF
361	#endif
362
363	IF (clef_files(iff)) THEN
364	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
365	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
366	IF (phys_out_regfkey(iff)) THEN
367	imin_ins=1
368	imax_ins=nbp_lon
369	jmin_ins=1
370	jmax_ins=jjmp1
371
372	! correction abderr
373	DO i=1,nbp_lon
374	WRITE(lunout,*)'io_lon(i)=',io_lon(i)
375	IF (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
376	IF (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
377	ENDDO
378
379	DO j=1,jjmp1
380	WRITE(lunout,*)'io_lat(j)=',io_lat(j)
381	IF (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
382	IF (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
383	ENDDO
384
385	WRITE(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
386	imin_ins,imax_ins,jmin_ins,jmax_ins
387	WRITE(lunout,*)'longitudes : ', &
388	io_lon(imin_ins),io_lon(imax_ins), &
389	'latitudes : ', &
390	io_lat(jmax_ins),io_lat(jmin_ins)
391
392	CALL histbeg(phys_out_filenames(iff),nbp_lon,io_lon,jjmp1,io_lat, &
393	imin_ins,imax_ins-imin_ins+1, &
394	jmin_ins,jmax_ins-jmin_ins+1, &
395	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
396	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
397	!IM fichiers stations
398	ELSE IF (clef_stations(iff)) THEN
399
400	IF (prt_level >= 10) THEN
401	WRITE(lunout,*)'phys_output_open: iff=',iff,' phys_out_filenames(iff)=',phys_out_filenames(iff)
402	ENDIF
403
404	CALL histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
405	phys_out_filenames(iff), &
406	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
407	ELSE
408
409	IF (prt_level >= 10) THEN
410	WRITE(lunout,*)'phys_output_open: iff=',iff,' phys_out_filenames(iff)=',phys_out_filenames(iff)
411	ENDIF
412
413	CALL histbeg_phy_all(phys_out_filenames(iff),itau_phy,zjulian,&
414	dtime,nhorim(iff),nid_files(iff))
415	ENDIF
416
417	#ifndef CPP_IOIPSL_NO_OUTPUT
418	IF (iff.LE.6.OR.iff.EQ.10) THEN
419	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &
420	levmax(iff) - levmin(iff) + 1, &
421	presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")
422	!!!! Composantes de la coordonnee sigma-hybride
423	CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
424	levmax(iff) - levmin(iff) + 1,aps,nvertap(iff))
425
426	CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
427	levmax(iff) - levmin(iff) + 1,bps,nvertbp(iff))
428
429	CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &
430	levmax(iff) - levmin(iff) + 1,pseudoalt,nvertAlt(iff))
431
432	ELSE
433	! NMC files
434	CALL histvert(nid_files(iff), "plev", "pressure", "Pa", &
435	levmax(iff) - levmin(iff) + 1, &
436	rlevSTD(levmin(iff):levmax(iff)), nvertm(iff), "down")
437	ENDIF
438	#endif
439
440	ENDIF ! clef_files
441
442	IF (nqtot>=nqo+1) THEN
443	!
444	DO iq=nqo+1,nqtot
445	iiq=niadv(iq)
446	o_trac(iq-nqo) = ctrl_out((/ 1, 5, 5, 5, 10, 10, 11, 11, 11, 11 /), &
447	tname(iiq),'Tracer '//ttext(iiq), "-", &
448	(/ '', '', '', '', '', '', '', '', '', '' /))
449	o_dtr_vdf(iq-nqo) = ctrl_out((/ 4, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
450	'd'//trim(tname(iq))//'_vdf', &
451	'Tendance tracer '//ttext(iiq), "-" , &
452	(/ '', '', '', '', '', '', '', '', '', '' /))
453
454	o_dtr_the(iq-nqo) = ctrl_out((/ 5, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
455	'd'//trim(tname(iq))//'_the', &
456	'Tendance tracer '//ttext(iiq), "-", &
457	(/ '', '', '', '', '', '', '', '', '', '' /))
458
459	o_dtr_con(iq-nqo) = ctrl_out((/ 5, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
460	'd'//trim(tname(iq))//'_con', &
461	'Tendance tracer '//ttext(iiq), "-", &
462	(/ '', '', '', '', '', '', '', '', '', '' /))
463
464	o_dtr_lessi_impa(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
465	'd'//trim(tname(iq))//'_lessi_impa', &
466	'Tendance tracer '//ttext(iiq), "-", &
467	(/ '', '', '', '', '', '', '', '', '', '' /))
468
469	o_dtr_lessi_nucl(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
470	'd'//trim(tname(iq))//'_lessi_nucl', &
471	'Tendance tracer '//ttext(iiq), "-", &
472	(/ '', '', '', '', '', '', '', '', '', '' /))
473
474	o_dtr_insc(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
475	'd'//trim(tname(iq))//'_insc', &
476	'Tendance tracer '//ttext(iiq), "-", &
477	(/ '', '', '', '', '', '', '', '', '', '' /))
478
479	o_dtr_bcscav(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
480	'd'//trim(tname(iq))//'_bcscav', &
481	'Tendance tracer '//ttext(iiq), "-", &
482	(/ '', '', '', '', '', '', '', '', '', '' /))
483
484	o_dtr_evapls(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
485	'd'//trim(tname(iq))//'_evapls', &
486	'Tendance tracer '//ttext(iiq), "-", &
487	(/ '', '', '', '', '', '', '', '', '', '' /))
488
489	o_dtr_ls(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
490	'd'//trim(tname(iq))//'_ls', &
491	'Tendance tracer '//ttext(iiq), "-", &
492	(/ '', '', '', '', '', '', '', '', '', '' /))
493
494	o_dtr_trsp(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
495	'd'//trim(tname(iq))//'_trsp', &
496	'Tendance tracer '//ttext(iiq), "-", &
497	(/ '', '', '', '', '', '', '', '', '', '' /))
498
499	o_dtr_sscav(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
500	'd'//trim(tname(iq))//'_sscav', &
501	'Tendance tracer '//ttext(iiq), "-", &
502	(/ '', '', '', '', '', '', '', '', '', '' /))
503
504	o_dtr_sat(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
505	'd'//trim(tname(iq))//'_sat', &
506	'Tendance tracer '//ttext(iiq), "-", &
507	(/ '', '', '', '', '', '', '', '', '', '' /))
508
509	o_dtr_uscav(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
510	'd'//trim(tname(iq))//'_uscav', &
511	'Tendance tracer '//ttext(iiq), "-", &
512	(/ '', '', '', '', '', '', '', '', '', '' /))
513
514	o_dtr_dry(iq-nqo) = ctrl_out((/ 7, 7, 7, 7, 10, 10, 11, 11, 11, 11 /), &
515	'cum'//'d'//trim(tname(iq))//'_dry', &
516	'tracer tendency dry deposition'//ttext(iiq), "-", &
517	(/ '', '', '', '', '', '', '', '', '', '' /))
518
519	o_trac_cum(iq-nqo) = ctrl_out((/ 1, 4, 10, 10, 10, 10, 11, 11, 11, 11 /), &
520	'cum'//tname(iiq),&
521	'Cumulated tracer '//ttext(iiq), "-", &
522	(/ '', '', '', '', '', '', '', '', '', '' /))
523	ENDDO
524	ENDIF
525
526	ENDDO ! iff
527
528	! Updated write frequencies due to phys_out_filetimesteps.
529	! Write frequencies are now in seconds.
530	ecrit_mth = ecrit_files(1)
531	ecrit_day = ecrit_files(2)
532	ecrit_hf = ecrit_files(3)
533	ecrit_ins = ecrit_files(4)
534	ecrit_LES = ecrit_files(5)
535	ecrit_ins = ecrit_files(6)
536
537	IF (prt_level >= 10) THEN
538	WRITE(lunout,*)'swaero_diag=',swaero_diag
539	WRITE(lunout,*)'dryaod_diag=',dryaod_diag
540	WRITE(lunout,*)'ok_4xCO2atm=',ok_4xCO2atm
541	WRITE(lunout,*)'phys_output_open: ends here'
542	ENDIF
543
544	END SUBROUTINE phys_output_open
545
546	SUBROUTINE convers_timesteps(str,dtime,timestep)
547
548	use ioipsl
549	USE phys_cal_mod
550	USE time_phylmdz_mod, ONLY: day_ref, annee_ref
551	USE print_control_mod, ONLY: lunout
552
553	IMPLICIT NONE
554
555	CHARACTER(LEN=20) :: str
556	CHARACTER(LEN=10) :: type
557	INTEGER :: ipos,il
558	real :: ttt,xxx,timestep,dayseconde,dtime
559	parameter (dayseconde=86400.)
560
561	ipos=scan(str,'0123456789.',.TRUE.)
562	!
563	il=len_trim(str)
564	WRITE(lunout,*) "ipos = ", ipos
565	WRITE(lunout,*) "il = ", il
566	IF (ipos == 0) CALL abort_physic("convers_timesteps", "bad str", 1)
567	read(str(1:ipos),*) ttt
568	WRITE(lunout,*)ttt
569	type=str(ipos+1:il)
570
571	IF ( il == ipos ) THEN
572	type='day'
573	ENDIF
574
575	IF ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde
576	IF ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) THEN
577	WRITE(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,mth_len
578	timestep = ttt * dayseconde * mth_len
579	ENDIF
580	IF ( type == 'hours'.or.type == 'hr'.or.type == 'heurs') timestep = ttt * dayseconde / 24.
581	IF ( type == 'mn'.or.type == 'minutes' ) timestep = ttt * 60.
582	IF ( type == 's'.or.type == 'sec'.or.type == 'secondes' ) timestep = ttt
583	IF ( type == 'TS' ) timestep = ttt * dtime
584
585	WRITE(lunout,*)'type = ',type
586	WRITE(lunout,*)'nb j/h/m = ',ttt
587	WRITE(lunout,*)'timestep(s)=',timestep
588
589	END SUBROUTINE convers_timesteps
590
591	END MODULE phys_output_mod

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

Download in other formats:

Original Format