Context Navigation

phys_output_mod.F90 @ 2031

Last change on this file since 2031 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

Rev	Line
[1279]	1	! $Id: phys_output_mod.F90 2002 2014-04-04 12:39:54Z 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, &
[1828]	28	jjmp1,nlevSTD,clevSTD,rlevSTD,nbteta, &
[1539]	29	ctetaSTD, dtime, ok_veget, &
[1279]	30	type_ocean, iflag_pbl,ok_mensuel,ok_journe, &
	31	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
[1539]	32	phys_out_filestations, &
[1807]	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)
[907]	36
[1562]	37	USE iophy
	38	USE dimphy
	39	USE infotrac
	40	USE ioipsl
[1742]	41	USE phys_cal_mod, only : hour
[1562]	42	USE mod_phys_lmdz_para
	43	USE aero_mod, only : naero_spc,name_aero
[1865]	44	!Martin
	45	USE surface_data, ONLY : ok_snow
[1807]	46	USE phys_output_ctrlout_mod
[1912]	47	USE mod_grid_phy_lmdz, only: klon_glo
[907]	48
[1825]	49	#ifdef CPP_XIOS
	50	! ug Pour les sorties XIOS
	51	USE wxios
	52	#endif
	53
[1562]	54	IMPLICIT NONE
	55	include "dimensions.h"
	56	include "temps.h"
	57	include "clesphys.h"
	58	include "thermcell.h"
	59	include "comvert.h"
[1575]	60	include "iniprint.h"
[907]	61
[1807]	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
[1539]	78	INTEGER, DIMENSION(pim) :: tabij
[1807]	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
[1539]	82
[1791]	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
[1807]	88	INTEGER, INTENT(IN):: read_climoz ! read ozone climatology
[1562]	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
[1279]	94
[1791]	95	REAL :: dtime
	96	INTEGER :: idayref
	97	REAL :: zjulian
	98	REAL, DIMENSION(klev) :: Ahyb, Bhyb, Alt
	99	CHARACTER(LEN=4), DIMENSION(nlevSTD) :: clevSTD
[1828]	100	REAL, DIMENSION(nlevSTD) :: rlevSTD
[1791]	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
[1884]	115	CHARACTER(LEN=20), DIMENSION(nfiles) :: chtimestep = (/ 'Default', 'Default', 'Default', 'Default', 'Default', &
	116	'Default', 'Default', 'Default', 'Default' /)
[1791]	117	LOGICAL, DIMENSION(nfiles) :: phys_out_filekeys
	118	LOGICAL, DIMENSION(nfiles) :: phys_out_filestations
[907]	119
[1065]	120	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	121	! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]
[1054]	122
[1828]	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. /)
[1825]	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
[1791]	139	WRITE(lunout,*) 'Debut phys_output_mod.F90'
[1562]	140	! Initialisations (Valeurs par defaut
[1403]	141
[1791]	142	IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
	143	IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot))
[1813]	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))
[1403]	149
[1813]	150
[1828]	151	levmax = (/ klev, klev, klev, klev, klev, klev, nlevSTD, nlevSTD, nlevSTD /)
[1065]	152
[1562]	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'
[1828]	159	phys_out_filenames(7) = 'histmthNMC'
	160	phys_out_filenames(8) = 'histdayNMC'
	161	phys_out_filenames(9) = 'histhfNMC.nc'
[907]	162
[1562]	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)'
[1828]	169	type_ecri(7) = 'inst(X)'
	170	type_ecri(8) = 'inst(X)'
	171	type_ecri(9) = 'inst(X)'
[907]	172
[1562]	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
[1828]	179	clef_files(7) = ok_histNMC(1)
	180	clef_files(8) = ok_histNMC(2)
	181	clef_files(9) = ok_histNMC(3)
[907]	182
[1562]	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.
[1828]	190	clef_stations(7) = .FALSE.
	191	clef_stations(8) = .FALSE.
	192	clef_stations(9) = .FALSE.
[1539]	193
[1562]	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
[1828]	200	lev_files(7) = levout_histNMC(1)
	201	lev_files(8) = levout_histNMC(2)
	202	lev_files(9) = levout_histNMC(3)
[907]	203
[1562]	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
[1912]	210	ecrit_files(7) = freq_outNMC(1)
	211	ecrit_files(8) = freq_outNMC(2)
	212	ecrit_files(9) = freq_outNMC(3)
