1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
import netCDF4 as nc
import numpy as np
#reccuperation des données depuis le premier fichier
part_ds = nc.Dataset("ECC-2018-08-01.nc4", "r", format="NETCDF4")
particles = part_ds.variables['part']
Times = part_ds.variables['time']
longit = part_ds.variables['lon']
Latit = part_ds.variables['lat']
#création du nouveau fichier
fn = 'Test2Pc.nc'
ds = nc.Dataset(fn, 'w', format='NETCDF4')
part = ds.createDimension('part', None)
lon = ds.createDimension('lon', 87)
lat = ds.createDimension('lat', 87)
parts = ds.createVariable('part', 'f4', ('part',))
lons = ds.createVariable('lon', 'f4', ('lon',))
lats = ds.createVariable('lat', 'f4', ('lat',))
time = ds.createVariable('time', 'i4', ('part','lon', 'lat',))
#time.units = 'Unknown'
#ecriture des données dans le nouveau fichier
Tabtime = np.arange(0, 87, 1) #dimension de time
for i in Tabtime:
lons[i] = np.append(lons, longit[i])
lats[i] = np.append(lats, Latit[i])
parts[i] = np.append(parts, particles[i])
# time[i, :, :] = np.append( time, ( parts[i], longit[i], Latit[i] ) ) : beaucoup de doutes ici, je ne sais pas comment ranger les valeurs dans time, surtout que parts n'a pas la même dimension que longit et Latit.
ds.close() |
