Context Navigation

phys_output_mod.F90 @ 5411

Last change on this file since 5411 was 5394, checked in by evignon, 12 days ago
ecriture de certaines declarations sur plusieurs lignes par souci de clarté
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 `Id`
File size: 30.4 KB

Rev	Line
[1279]	1	! $Id: phys_output_mod.F90 5394 2024-12-06 13:10:59Z evignon $
	2	!
[1813]	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
[1920]	8	REAL, DIMENSION(nfiles),SAVE :: ecrit_files
[1813]	9
[907]	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
[5357]	17	! AI. nov 2024 : Modifs pour rajouter plus de choix pour la frequence temporelle d'archivage
	18	! dans les fichiers de sorties (avec IOIPSL) :
	19	! month(s),m,mth,mois,m,day(s),d,jour(s),j, ...
[907]	20	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	21
[1562]	22	CONTAINS
[907]	23
	24	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	25	!!!!!!!!! Ouverture des fichier et definition des variable de sortie !!!!!!!!
[1562]	26	!! histbeg, histvert et histdef
[907]	27	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	28
[1539]	29	SUBROUTINE phys_output_open(rlon,rlat,pim,tabij,ipt,jpt,plon,plat, &
[2037]	30	jjmp1,nlevSTD,clevSTD,rlevSTD, dtime, ok_veget, &
[2159]	31	type_ocean, iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
[1279]	32	ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, read_climoz, &
[1539]	33	phys_out_filestations, &
[3630]	34	aerosol_couple, flag_aerosol_strat, &
[1807]	35	pdtphys, paprs, pphis, pplay, lmax_th, ptconv, ptconvth, ivap, &
[2665]	36	d_u, d_t, qx, d_qx, zmasse, ok_sync)
[907]	37
[1562]	38	USE iophy
	39	USE dimphy
[5056]	40	USE infotrac_phy, ONLY: nqtot, tracers, niso, ntraciso=>ntiso
[4120]	41	USE strings_mod, ONLY: maxlen
[1562]	42	USE ioipsl
[2358]	43	USE phys_cal_mod, only : hour, calend
[1562]	44	USE mod_phys_lmdz_para
[2704]	45	!Martin
[3792]	46	USE surface_data, ONLY : landice_opt
[1807]	47	USE phys_output_ctrlout_mod
[2344]	48	USE mod_grid_phy_lmdz, only: klon_glo,nbp_lon,nbp_lat
[2311]	49	USE print_control_mod, ONLY: prt_level,lunout
[4228]	50	USE vertical_layers_mod, ONLY: ap,bp,preff,presnivs, aps, bps, pseudoalt, presinter
[2344]	51	USE time_phylmdz_mod, ONLY: day_ini, itau_phy, start_time, annee_ref, day_ref
[1825]	52	! ug Pour les sorties XIOS
[5310]	53	use wxios_mod
[4998]	54	USE infotrac_phy, ONLY: nbtr_bin
[5056]	55	#ifdef ISO
	56	USE isotopes_mod, ONLY: isoName,iso_HTO
	57	#ifdef ISOTRAC
	58	use isotrac_mod, only: index_zone,index_iso,strtrac
	59	#endif
	60	#endif
	61
[5282]	62	USE clesphys_mod_h
[5252]	63	USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_STRATAER
	64
[5285]	65	USE yomcst_mod_h
[5274]	66	IMPLICIT NONE
[907]	67
[5274]	68
[2110]	69	! ug Nouveaux arguments n\'ecessaires au histwrite_mod:
[1807]	70	INTEGER, INTENT(IN) :: ivap
	71	INTEGER, DIMENSION(klon), INTENT(IN) :: lmax_th
	72	LOGICAL, INTENT(IN) :: ok_sync
	73	LOGICAL, DIMENSION(klon, klev), INTENT(IN) :: ptconv, ptconvth
	74	REAL, INTENT(IN) :: pdtphys
	75	REAL, DIMENSION(klon), INTENT(IN) :: pphis
[2665]	76	REAL, DIMENSION(klon, klev), INTENT(IN) :: pplay, d_u, d_t
[1807]	77	REAL, DIMENSION(klon, klev+1), INTENT(IN) :: paprs
	78	REAL, DIMENSION(klon,klev,nqtot), INTENT(IN):: qx, d_qx
[2344]	79	REAL, DIMENSION(klon, klev), INTENT(IN) :: zmasse
[1807]	80
	81
	82	REAL,DIMENSION(klon),INTENT(IN) :: rlon
	83	REAL,DIMENSION(klon),INTENT(IN) :: rlat
	84	INTEGER, INTENT(IN) :: pim
[1539]	85	INTEGER, DIMENSION(pim) :: tabij
[1807]	86	INTEGER,DIMENSION(pim), INTENT(IN) :: ipt, jpt
	87	REAL,DIMENSION(pim), INTENT(IN) :: plat, plon
	88	REAL,DIMENSION(pim,2) :: plat_bounds, plon_bounds
[1539]	89
[1791]	90	INTEGER :: jjmp1
[2037]	91	INTEGER :: nlevSTD, radpas
[1791]	92	LOGICAL :: ok_mensuel, ok_journe, ok_hf, ok_instan
[2530]	93	LOGICAL :: ok_LES,ok_ade,ok_aie
	94	INTEGER :: flag_aerosol_strat
[3630]	95	LOGICAL :: aerosol_couple
[1807]	96	INTEGER, INTENT(IN):: read_climoz ! read ozone climatology
[1562]	97	! Allowed values are 0, 1 and 2
	98	! 0: do not read an ozone climatology
	99	! 1: read a single ozone climatology that will be used day and night
	100	! 2: read two ozone climatologies, the average day and night
	101	! climatology and the daylight climatology
