Context Navigation

iostart_PEM.F90 @ 3086

Last change on this file since 3086 was 3050, checked in by jbclement, 22 months ago
Mars PEM: Minor changes concerning the form of the code in the PEM. JBC
File size: 35.8 KB

Rev	Line
[2794]	1	MODULE iostart_PEM
[2855]	2	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	3	!!!
	4	!!! Purpose: Read and write specific netcdf for the PEM
	5	!!!
	6	!!!
	7	!!! Author: LL, inspired by iostart from the GCM
	8	!!!
	9	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2794]	10
	11	IMPLICIT NONE
	12	PRIVATE
	13	INTEGER,SAVE :: nid_start ! NetCDF file identifier for startfi.nc file
	14	INTEGER,SAVE :: nid_restart ! NetCDF file identifier for restartfi.nc file
	15
	16	! restartfi.nc file dimension identifiers: (see open_restartphy())
	17	INTEGER,SAVE :: idim1 ! "index" dimension
	18	INTEGER,SAVE :: idim2 ! "physical_points" dimension
	19	INTEGER,SAVE :: idim3 ! "subsurface_layers" dimension
	20	INTEGER,SAVE :: idim4 ! "nlayer_plus_1" dimension
	21	INTEGER,SAVE :: idim5 ! "number_of_advected_fields" dimension
	22	INTEGER,SAVE :: idim6 ! "nlayer" dimension
	23	INTEGER,SAVE :: idim7 ! "Time" dimension
	24	INTEGER,SAVE :: idim8 ! "nslope" dimension
	25	INTEGER,SAVE :: idim9 ! "inter slope" dimension
	26	INTEGER,SAVE :: timeindex ! current time index (for time-dependent fields)
	27	INTEGER,PARAMETER :: length=100 ! size of tab_cntrl array
	28
	29	INTERFACE get_field
	30	MODULE PROCEDURE Get_field_r1,Get_field_r2,Get_field_r3
	31	END INTERFACE get_field
	32
	33	INTERFACE get_var
	34	MODULE PROCEDURE get_var_r0,Get_var_r1,Get_var_r2,Get_var_r3
	35	END INTERFACE get_var
	36
	37	INTERFACE put_field
	38	MODULE PROCEDURE put_field_r1,put_field_r2,put_field_r3
	39	END INTERFACE put_field
	40
	41	INTERFACE put_var
	42	MODULE PROCEDURE put_var_r0,put_var_r1,put_var_r2,put_var_r3
	43	END INTERFACE put_var
	44
	45	PUBLIC nid_start, length
	46	PUBLIC get_field,get_var,put_field,put_var
	47	PUBLIC inquire_dimension, inquire_dimension_length
	48	PUBLIC inquire_field, inquire_field_ndims
	49	PUBLIC open_startphy,close_startphy,open_restartphy,close_restartphy
	50
	51	CONTAINS
	52
	53	SUBROUTINE open_startphy(filename)
	54	USE netcdf, only: NF90_OPEN, NF90_NOERR, NF90_NOWRITE, nf90_strerror
	55	USE mod_phys_lmdz_para, only: is_master, bcast
	56	IMPLICIT NONE
	57	CHARACTER(LEN=*) :: filename
	58	INTEGER :: ierr
	59
	60	IF (is_master) THEN
	61	ierr = NF90_OPEN (filename, NF90_NOWRITE, nid_start)
	62	IF (ierr.NE.NF90_NOERR) THEN
	63	write(,)'open_startphy: problem opening file '//trim(filename)
	64	write(,)trim(nf90_strerror(ierr))
	65	CALL abort_physic("open_startphy","Cannot open file",1)
	66	ENDIF
	67	ENDIF
	68
	69	CALL bcast(nid_start) ! tell all procs about nid_start
	70
	71	END SUBROUTINE open_startphy
	72
	73	SUBROUTINE close_startphy
	74	USE netcdf, only: NF90_CLOSE
	75	USE mod_phys_lmdz_para, only: is_master
	76	IMPLICIT NONE
	77	INTEGER :: ierr
	78
	79	IF (is_master) THEN
	80	ierr = NF90_CLOSE (nid_start)
	81	ENDIF
	82
	83	END SUBROUTINE close_startphy
	84
	85
	86	FUNCTION inquire_field(Field_name)
	87	! check if a given field is present in the input file
	88	USE netcdf, only: NF90_INQ_VARID, NF90_NOERR
	89	USE mod_phys_lmdz_para, only: is_master, bcast
	90	IMPLICIT NONE
	91	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	92	LOGICAL :: inquire_field
	93	INTEGER :: varid
	94	INTEGER :: ierr
	95
	96	IF (is_master) THEN
	97	ierr=NF90_INQ_VARID(nid_start,Field_name,varid)
	98	IF (ierr==NF90_NOERR) THEN
	99	Inquire_field=.TRUE.
	100	ELSE
	101	Inquire_field=.FALSE.
	102	ENDIF
	103	ENDIF
	104
	105	CALL bcast(inquire_field)
	106
	107	END FUNCTION inquire_field
	108
	109
	110	FUNCTION inquire_field_ndims(Field_name)
	111	! give the number of dimensions of "Field_name" stored in the input file
	112	USE netcdf, only: nf90_inq_varid, nf90_inquire_variable, &
	113	NF90_NOERR, nf90_strerror
	114	USE mod_phys_lmdz_para, only: is_master, bcast
	115	IMPLICIT NONE
	116	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	117	INTEGER :: inquire_field_ndims
	118	INTEGER :: varid
	119	INTEGER :: ierr
	120
	121	IF (is_master) THEN
	122	ierr=nf90_inq_varid(nid_start,Field_name,varid)
	123	ierr=nf90_inquire_variable(nid_start,varid,&
	124	ndims=inquire_field_ndims)
	125	IF (ierr.NE.NF90_NOERR) THEN
	126	write(,)'inquire_field_ndims: problem obtaining ndims of '&
	127	//trim(field_name)
	128	write(,)trim(nf90_strerror(ierr))
	129	CALL abort_physic("inquire_field_ndims","Failed to get ndims",1)
	130	ENDIF
	131	ENDIF
	132
	133	CALL bcast(inquire_field_ndims)
	134
	135	END FUNCTION inquire_field_ndims
	136
	137
	138	FUNCTION inquire_dimension(Field_name)
	139	! check if a given dimension is present in the input file
	140	USE netcdf, only: nf90_inq_dimid, NF90_NOERR
	141	USE mod_phys_lmdz_para, only: is_master, bcast
	142	IMPLICIT NONE
	143	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	144	LOGICAL :: inquire_dimension
	145	INTEGER :: varid
	146	INTEGER :: ierr
	147
	148	IF (is_master) THEN
	149	ierr=NF90_INQ_DIMID(nid_start,Field_name,varid)
	150	IF (ierr==NF90_NOERR) THEN
	151	Inquire_dimension=.TRUE.
	152	ELSE
	153	Inquire_dimension=.FALSE.
	154	ENDIF
	155	ENDIF
	156
	157	CALL bcast(inquire_dimension)
	158
	159	END FUNCTION inquire_dimension
	160
	161	FUNCTION inquire_dimension_length(Field_name)
	162	! give the length of the "Field_name" dimension stored in the input file
	163	USE netcdf, only: nf90_inquire_dimension, nf90_inq_dimid, &
	164	NF90_NOERR, nf90_strerror
	165	USE mod_phys_lmdz_para, only: is_master, bcast
	166	IMPLICIT NONE
	167	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	168	INTEGER :: inquire_dimension_length
	169	INTEGER :: varid
	170	INTEGER :: ierr
	171
	172	IF (is_master) THEN
	173	ierr=nf90_inq_dimid(nid_start,Field_name,varid)
	174	ierr=nf90_inquire_dimension(nid_start,varid,&
	175	len=inquire_dimension_length)
	176	IF (ierr.NE.NF90_NOERR) THEN
	177	write(,)'inquire_field_length: problem obtaining length of '&
	178	//trim(field_name)
	179	write(,)trim(nf90_strerror(ierr))
	180	CALL abort_physic("inquire_field_ndims","Failed to get length",1)
	181	ENDIF
	182	ENDIF
	183
	184	CALL bcast(inquire_dimension_length)
	185
	186	END FUNCTION inquire_dimension_length
	187
	188
	189
	190	SUBROUTINE Get_Field_r1(field_name,field,found,timeindex)
	191	! For a surface field
	192	use mod_grid_phy_lmdz, only: klon_glo ! number of atmospheric columns (full grid)
	193	IMPLICIT NONE
	194	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	195	REAL,INTENT(INOUT) :: Field(:)
	196	LOGICAL,INTENT(OUT),OPTIONAL :: found
	197	INTEGER,INTENT(IN),OPTIONAL :: timeindex ! time index of sought data
	198
	199	integer :: edges(4), corners(4)
	200
	201	edges(1)=klon_glo
	202	edges(2:4)=1
	203	corners(1:4)=1
	204	if (PRESENT(timeindex)) then
	205	corners(2)=timeindex
	206	endif
	207
	208	IF (PRESENT(found)) THEN
	209	CALL Get_field_rgen(field_name,field,1,corners,edges,found)
	210	ELSE
	211	CALL Get_field_rgen(field_name,field,1,corners,edges)
	212	ENDIF
	213
	214	END SUBROUTINE Get_Field_r1
	215
	216	SUBROUTINE Get_Field_r2(field_name,field,found,timeindex)
	217	! For a "3D" horizontal-vertical field
	218	use mod_grid_phy_lmdz, only: klon_glo ! number of atmospheric columns (full grid)
	219	IMPLICIT NONE
	220	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	221	REAL,INTENT(INOUT) :: Field(:,:)
	222	LOGICAL,INTENT(OUT),OPTIONAL :: found
	223	INTEGER,INTENT(IN),OPTIONAL :: timeindex ! time index of sought data
	224
	225	integer :: edges(4), corners(4)
	226
	227	edges(1)=klon_glo
	228	edges(2)=size(field,2)
	229	edges(3:4)=1
	230	corners(1:4)=1
	231	if (PRESENT(timeindex)) then
	232	corners(3)=timeindex
	233	endif
	234
	235	IF (PRESENT(found)) THEN
	236	CALL Get_field_rgen(field_name,field,size(field,2),&
	237	corners,edges,found)
	238	ELSE
	239	CALL Get_field_rgen(field_name,field,size(field,2),&
	240	corners,edges)
	241	ENDIF
	242
	243
	244	END SUBROUTINE Get_Field_r2
	245
	246	SUBROUTINE Get_Field_r3(field_name,field,found,timeindex)
	247	! for a "4D" field surf/alt/??
	248	use mod_grid_phy_lmdz, only: klon_glo ! number of atmospheric columns (full grid)
	249	IMPLICIT NONE
	250	CHARACTER(LEN=*),INTENT(IN) :: Field_name
	251	REAL,INTENT(INOUT) :: Field(:,:,:)
	252	LOGICAL,INTENT(OUT),OPTIONAL :: found
	253	INTEGER,INTENT(IN),OPTIONAL :: timeindex ! time index of sought data
	254
	255	integer :: edges(4), corners(4)
	256
	257	edges(1)=klon_glo
	258	edges(2)=size(field,2)
	259	edges(3)=size(field,3)
	260	edges(4)=1
	261	corners(1:4)=1
	262	if (PRESENT(timeindex)) then
	263	corners(4)=timeindex
	264	endif
	265
	266	IF (PRESENT(found)) THEN
	267	CALL Get_field_rgen(field_name,field,size(field,2)*size(field,3),&
	268	corners,edges,found)
	269	ELSE
	270	CALL Get_field_rgen(field_name,field,size(field,2)*size(field,3),&
	271	corners,edges)
	272	ENDIF
	273
	274	END SUBROUTINE Get_Field_r3
	275
	276	SUBROUTINE Get_field_rgen(field_name,field,field_size, &
	277	corners,edges,found)
	278	USE netcdf
	279	USE dimphy
	280	USE mod_grid_phy_lmdz
	281	USE mod_phys_lmdz_para
	282	IMPLICIT NONE
	283	CHARACTER(LEN=*) :: Field_name
	284	INTEGER :: field_size
	285	REAL :: field(klon,field_size)
	286	INTEGER,INTENT(IN) :: corners(4)
	287	INTEGER,INTENT(IN) :: edges(4)
	288	LOGICAL,OPTIONAL :: found
	289
	290	REAL :: field_glo(klon_glo,field_size)
	291	LOGICAL :: tmp_found
	292	INTEGER :: varid
	293	INTEGER :: ierr
	294
	295	IF (is_master) THEN
	296
	297	ierr=NF90_INQ_VARID(nid_start,Field_name,varid)
	298
	299	IF (ierr==NF90_NOERR) THEN
	300	CALL body(field_glo)
	301	tmp_found=.TRUE.
	302	ELSE
	303	tmp_found=.FALSE.
	304	ENDIF
	305
	306	ENDIF
	307
	308	CALL bcast(tmp_found)
	309
	310	IF (tmp_found) THEN
	311	CALL scatter(field_glo,field)
	312	ENDIF
	313
	314	IF (PRESENT(found)) THEN
	315	found=tmp_found
	316	ELSE
	317	IF (.NOT. tmp_found) THEN
