Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/scripts_Laura/read_SSMIS_test.py

    r20 r23  
    77from mpl_toolkits.basemap import Basemap 
    88from mpl_toolkits.basemap import shiftgrid, cm 
    9 import netCDF4 
     9from netCDF4 import Dataset 
    1010import ffgrid2 
    1111 
    1212 
    13 ################ fichiers par mois - CH1 ################################################### 
     13 
     14################ fichiers par canal - mois de juin ################################################### 
     15 
    1416 
    1517f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH' 
    1618f3 = '_ANTARC_JUNE2010.DAT' 
    17 #date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    1819channel = np.array([12, 13, 15, 16, 17, 18]) 
    1920numlines = np.zeros([len(channel)],int) 
    … …  
    2627     for line in fichier: numlines[ich] += 1 
    2728 
    28  
    29      fichier.close() 
    30  
    31 ich = 2 # 37GHz, H polar 
     29      fichier.close() 
     30 
     31 
     32 
     33 
     34 
     35ich = 0 # 37GHz, H  
    3236fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    3337ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    159163 
    160164 
    161 ich = 5 # 91.66GHz, V polar 
     165ich = 5 # 91.66GHz, H polar 
    162166fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    163167ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    184188orog18_JUN=ssch18_JUN[11,:] 
    185189 
    186 ######################## 
    187 # EVOLUTION TEMPORELLE # 
    188 ######################## 
    189 ## Autour de "Dome C" (lon=123.23;lat=-75.06) - mois: JUNE ## 
    190 ## Autre zone de glace de mer / de continent (lon entre -40 et -60 // lat entre -75 et -85)##  
    191 ## ch17 ## 
    192 #bbzone_CH17_JUN = nonzero((lon17_JUN>=120.23)&(lon17_JUN<=126.23)&(lat17_JUN>=-78.06)&(lat17_JUN<=-72.06)) 
    193 bbzone_CH17_JUN = nonzero((lon17_JUN >= -60.) & (lon17_JUN <= -40.) & (lat17_JUN >= -85.) & (lat17_JUN <= -75.)) 
    194 emis17_JUN_moy = np.zeros([30],float) 
    195 tb17_JUN_moy = np.zeros([30],float) 
    196 ts17_JUN_moy = np.zeros([30],float) 
    197  
    198 ## ch18 ## 
    199 #bbzone_CH18_JUN = nonzero((lon18_JUN >= 120.23) & (lon18_JUN <= 126.23) & (lat18_JUN >= -78.06) & (lat18_JUN <= -72.06)) 
    200 bbzone_CH18_JUN = nonzero((lon18_JUN >= -60.) & (lon18_JUN <= -40.) & (lat18_JUN >= -85.) & (lat18_JUN <= -75.)) 
    201 emis18_JUN_moy = np.zeros([30],float) 
    202 tb18_JUN_moy = np.zeros([30],float) 
    203 ts18_JUN_moy = np.zeros([30],float) 
    204  
    205 ## ch 16 ## 
    206 #bbzone_CH16_JUN = nonzero((lon16_JUN >= 120.23) & (lon16_JUN <= 126.23) & (lat16_JUN >= -78.06) & (lat16_JUN <= -72.06)) 
    207 #emis16_JUN_moy = np.zeros([30],float) 
    208 #tb16_JUN_moy = np.zeros([30],float) 
    209 #ts16_JUN_moy = np.zeros([30],float) 
    210  
    211 for ijr in range (0,30): 
    212      print 'jour=', ijr+1 
    213      ## ch17 ## 
    214      ind_jr17 = np.where(jjr17_JUN[bbzone_CH17_JUN]==ijr+1)[0] 
    215      a = emis17_JUN[bbzone_CH17_JUN][ind_jr17] 
    216      b = tb17_JUN[bbzone_CH17_JUN][ind_jr17] 
    217      c = ts17_JUN[bbzone_CH17_JUN][ind_jr17] 
    218      emis17_JUN_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))]) 
    219      tb17_JUN_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))]) 
    220      ts17_JUN_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))]) 
    221      ## ch18 ## 
    222      ind_jr18 = np.where(jjr18_JUN[bbzone_CH18_JUN]==ijr+1)[0] 
    223      d = emis18_JUN[bbzone_CH18_JUN][ind_jr18] 
    224      e = tb18_JUN[bbzone_CH18_JUN][ind_jr18] 
    225      f = ts18_JUN[bbzone_CH18_JUN][ind_jr18] 
    226      emis18_JUN_moy[ijr] = mean(d[nonzero((d != -500.) & (d <= 1.))]) 
    227      tb18_JUN_moy[ijr] = mean(e[nonzero((e != -500.) & (e != 0.))]) 
    228      ts18_JUN_moy[ijr] = mean(f[nonzero((f != -500.) & (f != 0.))]) 
    229      ## ch16 ## 
    230 #     ind_jr16 = np.where(jjr16_JUN[bbzone_CH16_JUN]==ijr+1)[0] 
    231 #     g = emis16_JUN[bbzone_CH16_JUN][ind_jr16] 
    232 #     h = tb16_JUN[bbzone_CH16_JUN][ind_jr16] 
    233 #     l = ts16_JUN[bbzone_CH16_JUN][ind_jr16] 
    234 #     emis16_JUN_moy[ijr] = mean(g[nonzero((g != -500.) & (g <= 1.))]) 
    235 #     tb16_JUN_moy[ijr] = mean(h[nonzero((h != -500.) & (h != 0.))]) 
    236 #     ts16_JUN_moy[ijr] = mean(l[nonzero((l != -500.) & (l != 0.))]) 
    237  
    238  
    239 ## ch17 - ch18 ## 
    240 fig=plt.figure() 
    241 plt.subplot(3,1,1) 
    242 plt.plot(arange(0,30,1),emis17_JUN_moy,label='91.66GHz-V',c='r') 
    243 plt.plot(arange(0,30,1),emis18_JUN_moy,label='91.66GHz-H',c='b') 
    244 plt.xticks(arange(1,31,1)) 
    245 grid(True) 
    246 ylabel('emissivity') 
    247 legend(loc=7) 
    248 plt.title('SSMIS - JUNE 2010') 
    249 plt.subplot(3,1,2) 
    250 plt.plot(arange(0,30,1),tb17_JUN_moy,label='91.66GHz-V',c='r') 
    251 plt.plot(arange(0,30,1),tb18_JUN_moy,label='91.66GHz-H',c='b') 
    252 plt.xticks(arange(1,31,1)) 
    253 grid(True) 
    254 ylabel('brightness temperature') 
    255 legend(loc=7) 
    256 plt.subplot(3,1,3) 
    257 plt.plot(arange(0,30,1),ts17_JUN_moy,label='91.66GHz-V',c='r') 
    258 plt.plot(arange(0,30,1),ts18_JUN_moy,label='91.66GHz-H',c='b') 
    259 plt.xticks(arange(1,31,1)) 
    260 grid(True) 
    261 ylabel('surface temperature') 
    262 legend(loc=7) 
    263 fig.show() 
    264  
    265 ## ch16 ## 
    266 fig=plt.figure() 
    267 plt.subplot(3,1,1) 
    268 plt.plot(arange(0,30,1),emis16_JUN_moy,label='37GHz-V',c='r') 
    269 plt.plot(arange(0,30,1),emis17_JUN_moy,label='91.66GHz-V',c='g') 
    270 plt.xticks(arange(1,31,1)) 
    271 grid(True) 
    272 ylabel('emissivity') 
    273 legend(loc=7) 
    274 plt.title('SSMIS - JUNE 2010') 
    275 plt.subplot(3,1,2) 
    276 plt.plot(arange(0,30,1),tb16_JUN_moy,label='37GHz-V',c='r') 
    277 plt.plot(arange(0,30,1),tb17_JUN_moy,label='91.66GHz-V',c='g') 
    278 plt.xticks(arange(1,31,1)) 
    279 grid(True) 
    280 ylabel('brightness temperature') 
    281 legend(loc=7) 
    282 plt.subplot(3,1,3) 
    283 plt.plot(arange(0,30,1),ts16_JUN_moy,label='37GHz-V',c='r') 
    284 plt.plot(arange(0,30,1),ts17_JUN_moy,label='91.66GHz-V',c='g') 
    285 plt.xticks(arange(1,31,1)) 
    286 grid(True) 
    287 ylabel('surface temperature') 
    288 legend(loc=7) 
    289 fig.show() 
    290  
    291  
    292  
    293 tb17_JUN_anom = np.zeros([30],float) 
    294 tb18_JUN_anom = np.zeros([30],float) 
    295 for ijr in range (0,30): 
    296      print 'jour=', ijr + 1 
    297      tb17_JUN_anom[ijr]=tb17_JUN_moy[ijr]-tb17_JUN_moy.mean() 
    298      tb18_JUN_anom[ijr]=tb18_JUN_moy[ijr]-tb18_JUN_moy.mean() 
    299  
    300  
    301 fig=plt.figure() 
    302 plt.plot(arange(0,30,1),tb17_JUN_anom,label='91.66Ghz-V',c='r') 
    303 plt.plot(arange(0,30,1),tb18_JUN_anom,label='91.66GHz-H',c='b') 
    304 plt.plot(arange(0,31,1),np.zeros([31]),'--',c='k') 
    305 ylabel('Tb anomaly') 
    306 legend() 
    307 twinx().plot(arange(0,30,1),ts17_JUN_moy,label='surface temperature',c='g') 
    308 ylabel('surface temperature') 
    309 plt.xticks(arange(0, 30, 1), arange(1, 31, 1)) 
    310 grid(True) 
    311 plt.title('SSMIS - JUNE 2010') 
    312 fig.show() 
    313  
    314  
    315  
    316  
    317 ################ 
    318 # CARTOGRAPHIE # 
    319 ################ 
    320 ## Etude sur l'Antarctique ## 
    321 dx=5. 
    322 dy=5. 
    323 x0, x1 = -180, 180 
    324 y0, y1 = -90, -30 
    325  
    326 ## ch17 ## 
    327 bbemis_ch17_JUN=nonzero((emis17_JUN!=-500.)&(emis17_JUN<1.)&(emis17_JUN>0.)) 
    328 bbtb_ch17_JUN=nonzero((tb17_JUN!=-500.)&(tb17_JUN!=0.)) 
    329 OUTZCH17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float) 
    330 outzch17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float) 
    331 lonch17_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float) 
    332 latch17_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float) 
    333  
    334 ## ch18 ## 
    335 bbemis_ch18_JUN=nonzero((emis18_JUN!=-500.)&(emis18_JUN<1.)&(emis18_JUN>0.)) 
    336 bbtb_ch18_JUN=nonzero((tb18_JUN!=-500.)&(tb18_JUN!=0.)) 
    337 OUTZCH18_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float) 
    338 outzch18_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float) 
    339 lonch18_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float) 
    340 latch18_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float) 
    341  
    342 for ijr in range(0,30): 
    343     print 'jour=', ijr+1 
    344     ## ch17 ## 
    345     ind_jr17_JUN = np.where(jjr17_JUN[bbtb_ch17_JUN] == ijr+1)[0] 
    346     xx = lon17_JUN[bbtb_ch17_JUN][ind_jr17_JUN] 
    347     yy = lat17_JUN[bbtb_ch17_JUN][ind_jr17_JUN] 
    348     zz = tb17_JUN[bbtb_ch17_JUN][ind_jr17_JUN] 
    349     zz0 = min(zz) 
    350     zz1 = max(zz) 
    351     outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1) 
    352     outzch17_JUN = outz 
    353     lonch17_JUN = outx 
    354     latch17_JUN = outy 
    355     OUTZCH17_JUN[:,:,ijr] = outzch17_JUN[:,:] 
    356     ## ch18 ## 
    357     ind_jr18_JUN = np.where(jjr18_JUN[bbtb_ch18_JUN] == ijr+1)[0] 
    358     xx = lon18_JUN[bbtb_ch18_JUN][ind_jr18_JUN] 
    359     yy = lat18_JUN[bbtb_ch18_JUN][ind_jr18_JUN] 
    360     zz = tb18_JUN[bbtb_ch18_JUN][ind_jr18_JUN] 
    361     zz0 = min(zz) 
    362     zz1 = max(zz) 
    363     outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1) 
    364     outzch18_JUN = outz 
    365     lonch18_JUN = outx 
    366     latch18_JUN = outy 
    367     OUTZCH18_JUN[:,:,ijr] = outzch18_JUN[:,:] 
    368  
    369      
    370 ## ch17 ## 
    371 mean_outzch17_JUN = np.zeros([len(latch17_JUN), len(lonch17_JUN)], float) 
    372 for ilon in range (0,len(lonch17_JUN)): 
    373     for ilat in range (0,len(latch17_JUN)): 
    374         mean_outzch17_JUN[ilat,ilon] = OUTZCH17_JUN[ilat,ilon,:].mean() 
    375  
    376  
    377 ## ch18 ## 
    378 mean_outzch18_JUN = np.zeros([len(latch18_JUN), len(lonch18_JUN)], float) 
    379 for ilon in range (0,len(lonch18_JUN)): 
    380     for ilat in range (0,len(latch18_JUN)): 
    381         mean_outzch18_JUN[ilat,ilon] = OUTZCH18_JUN[ilat,ilon,:].mean() 
    382  
    383  
    384 ## ch17 ## 
    385 tbch17_anom_JUN = np.zeros([len(latch17_JUN),len(lonch17_JUN),30], float) 
    386 for ijr in range (0,30): 
    387     tbch17_anom_JUN[:,:,ijr] = OUTZCH17_JUN[:,:,ijr] - mean_outzch17_JUN[:,:] 
    388  
    389  
    390 ## ch18 ## 
    391 tbch18_anom_JUN = np.zeros([len(latch18_JUN),len(lonch18_JUN),30], float) 
    392 for ijr in range (0,30): 
    393     tbch18_anom_JUN[:,:,ijr] = OUTZCH18_JUN[:,:,ijr] - mean_outzch18_JUN[:,:] 
    394  
    395  
    396 for ijr in range (23, 30): 
    397      figure() 
    398      plt.ion() 
    399      m = Basemap(llcrnrlon=-60, urcrnrlon=-40, llcrnrlat=-85, urcrnrlat=-75, projection='cyl', resolution='c', fix_aspect=True) 
    400      m.drawcoastlines(linewidth = 1) 
    401      m.drawparallels(np.arange(-85., 75., 2)) 
    402      m.drawmeridians(np.arange(-60., -40., 2)) 
    403      #m.fillcontinents() 
    404      clevs = arange(-25., 5., 0.1) 
    405      xii,yii = m(*np.meshgrid(lonch17_JUN, latch17_JUN)) 
    406      cs = m.contourf(xii, yii, tbch17_anom_JUN[:,:,ijr], clevs, cmap=cm.s3pcpn_l_r) 
    407      cbar = colorbar(cs) 
    408      cbar.set_label('Tb anomaly SSMIS CH17 - JUNE') 
    409      xticks(np.arange(-60., -40., 2)) 
    410      yticks(np.arange(-85., 75., 2)) 
    411      plt.savefig('/usr/home/lahlod/twice_d/figures_output_ANTARC/SSMIS/mean_tb_anomaly_map_'+str(ijr+1)+'JUN_ch17_SSMIS_Zoom_zone2') 
    412  
    413  
    414 figure() 
    415 plt.ion() 
    416 m = Basemap(llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-90, urcrnrlat=-30, projection='cyl', resolution='c', fix_aspect=True) 
    417 m.drawcoastlines(linewidth = 1) 
    418 m.drawparallels(np.arange(-90., 90., 20)) 
    419 m.drawmeridians(np.arange(-180., 180., 20)) 
    420 #m.fillcontinents() 
    421 xii,yii = m(*np.meshgrid(lonch17_JUN, latch17_JUN)) 
    422  
    423 ## BIAIS ch17-ch18 ##  
    424 biais_JUN = mean_outzch17_JUN - mean_outzch18_JUN 
    425 clevs = arange(0., 45., 1.) 
    426 cs = m.contourf(xii, yii, biais_JUN, clevs, cmap=cm.s3pcpn_l_r) 
    427 cbar = colorbar(cs) 
    428 cbar.set_label('Bias [CH17-CH18] of Mean Tb JUNE - SSMIS') 
    429  
    430 ## VARIANCE ch17 - ch18 ## 
    431 std1 = np.zeros([len(latch17_JUN), len(lonch17_JUN)], float) 
    432 for ilat in range (0, len(latch17_JUN)): 
    433     for ilon in range (0, len(lonch17_JUN)): 
    434         std1[ilat, ilon] = (mean_outzch17_JUN[ilat, ilon]**2)-(mean_outzch18_JUN[ilat, ilon]**2) 
    435  
    436  
    437 N = len(lonch17_JUN)*len(latch17_JUN) 
    438 std_JUN = std1/N 
    439 clevs = arange(0., 25., 0.1) 
    440 cs = m.contourf(xii, yii, std_JUN, clevs, cmap=cm.s3pcpn_l_r) 
    441 cbar = colorbar(cs) 
    442 cbar.set_label('Stantard deviation [CH17-CH18] of Mean Tb JUNE - SSMIS') 
    443  
    444  
    445  
    446 ########################## 
    447 # DIAGRAMME DE HOVMOLLER # 
    448 ########################## 
    449  
    450 # shape(tbch17_anom_JUN) = [ilat, ilon, ijr] 
    451  
    452 ## tranche de latitude étudiée ## 
    453 #bbtranche17_JUN = nonzero((latch17_JUN == -75.)) 
    454 bbtranche17_JUN = nonzero((latch17_JUN >= -85.) & (latch17_JUN <= -75)) 
    455 #bbtranche17_JUN = nonzero((latch17_JUN >= -90.) & (latch17_JUN <= -85)) 
    456 mean_tbch17_anom_JUN = np.zeros([len(lonch17_JUN), 30], float) 
    457 for ilon in range (0, len(lonch17_JUN)): 
    458     for ijr in range (0,30): 
    459         mean_tbch17_anom_JUN[ilon,ijr] = tbch17_anom_JUN[bbtranche17_JUN][:,ilon,ijr].mean() 
    460  
    461  
    462 y_time, x_space = m(*np.meshgrid(arange(0,30,1), lonch17_JUN)) 
    463 fig = plt.figure() 
    464 plt.pcolor(x_space, y_time, mean_tbch17_anom_JUN, cmap=cm.s3pcpn_l_r, vmin = -30., vmax = 25.) 
    465 plt.axis([-180., 180., 0, 29]) 
    466 cb = plt.colorbar() 
    467 cb.set_label('Tb anomaly - SSMIS CH17') 
    468 plt.xticks(arange(-180.,220.,40)) 
    469 plt.yticks(arange(0, 30, 1), arange(1,31,1)) 
    470 plt.xlabel('longitude') 
    471 plt.ylabel('JUNE 2010') 
    472  
     190 
     191 
     192 
     193############### fichiers par mois pour deux canaux (polars H et V) ################################## 
     194 
     195######## 
     196# ch17 # 
     197######## 
     198 
     199f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH17_ANTARC_' 
     200f3 = '2010.DAT' 
     201date=np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     202numlines = np.zeros([len(date)],int) 
     203 
     204for imo in range (0, len(date)): 
     205     print date[imo]  
     206     f = f1 + str(date[imo]) + f3 
     207     fichier = open(f, 'r') 
     208     numlines[imo] = 0 
     209     for line in fichier: numlines[imo] += 1 
     210 
     211      fichier.close() 
     212 
     213 
     214imo = 0 # FEBUARY 
     215fichier = open(f1 + str(date[imo]) + f3, 'r') 
     216ssmis = np.zeros([18, numlines[imo]], float) 
     217for iligne in range (0,numlines[imo]): 
     218    line = fichier.readline() 
     219    liste = line.split() 
     220    for j in range(0,18): 
     221        ssmis[j,iligne] = float(liste[j]) 
     222 
     223    fichier.close 
     224 
     225ssch17_FEB=ssmis 
     226lon17_FEB=ssch17_FEB[0,:] 
     227lat17_FEB=ssch17_FEB[1,:] 
     228jjr17_FEB=ssch17_FEB[4,:] 
     229ts17_FEB=ssch17_FEB[8,:] 
     230emis17_FEB=ssch17_FEB[14,:] 
     231tb17_FEB=ssch17_FEB[13,:] 
     232tup17_FEB=ssch17_FEB[16,:] 
     233tdn17_FEB=ssch17_FEB[15,:] 
     234trans17_FEB=ssch17_FEB[17,:] 
     235orog17_FEB=ssch17_FEB[11,:] 
     236 
     237 
     238imo = 1 # APRIL 
     239fichier = open(f1 + str(date[imo]) + f3, 'r') 
     240ssmis = np.zeros([18, numlines[imo]], float) 
     241for iligne in range (0,numlines[imo]): 
     242    line = fichier.readline() 
     243    liste = line.split() 
     244    for j in range(0,18): 
     245        ssmis[j,iligne] = float(liste[j]) 
     246 
     247    fichier.close 
     248 
     249ssch17_APR=ssmis 
     250lon17_APR=ssch17_APR[0,:] 
     251lat17_APR=ssch17_APR[1,:] 
     252jjr17_APR=ssch17_APR[4,:] 
     253ts17_APR=ssch17_APR[8,:] 
     254emis17_APR=ssch17_APR[14,:] 
     255tb17_APR=ssch17_APR[13,:] 
     256tup17_APR=ssch17_APR[16,:] 
     257tdn17_APR=ssch17_APR[15,:] 
     258trans17_APR=ssch17_APR[17,:] 
     259orog17_APR=ssch17_APR[11,:] 
     260 
     261 
     262imo = 2 # MAY 
     263fichier = open(f1 + str(date[imo]) + f3, 'r') 
     264ssmis = np.zeros([18, numlines[imo]], float) 
     265for iligne in range (0,numlines[imo]): 
     266    line = fichier.readline() 
     267    liste = line.split() 
     268    for j in range(0,18): 
     269        ssmis[j,iligne] = float(liste[j]) 
     270 
     271    fichier.close 
     272 
     273ssch17_MAY=ssmis 
     274lon17_MAY=ssch17_MAY[0,:] 
     275lat17_MAY=ssch17_MAY[1,:] 
     276jjr17_MAY=ssch17_MAY[4,:] 
     277ts17_MAY=ssch17_MAY[8,:] 
     278emis17_MAY=ssch17_MAY[14,:] 
     279tb17_MAY=ssch17_MAY[13,:] 
     280tup17_MAY=ssch17_MAY[16,:] 
     281tdn17_MAY=ssch17_MAY[15,:] 
     282trans17_MAY=ssch17_MAY[17,:] 
     283orog17_MAY=ssch17_MAY[11,:] 
     284 
     285 
     286imo = 3 # JUNE 
     287fichier = open(f1 + str(date[imo]) + f3, 'r') 
     288ssmis = np.zeros([18, numlines[imo]], float) 
     289for iligne in range (0,numlines[imo]): 
     290    line = fichier.readline() 
     291    liste = line.split() 
     292    for j in range(0,18): 
     293        ssmis[j,iligne] = float(liste[j]) 
     294 
     295    fichier.close 
     296 
     297ssch17_JUN=ssmis 
     298lon17_JUN=ssch17_JUN[0,:] 
     299lat17_JUN=ssch17_JUN[1,:] 
     300jjr17_JUN=ssch17_JUN[4,:] 
     301ts17_JUN=ssch17_JUN[8,:] 
     302emis17_JUN=ssch17_JUN[14,:] 
     303tb17_JUN=ssch17_JUN[13,:] 
     304tup17_JUN=ssch17_JUN[16,:] 
     305tdn17_JUN=ssch17_JUN[15,:] 
     306trans17_JUN=ssch17_JUN[17,:] 
     307orog17_JUN=ssch17_JUN[11,:] 
     308 
     309 
     310imo = 4 # JULY 
     311fichier = open(f1 + str(date[imo]) + f3, 'r') 
     312ssmis = np.zeros([18, numlines[imo]], float) 
     313for iligne in range (0,numlines[imo]): 
     314    line = fichier.readline() 
     315    liste = line.split() 
     316    for j in range(0,18): 
     317        ssmis[j,iligne] = float(liste[j]) 
     318 
     319    fichier.close 
     320 
     321ssch17_JUL=ssmis 
     322lon17_JUL=ssch17_JUL[0,:] 
     323lat17_JUL=ssch17_JUL[1,:] 
     324jjr17_JUL=ssch17_JUL[4,:] 
     325ts17_JUL=ssch17_JUL[8,:] 
     326emis17_JUL=ssch17_JUL[14,:] 
     327tb17_JUL=ssch17_JUL[13,:] 
     328tup17_JUL=ssch17_JUL[16,:] 
     329tdn17_JUL=ssch17_JUL[15,:] 
     330trans17_JUL=ssch17_JUL[17,:] 
     331orog17_JUL=ssch17_JUL[11,:] 
     332 
     333 
     334######## 
     335# ch18 # 
     336######## 
     337 
     338f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH18_ANTARC_' 
     339f3 = '2010.DAT' 
     340date=np.array(['FEBRUARY', 'APRIL', 'JULY']) 
     341numlines = np.zeros([len(date)],int) 
     342 
     343for imo in range (0, len(date)): 
     344     print date[imo]  
     345     f = f1 + str(date[imo]) + f3 
     346     fichier = open(f, 'r') 
     347     numlines[imo] = 0 
     348     for line in fichier: numlines[imo] += 1 
     349 
     350      fichier.close() 
     351 
     352 
     353imo = 0 # FEBUARY 
     354fichier = open(f1 + str(date[imo]) + f3, 'r') 
     355ssmis = np.zeros([18, numlines[imo]], float) 
     356for iligne in range (0,numlines[imo]): 
     357    line = fichier.readline() 
     358    liste = line.split() 
     359    for j in range(0,18): 
     360        ssmis[j,iligne] = float(liste[j]) 
     361 
     362    fichier.close 
     363 
     364ssch18_FEB=ssmis 
     365lon18_FEB=ssch18_FEB[0,:] 
     366lat18_FEB=ssch18_FEB[1,:] 
     367jjr18_FEB=ssch18_FEB[4,:] 
     368ts18_FEB=ssch18_FEB[8,:] 
     369emis18_FEB=ssch18_FEB[14,:] 
     370tb18_FEB=ssch18_FEB[13,:] 
     371tup18_FEB=ssch18_FEB[16,:] 
     372tdn18_FEB=ssch18_FEB[15,:] 
     373trans18_FEB=ssch18_FEB[17,:] 
     374orog18_FEB=ssch18_FEB[11,:] 
     375 
     376 
     377imo = 1 # APRIL 
     378fichier = open(f1 + str(date[imo]) + f3, 'r') 
     379ssmis = np.zeros([18, numlines[imo]], float) 
     380for iligne in range (0,numlines[imo]): 
     381    line = fichier.readline() 
     382    liste = line.split() 
     383    for j in range(0,18): 
     384        ssmis[j,iligne] = float(liste[j]) 
     385 
     386    fichier.close 
     387 
     388ssch18_APR=ssmis 
     389lon18_APR=ssch18_APR[0,:] 
     390lat18_APR=ssch18_APR[1,:] 
     391jjr18_APR=ssch18_APR[4,:] 
     392ts18_APR=ssch18_APR[8,:] 
     393emis18_APR=ssch18_APR[14,:] 
     394tb18_APR=ssch18_APR[13,:] 
     395tup18_APR=ssch18_APR[16,:] 
     396tdn18_APR=ssch18_APR[15,:] 
     397trans18_APR=ssch18_APR[17,:] 
     398orog18_APR=ssch18_APR[11,:] 
     399 
     400 
     401 
     402imo = 2 # JULY 
     403fichier = open(f1 + str(date[imo]) + f3, 'r') 
     404ssmis = np.zeros([18, numlines[imo]], float) 
     405for iligne in range (0,numlines[imo]): 
     406    line = fichier.readline() 
     407    liste = line.split() 
     408    for j in range(0,18): 
     409        ssmis[j,iligne] = float(liste[j]) 
     410 
     411    fichier.close 
     412 
     413ssch18_JUL=ssmis 
     414lon18_JUL=ssch18_JUL[0,:] 
     415lat18_JUL=ssch18_JUL[1,:] 
     416jjr18_JUL=ssch18_JUL[4,:] 
     417ts18_JUL=ssch18_JUL[8,:] 
     418emis18_JUL=ssch18_JUL[14,:] 
     419tb18_JUL=ssch18_JUL[13,:] 
     420tup18_JUL=ssch18_JUL[16,:] 
     421tdn18_JUL=ssch18_JUL[15,:] 
     422trans18_JUL=ssch18_JUL[17,:] 
     423orog18_JUL=ssch18_JUL[11,:] 
     424 
     425 
     426########## 
     427## ch12 ## 
     428########## 
     429 
     430f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH12_ANTARC_' 
     431f3 = '2010.DAT' 
     432month = np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     433numlines = np.zeros([len(month)],int) 
     434 
     435for imo in range (0, len(month)): 
     436     print month[imo]  
     437     f = f1 + str(month[imo]) + f3 
     438     fichier = open(f, 'r') 
     439     numlines[imo] = 0 
     440     for line in fichier: numlines[imo] += 1 
     441 
     442      fichier.close() 
     443 
     444 
     445imo = 0 # FEB  
     446fichier = open(f1 + str(month[imo]) + f3, 'r') 
     447ssmis = np.zeros([18, numlines[imo]], float) 
     448for iligne in range (0,numlines[imo]): 
     449    line = fichier.readline() 
     450    liste = line.split() 
     451    for j in range(0,18): 
     452        ssmis[j,iligne] = float(liste[j]) 
     453 
     454 
     455    fichier.close 
     456 
     457 
     458ssch12_FEB=ssmis 
     459lon12_FEB=ssch12_FEB[0,:] 
     460lat12_FEB=ssch12_FEB[1,:] 
     461jjr12_FEB=ssch12_FEB[4,:] 
     462ts12_FEB=ssch12_FEB[8,:] 
     463emis12_FEB=ssch12_FEB[14,:] 
     464tb12_FEB=ssch12_FEB[13,:] 
     465tup12_FEB=ssch12_FEB[16,:] 
     466tdn12_FEB=ssch12_FEB[15,:] 
     467trans12_FEB=ssch12_FEB[17,:] 
     468orog12_FEB=ssch12_FEB[11,:] 
     469 
     470imo = 1 # APR  
     471fichier = open(f1 + str(month[imo]) + f3, 'r') 
     472ssmis = np.zeros([18, numlines[imo]], float) 
     473for iligne in range (0,numlines[imo]): 
     474    line = fichier.readline() 
     475    liste = line.split() 
     476    for j in range(0,18): 
     477        ssmis[j,iligne] = float(liste[j]) 
     478 
     479 
     480    fichier.close 
     481 
     482ssch12_APR=ssmis 
     483lon12_APR=ssch12_APR[0,:] 
     484lat12_APR=ssch12_APR[1,:] 
     485jjr12_APR=ssch12_APR[4,:] 
     486ts12_APR=ssch12_APR[8,:] 
     487emis12_APR=ssch12_APR[14,:] 
     488tb12_APR=ssch12_APR[13,:] 
     489tup12_APR=ssch12_APR[16,:] 
     490tdn12_APR=ssch12_APR[15,:] 
     491trans12_APR=ssch12_APR[17,:] 
     492orog12_APR=ssch12_APR[11,:] 
     493 
     494 
     495imo = 4 # JUL  
     496fichier = open(f1 + str(month[imo]) + f3, 'r') 
     497ssmis = np.zeros([18, numlines[imo]], float) 
     498for iligne in range (0,numlines[imo]): 
     499    line = fichier.readline() 
     500    liste = line.split() 
     501    for j in range(0,18): 
     502        ssmis[j,iligne] = float(liste[j]) 
     503 
     504 
     505    fichier.close 
     506 
     507ssch12_JUL=ssmis 
     508lon12_JUL=ssch12_JUL[0,:] 
     509lat12_JUL=ssch12_JUL[1,:] 
     510jjr12_JUL=ssch12_JUL[4,:] 
     511ts12_JUL=ssch12_JUL[8,:] 
     512emis12_JUL=ssch12_JUL[14,:] 
     513tb12_JUL=ssch12_JUL[13,:] 
     514tup12_JUL=ssch12_JUL[16,:] 
     515tdn12_JUL=ssch12_JUL[15,:] 
     516trans12_JUL=ssch12_JUL[17,:] 
     517orog12_JUL=ssch12_JUL[11,:] 
     518 
     519 
     520########## 
     521## ch13 ## 
     522########## 
     523 
     524f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH13_ANTARC_' 
     525f3 = '2010.DAT' 
     526month = np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     527numlines = np.zeros([len(month)],int) 
     528 
     529for imo in range (0, len(month)): 
     530     print month[imo]  
     531     f = f1 + str(month[imo]) + f3 
     532     fichier = open(f, 'r') 
     533     numlines[imo] = 0 
     534     for line in fichier: numlines[imo] += 1 
     535 
     536      fichier.close() 
     537 
     538 
     539imo = 0 # FEB  
     540fichier = open(f1 + str(month[imo]) + f3, 'r') 
     541ssmis = np.zeros([18, numlines[imo]], float) 
     542for iligne in range (0,numlines[imo]): 
     543    line = fichier.readline() 
     544    liste = line.split() 
     545    for j in range(0,18): 
     546        ssmis[j,iligne] = float(liste[j]) 
     547 
     548 
     549    fichier.close 
     550 
     551 
     552ssch13_FEB=ssmis 
     553lon13_FEB=ssch13_FEB[0,:] 
     554lat13_FEB=ssch13_FEB[1,:] 
     555jjr13_FEB=ssch13_FEB[4,:] 
     556ts13_FEB=ssch13_FEB[8,:] 
     557emis13_FEB=ssch13_FEB[14,:] 
     558tb13_FEB=ssch13_FEB[13,:] 
     559tup13_FEB=ssch13_FEB[16,:] 
     560tdn13_FEB=ssch13_FEB[15,:] 
     561trans13_FEB=ssch13_FEB[17,:] 
     562orog13_FEB=ssch13_FEB[11,:] 
     563 
     564imo = 1 # APR  
     565fichier = open(f1 + str(month[imo]) + f3, 'r') 
     566ssmis = np.zeros([18, numlines[imo]], float) 
     567for iligne in range (0,numlines[imo]): 
     568    line = fichier.readline() 
     569    liste = line.split() 
     570    for j in range(0,18): 
     571        ssmis[j,iligne] = float(liste[j]) 
     572 
     573 
     574    fichier.close 
     575 
     576ssch13_APR=ssmis 
     577lon13_APR=ssch13_APR[0,:] 
     578lat13_APR=ssch13_APR[1,:] 
     579jjr13_APR=ssch13_APR[4,:] 
     580ts13_APR=ssch13_APR[8,:] 
     581emis13_APR=ssch13_APR[14,:] 
     582tb13_APR=ssch13_APR[13,:] 
     583tup13_APR=ssch13_APR[16,:] 
     584tdn13_APR=ssch13_APR[15,:] 
     585trans13_APR=ssch13_APR[17,:] 
     586orog13_APR=ssch13_APR[11,:] 
     587 
     588 
     589imo = 4 # JUL  
     590fichier = open(f1 + str(month[imo]) + f3, 'r') 
     591ssmis = np.zeros([18, numlines[imo]], float) 
     592for iligne in range (0,numlines[imo]): 
     593    line = fichier.readline() 
     594    liste = line.split() 
     595    for j in range(0,18): 
     596        ssmis[j,iligne] = float(liste[j]) 
     597 
     598 
     599    fichier.close 
     600 
     601ssch13_JUL=ssmis 
     602lon13_JUL=ssch13_JUL[0,:] 
     603lat13_JUL=ssch13_JUL[1,:] 
     604jjr13_JUL=ssch13_JUL[4,:] 
     605ts13_JUL=ssch13_JUL[8,:] 
     606emis13_JUL=ssch13_JUL[14,:] 
     607tb13_JUL=ssch13_JUL[13,:] 
     608tup13_JUL=ssch13_JUL[16,:] 
     609tdn13_JUL=ssch13_JUL[15,:] 
     610trans13_JUL=ssch13_JUL[17,:] 
     611orog13_JUL=ssch13_JUL[11,:] 
Note: See TracChangeset for help on using the changeset viewer.