[1279]	102
[1791]	103	REAL :: dtime
	104	INTEGER :: idayref
[2095]	105	REAL :: zjulian_start, zjulian
[1791]	106	CHARACTER(LEN=4), DIMENSION(nlevSTD) :: clevSTD
[1828]	107	REAL, DIMENSION(nlevSTD) :: rlevSTD
[4998]	108	INTEGER :: nsrf, k, iq, iff, i, j, ilev, itr, itrb, ixt, iiso, izone
[1791]	109	INTEGER :: naero
	110	LOGICAL :: ok_veget
	111	INTEGER :: iflag_pbl
[2159]	112	INTEGER :: iflag_pbl_split
[1791]	113	CHARACTER(LEN=4) :: bb2
	114	CHARACTER(LEN=2) :: bb3
	115	CHARACTER(LEN=6) :: type_ocean
[2344]	116	INTEGER, DIMENSION(nbp_lon*jjmp1) :: ndex2d
	117	INTEGER, DIMENSION(nbp_lonjjmp1klev) :: ndex3d
[1791]	118	INTEGER :: imin_ins, imax_ins
	119	INTEGER :: jmin_ins, jmax_ins
	120	INTEGER, DIMENSION(nfiles) :: phys_out_levmin, phys_out_levmax
	121	INTEGER, DIMENSION(nfiles) :: phys_out_filelevels
[1884]	122	CHARACTER(LEN=20), DIMENSION(nfiles) :: chtimestep = (/ 'Default', 'Default', 'Default', 'Default', 'Default', &
[2752]	123	'Default', 'Default', 'Default', 'Default', 'Default' /)
[1791]	124	LOGICAL, DIMENSION(nfiles) :: phys_out_filekeys
	125	LOGICAL, DIMENSION(nfiles) :: phys_out_filestations
[907]	126
[5056]	127	#ifdef ISO
	128	CHARACTER(LEN=maxlen) :: outiso
	129	CHARACTER(LEN=20) :: unit
	130	#endif
[4046]	131	CHARACTER(LEN=maxlen) :: tnam, lnam, dn
	132	INTEGER :: flag(nfiles)
	133
[1065]	134	!!!!!!!!!! stockage dans une region limitee pour chaque fichier !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	135	! entre [phys_out_lonmin,phys_out_lonmax] et [phys_out_latmin,phys_out_latmax]
[5394]	136	LOGICAL, DIMENSION(nfiles), SAVE :: phys_out_regfkey = [.FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE., .FALSE.]
	137	REAL, DIMENSION(nfiles), SAVE :: phys_out_lonmin = [ -180., -180., -180., -180., -180., -180., -180., -180., -180., -180.]
	138	REAL, DIMENSION(nfiles), SAVE :: phys_out_lonmax = [ 180., 180., 180., 180., 180., 180., 180., 180., 180., 180.]
	139	REAL, DIMENSION(nfiles), SAVE :: phys_out_latmin = [ -90., -90., -90., -90., -90., -90., -90., -90., -90., -90.]
	140	REAL, DIMENSION(nfiles), SAVE :: phys_out_latmax = [ 90., 90., 90., 90., 90., 90., 90., 90., 90., 90.]
[3003]	141	REAL, DIMENSION(klev,2) :: Ahyb_bounds, Bhyb_bounds
[3474]	142	REAL, DIMENSION(klev+1) :: lev_index
[3003]	143
[2110]	144	! ug Variables utilis\'ees pour r\'ecup\'erer le calendrier pour xios
[1825]	145	INTEGER :: x_an, x_mois, x_jour
	146	REAL :: x_heure
[2095]	147	INTEGER :: ini_an, ini_mois, ini_jour
	148	REAL :: ini_heure
[4619]	149
[3048]	150	INTEGER :: ISW
	151	REAL, DIMENSION(NSW) :: wl1_sun, wl2_sun !wavelength bounds (in um) for SW
	152	REAL, DIMENSION(NSW) :: wn1_sun, wn2_sun !wavenumber bounds (in m-1) for SW
	153	REAL, DIMENSION(NSW) :: spectband !mean wavenumb. of each sp.band
	154	REAL, DIMENSION(NSW,2) :: spbnds_sun !bounds of spectband