[3050]	318	write(,) 'get_field_rgen: Field <'//field_name//'> not found'
[2794]	319	CALL abort_physic("get_field_rgen","Failed to get field",1)
	320	ENDIF
	321	ENDIF
	322
	323
	324	CONTAINS
	325
	326	SUBROUTINE body(field_glo)
	327	REAL :: field_glo(klon_glo*field_size)
	328	ierr=NF90_GET_VAR(nid_start,varid,field_glo,corners,edges)
	329	IF (ierr/=NF90_NOERR) THEN
	330	! La variable exist dans le fichier mais la lecture a echouee.
[3050]	331	write(,) 'get_field_rgen: Failed reading <'//field_name//'>'
[2794]	332
	333	! IF (field_name=='CLWCON' .OR. field_name=='RNEBCON' .OR. field_name=='RATQS') THEN
	334	! ! Essaye de lire le variable sur surface uniqument, comme fait avant
	335	! field_glo(:)=0.
	336	! ierr=NF90_GET_VAR(nid_start,varid,field_glo(1:klon_glo))
	337	! IF (ierr/=NF90_NOERR) THEN
[3050]	338	! write(,) 'phyetat0: Lecture echouee aussi en 2D pour <'//field_name//'>'
[2794]	339	! CALL abort
	340	! ELSE
