Changeset 9 for trunk/src/read_antarc.py


Ignore:
Timestamp:
02/18/14 16:43:11 (10 years ago)
Author:
pinsard
Message:

fix thanks to coding rules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/read_antarc.py

    r4 r9  
    3838 
    3939while (iligne < nbtotal-1) : 
    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 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+1 
    52           
     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 
    5353 
    5454fichier.close 
    5555 
    56   
     56 
    5757 
    5858dx=1.0 
     
    7070zz1= 300 
    7171outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,zz0, zz1) 
    72 monthly_outz=outz     
     72monthly_outz=outz 
    7373monthly_lon=outx 
    7474monthly_lat=outy 
     
    8282zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    8383outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    84 monthly_outz_ssmis_polar[0,:,:]=outz     
     84monthly_outz_ssmis_polar[0,:,:]=outz 
    8585del outz, outx, outy, zz 
    8686 
    8787zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    8888outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    89 monthly_outz_ssmis_polar[1,:,:]=outz     
     89monthly_outz_ssmis_polar[1,:,:]=outz 
    9090del outz, outx, outy, zz 
    9191 
    9292zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    9393outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    94 monthly_outz_ssmis_polar[2,:,:]=outz     
     94monthly_outz_ssmis_polar[2,:,:]=outz 
    9595del outz, outx, outy, zz, xx, yy 
    9696 
     
    100100zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    101101outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    102 monthly_outz_ssmis_diff[0,:,:]=outz     
     102monthly_outz_ssmis_diff[0,:,:]=outz 
    103103del outz, outx, outy, zz 
    104104 
    105105zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    106106outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    107 monthly_outz_ssmis_diff[1,:,:]=outz     
     107monthly_outz_ssmis_diff[1,:,:]=outz 
    108108del outz, outx, outy, zz 
    109109 
    110110zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    111111outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    112 monthly_outz_ssmis_diff[2,:,:]=outz     
     112monthly_outz_ssmis_diff[2,:,:]=outz 
    113113del outz, outx, outy, zz, xx, yy 
    114114 
     
    132132 
    133133bins=arange(0.3,1,0.001) 
    134 bb=(lat_ssmis >= 0)  
    135   
     134bb=(lat_ssmis >= 0) 
     135 
    136136plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    137137plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
     
    222222from netCDF4 import Dataset 
    223223rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    224   
     224 
    225225rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    226226rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 
     
    245245 
    246246rootgrp.close() 
    247  
    248  
Note: See TracChangeset for help on using the changeset viewer.