Context Navigation

phys_output_mod.F90 @ 2236

Last change on this file since 2236 was 2167, checked in by musat, 10 years ago
Correction sorties fichier 6heures en "inst" au lieu de "ave". Changement noms par defaut des fichiers de sorties: histmth, histday, histhf6h, histhf3h, histhf3hm FH/IM
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.5 KB

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

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

Download in other formats:

Original Format