[3050]	341	! write(,) 'phyetat0: La variable <'//field_name//'> lu sur surface seulement'!, selon ancien format, le reste mis a zero'
[2794]	342	! END IF
	343	! ELSE
	344	CALL abort_physic("get_field_rgen","Failed to read field",1)
	345	! ENDIF
	346	ENDIF
	347
	348	END SUBROUTINE body
	349
	350	END SUBROUTINE Get_field_rgen
	351
	352
	353	SUBROUTINE get_var_r0(var_name,var,found)
	354	! Get a scalar from input file
	355	IMPLICIT NONE
	356	CHARACTER(LEN=*),INTENT(IN) :: var_name
	357	REAL,INTENT(INOUT) :: var
	358	LOGICAL,OPTIONAL,INTENT(OUT) :: found
	359
	360	REAL :: varout(1)
	361
	362	IF (PRESENT(found)) THEN
	363	CALL Get_var_rgen(var_name,varout,size(varout),found)
	364	ELSE
	365	CALL Get_var_rgen(var_name,varout,size(varout))
	366	ENDIF
	367	var=varout(1)
	368
	369	END SUBROUTINE get_var_r0
	370
	371	SUBROUTINE get_var_r1(var_name,var,found)
	372	! Get a vector from input file
	373	IMPLICIT NONE
	374	CHARACTER(LEN=*),INTENT(IN) :: var_name
	375	REAL,INTENT(INOUT) :: var(:)
	376	LOGICAL,OPTIONAL,INTENT(OUT) :: found
	377
	378	IF (PRESENT(found)) THEN
	379	CALL Get_var_rgen(var_name,var,size(var),found)
	380	ELSE
	381	CALL Get_var_rgen(var_name,var,size(var))
	382	ENDIF
	383
	384	END SUBROUTINE get_var_r1
	385
	386	SUBROUTINE get_var_r2(var_name,var,found)
	387	! Get a 2D field from input file
	388	IMPLICIT NONE
	389	CHARACTER(LEN=*),INTENT(IN) :: var_name
	390	REAL,INTENT(OUT) :: var(:,:)
	391	LOGICAL,OPTIONAL,INTENT(OUT) :: found
	392
	393	IF (PRESENT(found)) THEN
	394	CALL Get_var_rgen(var_name,var,size(var),found)
	395	ELSE
	396	CALL Get_var_rgen(var_name,var,size(var))
	397	ENDIF
	398
	399	END SUBROUTINE get_var_r2
	400
	401	SUBROUTINE get_var_r3(var_name,var,found)
	402	! Get a 3D field frominput file
	403	IMPLICIT NONE
	404	CHARACTER(LEN=*),INTENT(IN) :: var_name
	405	REAL,INTENT(INOUT) :: var(:,:,:)
	406	LOGICAL,OPTIONAL,INTENT(OUT) :: found
	407
	408	IF (PRESENT(found)) THEN
	409	CALL Get_var_rgen(var_name,var,size(var),found)
	410	ELSE
	411	CALL Get_var_rgen(var_name,var,size(var))
	412	ENDIF
	413
	414	END SUBROUTINE get_var_r3
	415
	416	SUBROUTINE Get_var_rgen(var_name,var,var_size,found)
	417	USE netcdf
	418	USE dimphy
	419	USE mod_grid_phy_lmdz
	420	USE mod_phys_lmdz_para
	421	IMPLICIT NONE
	422	CHARACTER(LEN=*) :: var_name
	423	INTEGER :: var_size
	424	REAL :: var(var_size)
	425	LOGICAL,OPTIONAL :: found
	426
	427	LOGICAL :: tmp_found
	428	INTEGER :: varid
	429	INTEGER :: ierr
	430
	431	IF (is_mpi_root .AND. is_omp_root) THEN
	432
	433	ierr=NF90_INQ_VARID(nid_start,var_name,varid)
	434
	435	IF (ierr==NF90_NOERR) THEN
	436	ierr=NF90_GET_VAR(nid_start,varid,var)
	437	IF (ierr/=NF90_NOERR) THEN
