Changeset 9 for trunk/src/read_antarc.py
- Timestamp:
- 02/18/14 16:43:11 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/read_antarc.py
r4 r9 38 38 39 39 while (iligne < nbtotal-1) : 40 41 42 43 # exemple : listeCoord ['0.22', '2.3', '5.0', '6'] (liste de chaine de caract?es)44 lat[iligne] = float(liste[1])45 lon[iligne] = float(liste[0])46 jjr[iligne] = float(liste[4])47 ts[iligne] = float(liste[10])48 tb[iligne] = float(liste[13])49 emis[iligne] = float(liste[16])50 orog[iligne] = float(liste[13])51 iligne=iligne+152 40 line=fichier.readline() 41 # exemple : line = "0.22 2.3 5.0 6" 42 liste = line.split() 43 # exemple : listeCoord ['0.22', '2.3', '5.0', '6'] (liste de chaine de 44 # caractÚres) 45 lat[iligne] = float(liste[1]) 46 lon[iligne] = float(liste[0]) 47 jjr[iligne] = float(liste[4]) 48 ts[iligne] = float(liste[10]) 49 tb[iligne] = float(liste[13]) 50 emis[iligne] = float(liste[16]) 51 orog[iligne] = float(liste[13]) 52 iligne=iligne+1 53 53 54 54 fichier.close 55 55 56 56 57 57 58 58 dx=1.0 … … 70 70 zz1= 300 71 71 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,zz0, zz1) 72 monthly_outz=outz 72 monthly_outz=outz 73 73 monthly_lon=outx 74 74 monthly_lat=outy … … 82 82 zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 83 83 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 84 monthly_outz_ssmis_polar[0,:,:]=outz 84 monthly_outz_ssmis_polar[0,:,:]=outz 85 85 del outz, outx, outy, zz 86 86 87 87 zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 88 88 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 89 monthly_outz_ssmis_polar[1,:,:]=outz 89 monthly_outz_ssmis_polar[1,:,:]=outz 90 90 del outz, outx, outy, zz 91 91 92 92 zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 93 93 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 94 monthly_outz_ssmis_polar[2,:,:]=outz 94 monthly_outz_ssmis_polar[2,:,:]=outz 95 95 del outz, outx, outy, zz, xx, yy 96 96 … … 100 100 zz = emis_ssmis[1,:]-emis_ssmis[2,:] 101 101 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 102 monthly_outz_ssmis_diff[0,:,:]=outz 102 monthly_outz_ssmis_diff[0,:,:]=outz 103 103 del outz, outx, outy, zz 104 104 105 105 zz = emis_ssmis[4,:]-emis_ssmis[5,:] 106 106 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 107 monthly_outz_ssmis_diff[1,:,:]=outz 107 monthly_outz_ssmis_diff[1,:,:]=outz 108 108 del outz, outx, outy, zz 109 109 110 110 zz = emis_ssmis[6,:]-emis_ssmis[7,:] 111 111 outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 112 monthly_outz_ssmis_diff[2,:,:]=outz 112 monthly_outz_ssmis_diff[2,:,:]=outz 113 113 del outz, outx, outy, zz, xx, yy 114 114 … … 132 132 133 133 bins=arange(0.3,1,0.001) 134 bb=(lat_ssmis >= 0) 135 134 bb=(lat_ssmis >= 0) 135 136 136 plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 137 137 plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') … … 222 222 from netCDF4 import Dataset 223 223 rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 224 224 225 225 rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 226 226 rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) … … 245 245 246 246 rootgrp.close() 247 248
Note: See TracChangeset
for help on using the changeset viewer.