[1825]	155
[1791]	156	WRITE(lunout,*) 'Debut phys_output_mod.F90'
[1562]	157	! Initialisations (Valeurs par defaut
[1403]	158
[3003]	159	DO ilev=1,klev
	160	Ahyb_bounds(ilev,1) = ap(ilev)
	161	Ahyb_bounds(ilev,2) = ap(ilev+1)
	162	Bhyb_bounds(ilev,1) = bp(ilev)
	163	Bhyb_bounds(ilev,2) = bp(ilev+1)
	164	lev_index(ilev) = REAL(ilev)
	165	END DO
[3474]	166	lev_index(klev+1) = REAL(klev+1)
[3003]	167
[1791]	168	IF (.NOT. ALLOCATED(o_trac)) ALLOCATE(o_trac(nqtot))
	169	IF (.NOT. ALLOCATED(o_trac_cum)) ALLOCATE(o_trac_cum(nqtot))
[1813]	170	ALLOCATE(o_dtr_the(nqtot),o_dtr_con(nqtot),o_dtr_lessi_impa(nqtot))
	171	ALLOCATE(o_dtr_lessi_nucl(nqtot),o_dtr_insc(nqtot),o_dtr_bcscav(nqtot))
	172	ALLOCATE(o_dtr_evapls(nqtot),o_dtr_ls(nqtot),o_dtr_trsp(nqtot))
	173	ALLOCATE(o_dtr_sscav(nqtot),o_dtr_sat(nqtot),o_dtr_uscav(nqtot))
	174	ALLOCATE(o_dtr_dry(nqtot),o_dtr_vdf(nqtot))
[5252]	175	IF (CPPKEY_STRATAER) THEN
[4998]	176	ALLOCATE(o_nd_mode(nbtr_bin),o_sulfmmr_mode(nbtr_bin))
[5252]	177	END IF
[5056]	178	#ifdef ISO
	179	ALLOCATE(o_xtprecip(ntraciso))
	180	ALLOCATE(o_xtplul(ntraciso))
	181	ALLOCATE(o_xtpluc(ntraciso))
	182	ALLOCATE(o_xtevap(ntraciso))
	183	ALLOCATE(o_xtevap_srf(ntraciso,4))
	184	ALLOCATE(o_xtovap(ntraciso))
	185	ALLOCATE(o_xtoliq(ntraciso))
	186	ALLOCATE(o_xtcond(ntraciso))
	187	ALLOCATE(o_xtrunoff_diag(ntraciso))
	188	ALLOCATE(o_dxtdyn(ntraciso))
	189	ALLOCATE(o_dxtldyn(ntraciso))
	190	ALLOCATE(o_dxtcon(ntraciso))
	191	ALLOCATE(o_dxtlsc(ntraciso))
	192	ALLOCATE(o_dxteva(ntraciso))
	193	ALLOCATE(o_dxtajs(ntraciso))
	194	ALLOCATE(o_dxtvdf(ntraciso))
	195	ALLOCATE(o_dxtthe(ntraciso))
	196	ALLOCATE(o_dxtch4(ntraciso))
	197	if (iso_HTO.gt.0) then
	198	ALLOCATE(o_dxtprod_nucl(ntraciso))
	199	ALLOCATE(o_dxtcosmo(ntraciso))
	200	ALLOCATE(o_dxtdecroiss(ntraciso))
	201	endif
	202	#endif
	203
[4056]	204	levmax = [klev, klev, klev, klev, klev, klev, nlevSTD, nlevSTD, nlevSTD, klev]
[1813]	205
[1562]	206	phys_out_filenames(1) = 'histmth'
	207	phys_out_filenames(2) = 'histday'
[2167]	208	phys_out_filenames(3) = 'histhf6h'
	209	phys_out_filenames(4) = 'histhf3h'
	210	phys_out_filenames(5) = 'histhf3hm'
[1562]	211	phys_out_filenames(6) = 'histstn'
[1828]	212	phys_out_filenames(7) = 'histmthNMC'
	213	phys_out_filenames(8) = 'histdayNMC'
[2752]	214	phys_out_filenames(9) = 'histhfNMC'
	215	phys_out_filenames(10)= 'histstrataer'
[907]	216
[1562]	217	type_ecri(1) = 'ave(X)'
	218	type_ecri(2) = 'ave(X)'
[2167]	219	type_ecri(3) = 'inst(X)'
[1562]	220	type_ecri(4) = 'inst(X)'
	221	type_ecri(5) = 'ave(X)'
	222	type_ecri(6) = 'inst(X)'
[1828]	223	type_ecri(7) = 'inst(X)'
	224	type_ecri(8) = 'inst(X)'
	225	type_ecri(9) = 'inst(X)'
[2752]	226	type_ecri(10)= 'ave(X)'
[907]	227
[1562]	228	clef_files(1) = ok_mensuel
	229	clef_files(2) = ok_journe
	230	clef_files(3) = ok_hf
	231	clef_files(4) = ok_instan
	232	clef_files(5) = ok_LES
	233	clef_files(6) = ok_instan
[1828]	234	clef_files(7) = ok_histNMC(1)
	235	clef_files(8) = ok_histNMC(2)
	236	clef_files(9) = ok_histNMC(3)
[5252]	237	IF (CPPKEY_STRATAER) THEN
	238	clef_files(10)= .TRUE.
	239	ELSE
	240	clef_files(10)= .FALSE.
	241	END IF
[907]	242
[1562]	243	!sortir des fichiers "stations" si clef_stations(:)=.TRUE.
	244	clef_stations(1) = .FALSE.
	245	clef_stations(2) = .FALSE.
	246	clef_stations(3) = .FALSE.
	247	clef_stations(4) = .FALSE.
	248	clef_stations(5) = .FALSE.
	249	clef_stations(6) = .FALSE.
[1828]	250	clef_stations(7) = .FALSE.
	251	clef_stations(8) = .FALSE.
	252	clef_stations(9) = .FALSE.
[2752]	253	clef_stations(10)= .FALSE.
[1539]	254
[1562]	255	lev_files(1) = lev_histmth
	256	lev_files(2) = lev_histday
	257	lev_files(3) = lev_histhf
	258	lev_files(4) = lev_histins
	259	lev_files(5) = lev_histLES
	260	lev_files(6) = lev_histins
[1828]	261	lev_files(7) = levout_histNMC(1)
	262	lev_files(8) = levout_histNMC(2)
	263	lev_files(9) = levout_histNMC(3)
[2752]	264	lev_files(10)= 5
[907]	265
[1562]	266	ecrit_files(1) = ecrit_mth
	267	ecrit_files(2) = ecrit_day
	268	ecrit_files(3) = ecrit_hf
	269	ecrit_files(4) = ecrit_ins
	270	ecrit_files(5) = ecrit_LES
	271	ecrit_files(6) = ecrit_ins
[1912]	272	ecrit_files(7) = freq_outNMC(1)
	273	ecrit_files(8) = freq_outNMC(2)
	274	ecrit_files(9) = freq_outNMC(3)
[2752]	275	ecrit_files(10)= ecrit_mth
[1279]	276
[1562]	277	!! Lectures des parametres de sorties dans physiq.def
[1279]	278
[1791]	279	CALL getin('phys_out_regfkey',phys_out_regfkey)
	280	CALL getin('phys_out_lonmin',phys_out_lonmin)
	281	CALL getin('phys_out_lonmax',phys_out_lonmax)
	282	CALL getin('phys_out_latmin',phys_out_latmin)
	283	CALL getin('phys_out_latmax',phys_out_latmax)
[1562]	284	phys_out_levmin(:)=levmin(:)
[1791]	285	CALL getin('phys_out_levmin',levmin)
[1562]	286	phys_out_levmax(:)=levmax(:)
[1791]	287	CALL getin('phys_out_levmax',levmax)
	288	CALL getin('phys_out_filenames',phys_out_filenames)
[1562]	289	phys_out_filekeys(:)=clef_files(:)
[1791]	290	CALL getin('phys_out_filekeys',clef_files)
[1562]	291	phys_out_filestations(:)=clef_stations(:)
[1791]	292	CALL getin('phys_out_filestations',clef_stations)
[1562]	293	phys_out_filelevels(:)=lev_files(:)
[1791]	294	CALL getin('phys_out_filelevels',lev_files)
	295	CALL getin('phys_out_filetimesteps',chtimestep)
[1562]	296	phys_out_filetypes(:)=type_ecri(:)
[1791]	297	CALL getin('phys_out_filetypes',type_ecri)
[1279]	298
[1562]	299	type_ecri_files(:)=type_ecri(:)
[1279]	300
[2114]	301	! if (ok_all_xml) phys_out_filelevels = 999
	302
[1791]	303	WRITE(lunout,*)'phys_out_lonmin=',phys_out_lonmin
	304	WRITE(lunout,*)'phys_out_lonmax=',phys_out_lonmax
	305	WRITE(lunout,*)'phys_out_latmin=',phys_out_latmin
	306	WRITE(lunout,*)'phys_out_latmax=',phys_out_latmax
	307	WRITE(lunout,*)'phys_out_filenames=',phys_out_filenames
	308	WRITE(lunout,*)'phys_out_filetypes=',type_ecri
	309	WRITE(lunout,*)'phys_out_filekeys=',clef_files
	310	WRITE(lunout,*)'phys_out_filestations=',clef_stations
	311	WRITE(lunout,*)'phys_out_filelevels=',lev_files
[2752]	312	WRITE(lunout,*)'phys_out_regfkey=',phys_out_regfkey
[1562]	313
[1971]	314	! A noter pour
	315	! l heure initiale - dans les fichiers histoire hist* - on met comme
	316	! heure de debut soit la vraie heure (pour le 1D) soit 0h (pour le 3D)
	317	! afin d avoir une seule sortie mensuelle par mois lorsque l on tourne
	318	! par annee (IM).
	319	!
	320	idayref = day_ref
	321	IF (klon_glo==1) THEN
[2491]	322	! current_time (used to compute hour) is updated at the begining of
	323	! the physics; to set the correct outputs "initial time" we thus
	324	! have to use (hour-dtphys).
	325	CALL ymds2ju(annee_ref, 1, idayref, hour-pdtphys, zjulian)
	326	print *,'phys_output_mod: annee,iday,hour,zjulian=',annee_ref,idayref, hour, zjulian
[1971]	327	ELSE
	328	CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
[2551]	329	CALL ymds2ju(annee_ref, 1, day_ini, start_time*rday, zjulian_start)
[2704]	330	ENDIF
[1971]	331
[4619]	332	IF (using_xios) THEN
	333	! ug R\'eglage du calendrier xios
	334	!Temps julian => an, mois, jour, heure
	335	CALL ju2ymds(zjulian, x_an, x_mois, x_jour, x_heure)
	336	CALL ju2ymds(zjulian_start, ini_an, ini_mois, ini_jour, ini_heure)
	337	CALL wxios_set_cal(dtime, calend, x_an, x_mois, x_jour, x_heure, ini_an, &
	338	ini_mois, ini_jour, ini_heure )
	339	ENDIF
[1825]	340
[907]	341	!!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	342	! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
	343	! Appel des histbeg pour definir les variables (nom, moy ou inst, freq de sortie ..
[907]	344	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	345
[1807]	346	zdtime_moy = dtime ! Frequence ou l on moyenne
[1279]	347
[2380]	348
	349	ecrit_files(7) = ecrit_files(1)
	350	ecrit_files(8) = ecrit_files(2)
	351	ecrit_files(9) = ecrit_files(3)
	352
[1971]	353	DO iff=1,nfiles
[929]	354
[1641]	355	! Calculate ecrit_files for all files
[2704]	356	IF ( chtimestep(iff).eq.'Default' ) THEN
[1884]	357	! Par defaut ecrit_files = (ecrit_mensuel ecrit_jour ecrit_hf
	358	! ...)*86400.
[1641]	359	ecrit_files(iff)=ecrit_files(iff)*86400.
[2704]	360	ELSE IF (chtimestep(iff).eq.'-1') THEN
[1912]	361	PRINT*,'ecrit_files(',iff,') < 0 so IOIPSL work on different'
	362	PRINT*,'months length'
	363	ecrit_files(iff)=-1.
[2704]	364	ELSE
[1884]	365	CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff))
	366	ENDIF
	367
