Changeset 7


Ignore:
Timestamp:
02/11/14 15:15:35 (10 years ago)
Author:
pinsard
Message:

no more white space in file names

Location:
trunk/src/python_script
Files:
37 edited
4 moved

Legend:

Unmodified
Added
Removed

  • trunk/src/python_script/bina_gc.py

    r6 r7  
    66    # [i, nbins] = bin(x, dx, x0, x1); 
    77    # 
    8     # Returns the vector of indices, starting from 1,  
    9     # corresponding to the chosen bin size, dx,  
    10     # start x0 and end x1. If x1 is omitted, x1 = max(x) - dx/2.  
    11     # If x0 is omitted, x0 = min(x) + dx/2. If dx is omitted, the data  
     8    # Returns the vector of indices, starting from 1, 
     9    # corresponding to the chosen bin size, dx, 
     10    # start x0 and end x1. If x1 is omitted, x1 = max(x) - dx/2. 
     11    # If x0 is omitted, x0 = min(x) + dx/2. If dx is omitted, the data 
    1212    # are divided into 10 classes. Note that outliers are not removed. 
    1313    # 
    … …  
    2525    #end 
    2626    #if nargin < 3: 
    27     #x0 = min(x)  
     27    #x0 = min(x) 
    2828    #end 
    2929    #if nargin < 4: 
    3030    #x1 = max(x) + dx 
    3131    #end 
    32      
     32 
    3333    nn=size(x) 
    3434    indic = zeros(nn, float) 
    35      
     35 
    3636    nbins= round((x1 - x0)/dx) 
    37      
     37 
    3838    for kk in range(0,nn): indic[kk] = math.ceil((x[kk] - x0)/dx)-1 
    39     #indicp=((x - x0)/dx) + 1  
     39    #indicp=((x - x0)/dx) + 1 
    4040    #indic = indicp.round(); 
    4141 

  • trunk/src/python_script/draw_carte_terre.py

    r6 r7  
    3434llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    3535xii, yii = map(*np.meshgrid(xvec,yvec)) 
    36 #, clevs, cmap=my_cmap)  
     36#, clevs, cmap=my_cmap) 
    3737clevs=arange(0,1,0.001)#star, stop, step 
    3838cs=map.pcolor(xii,yii,zz1) 

  • trunk/src/python_script/draw_map_emis_amsuach2.py

    r6 r7  
    3636llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    3737xii, yii = map(*np.meshgrid(xvec,yvec)) 
    38 #, clevs, cmap=my_cmap)  
     38#, clevs, cmap=my_cmap) 
    3939clevs=arange(0.5,1.1,0.001)#star, stop, step 
    4040cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 

  • trunk/src/python_script/draw_map_mod_gc1.py

    r6 r7  
    1818        m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-90,urcrnrlat=90,projection='cyl',resolution='c',fix_aspect=True) 
    1919        #m = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='c') 
    20         #m = Basemap(resolution='c',projection='sinu',lon_0=60) 
    21         #m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-90,urcrnrlat=90,projection='cyl',resolution='i',fix_aspect=True) 
     20        #m = Basemap(resolution='c',projection='sinu',lon_0=60) 
     21        #m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-90,urcrnrlat=90,projection='cyl',resolution='i',fix_aspect=True) 
    2222        #m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-60,urcrnrlat=60,projection='cyl',resolution='c',fix_aspect=True) 
    23          
     23 
    2424        #ax = plt.gca() 
    25         xii, yii = m(*np.meshgrid(x,y)) 
    26         #ax = plt.gca() 
    27         #clevs=arange(0,0.8,0.005) 
    28         clevs=arange(cax1,cax2,cax3) 
    29         #boxColors= ['darkkhaki','royalblue'] 
     25        xii, yii = m(*np.meshgrid(x,y)) 
     26        #ax = plt.gca() 
     27        #clevs=arange(0,0.8,0.005) 
     28        clevs=arange(cax1,cax2,cax3) 
     29        #boxColors= ['darkkhaki','royalblue'] 
    3030 
    31         m.drawcoastlines(linewidth=1) 
    32         #m.drawlsmask(land_color='white',ocean_color='#F4EAF1', lakes=False) 
     31        m.drawcoastlines(linewidth=1) 
     32        #m.drawlsmask(land_color='white',ocean_color='#F4EAF1', lakes=False) 
    3333        cs=m.contourf(xii,yii,z,20) # emissivity 
    34         cbar =colorbar(cs) 
     34        cbar =colorbar(cs) 
    3535        cbar.set_label(txt) 
    3636        #cbar.set_clim(0.65,1) 
    37         #m.drawmapboundary(fill_color='aqua') 
    38         #m.drawparallels(np.arange(-90.,90.,10.), labels=[1,0,0,0]) 
    39         #m.drawmeridians(np.arange(-180.,180.,30.), labels=[0,0,0,1]) 
    40         m.drawparallels(np.arange(-90.,90.,45.), color='#909090', linewidth=1.0, zorder=None, dashes=[1, 1]) 
    41         #m.drawmeridians(np.arange(-180.,180.,50.)) 
    42         # stations: Thulé, Resolute, Longyearbyen, lqaluit, murmonsk, arkangelek 
     37        #m.drawmapboundary(fill_color='aqua') 
     38        #m.drawparallels(np.arange(-90.,90.,10.), labels=[1,0,0,0]) 
     39        #m.drawmeridians(np.arange(-180.,180.,30.), labels=[0,0,0,1]) 
     40        m.drawparallels(np.arange(-90.,90.,45.), color='#909090', linewidth=1.0, zorder=None, dashes=[1, 1]) 
     41        #m.drawmeridians(np.arange(-180.,180.,50.)) 
     42        # stations: Thulé, Resolute, Longyearbyen, lqaluit, murmonsk, arkangelek 
    4343 
    44         plt.savefig(outfile) 
    45         #close() 
     44        plt.savefig(outfile) 
     45        #close() 

  • trunk/src/python_script/draw_map_ssmisch15.py

    r6 r7  
    3434llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    3535xii, yii = map(*np.meshgrid(xvec,yvec)) 
    36 #, clevs, cmap=my_cmap)  
     36#, clevs, cmap=my_cmap) 
    3737clevs=arange(100,300,1)#star, stop, step 
    3838cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    8080llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    8181xii, yii = map(*np.meshgrid(xvec,yvec)) 
    82 #, clevs, cmap=my_cmap)  
     82#, clevs, cmap=my_cmap) 
    8383clevs=arange(0,40,0.1)#star, stop, step 
    8484cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    125125llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    126126xii, yii = map(*np.meshgrid(xvec,yvec)) 
    127 #, clevs, cmap=my_cmap)  
     127#, clevs, cmap=my_cmap) 
    128128clevs=arange(0,300,1)#star, stop, step 
    129129cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    167167llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    168168xii, yii = map(*np.meshgrid(xvec,yvec)) 
    169 #, clevs, cmap=my_cmap)  
     169#, clevs, cmap=my_cmap) 
    170170clevs=arange(0,300,1)#star, stop, step 
    171171cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 

  • trunk/src/python_script/draw_map_ssmisch16.py

    r6 r7  
    3434llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    3535xii, yii = map(*np.meshgrid(xvec,yvec)) 
    36 #, clevs, cmap=my_cmap)  
     36#, clevs, cmap=my_cmap) 
    3737clevs=arange(100,300,1)#star, stop, step 
    3838cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    8080llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    8181xii, yii = map(*np.meshgrid(xvec,yvec)) 
    82 #, clevs, cmap=my_cmap)  
     82#, clevs, cmap=my_cmap) 
    8383clevs=arange(0,40,0.1)#star, stop, step 
    8484cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    125125llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    126126xii, yii = map(*np.meshgrid(xvec,yvec)) 
    127 #, clevs, cmap=my_cmap)  
     127#, clevs, cmap=my_cmap) 
    128128clevs=arange(100,600,1)#star, stop, step 
    129129cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 
    … …  
    167167llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    168168xii, yii = map(*np.meshgrid(xvec,yvec)) 
    169 #, clevs, cmap=my_cmap)  
     169#, clevs, cmap=my_cmap) 
    170170clevs=arange(0,300,1)#star, stop, step 
    171171cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 

  • trunk/src/python_script/draw_map_tr_amsuach2.py

    r6 r7  
    3434llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    3535xii, yii = map(*np.meshgrid(xvec,yvec)) 
    36 #, clevs, cmap=my_cmap)  
     36#, clevs, cmap=my_cmap) 
    3737clevs=arange(0,1,0.001)#star, stop, step 
    3838cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 

  • trunk/src/python_script/draw_map_tup.py

    r6 r7  
    298298llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) 
    299299xii, yii = map(*np.meshgrid(xvec,yvec)) 
    300 #, clevs, cmap=my_cmap)  
     300#, clevs, cmap=my_cmap) 
    301301clevs=arange(a1,a2,a3)#star, stop, step 
    302302cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) 

  • trunk/src/python_script/ffgrid2_gc.py

    r6 r7  
    1010    z=zz 
    1111    dx=0.1 
    12     dy=1         
     12    dy=1 
    1313    dparx=0.05 
    1414    dpary=0.5 
    … …  
    2121    dx=0.1 
    2222    dy=1.0 
    23     n=size(x)     
     23    n=size(x) 
    2424    xvec = arange(x0,x1,dx) 
    2525    yvec = arange(y0,y11,dy) 
    2626    nx = size(xvec) 
    2727    ny = size(yvec) 
    28     ix = zeros(n, float)     
    29     iy = zeros(n, float)     
     28    ix = zeros(n, float) 
     29    iy = zeros(n, float) 
    3030 
    3131    for kk in range(0,n): 
    … …  
    4141           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4242 
    43     inx = (ix >= 0) & (ix < nx)  
     43    inx = (ix >= 0) & (ix < nx) 
    4444    iny = (iy >= 0) & (iy < ny) 
    45     inz = (z >= z0) & (z <= z1)  
     45    inz = (z >= z0) & (z <= z1) 
    4646    inn = inx & iny & inz 
    4747    iix = ix[inn] 
    … …  
    6262        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    6363        z2grid[iix[i], iiy[i]] = z2grid[iix[i],iiy[i]]+zz[i]*zz[i] 
    64         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    65          
     64        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     65 
    6666    for i in range(0,ny): 
    6767        jj=round(3600/dxx[i]) 
    … …  
    7171        e=d.sum(axis=1) 
    7272        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    73         g=f.sum(axis=1)   
     73        g=f.sum(axis=1) 
    7474        for j in range(0, nx): 
    7575            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7878            zzdgrid[j,i]=g[ii] 
    7979 
    80     zzpgrid = zeros([nx, ny], float)     
     80    zzpgrid = zeros([nx, ny], float) 
    8181    zzpgrid = zzgrid/nngrid 
    8282    sigma_grid=sqrt(zzdgrid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/ffgrid2_gc_nof.py

    r6 r7  
    44 
    55 
    6    
     6 
    77x=xx 
    88y=yy 
    99z=zz 
    1010dx=0.1 
    11 dy=1     
     11dy=1 
    1212dparx=0.05 
    1313dpary=0.5 
    … …  
    2222dy=1.0 
    2323 
    24 n=size(lon)     
     24n=size(lon) 
    2525xvec = arange(x0,x1,dx) 
    2626yvec = arange(y0,y11,dy) 
    … …  
    2828ny = size(yvec) 
    2929 
    30 ix = zeros(n, float)     
    31 iy = zeros(n, float)     
     30ix = zeros(n, float) 
     31iy = zeros(n, float) 
    3232 
    3333for kk in range(0,n): 
    … …  
    4444 
    4545#save txt file 
    46 #np.savetxt('ix_py.txt', ix, fmt='%.18e', delimiter=';')     
     46#np.savetxt('ix_py.txt', ix, fmt='%.18e', delimiter=';') 
    4747 
    4848 
    49 inx = (ix >= 0) & (ix < nx)  
     49inx = (ix >= 0) & (ix < nx) 
    5050iny = (iy >= 0) & (iy < ny) 
    51 inz = (z >= z0) & (z <= z1)  
     51inz = (z >= z0) & (z <= z1) 
    5252inn = inx & iny & inz 
    5353iix = ix[inn] 
    … …  
    6060Rt=6371 
    6161beta=math.pi/180 
    62 R=beta*Rt  
     62R=beta*Rt 
    6363L=zeros(ny,float)#arc vector 
    6464L[0]=2*math.pi*R 
    … …  
    9292zgrid = zeros([nx, ny], float)# z-coordinate 
    9393zzgrid = zeros([nx, ny], float)# z-coordinate grouped 
    94 z2grid = zeros([nx, ny], float)# z-coordinate  
     94z2grid = zeros([nx, ny], float)# z-coordinate 
    9595zz2grid = zeros([nx, ny], float)# z-coordinate  grouped 
    9696 
    … …  
    9898    zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    9999    z2grid[iix[i], iiy[i]] = z2grid[iix[i],iiy[i]]+zz[i]*zz[i] 
    100     ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    101          
     100    ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     101 
    102102for i in range(0,ny): 
    103103    jj=round(3600/dxx[i]) 
    … …  
    107107    e=d.sum(axis=1) 
    108108    f=reshape(z2grid[:,i],(jj,dxx[i])) 
    109     g=f.sum(axis=1)   
     109    g=f.sum(axis=1) 
    110110    for j in range(0, nx): 
    111111        ii=ceil((j+1)/dxx[i])-1 
    … …  
    114114        zz2grid[j,i]=g[ii] 
    115115 
    116 zzpgrid = zeros([nx, ny], float)     
     116zzpgrid = zeros([nx, ny], float) 
    117117zzpgrid = zzgrid/nngrid 
    118118sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 
    119      
     119 
    120120 
    121121zgrid = np.transpose(zgrid) 

  • trunk/src/python_script/ffgrid2_gc_nof_pos.py

    r6 r7  
    44 
    55 
    6    
     6 
    77x=xx 
    88y=yy 
    99z=zz 
    1010dx=0.1 
    11 dy=1     
     11dy=1 
    1212dparx=0.05 
    1313dpary=0.5 
    … …  
    2121dy=1.0 
    2222 
    23 n=size(lon)     
     23n=size(lon) 
    2424xvec = arange(x0,x1,dx) 
    2525yvec = arange(y0,y11,dy) 
    … …  
    2727ny = size(yvec) 
    2828 
    29 ix = zeros(n, float)     
    30 iy = zeros(n, float)     
     29ix = zeros(n, float) 
     30iy = zeros(n, float) 
    3131 
    3232for kk in range(0,n): 
    … …  
    4343 
    4444#save txt file 
    45 #np.savetxt('ix_py.txt', ix, fmt='%.18e', delimiter=';')     
     45#np.savetxt('ix_py.txt', ix, fmt='%.18e', delimiter=';') 
    4646 
    4747 
    48 inx = (ix >= 0) & (ix < nx)  
     48inx = (ix >= 0) & (ix < nx) 
    4949iny = (iy >= 0) & (iy < ny) 
    50 inz = (z >= z0) & (z <= z1)  
    51 ipos = ((pos >=11) & (pos<=20)) #| ((pos >=21) & (pos <=25))  
     50inz = (z >= z0) & (z <= z1) 
     51ipos = ((pos >=11) & (pos<=20)) #| ((pos >=21) & (pos <=25)) 
    5252inn = inx & iny & inz & ipos 
    5353iix = ix[inn] 
    … …  
    6060Rt=6371 
    6161beta=math.pi/180 
    62 R=beta*Rt  
     62R=beta*Rt 
    6363L=zeros(ny,float)#arc vector 
    6464L[0]=2*math.pi*R 
    … …  
    9898    zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    9999    z2grid[iix[i], iiy[i]] = z2grid[iix[i],iiy[i]]+zz[i]*zz[i] 
    100     ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    101          
     100    ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     101 
    102102for i in range(0,ny): 
    103103    jj=round(3600/dxx[i]) 
    … …  
    107107    e=d.sum(axis=1) 
    108108    f=reshape(z2grid[:,i],(jj,dxx[i])) 
    109     g=f.sum(axis=1)   
     109    g=f.sum(axis=1) 
    110110    for j in range(0, nx): 
    111111        ii=ceil((j+1)/dxx[i])-1 
    … …  
    114114        zz2grid[j,i]=g[ii] 
    115115 
    116 zzpgrid = zeros([nx, ny], float)     
     116zzpgrid = zeros([nx, ny], float) 
    117117zzpgrid = zzgrid/nngrid 
    118118sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 
    119      
     119 
    120120 
    121121zgrid = np.transpose(zgrid) 

  • trunk/src/python_script/hist_tb_amsua_ch2.py

    r6 r7  
    119119close() 
    120120 
    121 #ឧist tb cont pos1  
     121#ឧist tb cont pos1 
    122122 
    123123tbp_1nc=np.reshape(tbpgrid_1nc,-1) 
    … …  
    130130close() 
    131131 
    132 #ឧist tb cont pos2  
     132#ឧist tb cont pos2 
    133133 
    134134tbp_2nc=np.reshape(tbpgrid_2nc,-1) 
    … …  
    141141close() 
    142142 
    143 #ឧist tb cont pos3  
     143#ឧist tb cont pos3 
    144144 
    145145tbp_3nc=np.reshape(tbpgrid_3nc,-1) 
    … …  
    157157mertb_1n 
    158158 
    159 #ឧist tb mer pos1  
     159#ឧist tb mer pos1 
    160160 
    161161tbmer_1n=np.reshape(mertb_1n,-1) 
    … …  
    208208 
    209209 
    210 #hist tb mer pas moyenné  
     210#hist tb mer pas moyenné 
    211211imer =  (mm_t50 < 0.5) 
    212212tb_tmer=tb_t50[imer] 

  • trunk/src/python_script/newgrid.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2222    dx=0.1 
    2323    dy=1.0 
    24     n=size(x)     
     24    n=size(x) 
    2525    xvec = arange(x0,x1,dx) 
    2626    yvec = arange(y0,y11,dy) 
    2727    nx = size(xvec) 
    2828    ny = size(yvec) 
    29     ix = zeros(n, float)     
    30     iy = zeros(n, float)     
     29    ix = zeros(n, float) 
     30    iy = zeros(n, float) 
    3131 
    3232    for kk in range(0,n): 
    … …  
    4242           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4343 
    44     inx = (ix >= 0) & (ix < nx)  
     44    inx = (ix >= 0) & (ix < nx) 
    4545    iny = (iy >= 0) & (iy < ny) 
    46     inz = (z >= z0) & (z <= z1)  
     46    inz = (z >= z0) & (z <= z1) 
    4747    inn = inx & iny & inz 
    4848    iix = ix[inn] 
    … …  
    6363        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    6464        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    65         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    66          
     65        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     66 
    6767    for i in range(0,ny): 
    6868        jj=round(3600/dxx[i]) 
    … …  
    7272        e=d.sum(axis=1) 
    7373        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    74         g=f.sum(axis=1)   
     74        g=f.sum(axis=1) 
    7575        for j in range(0, nx): 
    7676            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7979            zz2grid[j,i]=g[ii] 
    8080 
    81     zzpgrid = zeros([nx, ny], float)     
     81    zzpgrid = zeros([nx, ny], float) 
    8282    zzpgrid = zzgrid/nngrid 
    8383    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgrid1.py

    r6 r7  
    99 
    1010    dx=0.1 
    11     dy=1         
     11    dy=1 
    1212    dparx=0.05 
    1313    dpary=0.5 
    … …  
    2020    dx=0.1 
    2121    dy=1.0 
    22     n=size(x)     
     22    n=size(x) 
    2323    xvec = arange(x0,x1,dx) 
    2424    yvec = arange(y0,y11,dy) 
    2525    nx = size(xvec) 
    2626    ny = size(yvec) 
    27     ix = zeros(n, float)     
    28     iy = zeros(n, float)     
     27    ix = zeros(n, float) 
     28    iy = zeros(n, float) 
    2929 
    3030    for kk in range(0,n): 
    … …  
    4040           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4141 
    42     inx = (ix >= 0) & (ix < nx)  
     42    inx = (ix >= 0) & (ix < nx) 
    4343    iny = (iy >= 0) & (iy < ny) 
    44     inz = (z >= z0) & (z <= z1)  
     44    inz = (z >= z0) & (z <= z1) 
    4545    inn = inx & iny & inz 
    4646    iix = ix[inn] 
    … …  
    6161        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    6262        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    63         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    64          
     63        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     64 
    6565    for i in range(0,ny): 
    6666        jj=round(3600/dxx[i]) 
    … …  
    7070        e=d.sum(axis=1) 
    7171        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    72         g=f.sum(axis=1)   
     72        g=f.sum(axis=1) 
    7373        for j in range(0, nx): 
    7474            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7777            zz2grid[j,i]=g[ii] 
    7878 
    79     zzpgrid = zeros([nx, ny], float)     
     79    zzpgrid = zeros([nx, ny], float) 
    8080    zzpgrid = zzgrid/nngrid 
    8181    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgrid_z.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    5858        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    5959        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    60         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    61          
     60        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     61 
    6262    for i in range(0,ny): 
    6363        jj=round(3600/dxx[i]) 
    … …  
    6767        e=d.sum(axis=1) 
    6868        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    69         g=f.sum(axis=1)   
     69        g=f.sum(axis=1) 
    7070        for j in range(0, nx): 
    7171            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7474            zz2grid[j,i]=g[ii] 
    7575 
    76     zzpgrid = zeros([nx, ny], float)     
     76    zzpgrid = zeros([nx, ny], float) 
    7777    zzpgrid = zzgrid/nngrid 
    7878    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgrid_zm.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    6060        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    6161        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    62         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    63          
     62        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     63 
    6464    for i in range(0,ny): 
    6565        jj=round(3600/dxx[i]) 
    … …  
    6969        e=d.sum(axis=1) 
    7070        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    71         g=f.sum(axis=1)   
     71        g=f.sum(axis=1) 
    7272        for j in range(0, nx): 
    7373            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7676            zz2grid[j,i]=g[ii] 
    7777 
    78     zzpgrid = zeros([nx, ny], float)     
     78    zzpgrid = zeros([nx, ny], float) 
    7979    zzpgrid = zzgrid/nngrid 
    8080    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgrid_zml.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    6060        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    6161        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    62         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    63          
     62        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     63 
    6464    for i in range(0,ny): 
    6565        jj=round(3600/dxx[i]) 
    … …  
    6969        e=d.sum(axis=1) 
    7070        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    71         g=f.sum(axis=1)   
     71        g=f.sum(axis=1) 
    7272        for j in range(0, nx): 
    7373            ii=ceil((j+1)/dxx[i])-1 
    … …  
    7676            zz2grid[j,i]=g[ii] 
    7777 
    78     zzpgrid = zeros([nx, ny], float)     
     78    zzpgrid = zeros([nx, ny], float) 
    7979    zzpgrid = zzgrid/nngrid 
    8080    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgridn.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2222    dx=0.1 
    2323    dy=1.0 
    24     n=size(x)     
     24    n=size(x) 
    2525    xvec = arange(x0,x1,dx) 
    2626    yvec = arange(y0,y11,dy) 
    2727    nx = size(xvec) 
    2828    ny = size(yvec) 
    29     ix = zeros(n, float)     
    30     iy = zeros(n, float)     
     29    ix = zeros(n, float) 
     30    iy = zeros(n, float) 
    3131 
    3232    for kk in range(0,n): 
    … …  
    4242           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4343 
    44      
     44 
    4545    if nscan==1: 
    4646       ipos = ((pos <=5) | (pos >=26)) 
    47     if nscan==2:     
     47    if nscan==2: 
    4848       ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) 
    4949    if nscan==3: 
    5050       ipos = ((pos >=11) & (pos<=20)) 
    51      
    52     inx = (ix >= 0) & (ix < nx)  
     51 
     52    inx = (ix >= 0) & (ix < nx) 
    5353    iny = (iy >= 0) & (iy < ny) 
    54     inz = (z >= z0) & (z <= z1)  
     54    inz = (z >= z0) & (z <= z1) 
    5555    inn = inx & iny & inz & ipos 
    5656    iix = ix[inn] 
    … …  
    7171        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    7272        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    73         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    74          
     73        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     74 
    7575    for i in range(0,ny): 
    7676        jj=round(3600/dxx[i]) 
    … …  
    8080        e=d.sum(axis=1) 
    8181        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    82         g=f.sum(axis=1)   
     82        g=f.sum(axis=1) 
    8383        for j in range(0, nx): 
    8484            ii=ceil((j+1)/dxx[i])-1 
    … …  
    8787            zz2grid[j,i]=g[ii] 
    8888 
    89     zzpgrid = zeros([nx, ny], float)     
     89    zzpgrid = zeros([nx, ny], float) 
    9090    zzpgrid = zzgrid/nngrid 
    9191    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgridns.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2222    dx=0.1 
    2323    dy=1.0 
    24     n=size(x)     
     24    n=size(x) 
    2525    xvec = arange(x0,x1,dx) 
    2626    yvec = arange(y0,y11,dy) 
    2727    nx = size(xvec) 
    2828    ny = size(yvec) 
    29     ix = zeros(n, float)     
    30     iy = zeros(n, float)     
     29    ix = zeros(n, float) 
     30    iy = zeros(n, float) 
    3131 
    3232    for kk in range(0,n): 
    … …  
    4242           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4343 
    44      
     44 
    4545    if nscan==1: 
    4646       ipos = ((pos <=5) | (pos >=26)) 
    47     if nscan==2:     
     47    if nscan==2: 
    4848       ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) 
    4949    if nscan==3: 
    5050       ipos = ((pos >=11) & (pos<=20)) 
    51      
    52     inx = (ix >= 0) & (ix < nx)  
     51 
     52    inx = (ix >= 0) & (ix < nx) 
    5353    iny = (iy >= 0) & (iy < ny) 
    54     inz = (z >= z0) & (z <= z1)  
     54    inz = (z >= z0) & (z <= z1) 
    5555    inn = inx & iny & inz & ipos 
    5656    iix = ix[inn] 
    … …  
    7171        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    7272        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    73         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    74          
     73        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     74 
    7575    for i in range(0,ny): 
    7676        jj=round(3600/dxx[i]) 
    … …  
    8080        e=d.sum(axis=1) 
    8181        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    82         g=f.sum(axis=1)   
     82        g=f.sum(axis=1) 
    8383        for j in range(0, nx): 
    8484            ii=ceil((j+1)/dxx[i])-1 
    … …  
    8787            zz2grid[j,i]=g[ii] 
    8888 
    89     zzpgrid = zeros([nx, ny], float)     
     89    zzpgrid = zeros([nx, ny], float) 
    9090    zzpgrid = zzgrid/nngrid 
    9191    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgridns1.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2222    dx=0.1 
    2323    dy=1.0 
    24     n=size(x)     
     24    n=size(x) 
    2525    xvec = arange(x0,x1,dx) 
    2626    yvec = arange(y0,y11,dy) 
    2727    nx = size(xvec) 
    2828    ny = size(yvec) 
    29     ix = zeros(n, float)     
    30     iy = zeros(n, float)     
     29    ix = zeros(n, float) 
     30    iy = zeros(n, float) 
    3131 
    3232    for kk in range(0,n): 
    … …  
    4242           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4343 
    44      
     44 
    4545    if nscan==1: 
    4646       ipos = ((pos <=5) | (pos >=26)) 
    47     if nscan==2:     
     47    if nscan==2: 
    4848       ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) 
    4949    if nscan==3: 
    5050       ipos = ((pos >=11) & (pos<=20)) 
    51      
    52     inx = (ix >= 0) & (ix < nx)  
     51 
     52    inx = (ix >= 0) & (ix < nx) 
    5353    iny = (iy >= 0) & (iy < ny) 
    54     inz = (z >= z0) & (z <= z1)  
     54    inz = (z >= z0) & (z <= z1) 
    5555    inn = inx & iny & inz & ipos 
    5656    iix = ix[inn] 
    … …  
    7171        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    7272        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    73         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    74          
     73        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     74 
    7575    for i in range(0,ny): 
    7676        jj=round(3600/dxx[i]) 
    … …  
    8080        e=d.sum(axis=1) 
    8181        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    82         g=f.sum(axis=1)   
     82        g=f.sum(axis=1) 
    8383        for j in range(0, nx): 
    8484            ii=ceil((j+1)/dxx[i])-1 
    … …  
    8787            zz2grid[j,i]=g[ii] 
    8888 
    89     zzpgrid = zeros([nx, ny], float)     
     89    zzpgrid = zeros([nx, ny], float) 
    9090    zzpgrid = zzgrid/nngrid 
    9191    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgridns2.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2222    dx=0.1 
    2323    dy=1.0 
    24     n=size(x)     
     24    n=size(x) 
    2525    xvec = arange(x0,x1,dx) 
    2626    yvec = arange(y0,y11,dy) 
    2727    nx = size(xvec) 
    2828    ny = size(yvec) 
    29     ix = zeros(n, float)     
    30     iy = zeros(n, float)     
     29    ix = zeros(n, float) 
     30    iy = zeros(n, float) 
    3131 
    3232    for kk in range(0,n): 
    … …  
    4242           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4343 
    44      
     44 
    4545    if nscan==1: 
    4646       ipos = ((pos <=5) | (pos >=26)) 
    47     if nscan==2:     
     47    if nscan==2: 
    4848       ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) 
    4949    if nscan==3: 
    5050       ipos = ((pos >=11) & (pos<=20)) 
    51      
    52     inx = (ix >= 0) & (ix < nx)  
     51 
     52    inx = (ix >= 0) & (ix < nx) 
    5353    iny = (iy >= 0) & (iy < ny) 
    54     inz = (z >= z0) & (z <= z1)  
     54    inz = (z >= z0) & (z <= z1) 
    5555    inn = inx & iny & inz & ipos 
    5656    iix = ix[inn] 
    … …  
    7171        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    7272        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    73         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    74          
     73        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     74 
    7575    for i in range(0,ny): 
    7676        jj=round(3600/dxx[i]) 
    … …  
    8080        e=d.sum(axis=1) 
    8181        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    82         g=f.sum(axis=1)   
     82        g=f.sum(axis=1) 
    8383        for j in range(0, nx): 
    8484            ii=ceil((j+1)/dxx[i])-1 
    … …  
    8787            zz2grid[j,i]=g[ii] 
    8888 
    89     zzpgrid = zeros([nx, ny], float)     
     89    zzpgrid = zeros([nx, ny], float) 
    9090    zzpgrid = zzgrid/nngrid 
    9191    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newgridns_z.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    7070        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    7171        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    72          
     72 
    7373    for i in range(0,ny): 
    7474        jj=round(3600/dxx[i]) 

  • trunk/src/python_script/newgridns_zm.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 
    … …  
    2020    dx=0.1 
    2121    dy=1.0 
    22     n=size(x)     
     22    n=size(x) 
    2323    xvec = arange(x0,x1,dx) 
    2424    yvec = arange(y0,y11,dy) 
    2525    nx = size(xvec) 
    2626    ny = size(yvec) 
    27     ix = zeros(n, float)     
    28     iy = zeros(n, float)     
     27    ix = zeros(n, float) 
     28    iy = zeros(n, float) 
    2929 
    3030    for kk in range(0,n): 
    … …  
    4040           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    4141 
    42      
     42 
    4343    if nscan==1: 
    4444       ipos = ((pos <=5) | (pos >=26)) 
    45     if nscan==2:     
     45    if nscan==2: 
    4646       ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) 
    4747    if nscan==3: 
    4848       ipos = ((pos >=11) & (pos<=20)) 
    49      
    50     inx = (ix >= 0) & (ix < nx)  
     49 
     50    inx = (ix >= 0) & (ix < nx) 
    5151    iny = (iy >= 0) & (iy < ny) 
    52     inz = (z >= z0) & (z <= z1)  
     52    inz = (z >= z0) & (z <= z1) 
    5353    inn = inx & iny & inz & ipos 
    5454    iix = ix[inn] 
    … …  
    7070        zgrid[iix[i], iiy[i]] = zgrid[iix[i],iiy[i]]+zz[i] 
    7171        zdgrid[iix[i], iiy[i]] = zdgrid[iix[i],iiy[i]]+zz[i]*zz[i] 
    72         ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1  
    73          
     72        ngrid[iix[i], iiy[i]] = ngrid[iix[i],iiy[i]]+ 1 
     73 
    7474    for i in range(0,ny): 
    7575        jj=round(3600/dxx[i]) 
    … …  
    7979        e=d.sum(axis=1) 
    8080        f=reshape(zdgrid[:,i],(jj,dxx[i])) 
    81         g=f.sum(axis=1)   
     81        g=f.sum(axis=1) 
    8282        for j in range(0, nx): 
    8383            ii=ceil((j+1)/dxx[i])-1 
    … …  
    8686            zz2grid[j,i]=g[ii] 
    8787 
    88     zzpgrid = zeros([nx, ny], float)     
     88    zzpgrid = zeros([nx, ny], float) 
    8989    zzpgrid = zzgrid/nngrid 
    9090    sigma_grid=sqrt(zz2grid/nngrid-zzpgrid*zzpgrid) 

  • trunk/src/python_script/newvec_x.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    5050 
    5151    N = size(iix)# how many datapoints are left now? 
    52      
     52 
    5353    return zz, alt, mm, lo, la 

  • trunk/src/python_script/newvec_zml.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    4747 
    4848    N = size(iix)# how many datapoints are left now? 
    49      
     49 
    5050    return zz, mm 

  • trunk/src/python_script/newvec_zmla.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    4848 
    4949    N = size(iix)# how many datapoints are left now? 
    50      
     50 
    5151    return zz, alt, mm 

  • trunk/src/python_script/newvec_zmll.py

    r6 r7  
    88 
    99    dx=0.1 
    10     dy=1         
     10    dy=1 
    1111    dparx=0.05 
    1212    dpary=0.5 
    … …  
    1717    dx=0.1 
    1818    dy=1.0 
    19     n=size(x)     
     19    n=size(x) 
    2020    xvec = arange(x0,x1,dx) 
    2121    yvec = arange(y0,y11,dy) 
    2222    nx = size(xvec) 
    2323    ny = size(yvec) 
    24     ix = zeros(n, float)     
    25     iy = zeros(n, float)     
     24    ix = zeros(n, float) 
     25    iy = zeros(n, float) 
    2626 
    2727    for kk in range(0,n): 
    … …  
    3737           iy[kk] = math.ceil((y[kk] - y0)/dy)-1 
    3838 
    39     inx = (ix >= 0) & (ix < nx)  
     39    inx = (ix >= 0) & (ix < nx) 
    4040    iny = (iy >= 0) & (iy < ny) 
    41     inz = (z >= z0) & (z <= z1)  
     41    inz = (z >= z0) & (z <= z1) 
    4242    inn = inx & iny & inz 
    4343    iix = ix[inn] 
    … …  
    4949 
    5050    N = size(iix)# how many datapoints are left now? 
    51      
     51 
    5252    return zz, mm, llat, llon 

  • trunk/src/python_script/newvecns_zml.py

    r6 r7  
    1111    z=zz 
    1212    dx=0.1 
    13     dy=1         
     13    dy=1 
    1414    dparx=0.05 
    1515    dpary=0.5 

  • trunk/src/python_script/read_AMSUA.py

    r6 r7  
    4545    for j in range(0,20): 
    4646        amch1[j,iligne] = float(liste[j]) 
    47          
     47 
    4848 
    4949fichier.close 

  • trunk/src/python_script/read_AMSUB.py

    r6 r7  
    4343    for j in range(0,19): 
    4444        ambch1[j,iligne] = float(liste[j]) 
    45          
     45 
    4646 
    4747fichier.close 
    … …  
    5959    for j in range(0,19): 
    6060        ambch2[j,iligne] = float(liste[j]) 
    61          
     61 
    6262 
    6363fichier.close 
    … …  
    7474    for j in range(0,19): 
    7575        ambch3[j,iligne] = float(liste[j]) 
    76          
     76 
    7777 
    7878fichier.close 
    … …  
    8989    for j in range(0,19): 
    9090        ambch4[j,iligne] = float(liste[j]) 
    91          
     91 
    9292 
    9393fichier.close 
    … …  
    104104    for j in range(0,19): 
    105105        ambch5[j,iligne] = float(liste[j]) 
    106          
     106 
    107107 
    108108fichier.close 

  • trunk/src/python_script/read_SSMIS_CH15_june.py

    r6 r7  
    5353         trans[iligne] = float(liste[17]) 
    5454         iligne=iligne+1 
    55           
     55 
    5656 
    5757fichier.close 
    … …  
    171171zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    172172outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    173 monthly_outz_ssmis_polar[0,:,:]=outz     
     173monthly_outz_ssmis_polar[0,:,:]=outz 
    174174del outz, outx, outy, zz 
    175175 
    176176zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    177177outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    178 monthly_outz_ssmis_polar[1,:,:]=outz     
     178monthly_outz_ssmis_polar[1,:,:]=outz 
    179179del outz, outx, outy, zz 
    180180 
    181181zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    182182outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    183 monthly_outz_ssmis_polar[2,:,:]=outz     
     183monthly_outz_ssmis_polar[2,:,:]=outz 
    184184del outz, outx, outy, zz, xx, yy 
    185185 
    … …  
    189189zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    190190outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    191 monthly_outz_ssmis_diff[0,:,:]=outz     
     191monthly_outz_ssmis_diff[0,:,:]=outz 
    192192del outz, outx, outy, zz 
    193193 
    194194zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    195195outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    196 monthly_outz_ssmis_diff[1,:,:]=outz     
     196monthly_outz_ssmis_diff[1,:,:]=outz 
    197197del outz, outx, outy, zz 
    198198 
    199199zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    200200outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    201 monthly_outz_ssmis_diff[2,:,:]=outz     
     201monthly_outz_ssmis_diff[2,:,:]=outz 
    202202del outz, outx, outy, zz, xx, yy 
    203203 
    … …  
    221221 
    222222bins=arange(0.3,1,0.001) 
    223 bb=(lat_ssmis >= 0)  
    224   
     223bb=(lat_ssmis >= 0) 
     224 
    225225plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    226226plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
    … …  
    309309from netCDF4 import Dataset 
    310310rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    311   
     311 
    312312rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    313313rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 

  • trunk/src/python_script/read_SSMIS_CH16_june.py

    r6 r7  
    5353         trans16[iligne] = float(liste[17]) 
    5454         iligne=iligne+1 
    55           
     55 
    5656 
    5757fichier.close 
    … …  
    171171zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    172172outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    173 monthly_outz_ssmis_polar[0,:,:]=outz     
     173monthly_outz_ssmis_polar[0,:,:]=outz 
    174174del outz, outx, outy, zz 
    175175 
    176176zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    177177outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    178 monthly_outz_ssmis_polar[1,:,:]=outz     
     178monthly_outz_ssmis_polar[1,:,:]=outz 
    179179del outz, outx, outy, zz 
    180180 
    181181zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    182182outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    183 monthly_outz_ssmis_polar[2,:,:]=outz     
     183monthly_outz_ssmis_polar[2,:,:]=outz 
    184184del outz, outx, outy, zz, xx, yy 
    185185 
    … …  
    189189zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    190190outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    191 monthly_outz_ssmis_diff[0,:,:]=outz     
     191monthly_outz_ssmis_diff[0,:,:]=outz 
    192192del outz, outx, outy, zz 
    193193 
    194194zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    195195outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    196 monthly_outz_ssmis_diff[1,:,:]=outz     
     196monthly_outz_ssmis_diff[1,:,:]=outz 
    197197del outz, outx, outy, zz 
    198198 
    199199zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    200200outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    201 monthly_outz_ssmis_diff[2,:,:]=outz     
     201monthly_outz_ssmis_diff[2,:,:]=outz 
    202202del outz, outx, outy, zz, xx, yy 
    203203 
    … …  
    221221 
    222222bins=arange(0.3,1,0.001) 
    223 bb=(lat_ssmis >= 0)  
    224   
     223bb=(lat_ssmis >= 0) 
     224 
    225225plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    226226plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
    … …  
    309309from netCDF4 import Dataset 
    310310rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    311   
     311 
    312312rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    313313rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 

  • trunk/src/python_script/read_SSMIS_CH2_ANTARC_JUNE2010.DAT_gc.py

    r6 r7  
    4848         orog[iligne] = float(liste[11]) 
    4949         iligne=iligne+1 
    50           
     50 
    5151 
    5252fichier.close 
    … …  
    9090zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    9191outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    92 monthly_outz_ssmis_polar[0,:,:]=outz     
     92monthly_outz_ssmis_polar[0,:,:]=outz 
    9393del outz, outx, outy, zz 
    9494 
    9595zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    9696outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    97 monthly_outz_ssmis_polar[1,:,:]=outz     
     97monthly_outz_ssmis_polar[1,:,:]=outz 
    9898del outz, outx, outy, zz 
    9999 
    100100zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    101101outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    102 monthly_outz_ssmis_polar[2,:,:]=outz     
     102monthly_outz_ssmis_polar[2,:,:]=outz 
    103103del outz, outx, outy, zz, xx, yy 
    104104 
    … …  
    108108zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    109109outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    110 monthly_outz_ssmis_diff[0,:,:]=outz     
     110monthly_outz_ssmis_diff[0,:,:]=outz 
    111111del outz, outx, outy, zz 
    112112 
    113113zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    114114outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    115 monthly_outz_ssmis_diff[1,:,:]=outz     
     115monthly_outz_ssmis_diff[1,:,:]=outz 
    116116del outz, outx, outy, zz 
    117117 
    118118zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    119119outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    120 monthly_outz_ssmis_diff[2,:,:]=outz     
     120monthly_outz_ssmis_diff[2,:,:]=outz 
    121121del outz, outx, outy, zz, xx, yy 
    122122 
    … …  
    140140 
    141141bins=arange(0.3,1,0.001) 
    142 bb=(lat_ssmis >= 0)  
    143   
     142bb=(lat_ssmis >= 0) 
     143 
    144144plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    145145plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
    … …  
    228228from netCDF4 import Dataset 
    229229rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    230   
     230 
    231231rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    232232rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 

  • trunk/src/python_script/read_amsua_ch2_gc.py

    r6 r7  
    5555         trans_a2[iligne] = float(liste[19]) 
    5656         iligne=iligne+1 
    57           
     57 
    5858 
    5959fichier.close 
    … …  
    251251 
    252252#################################################" 
    253 #vec  
     253#vec 
    254254import newvec_zml 
    255255x=lon 
    … …  
    331331zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    332332outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    333 monthly_outz_ssmis_polar[0,:,:]=outz     
     333monthly_outz_ssmis_polar[0,:,:]=outz 
    334334del outz, outx, outy, zz 
    335335 
    336336zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    337337outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    338 monthly_outz_ssmis_polar[1,:,:]=outz     
     338monthly_outz_ssmis_polar[1,:,:]=outz 
    339339del outz, outx, outy, zz 
    340340 
    341341zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    342342outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    343 monthly_outz_ssmis_polar[2,:,:]=outz     
     343monthly_outz_ssmis_polar[2,:,:]=outz 
    344344del outz, outx, outy, zz, xx, yy 
    345345 
    … …  
    349349zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    350350outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    351 monthly_outz_ssmis_diff[0,:,:]=outz     
     351monthly_outz_ssmis_diff[0,:,:]=outz 
    352352del outz, outx, outy, zz 
    353353 
    354354zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    355355outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    356 monthly_outz_ssmis_diff[1,:,:]=outz     
     356monthly_outz_ssmis_diff[1,:,:]=outz 
    357357del outz, outx, outy, zz 
    358358 
    359359zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    360360outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    361 monthly_outz_ssmis_diff[2,:,:]=outz     
     361monthly_outz_ssmis_diff[2,:,:]=outz 
    362362del outz, outx, outy, zz, xx, yy 
    363363 
    … …  
    381381 
    382382bins=arange(0.3,1,0.001) 
    383 bb=(lat_ssmis >= 0)  
    384   
     383bb=(lat_ssmis >= 0) 
     384 
    385385plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    386386plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
    … …  
    469469from netCDF4 import Dataset 
    470470rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    471   
     471 
    472472rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    473473rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 

  • trunk/src/python_script/read_amsua_ch2_n.py

    r6 r7  
    4949         pos[iligne] = float(liste[7]) 
    5050         iligne=iligne+1 
    51           
     51 
    5252 
    5353fichier.close 
    … …  
    100100zz = 0.5*(emis_ssmis[1,:]+emis_ssmis[2,:]) 
    101101outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    102 monthly_outz_ssmis_polar[0,:,:]=outz     
     102monthly_outz_ssmis_polar[0,:,:]=outz 
    103103del outz, outx, outy, zz 
    104104 
    105105zz = 0.5*(emis_ssmis[4,:]+emis_ssmis[5,:]) 
    106106outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    107 monthly_outz_ssmis_polar[1,:,:]=outz     
     107monthly_outz_ssmis_polar[1,:,:]=outz 
    108108del outz, outx, outy, zz 
    109109 
    110110zz = 0.5*(emis_ssmis[6,:]+emis_ssmis[7,:]) 
    111111outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,0.1, 1) 
    112 monthly_outz_ssmis_polar[2,:,:]=outz     
     112monthly_outz_ssmis_polar[2,:,:]=outz 
    113113del outz, outx, outy, zz, xx, yy 
    114114 
    … …  
    118118zz = emis_ssmis[1,:]-emis_ssmis[2,:] 
    119119outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    120 monthly_outz_ssmis_diff[0,:,:]=outz     
     120monthly_outz_ssmis_diff[0,:,:]=outz 
    121121del outz, outx, outy, zz 
    122122 
    123123zz = emis_ssmis[4,:]-emis_ssmis[5,:] 
    124124outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    125 monthly_outz_ssmis_diff[1,:,:]=outz     
     125monthly_outz_ssmis_diff[1,:,:]=outz 
    126126del outz, outx, outy, zz 
    127127 
    128128zz = emis_ssmis[6,:]-emis_ssmis[7,:] 
    129129outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,-0.05, 0.2) 
    130 monthly_outz_ssmis_diff[2,:,:]=outz     
     130monthly_outz_ssmis_diff[2,:,:]=outz 
    131131del outz, outx, outy, zz, xx, yy 
    132132 
    … …  
    150150 
    151151bins=arange(0.3,1,0.001) 
    152 bb=(lat_ssmis >= 0)  
    153   
     152bb=(lat_ssmis >= 0) 
     153 
    154154plt.hist(emis_ssmis[0,nonzero(bb)[0]], bins=bins,histtype='step', label='e50V',normed='True',color='#4BB5C1') 
    155155plt.hist(emis_ssmis[1,nonzero(bb)[0]], bins=bins,histtype='step', label='e19V',normed='True',color='black') 
    … …  
    238238from netCDF4 import Dataset 
    239239rootgrp = Dataset('..\EMIS\EMIS_SSMIS_'+le_mois+'.nc', 'w', format='NETCDF4') 
    240   
     240 
    241241rootgrp.createDimension('longitude', len(monthly_lon_ssmis)) 
    242242rootgrp.createDimension('latitude', len(monthly_lat_ssmis)) 

  • trunk/src/python_script/vnx.py

    r6 r7  
    88Rt=6371 
    99beta=math.pi/180 
    10 R=beta*Rt  
     10R=beta*Rt 
    1111L=zeros(ny,float)#arc vector 
    1212L[0]=2*math.pi*R 

  • trunk/src/python_script/vnx1.py

    r6 r7  
    2323 
    2424########################## 
    25 inlat = (y >= -90) & (y < -88)  
    26 inz = (z >= 100) & (z <= 300)  
     25inlat = (y >= -90) & (y < -88) 
     26inz = (z >= 100) & (z <= 300) 
    2727innn = inlat & inz 
    2828 
Note: See TracChangeset for help on using the changeset viewer.