[3050]	438	write(,) 'phyetat0: Failed loading <'//trim(var_name)//'>'
[2794]	439	CALL abort_physic("get_var_rgen","Failed to read variable",1)
	440	ENDIF
	441	tmp_found=.TRUE.
	442	ELSE
	443	tmp_found=.FALSE.
	444	ENDIF
	445
	446	ENDIF
	447
	448	CALL bcast(tmp_found)
	449
	450	IF (tmp_found) THEN
	451	CALL bcast(var)
	452	ENDIF
	453
	454	IF (PRESENT(found)) THEN
	455	found=tmp_found
	456	ELSE
	457	IF (.NOT. tmp_found) THEN
[3050]	458	write(,) 'phyetat0: Variable <'//trim(var_name)//'> not found'
[2794]	459	CALL abort_physic("get_var_rgen","Failed to read variable",1)
	460	ENDIF
	461	ENDIF
	462
	463	END SUBROUTINE Get_var_rgen
	464
	465	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	466
	467	SUBROUTINE open_restartphy(filename)
	468	USE netcdf, only: NF90_CREATE, NF90_CLOBBER, NF90_64BIT_OFFSET, &
	469	NF90_NOERR, nf90_strerror, &
	470	NF90_PUT_ATT, NF90_GLOBAL, NF90_DEF_DIM, &
	471	NF90_UNLIMITED, NF90_ENDDEF, &
	472	NF90_WRITE, NF90_OPEN
	473	USE mod_phys_lmdz_para, only: is_master
	474	USE mod_grid_phy_lmdz, only: klon_glo
	475	USE dimphy, only: klev, klevp1
[2842]	476	#ifndef CPP_STD
[2794]	477	USE tracer_mod, only: nqmx
[2842]	478	#else
	479	use tracer_h, only: nqtot
	480	#endif
[2794]	481	USE comsoil_h_PEM, only: nsoilmx_PEM
	482	USE comslope_mod, only: nslope
	483	IMPLICIT NONE
	484	CHARACTER(LEN=*),INTENT(IN) :: filename
	485	INTEGER :: ierr
	486	LOGICAL,SAVE :: already_created=.false.
[2842]	487
	488	#ifdef CPP_STD
	489	INTEGER :: nqmx
	490	nqmx=nqtot
	491	#endif