[1279]	213
[1562]	214	!! Lectures des parametres de sorties dans physiq.def
[1279]	215
[1791]	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)
[1562]	221	phys_out_levmin(:)=levmin(:)
[1791]	222	CALL getin('phys_out_levmin',levmin)
[1562]	223	phys_out_levmax(:)=levmax(:)
[1791]	224	CALL getin('phys_out_levmax',levmax)
	225	CALL getin('phys_out_filenames',phys_out_filenames)
[1562]	226	phys_out_filekeys(:)=clef_files(:)
[1791]	227	CALL getin('phys_out_filekeys',clef_files)
[1562]	228	phys_out_filestations(:)=clef_stations(:)
[1791]	229	CALL getin('phys_out_filestations',clef_stations)
[1562]	230	phys_out_filelevels(:)=lev_files(:)
[1791]	231	CALL getin('phys_out_filelevels',lev_files)
	232	CALL getin('phys_out_filetimesteps',chtimestep)
[1562]	233	phys_out_filetypes(:)=type_ecri(:)
[1791]	234	CALL getin('phys_out_filetypes',type_ecri)
[1279]	235
[1562]	236	type_ecri_files(:)=type_ecri(:)
[1279]	237
[1791]	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
[1562]	247
[1971]	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
[1825]	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
[907]	268	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	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 ..
[907]	271	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	272
[1807]	273	zdtime_moy = dtime ! Frequence ou l on moyenne
[1279]	274
[1562]	275	! Calcul des Ahyb, Bhyb et Alt
[1791]	276	DO k=1,klev
[1562]	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.
[1791]	280	ENDDO
[1562]	281	! if(prt_level.ge.1) then
[1791]	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)
[1562]	285	! endif
[1971]	286	DO iff=1,nfiles
[929]	287
[1641]	288	! Calculate ecrit_files for all files
[1884]	289	IF ( chtimestep(iff).eq.'Default' ) then
	290	! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf
	291	! ...)*86400.
[1641]	292	ecrit_files(iff)=ecrit_files(iff)*86400.
[1912]	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.
[1884]	297	else
	298	CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
	299	ENDIF
	300
[1791]	301	WRITE(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
[1641]	302	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
[1279]	303
[1971]	304	IF (clef_files(iff)) THEN
	305
[1852]	306	#ifdef CPP_XIOS
[1971]	307	!!! Ouverture de chaque fichier XIOS !!!!!!!!!!!
[2001]	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
[1971]	311	CALL wxios_add_file(phys_out_filenames(iff),chtimestep(iff),lev_files(iff))
	312
	313	!!! Declaration des axes verticaux de chaque fichier:
[2001]	314	if (prt_level >= 10) then
	315	print*,'phys_output_open: Declare vertical axes for each file'
	316	endif
[1971]	317	if (iff.le.6) then
[2002]	318	CALL wxios_add_vaxis("presnivs", &
[1852]	319	levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
[2002]	320	CALL wxios_add_vaxis("Ahyb", &
[1852]	321	levmax(iff) - levmin(iff) + 1, Ahyb)
[2002]	322	CALL wxios_add_vaxis("Bhyb", &
[1852]	323	levmax(iff) - levmin(iff) + 1, Bhyb)
[2002]	324	CALL wxios_add_vaxis("Alt", &
[1852]	325	levmax(iff) - levmin(iff) + 1, Alt)
[1971]	326	else
[2001]	327	! NMC files
[2002]	328	CALL wxios_add_vaxis("plev", &
[1971]	329	levmax(iff) - levmin(iff) + 1, rlevSTD(levmin(iff):levmax(iff)))
	330	endif
[1852]	331	#endif
	332
[907]	333	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
	334	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1791]	335	IF (phys_out_regfkey(iff)) then
[1562]	336	imin_ins=1
	337	imax_ins=iim
	338	jmin_ins=1
	339	jmax_ins=jjmp1
[907]	340
[1562]	341	! correction abderr
	342	do i=1,iim
[1791]	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
[1562]	346	enddo
[907]	347
[1562]	348	do j=1,jjmp1
[1791]	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
[1562]	352	enddo
[907]	353
[1791]	354	WRITE(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
[1562]	355	imin_ins,imax_ins,jmin_ins,jmax_ins
[1791]	356	WRITE(lunout,*)'longitudes : ', &
[1562]	357	io_lon(imin_ins),io_lon(imax_ins), &
	358	'latitudes : ', &
	359	io_lat(jmax_ins),io_lat(jmin_ins)
[907]	360
[1562]	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))
[907]	365	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	366	!IM fichiers stations
[1791]	367	else IF (clef_stations(iff)) THEN
[1539]	368
[2001]	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
[1791]	373	CALL histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
[1562]	374	phys_out_filenames(iff), &
	375	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
	376	else
