Context Navigation

phys_output_mod.F90 @ 2002

Last change on this file since 2002 was 2002, checked in by Ehouarn Millour, 10 years ago
Further cleanup concerning XIOS (mainly about axes being defined as axes and not as groups of axes). EM
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: 23.1 KB

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

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

Download in other formats:

Original Format