[2794]	492
	493	IF (is_master) THEN
	494	if (.not.already_created) then
	495	! At the very first call, create the file
	496	ierr=NF90_CREATE(filename,IOR(NF90_CLOBBER,NF90_64BIT_OFFSET), &
	497	nid_restart)
	498	IF (ierr/=NF90_NOERR) THEN
	499	write(,)'open_restartphy: problem creating file '//trim(filename)
	500	write(,)trim(nf90_strerror(ierr))
	501	CALL abort_physic("open_restartphy","Failed creating file",1)
	502	ENDIF
	503	already_created=.true.
	504	else
	505	! Just open the file
	506	ierr=NF90_OPEN(filename,NF90_WRITE,nid_restart)
	507	IF (ierr/=NF90_NOERR) THEN
	508	write(,)'open_restartphy: problem opening file '//trim(filename)
	509	write(,)trim(nf90_strerror(ierr))
	510	CALL abort_physic("open_restartphy","Failed opening file",1)
	511	ENDIF
	512	return
	513	endif ! of if (.not.already_created)
	514
	515	ierr=NF90_PUT_ATT(nid_restart,NF90_GLOBAL,"title",&
	516	"Physics start file")
	517	IF (ierr/=NF90_NOERR) THEN
	518	write(,)'open_restartphy: problem writing title '
	519	write(,)trim(nf90_strerror(ierr))
	520	ENDIF
	521
	522	ierr=NF90_DEF_DIM(nid_restart,"index",length,idim1)
	523	IF (ierr/=NF90_NOERR) THEN
	524	write(,)'open_restartphy: problem defining index dimension '
	525	write(,)trim(nf90_strerror(ierr))
	526	CALL abort_physic("open_restartphy","Failed defining index",1)
	527	ENDIF
	528
	529	ierr=NF90_DEF_DIM(nid_restart,"physical_points",klon_glo,idim2)
	530	IF (ierr/=NF90_NOERR) THEN
	531	write(,)'open_restartphy: problem defining physical_points dimension '
	532	write(,)trim(nf90_strerror(ierr))
	533	CALL abort_physic("open_restartphy","Failed defining physical_points",1)
	534	ENDIF
	535
	536	ierr=NF90_DEF_DIM(nid_restart,"subsurface_layers",nsoilmx_PEM,idim3)
	537	IF (ierr/=NF90_NOERR) THEN
	538	write(,)'open_restartphy: problem defining subsurface_layers dimension '
	539	write(,)trim(nf90_strerror(ierr))
	540	CALL abort_physic("open_restartphy","Failed defining subsurface_layers",1)
	541	ENDIF
	542
	543	ierr=NF90_DEF_DIM(nid_restart,"nlayer_plus_1",klevp1,idim4)
	544	IF (ierr/=NF90_NOERR) THEN
	545	write(,)'open_restartphy: problem defining nlayer_plus_1 dimension '
	546	write(,)trim(nf90_strerror(ierr))
	547	CALL abort_physic("open_restartphy","Failed defining nlayer_plus_1",1)
	548	ENDIF
	549
	550	if (nqmx.ne.0) then
	551	ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",nqmx,idim5)
	552	else
	553	! pretend nqmx=1 because 0 size implies unlimited dimension for NetCDF
	554	ierr=NF90_DEF_DIM(nid_restart,"number_of_advected_fields",1,idim5)
	555	endif
	556	IF (ierr/=NF90_NOERR) THEN
	557	write(,)'open_restartphy: problem defining number_of_advected_fields dimension '
	558	write(,)trim(nf90_strerror(ierr))
	559	CALL abort_physic("open_restartphy","Failed defining number_of_advected_fields",1)
	560	ENDIF
	561
	562	ierr=NF90_DEF_DIM(nid_restart,"nlayer",klev,idim6)
	563	IF (ierr/=NF90_NOERR) THEN
	564	write(,)'open_restartphy: problem defining nlayer dimension '
	565	write(,)trim(nf90_strerror(ierr))
	566	CALL abort_physic("open_restartphy","Failed defining nlayer",1)
	567	ENDIF
	568
	569	ierr=NF90_DEF_DIM(nid_restart,"Time",NF90_UNLIMITED,idim7)
	570	IF (ierr/=NF90_NOERR) THEN
	571	write(,)'open_restartphy: problem defining Time dimension '
	572	write(,)trim(nf90_strerror(ierr))
	573	CALL abort_physic("open_restartphy","Failed defining Time",1)
	574	ENDIF
	575
	576	ierr=NF90_DEF_DIM(nid_restart,"nslope",nslope,idim8)
	577	IF (ierr/=NF90_NOERR) THEN
	578	write(,)'phyredem: problem defining nslope dimension'
	579	write(,)trim(nf90_strerror(ierr))
	580	CALL ABORT
	581	ENDIF
	582
	583	ierr=NF90_DEF_DIM(nid_restart,"inter slope",nslope+1,idim9)
	584	IF (ierr/=NF90_NOERR) THEN
	585	write(,)'phyredem: problem defining inter slope dimension'
	586	write(,)trim(nf90_strerror(ierr))
	587	CALL ABORT
	588	ENDIF
	589
	590
	591	ierr=NF90_ENDDEF(nid_restart)
	592	IF (ierr/=NF90_NOERR) THEN
	593	write(,)'open_restartphy: problem ending definition mode '
	594	write(,)trim(nf90_strerror(ierr))
	595	CALL abort_physic("open_restartphy","Failed ending definition mode",1)
	596	ENDIF
	597	ENDIF
	598
	599	END SUBROUTINE open_restartphy
	600
	601	SUBROUTINE close_restartphy
	602	USE netcdf, only: NF90_CLOSE
	603	USE mod_phys_lmdz_para, only: is_master
	604	IMPLICIT NONE
	605	INTEGER :: ierr
	606
	607	IF (is_master) THEN
	608	ierr = NF90_CLOSE (nid_restart)
	609	ENDIF
	610
	611	END SUBROUTINE close_restartphy
	612
	613	SUBROUTINE put_field_r1(field_name,title,field,time)
	614	! For a surface field
	615	IMPLICIT NONE
	616	CHARACTER(LEN=*),INTENT(IN) :: field_name
	617	CHARACTER(LEN=*),INTENT(IN) :: title
	618	REAL,INTENT(IN) :: field(:)
	619	REAL,OPTIONAL,INTENT(IN) :: time
	620
	621	IF (present(time)) THEN
	622	! if timeindex is present, it is a time-dependent variable
	623	CALL put_field_rgen(field_name,title,field,1,time)
	624	ELSE
	625	CALL put_field_rgen(field_name,title,field,1)
	626	ENDIF
	627
	628	END SUBROUTINE put_field_r1
	629
	630	SUBROUTINE put_field_r2(field_name,title,field,time)
	631	! For a "3D" horizontal-vertical field
	632	IMPLICIT NONE
	633	CHARACTER(LEN=*),INTENT(IN) :: field_name
	634	CHARACTER(LEN=*),INTENT(IN) :: title
	635	REAL,INTENT(IN) :: field(:,:)
	636	REAL,OPTIONAL,INTENT(IN) :: time
	637
	638	IF (present(time)) THEN
	639	! if timeindex is present, it is a time-dependent variable
	640	CALL put_field_rgen(field_name,title,field,size(field,2),time)
	641	ELSE
	642	CALL put_field_rgen(field_name,title,field,size(field,2))
	643	ENDIF
	644
	645	END SUBROUTINE put_field_r2
	646
	647	SUBROUTINE put_field_r3(field_name,title,field,time)
	648	! For a "4D" field surf/alt/??
	649	IMPLICIT NONE
	650	CHARACTER(LEN=*),INTENT(IN) :: field_name
	651	CHARACTER(LEN=*),INTENT(IN) :: title
	652	REAL,INTENT(IN) :: field(:,:,:)
	653	REAL,OPTIONAL,INTENT(IN) :: time
	654
	655	IF (present(time)) THEN
	656	! if timeindex is present, it is a time-dependent variable
	657	CALL put_field_rgen(field_name,title,field,size(field,2)*size(field,3),&
	658	time)
	659	ELSE
	660	CALL put_field_rgen(field_name,title,field,size(field,2)*size(field,3))
	661	ENDIF
	662
	663	END SUBROUTINE put_field_r3
	664
	665	SUBROUTINE put_field_rgen(field_name,title,field,field_size,time)
	666	USE netcdf
	667	USE dimphy
	668	USE comsoil_h_PEM, only: nsoilmx_PEM
	669	USE comslope_mod, ONLY: nslope
	670	USE mod_grid_phy_lmdz
	671	USE mod_phys_lmdz_para
	672	IMPLICIT NONE
	673	CHARACTER(LEN=*),INTENT(IN) :: field_name
	674	CHARACTER(LEN=*),INTENT(IN) :: title
	675	INTEGER,INTENT(IN) :: field_size
	676	REAL,INTENT(IN) :: field(klon,field_size)
	677	REAL,OPTIONAL,INTENT(IN) :: time
	678
	679	REAL :: field_glo(klon_glo,field_size)
	680	INTEGER :: ierr
	681	INTEGER :: nvarid
	682	INTEGER :: idim
	683
	684	CALL gather(field,field_glo)
	685
	686	IF (is_master) THEN
	687
	688	IF (field_size==1) THEN
	689	! input is a 1D "surface field" array
	690	if (.not.present(time)) then ! for a time-independent field
	691	ierr=NF90_REDEF(nid_restart)
	692	#ifdef NC_DOUBLE
	693	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	694	(/idim2/),nvarid)
	695	#else
	696	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	697	(/idim2/),nvarid)
	698	#endif
	699	if (ierr.ne.NF90_NOERR) then
	700	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	701	write(,)trim(nf90_strerror(ierr))
	702	endif
	703	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	704	ierr=NF90_ENDDEF(nid_restart)
	705	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo)
	706	else
	707	! check if the variable has already been defined:
	708	ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
	709	if (ierr/=NF90_NOERR) then ! variable not found, define it
	710	ierr=NF90_REDEF(nid_restart)
	711	#ifdef NC_DOUBLE
	712	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	713	(/idim2,idim7/),nvarid)
	714	#else
	715	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	716	(/idim2,idim7/),nvarid)
	717	#endif
	718	if (ierr.ne.NF90_NOERR) then
	719	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	720	write(,)trim(nf90_strerror(ierr))
	721	endif
	722	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	723	ierr=NF90_ENDDEF(nid_restart)
	724	endif
	725	! Write the variable
	726	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
	727	start=(/1,timeindex/))
	728	endif ! of if (.not.present(timeindex))
	729
	730	ELSE IF (field_size==klev) THEN
	731	! input is a 2D "atmospheric field" array
	732	if (.not.present(time)) then ! for a time-independent field
	733	ierr=NF90_REDEF(nid_restart)
	734	#ifdef NC_DOUBLE
	735	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	736	(/idim2,idim6/),nvarid)
	737	#else
	738	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	739	(/idim2,idim6/),nvarid)
	740	#endif
	741	if (ierr.ne.NF90_NOERR) then
	742	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	743	write(,)trim(nf90_strerror(ierr))
	744	endif
	745	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	746	ierr=NF90_ENDDEF(nid_restart)
	747	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo)
	748	else
	749	! check if the variable has already been defined:
	750	ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
	751	if (ierr/=NF90_NOERR) then ! variable not found, define it
	752	ierr=NF90_REDEF(nid_restart)
	753	#ifdef NC_DOUBLE
	754	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	755	(/idim2,idim6,idim7/),nvarid)
	756	#else
	757	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	758	(/idim2,idim6,idim7/),nvarid)
	759	#endif
	760	if (ierr.ne.NF90_NOERR) then
	761	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	762	write(,)trim(nf90_strerror(ierr))
	763	endif
	764	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	765	ierr=NF90_ENDDEF(nid_restart)
	766	endif
	767	! Write the variable
	768	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
	769	start=(/1,1,timeindex/))
	770	endif ! of if (.not.present(time))
	771
	772	ELSE IF (field_size==klevp1) THEN
	773	! input is a 2D "interlayer atmospheric field" array
	774	if (.not.present(time)) then ! for a time-independent field
	775	ierr=NF90_REDEF(nid_restart)
	776	#ifdef NC_DOUBLE
	777	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	778	(/idim2,idim4/),nvarid)
	779	#else
	780	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	781	(/idim2,idim4/),nvarid)
	782	#endif
	783	if (ierr.ne.NF90_NOERR) then
	784	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	785	write(,)trim(nf90_strerror(ierr))
	786	endif
	787	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	788	ierr = NF90_ENDDEF(nid_restart)
	789	ierr = NF90_PUT_VAR(nid_restart,nvarid,field_glo)
	790	else
	791	! check if the variable has already been defined:
	792	ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
	793	if (ierr/=NF90_NOERR) then ! variable not found, define it
	794	ierr=NF90_REDEF(nid_restart)
	795	#ifdef NC_DOUBLE
	796	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	797	(/idim2,idim4,idim7/),nvarid)
	798	#else
	799	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	800	(/idim2,idim4,idim7/),nvarid)
	801	#endif
	802	if (ierr.ne.NF90_NOERR) then
	803	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	804	write(,)trim(nf90_strerror(ierr))
	805	endif
	806	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	807	ierr=NF90_ENDDEF(nid_restart)
	808	endif
	809	! Write the variable
	810	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
	811	start=(/1,1,timeindex/))
	812	endif ! of if (.not.present(timeindex))
	813
	814	ELSE IF (field_size==nsoilmx_PEM) THEN
	815	! input is a 2D "subsurface field" array
	816	if (.not.present(time)) then ! for a time-independent field
	817	ierr = NF90_REDEF(nid_restart)
	818	#ifdef NC_DOUBLE
	819	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	820	(/idim2,idim3/),nvarid)
	821	#else
	822	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	823	(/idim2,idim3/),nvarid)
	824	#endif
	825	if (ierr.ne.NF90_NOERR) then
	826	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	827	write(,)trim(nf90_strerror(ierr))
	828	endif
	829	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	830	ierr = NF90_ENDDEF(nid_restart)
	831	ierr = NF90_PUT_VAR(nid_restart,nvarid,field_glo)
	832	else
	833	! check if the variable has already been defined:
	834	ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
	835	if (ierr/=NF90_NOERR) then ! variable not found, define it
	836	ierr=NF90_REDEF(nid_restart)
	837	#ifdef NC_DOUBLE
	838	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	839	(/idim2,idim3,idim7/),nvarid)
	840	#else
	841	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	842	(/idim2,idim3,idim7/),nvarid)
	843	#endif
	844	if (ierr.ne.NF90_NOERR) then
	845	write(,)"put_field_rgen error: failed to define "//trim(field_name)
	846	write(,)trim(nf90_strerror(ierr))
	847	endif
	848	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	849	ierr=NF90_ENDDEF(nid_restart)
	850	endif
	851	! Write the variable
	852	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
	853	start=(/1,1,timeindex/))
	854
	855	endif ! of if (.not.present(time))
	856
	857
	858	ELSE IF (field_size==nslope) THEN
	859	! input is a 2D "subsurface field" array
	860	if (.not.present(time)) then ! for a time-independent field
	861	ierr = NF90_REDEF(nid_restart)
	862	#ifdef NC_DOUBLE
	863	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	864	(/idim2,idim8/),nvarid)
	865	#else
	866	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	867	(/idim2,idim8/),nvarid)
	868	#endif
	869	if (ierr.ne.NF90_NOERR) then
	870	write(,)"put_field_rgen error: failed to define"//trim(field_name)
	871	write(,)trim(nf90_strerror(ierr))
	872	endif
	873	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	874	ierr = NF90_ENDDEF(nid_restart)
	875	ierr = NF90_PUT_VAR(nid_restart,nvarid,field_glo)
	876	else
	877	! check if the variable has already been defined:
	878	ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
	879	if (ierr/=NF90_NOERR) then ! variable not found, define it
	880	ierr=NF90_REDEF(nid_restart)
	881	#ifdef NC_DOUBLE
	882	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
	883	(/idim2,idim8,idim7/),nvarid)
	884	#else
	885	ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
	886	(/idim2,idim8,idim7/),nvarid)
	887	#endif
	888	if (ierr.ne.NF90_NOERR) then
	889	write(,)"put_field_rgen error: failed to define"//trim(field_name)
	890	write(,)trim(nf90_strerror(ierr))
	891	endif
	892	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	893	ierr=NF90_ENDDEF(nid_restart)
	894	endif
	895	! Write the variable
	896	ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
	897	start=(/1,1,timeindex/))
	898
	899	endif ! of if (.not.present(time))
	900
	901	ELSE