[1825]	377	CALL histbeg_phy_all(phys_out_filenames(iff),itau_phy,zjulian,&
[1852]	378	dtime,nhorim(iff),nid_files(iff))
[1562]	379	endif
[907]	380
[1852]	381	#ifndef CPP_NO_IOIPSL
[1828]	382	if (iff.le.6) then
	383	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &
[1562]	384	levmax(iff) - levmin(iff) + 1, &
	385	presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")
[1807]	386	!!!! Composantes de la coordonnee sigma-hybride
[1562]	387	CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
	388	levmax(iff) - levmin(iff) + 1,Ahyb,nvertap(iff))
[907]	389
[1562]	390	CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
	391	levmax(iff) - levmin(iff) + 1,Bhyb,nvertbp(iff))
[1279]	392
[1828]	393	CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &
[1562]	394	levmax(iff) - levmin(iff) + 1,Alt,nvertAlt(iff))
[1279]	395
[1828]	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
[907]	402
[1791]	403	IF (nqtot>=3) THEN
	404	DO iq=3,nqtot
	405	iiq=niadv(iq)
[1828]	406	o_trac(iq-2) = ctrl_out((/ 4, 5, 1, 1, 1, 10, 11, 11, 11 /), &
	407	tname(iiq),'Tracer '//ttext(iiq), "-", &
	408	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	409
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	414
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	419
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	424
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	429
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	434
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	439
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	444
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	449
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	454
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	459
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	464
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	469
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	474
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1813]	479
[1828]	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	(/ '', '', '', '', '', '', '', '', '' /))
[1791]	484	ENDDO
	485	ENDIF
[1539]	486
[1971]	487	ENDIF ! clef_files
[1539]	488
[1971]	489	ENDDO ! iff
[1539]	490
[1641]	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
[2001]	500	if (prt_level >= 10) then
	501	WRITE(lunout,*)'swaero_diag=',swaero_diag
	502	WRITE(lunout,*)'phys_output_open: ends here'
	503	endif
[1920]	504
[1791]	505	end SUBROUTINE phys_output_open
[907]	506
	507
	508
[1562]	509	SUBROUTINE convers_timesteps(str,dtime,timestep)
[1279]	510
[1562]	511	use ioipsl
	512	USE phys_cal_mod
[1279]	513
[1562]	514	IMPLICIT NONE
[1279]	515
[1791]	516	CHARACTER(LEN=20) :: str
	517	CHARACTER(LEN=10) :: type
	518	INTEGER :: ipos,il
[1562]	519	real :: ttt,xxx,timestep,dayseconde,dtime
	520	parameter (dayseconde=86400.)
	521	include "temps.h"
	522	include "comconst.h"
[1575]	523	include "iniprint.h"
[1279]	524
[1791]	525	ipos=scan(str,'0123456789.',.TRUE.)
[1562]	526	!
	527	il=len_trim(str)
[1937]	528	WRITE(lunout,*) "ipos = ", ipos
	529	WRITE(lunout,*) "il = ", il
	530	if (ipos == 0) call abort_gcm("convers_timesteps", "bad str", 1)
[1562]	531	read(str(1:ipos),*) ttt
[1791]	532	WRITE(lunout,*)ttt
[1562]	533	type=str(ipos+1:il)
[1279]	534
	535
[1791]	536	IF ( il == ipos ) then
[1562]	537	type='day'
	538	endif
[1279]	539
[1791]	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
[1562]	543	timestep = ttt * dayseconde * mth_len
	544	endif
[1791]	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
[1279]	549
[1791]	550	WRITE(lunout,*)'type = ',type
	551	WRITE(lunout,*)'nb j/h/m = ',ttt
	552	WRITE(lunout,*)'timestep(s)=',timestep
[1279]	553
[1562]	554	END SUBROUTINE convers_timesteps
[1279]	555
[907]	556	END MODULE phys_output_mod
	557

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

Download in other formats:

Original Format