[1791]	368	WRITE(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
[1641]	369	zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
[1279]	370
[1971]	371
[4619]	372	IF (using_xios) THEN
	373	!!! Ouverture de chaque fichier XIOS !!!!!!!!!!!
	374	IF (.not. ok_all_xml) THEN
	375	IF (prt_level >= 10) THEN
	376	print*,'phys_output_open: call wxios_add_file with phys_out_filenames(iff)=',trim(phys_out_filenames(iff))
	377	ENDIF
	378	CALL wxios_add_file(phys_out_filenames(iff),chtimestep(iff),lev_files(iff))
	379	ENDIF
	380
	381	!!! Declaration des axes verticaux de chaque fichier:
[2704]	382	IF (prt_level >= 10) THEN
[4619]	383	print*,'phys_output_open: Declare vertical axes for each file'
[2704]	384	ENDIF
[1971]	385
[4619]	386	IF (iff.LE.6.OR.iff.EQ.10) THEN
	387	CALL wxios_add_vaxis("presnivs", &
	388	levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
	389	CALL wxios_add_vaxis("presinter", &
	390	klev + 1, presinter(1:klev+1))
	391	CALL wxios_add_vaxis("Ahyb", &
	392	levmax(iff) - levmin(iff) + 1, aps(levmin(iff):levmax(iff)), positif='down', &
	393	bnds=Ahyb_bounds(levmin(iff):levmax(iff),:))
	394	CALL wxios_add_vaxis("Bhyb", &
	395	levmax(iff) - levmin(iff) + 1, bps(levmin(iff):levmax(iff)), positif='down', &
	396	bnds=Bhyb_bounds(levmin(iff):levmax(iff),:))
	397	CALL wxios_add_vaxis("klev", levmax(iff) - levmin(iff) + 1, &
	398	lev_index(levmin(iff):levmax(iff)))
	399	CALL wxios_add_vaxis("klevp1", klev+1, &
	400	lev_index(1:klev+1))
	401	CALL wxios_add_vaxis("bnds", 2, (/1.,2./))
	402
	403	CALL wxios_add_vaxis("Alt", &
	404	levmax(iff) - levmin(iff) + 1, pseudoalt)
[3003]	405
[4619]	406	! wl1_sun/wl2_sun: minimum/maximum bound of wavelength (in um)
	407	SELECT CASE(NSW)
	408	CASE(6)
	409	wl1_sun(1:6) = [0.180, 0.250, 0.440, 0.690, 1.190, 2.380]
	410	wl2_sun(1:6) = [0.250, 0.440, 0.690, 1.190, 2.380, 4.000]
	411	CASE(2)
	412	wl1_sun(1:2) = [0.250, 0.690]
	413	wl2_sun(1:2) = [0.690, 4.000]
	414	END SELECT
[3003]	415
[4619]	416	DO ISW=1, NSW
	417	wn1_sun(ISW)=1.e+6/wl1_sun(ISW)
	418	wn2_sun(ISW)=1.e+6/wl2_sun(ISW)
	419	spbnds_sun(ISW,1)=wn2_sun(ISW)
	420	spbnds_sun(ISW,2)=wn1_sun(ISW)
	421	spectband(ISW)=(wn1_sun(ISW)+wn2_sun(ISW))/2
	422	ENDDO
[3048]	423	!
	424	!!! ajout axe vertical spectband : solar band number
[4619]	425	CALL wxios_add_vaxis("spectband", NSW, spectband, positif='down')
	426	ELSE
	427	! NMC files
	428	CALL wxios_add_vaxis("plev", &
	429	levmax(iff) - levmin(iff) + 1, rlevSTD(levmin(iff):levmax(iff)))
	430	ENDIF
	431	ENDIF !using_xios
[1852]	432
[2114]	433	IF (clef_files(iff)) THEN
[907]	434	!!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
	435	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2704]	436	IF (phys_out_regfkey(iff)) THEN
[1562]	437	imin_ins=1
[2344]	438	imax_ins=nbp_lon
[1562]	439	jmin_ins=1
	440	jmax_ins=jjmp1
[907]	441
[1562]	442	! correction abderr
[2704]	443	DO i=1,nbp_lon
[1791]	444	WRITE(lunout,*)'io_lon(i)=',io_lon(i)
	445	IF (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
	446	IF (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
[2704]	447	ENDDO
[907]	448
[2704]	449	DO j=1,jjmp1
[1791]	450	WRITE(lunout,*)'io_lat(j)=',io_lat(j)
	451	IF (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
	452	IF (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
[2704]	453	ENDDO
[907]	454
[1791]	455	WRITE(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
[1562]	456	imin_ins,imax_ins,jmin_ins,jmax_ins
[1791]	457	WRITE(lunout,*)'longitudes : ', &
[1562]	458	io_lon(imin_ins),io_lon(imax_ins), &
	459	'latitudes : ', &
	460	io_lat(jmax_ins),io_lat(jmin_ins)
[907]	461
[2344]	462	CALL histbeg(phys_out_filenames(iff),nbp_lon,io_lon,jjmp1,io_lat, &
[1562]	463	imin_ins,imax_ins-imin_ins+1, &
	464	jmin_ins,jmax_ins-jmin_ins+1, &
	465	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
[907]	466	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1562]	467	!IM fichiers stations
[2704]	468	ELSE IF (clef_stations(iff)) THEN
[1539]	469
[2704]	470	IF (prt_level >= 10) THEN
[2001]	471	WRITE(lunout,*)'phys_output_open: iff=',iff,' phys_out_filenames(iff)=',phys_out_filenames(iff)
[2704]	472	ENDIF
[2001]	473
[1791]	474	CALL histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
[1562]	475	phys_out_filenames(iff), &
	476	itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
[2704]	477	ELSE
[2752]	478
	479	IF (prt_level >= 10) THEN
	480	WRITE(lunout,*)'phys_output_open: iff=',iff,' phys_out_filenames(iff)=',phys_out_filenames(iff)
	481	ENDIF
	482
[1825]	483	CALL histbeg_phy_all(phys_out_filenames(iff),itau_phy,zjulian,&
[1852]	484	dtime,nhorim(iff),nid_files(iff))
[2704]	485	ENDIF
[907]	486
[2097]	487	#ifndef CPP_IOIPSL_NO_OUTPUT
[2752]	488	IF (iff.LE.6.OR.iff.EQ.10) THEN
[1828]	489	CALL histvert(nid_files(iff), "presnivs", "Vertical levels", "Pa", &
[1562]	490	levmax(iff) - levmin(iff) + 1, &
	491	presnivs(levmin(iff):levmax(iff)), nvertm(iff),"down")
[1807]	492	!!!! Composantes de la coordonnee sigma-hybride
[1562]	493	CALL histvert(nid_files(iff), "Ahyb","Ahyb comp of Hyb Cord ", "Pa", &
[2832]	494	levmax(iff) - levmin(iff) + 1,aps,nvertap(iff))
[907]	495
[1562]	496	CALL histvert(nid_files(iff), "Bhyb","Bhyb comp of Hyb Cord", " ", &
[2832]	497	levmax(iff) - levmin(iff) + 1,bps,nvertbp(iff))
[1279]	498
[1828]	499	CALL histvert(nid_files(iff), "Alt","Height approx for scale heigh of 8km at levels", "Km", &
[2833]	500	levmax(iff) - levmin(iff) + 1,pseudoalt,nvertAlt(iff))
[1279]	501
[2704]	502	ELSE
[2752]	503	! NMC files
[1828]	504	CALL histvert(nid_files(iff), "plev", "pressure", "Pa", &
	505	levmax(iff) - levmin(iff) + 1, &
	506	rlevSTD(levmin(iff):levmax(iff)), nvertm(iff), "down")
[2704]	507	ENDIF
[1828]	508	#endif
[907]	509
[2114]	510	ENDIF ! clef_files
	511
[4998]	512	itr = 0; itrb = 0
[4056]	513	DO iq = 1, nqtot
[4071]	514	IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE
[4056]	515	itr = itr + 1
	516	dn = 'd'//TRIM(tracers(iq)%name)//'_'
[1813]	517
[4046]	518	flag = [1, 5, 5, 5, 10, 10, 11, 11, 11, 11]
[4056]	519	lnam = 'Tracer '//TRIM(tracers(iq)%longName)
	520	tnam = TRIM(tracers(iq)%name); o_trac (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	521
[4046]	522	flag = [4, 7, 7, 7, 10, 10, 11, 11, 11, 11]
[4056]	523	lnam = 'Tendance tracer '//TRIM(tracers(iq)%longName)
	524	tnam = TRIM(dn)//'vdf'; o_dtr_vdf (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
[1813]	525
[4046]	526	flag = [5, 7, 7, 7, 10, 10, 11, 11, 11, 11]
[4056]	527	tnam = TRIM(dn)//'the'; o_dtr_the (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	528	tnam = TRIM(dn)//'con'; o_dtr_con (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
[1813]	529
[4046]	530	flag = [7, 7, 7, 7, 10, 10, 11, 11, 11, 11]
[4056]	531	tnam = TRIM(dn)//'lessi_impa'; o_dtr_lessi_impa(itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	532	tnam = TRIM(dn)//'lessi_nucl'; o_dtr_lessi_nucl(itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	533	tnam = TRIM(dn)//'insc'; o_dtr_insc (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	534	tnam = TRIM(dn)//'bcscav'; o_dtr_bcscav (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	535	tnam = TRIM(dn)//'evapls'; o_dtr_evapls (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	536	tnam = TRIM(dn)//'ls'; o_dtr_ls (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	537	tnam = TRIM(dn)//'trsp'; o_dtr_trsp (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	538	tnam = TRIM(dn)//'sscav'; o_dtr_sscav (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	539	tnam = TRIM(dn)//'sat'; o_dtr_sat (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
	540	tnam = TRIM(dn)//'uscav'; o_dtr_uscav (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
[1813]	541
[4056]	542	lnam = 'tracer tendency dry deposition'//TRIM(tracers(iq)%longName)
	543	tnam = 'cum'//TRIM(dn)//'dry'; o_dtr_dry (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
[1813]	544
[4046]	545	flag = [1, 4, 10, 10, 10, 10, 11, 11, 11, 11]
[4056]	546	lnam = 'Cumulated tracer '//TRIM(tracers(iq)%longName)
	547	tnam = 'cum'//TRIM(tracers(iq)%name); o_trac_cum(itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
[4998]	548
[5252]	549	IF (CPPKEY_STRATAER) THEN
[4998]	550	if(tracers(iq)%name(1:3)=='BIN') then
	551	itrb = itrb + 1
	552	flag = [11, 11, 11, 11, 11, 11, 11, 11, 11, 1]
	553	lnam = 'Dry particle concentration in '//TRIM(tracers(iq)%longName)
	554	tnam = TRIM(tracers(iq)%name)//'_nd_mode'; o_nd_mode (itrb) = ctrl_out(flag, tnam, lnam, "part/m3", [('',i=1,nfiles)])
	555	lnam = 'Sulfate MMR in '//TRIM(tracers(iq)%longName)
	556	tnam = TRIM(tracers(iq)%name)//'_sulfmmr_mode';o_sulfmmr_mode (itrb) = ctrl_out(flag, tnam, lnam, "kg(H2SO4)/kg(air)", [('',i=1,nfiles)])
	557	endif
[5252]	558	END IF
[4998]	559	ENDDO
[1813]	560
[1971]	561	ENDDO ! iff
[1539]	562
[5056]	563	#ifdef ISO
	564	write(,) 'phys_output_mid 589'
	565	do ixt=1,ntraciso
	566	outiso = TRIM(isoName(ixt))
	567	i = INDEX(outiso, '_', .TRUE.)
	568	outiso = outiso(1:i-1)//outiso(i+1:LEN_TRIM(outiso))
	569
	570	flag = [1, 1, 1, 10, 5, 10, 11, 11, 11, 11]; unit = 'kg/(s*m2)'
	571	o_xtprecip(ixt)=ctrl_out(flag, 'precip'//TRIM(outiso), 'Precip Totale liq+sol', unit, [('',i=1,nfiles)])
	572	o_xtpluc (ixt)=ctrl_out(flag, 'pluc'//TRIM(outiso), 'Convective Precip.', unit, [('',i=1,nfiles)])
	573
	574	flag = [1, 1, 1, 10, 10, 10, 11, 11, 11, 11]
	575	o_xtplul (ixt)=ctrl_out(flag, 'plul'//TRIM(outiso), 'Large-scale Precip.', unit, [('',i=1,nfiles)])
	576	o_xtevap (ixt)=ctrl_out(flag, 'evap'//TRIM(outiso), 'Evaporat.', unit, [('',i=1,nfiles)])
	577
	578	! ajout Camille 8 mai 2023
	579	flag = [1, 6, 10, 10, 10, 10, 11, 11, 11, 11]
	580	o_xtevap_srf (ixt,1)=ctrl_out(flag, 'evap_ter'//TRIM(outiso), 'Evap sfc'//clnsurf(1), unit, [('',i=1,nfiles)])
	581	o_xtevap_srf (ixt,2)=ctrl_out(flag, 'evap_lic'//TRIM(outiso), 'Evap sfc'//clnsurf(2), unit, [('',i=1,nfiles)])
	582	o_xtevap_srf (ixt,3)=ctrl_out(flag, 'evap_oce'//TRIM(outiso), 'Evap sfc'//clnsurf(3), unit, [('',i=1,nfiles)])
	583	o_xtevap_srf (ixt,4)=ctrl_out(flag, 'evap_sic'//TRIM(outiso), 'Evap sfc'//clnsurf(4), unit, [('',i=1,nfiles)])
	584
	585	flag = [2, 3, 4, 10, 10, 10, 11, 11, 11, 11]; unit = 'kg/kg'
	586	o_xtovap (ixt)=ctrl_out(flag, 'ovap'//TRIM(outiso), 'Specific humidity', unit, [('',i=1,nfiles)])
	587	o_xtoliq (ixt)=ctrl_out(flag, 'oliq'//TRIM(outiso), 'Liquid water', unit, [('',i=1,nfiles)])
	588	o_xtcond (ixt)=ctrl_out(flag, 'ocond'//TRIM(outiso), 'Condensed water', unit, [('',i=1,nfiles)])
	589
	590	flag = [1, 1, 1, 10, 5, 10, 11, 11, 11, 11]; unit = 'kg/m2/s'
	591	o_xtrunoff_diag (ixt)=ctrl_out(flag, 'runoffland'//TRIM(outiso), 'Run-off rate land for bucket', unit, [('',i=1,nfiles)])
	592
	593	flag = [4, 10, 10, 10, 10, 10, 11, 11, 11, 11]; unit = '(kg/kg)/s'
	594	o_dxtdyn (ixt)=ctrl_out(flag, 'dqdyn'//TRIM(outiso), 'Dynamics dQ', unit, [('',i=1,nfiles)])
	595	o_dxtldyn (ixt)=ctrl_out(flag, 'dqldyn'//TRIM(outiso), 'Dynamics dQL', unit, [('',i=1,nfiles)])
	596	o_dxtcon (ixt)=ctrl_out(flag, 'dqcon'//TRIM(outiso), 'Convection dQ', unit, [('',i=1,nfiles)])
	597	o_dxteva (ixt)=ctrl_out(flag, 'dqeva'//TRIM(outiso), 'Reevaporation dQ', unit, [('',i=1,nfiles)])
	598	o_dxtlsc (ixt)=ctrl_out(flag, 'dqlsc'//TRIM(outiso), 'Condensation dQ', unit, [('',i=1,nfiles)])
	599	o_dxtajs (ixt)=ctrl_out(flag, 'dqajs'//TRIM(outiso), 'Dry adjust. dQ', unit, [('',i=1,nfiles)])
	600	o_dxtvdf (ixt)=ctrl_out(flag, 'dqvdf'//TRIM(outiso), 'Boundary-layer dQ', unit, [('',i=1,nfiles)])
	601	o_dxtthe (ixt)=ctrl_out(flag, 'dqthe'//TRIM(outiso), 'Thermal dQ', unit, [('',i=1,nfiles)])
	602
	603	IF(ok_qch4) o_dxtch4(ixt)=ctrl_out(flag, 'dqch4'//TRIM(outiso), 'H2O due to CH4 oxidation & photolysis', &
	604	unit, [('',i=1,nfiles)])
	605	IF(ixt == iso_HTO) THEN
	606	o_dxtprod_nucl(ixt)=ctrl_out(flag, 'dqprodnucl'//TRIM(outiso), 'dHTO/dt due to nuclear production', &
	607	unit, [('',i=1,nfiles)])
	608	o_dxtcosmo (ixt)=ctrl_out(flag, 'dqcosmo'//TRIM(outiso), 'dHTO/dt due to cosmogenic production', &
	609	unit, [('',i=1,nfiles)])
	610	o_dxtdecroiss (ixt)=ctrl_out(flag, 'dqdecroiss'//TRIM(outiso), 'dHTO/dt due to radiative destruction', &
	611	unit, [('',i=1,nfiles)])
	612	END IF
	613	enddo !do ixt=1,niso
	614	write(,) 'phys_output_mid 596'
	615	#endif
	616
	617	! Updated write frequencies due to phys_out_filetimesteps.
[1641]	618	! Write frequencies are now in seconds.
	619	ecrit_mth = ecrit_files(1)
	620	ecrit_day = ecrit_files(2)
	621	ecrit_hf = ecrit_files(3)
	622	ecrit_ins = ecrit_files(4)
	623	ecrit_LES = ecrit_files(5)
	624	ecrit_ins = ecrit_files(6)
	625
[2704]	626	IF (prt_level >= 10) THEN
[3082]	627	WRITE(lunout,*)'swaerofree_diag=',swaerofree_diag
[2001]	628	WRITE(lunout,*)'swaero_diag=',swaero_diag
[2854]	629	WRITE(lunout,*)'dryaod_diag=',dryaod_diag
[2989]	630	WRITE(lunout,*)'ok_4xCO2atm=',ok_4xCO2atm
[2001]	631	WRITE(lunout,*)'phys_output_open: ends here'
[2704]	632	ENDIF
[1920]	633
[4056]	634	! DO iq=1,nqtot
[4071]	635	! IF(.NOT.(tracers(iq)%isAdvected .AND. tracers(iq)%isInPhysics)) CYCLE
[4056]	636	! WRITE(*,'(a,i1,a,10i3)')'trac(',iq,')%flag = ',o_trac(iq)%flag
	637	! WRITE(*,'(a,i1,a)')'trac(',iq,')%name = '//TRIM(o_trac(iq)%name)
	638	! WRITE(*,'(a,i1,a)')'trac(',iq,')%description = '//TRIM(o_trac(iq)%description)
[4046]	639	! END DO
	640
[2704]	641	END SUBROUTINE phys_output_open
[907]	642
[1562]	643	SUBROUTINE convers_timesteps(str,dtime,timestep)
[1279]	644
[1562]	645	use ioipsl
	646	USE phys_cal_mod
[2344]	647	USE time_phylmdz_mod, ONLY: day_ref, annee_ref
[2311]	648	USE print_control_mod, ONLY: lunout
[1279]	649
[1562]	650	IMPLICIT NONE
[1279]	651
[1791]	652	CHARACTER(LEN=20) :: str
	653	CHARACTER(LEN=10) :: type
	654	INTEGER :: ipos,il
[1562]	655	real :: ttt,xxx,timestep,dayseconde,dtime
	656	parameter (dayseconde=86400.)
[1279]	657
[1791]	658	ipos=scan(str,'0123456789.',.TRUE.)
[1562]	659	!
	660	il=len_trim(str)
[1937]	661	WRITE(lunout,*) "ipos = ", ipos
	662	WRITE(lunout,*) "il = ", il
[2704]	663	IF (ipos == 0) CALL abort_physic("convers_timesteps", "bad str", 1)
[1562]	664	read(str(1:ipos),*) ttt
[1791]	665	WRITE(lunout,*)ttt
[1562]	666	type=str(ipos+1:il)
[1279]	667
[2704]	668	IF ( il == ipos ) THEN
[1562]	669	type='day'
[2704]	670	ENDIF
[1279]	671
[5357]	672	IF ( type == 'day'.or.type == 'days'.or.type == 'd'.or.type == 'jours'.or.type == 'jour'.or.type == 'j' )&
	673	& timestep = ttt * dayseconde
	674	IF ( type == 'months'.or.type == 'month'.or.type == 'mth'.or.type == 'mois'.or.type == 'm' ) THEN
[1791]	675	WRITE(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,mth_len
[1562]	676	timestep = ttt * dayseconde * mth_len
[2704]	677	ENDIF
[5357]	678	IF ( type == 'hours'.or.type == 'hour'.or.type == 'hr'.or.type == 'heures'.or.type == 'heure'.or.type =='h' )&
	679	& timestep = ttt * dayseconde / 24.
	680	IF ( type == 'mn'.or.type == 'minutes'.or.type == 'minute'.or.type == 'm' ) timestep = ttt * 60.
	681	IF ( type == 's'.or.type == 'sec'.or.type == 'secondes'.or.type =='seconde' ) timestep = ttt
[1791]	682	IF ( type == 'TS' ) timestep = ttt * dtime
[1279]	683
[1791]	684	WRITE(lunout,*)'type = ',type
	685	WRITE(lunout,*)'nb j/h/m = ',ttt
	686	WRITE(lunout,*)'timestep(s)=',timestep
[1279]	687
[1562]	688	END SUBROUTINE convers_timesteps
[1279]	689
[907]	690	END MODULE phys_output_mod

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

Download in other formats:

Original Format