[3050]	902	write(,) "Error phyredem(put_field_rgen) : wrong dimension for ",trim(field_name)
[2794]	903	write(,) " field_size =",field_size
	904	CALL abort_physic("put_field_rgen","wrong field dimension",1)
	905	ENDIF
	906
	907	! Check the writting of field to file went OK
	908	if (ierr.ne.NF90_NOERR) then
	909	write(,) " Error phyredem(put_field_rgen) : failed writing ",trim(field_name)
	910	write(,)trim(nf90_strerror(ierr))
	911	CALL abort_physic("put_field_rgen","Failed writing field",1)
	912	endif
	913
	914	ENDIF ! of IF (is_master)
	915
	916	END SUBROUTINE put_field_rgen
	917
	918	SUBROUTINE put_var_r0(var_name,title,var)
	919	! Put a scalar in file
	920	IMPLICIT NONE
	921	CHARACTER(LEN=*),INTENT(IN) :: var_name
	922	CHARACTER(LEN=*),INTENT(IN) :: title
	923	REAL,INTENT(IN) :: var
	924	REAL :: varin(1)
	925
	926	varin(1)=var
	927
	928	CALL put_var_rgen(var_name,title,varin,size(varin))
	929
	930	END SUBROUTINE put_var_r0
	931
	932
	933	SUBROUTINE put_var_r1(var_name,title,var)
	934	! Put a vector in file
	935	IMPLICIT NONE
	936	CHARACTER(LEN=*),INTENT(IN) :: var_name
	937	CHARACTER(LEN=*),INTENT(IN) :: title
	938	REAL,INTENT(IN) :: var(:)
	939
	940	CALL put_var_rgen(var_name,title,var,size(var))
	941
	942	END SUBROUTINE put_var_r1
	943
	944	SUBROUTINE put_var_r2(var_name,title,var)
	945	! Put a 2D field in file
	946	IMPLICIT NONE
	947	CHARACTER(LEN=*),INTENT(IN) :: var_name
	948	CHARACTER(LEN=*),INTENT(IN) :: title
	949	REAL,INTENT(IN) :: var(:,:)
	950
	951	CALL put_var_rgen(var_name,title,var,size(var))
	952
	953	END SUBROUTINE put_var_r2
	954
	955	SUBROUTINE put_var_r3(var_name,title,var)
	956	! Put a 3D field in file
	957	IMPLICIT NONE
	958	CHARACTER(LEN=*),INTENT(IN) :: var_name
	959	CHARACTER(LEN=*),INTENT(IN) :: title
	960	REAL,INTENT(IN) :: var(:,:,:)
	961
	962	CALL put_var_rgen(var_name,title,var,size(var))
	963
	964	END SUBROUTINE put_var_r3
	965
	966	SUBROUTINE put_var_rgen(var_name,title,var,var_size)
	967	USE netcdf, only: NF90_REDEF, NF90_DEF_VAR, NF90_ENDDEF, NF90_PUT_VAR, &
	968	NF90_FLOAT, NF90_DOUBLE, &
	969	NF90_PUT_ATT, NF90_NOERR, nf90_strerror, &
	970	nf90_inq_dimid, nf90_inquire_dimension, NF90_INQ_VARID
	971	USE comsoil_h_PEM, only: nsoilmx_PEM
	972	USE comslope_mod, only: nslope
	973	USE mod_phys_lmdz_para, only: is_master
	974	IMPLICIT NONE
	975	CHARACTER(LEN=*),INTENT(IN) :: var_name
	976	CHARACTER(LEN=*),INTENT(IN) :: title
	977	INTEGER,INTENT(IN) :: var_size
	978	REAL,INTENT(IN) :: var(var_size)
	979
	980	INTEGER :: ierr
	981	INTEGER :: nvarid
	982	INTEGER :: idim1d
	983	logical,save :: firsttime=.true.
	984
	985	IF (is_master) THEN
	986
	987	IF (var_name=="Time") THEN
	988	! Very specific case of "Time" variable
	989	if (firsttime) then
	990	! Create the "Time variable"
	991	ierr=NF90_REDEF(nid_restart)
	992	#ifdef NC_DOUBLE
	993	ierr=NF90_DEF_VAR(nid_restart,var_name,NF90_DOUBLE,&
	994	(/idim7/),nvarid)
	995	#else
	996	ierr=NF90_DEF_VAR(nid_restart,var_name,NF90_FLOAT,&
	997	(/idim7/),nvarid)
	998	#endif
	999	IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	1000	ierr=NF90_ENDDEF(nid_restart)
	1001
	1002	firsttime=.false.
	1003	endif
	1004	! Append "time" value
	1005	! get current length of "Time" dimension
	1006	ierr=nf90_inq_dimid(nid_restart,var_name,idim1d)
	1007	ierr=nf90_inquire_dimension(nid_restart,idim1d,len=timeindex)
	1008	timeindex=timeindex+1
	1009	ierr=NF90_INQ_VARID(nid_restart,var_name,nvarid)
	1010	ierr=NF90_PUT_VAR(nid_restart,nvarid,var,&
	1011	start=(/timeindex/))
	1012	IF (ierr/=NF90_NOERR) THEN
	1013	write(,)'put_var_rgen: problem writing Time'
	1014	write(,)trim(nf90_strerror(ierr))
	1015	CALL abort_physic("get_var_rgen","Failed to write Time",1)
	1016	ENDIF
	1017	return ! nothing left to do
	1018	ELSEIF (var_size==length) THEN
	1019	! We know it is a "controle" kind of 1D array
	1020	idim1d=idim1
	1021	ELSEIF (var_size==nsoilmx_PEM) THEN
	1022	! We know it is an "mlayer" kind of 1D array
	1023	idim1d=idim3
	1024	ELSEIF (var_size==nslope+1) THEN
	1025	! We know it is an "inter slope" kind of 1D array
	1026	idim1d=idim9
	1027	ELSE
[3050]	1028	write(,) "put_var_rgen error : wrong dimension"
[2794]	1029	write(,) " var_size =",var_size
	1030	CALL abort_physic("get_var_rgen","Wrong variable dimension",1)
	1031
	1032	ENDIF ! of IF (var_size==length) THEN
	1033
	1034	! Swich to NetCDF define mode
	1035	ierr=NF90_REDEF (nid_restart)
	1036	! Define the variable
	1037	#ifdef NC_DOUBLE
	1038	ierr=NF90_DEF_VAR(nid_restart,var_name,NF90_DOUBLE,(/idim1d/),nvarid)
	1039	#else
	1040	ierr=NF90_DEF_VAR(nid_restart,var_name,NF90_FLOAT,(/idim1d/),nvarid)
	1041	#endif
	1042	! Add a "title" attribute
	1043	IF (LEN_TRIM(title)>0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
	1044	! Swich out of define mode
	1045	ierr=NF90_ENDDEF(nid_restart)
	1046	! Write variable to file
	1047	ierr=NF90_PUT_VAR(nid_restart,nvarid,var)
	1048	IF (ierr/=NF90_NOERR) THEN
	1049	write(,)'put_var_rgen: problem writing '//trim(var_name)
	1050	write(,)trim(nf90_strerror(ierr))
	1051	CALL abort_physic("get_var_rgen","Failed writing variable",1)
	1052	ENDIF
	1053	ENDIF ! of IF (is_master)
	1054
	1055	END SUBROUTINE put_var_rgen
	1056
	1057	END MODULE iostart_PEM

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

Download in other formats:

Original Format