Changeset 21

Timestamp:

05/27/14 18:52:15 (10 years ago)

Author:

lahlod

Message:

modifs laura

Location:

trunk/src

Files:

• output_Laura (added)
• scripts_Laura/read_SSMIS_test.py (modified) (15 diffs)

trunk/src/scripts_Laura/read_SSMIS_test.py

@@ -7 +7 @@
 from mpl_toolkits.basemap import Basemap
 from mpl_toolkits.basemap import shiftgrid, cm
-import netCDF4
+from netCDF4 import Dataset
 import ffgrid2
 
-
-################ fichiers par mois - CH1 ###################################################
 
 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH'
 f3 = '_ANTARC_JUNE2010.DAT'
-#date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
 channel = np.array([12, 13, 15, 16, 17, 18])
 numlines = np.zeros([len(channel)],int)
@@ -26 +23 @@
      for line in fichier: numlines[ich] += 1
 
-
-     fichier.close()
-
-ich = 2 # 37GHz, H polar
+      fichier.close()
+
+
+################ fichiers par canal - mois de juin ###################################################
+
+
+ich = 0 # 37GHz, H 
 fichier = open(f1 + str(channel[ich]) + f3, 'r')
 ssmis = np.zeros([18, numlines[ich]], float)
@@ -159 +159 @@
 
 
-ich = 5 # 91.66GHz, V polar
+ich = 5 # 91.66GHz, H polar
 fichier = open(f1 + str(channel[ich]) + f3, 'r')
 ssmis = np.zeros([18, numlines[ich]], float)
@@ -184 +184 @@
 orog18_JUN=ssch18_JUN[11,:]
 
+
+
+################ fichiers par mois - canal 17 ###################################################
+
+f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH17_ANTARC_'
+f3 = '2010.DAT'
+date=np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY'])
+numlines = np.zeros([len(date)],int)
+
+for imo in range (0, len(date)):
+     print date[imo] 
+     f = f1 + str(date[imo]) + f3
+     fichier = open(f, 'r')
+     numlines[imo] = 0
+     for line in fichier: numlines[imo] += 1
+
+      fichier.close()
+
+
+imo = 0 # FEBUARY
+fichier = open(f1 + str(date[imo]) + f3, 'r')
+ssmis = np.zeros([18, numlines[imo]], float)
+for iligne in range (0,numlines[imo]):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+    fichier.close
+
+ssch17_FEB=ssmis
+lon17_FEB=ssch17_FEB[0,:]
+lat17_FEB=ssch17_FEB[1,:]
+jjr17_FEB=ssch17_FEB[4,:]
+ts17_FEB=ssch17_FEB[8,:]
+emis17_FEB=ssch17_FEB[14,:]
+tb17_FEB=ssch17_FEB[13,:]
+tup17_FEB=ssch17_FEB[16,:]
+tdn17_FEB=ssch17_FEB[15,:]
+trans17_FEB=ssch17_FEB[17,:]
+orog17_FEB=ssch17_FEB[11,:]
+
+
+imo = 1 # APRIL
+fichier = open(f1 + str(date[imo]) + f3, 'r')
+ssmis = np.zeros([18, numlines[imo]], float)
+for iligne in range (0,numlines[imo]):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+    fichier.close
+
+ssch17_APR=ssmis
+lon17_APR=ssch17_APR[0,:]
+lat17_APR=ssch17_APR[1,:]
+jjr17_APR=ssch17_APR[4,:]
+ts17_APR=ssch17_APR[8,:]
+emis17_APR=ssch17_APR[14,:]
+tb17_APR=ssch17_APR[13,:]
+tup17_APR=ssch17_APR[16,:]
+tdn17_APR=ssch17_APR[15,:]
+trans17_APR=ssch17_APR[17,:]
+orog17_APR=ssch17_APR[11,:]
+
+
+imo = 2 # MAY
+fichier = open(f1 + str(date[imo]) + f3, 'r')
+ssmis = np.zeros([18, numlines[imo]], float)
+for iligne in range (0,numlines[imo]):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+    fichier.close
+
+ssch17_MAY=ssmis
+lon17_MAY=ssch17_MAY[0,:]
+lat17_MAY=ssch17_MAY[1,:]
+jjr17_MAY=ssch17_MAY[4,:]
+ts17_MAY=ssch17_MAY[8,:]
+emis17_MAY=ssch17_MAY[14,:]
+tb17_MAY=ssch17_MAY[13,:]
+tup17_MAY=ssch17_MAY[16,:]
+tdn17_MAY=ssch17_MAY[15,:]
+trans17_MAY=ssch17_MAY[17,:]
+orog17_MAY=ssch17_MAY[11,:]
+
+
+imo = 3 # JUNE
+fichier = open(f1 + str(date[imo]) + f3, 'r')
+ssmis = np.zeros([18, numlines[imo]], float)
+for iligne in range (0,numlines[imo]):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+    fichier.close
+
+ssch17_JUN=ssmis
+lon17_JUN=ssch17_JUN[0,:]
+lat17_JUN=ssch17_JUN[1,:]
+jjr17_JUN=ssch17_JUN[4,:]
+ts17_JUN=ssch17_JUN[8,:]
+emis17_JUN=ssch17_JUN[14,:]
+tb17_JUN=ssch17_JUN[13,:]
+tup17_JUN=ssch17_JUN[16,:]
+tdn17_JUN=ssch17_JUN[15,:]
+trans17_JUN=ssch17_JUN[17,:]
+orog17_JUN=ssch17_JUN[11,:]
+
+
+imo = 4 # JULY
+fichier = open(f1 + str(date[imo]) + f3, 'r')
+ssmis = np.zeros([18, numlines[imo]], float)
+for iligne in range (0,numlines[imo]):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+    fichier.close
+
+ssch17_JUL=ssmis
+lon17_JUL=ssch17_JUL[0,:]
+lat17_JUL=ssch17_JUL[1,:]
+jjr17_JUL=ssch17_JUL[4,:]
+ts17_JUL=ssch17_JUL[8,:]
+emis17_JUL=ssch17_JUL[14,:]
+tb17_JUL=ssch17_JUL[13,:]
+tup17_JUL=ssch17_JUL[16,:]
+tdn17_JUL=ssch17_JUL[15,:]
+trans17_JUL=ssch17_JUL[17,:]
+orog17_JUL=ssch17_JUL[11,:]
+
+
 ########################
 # EVOLUTION TEMPORELLE #
 ########################
-## Autour de "Dome C" (lon=123.23;lat=-75.06) - mois: JUNE ##
-## Autre zone de glace de mer / de continent (lon entre -40 et -60 // lat entre -75 et -85)## 
+
+## ZONE1 = Autour de "Dome C" (lon=123.23;lat=-75.06) - mois: JUNE ##
 ## ch17 ##
-#bbzone_CH17_JUN = nonzero((lon17_JUN>=120.23)&(lon17_JUN<=126.23)&(lat17_JUN>=-78.06)&(lat17_JUN<=-72.06))
+#bbzone_CH17_JUN = nonzero((lon17_JUN >= 120.23) & (lon17_JUN <= 126.23) & (lat17_JUN >= -78.06) & (lat17_JUN <= -72.06))
+
+## ch18 ##
+#bbzone_CH18_JUN = nonzero((lon18_JUN >= 120.23) & (lon18_JUN <= 126.23) & (lat18_JUN >= -78.06) & (lat18_JUN <= -72.06))
+
+
+## ZONE2 = Autre zone de glace de mer / de continent (lon entre -40 et -60 // lat entre -75 et -85) ##
+## ch17 ##
+bbzone_CH17_FEB = nonzero((lon17_FEB >= -60.) & (lon17_FEB <= -40.) & (lat17_FEB >= -85.) & (lat17_FEB <= -75.))
+bbzone_CH17_APR = nonzero((lon17_APR >= -60.) & (lon17_APR <= -40.) & (lat17_APR >= -85.) & (lat17_APR <= -75.))
+bbzone_CH17_MAY = nonzero((lon17_MAY >= -60.) & (lon17_MAY <= -40.) & (lat17_MAY >= -85.) & (lat17_MAY <= -75.))
 bbzone_CH17_JUN = nonzero((lon17_JUN >= -60.) & (lon17_JUN <= -40.) & (lat17_JUN >= -85.) & (lat17_JUN <= -75.))
+bbzone_CH17_JUL = nonzero((lon17_JUL >= -60.) & (lon17_JUL <= -40.) & (lat17_JUL >= -85.) & (lat17_JUL <= -75.))
+
+## ch18 ##
+#bbzone_CH18_JUN = nonzero((lon18_JUN >= -60.) & (lon18_JUN <= -40.) & (lat18_JUN >= -85.) & (lat18_JUN <= -75.))
+
+## ch 16 ##
+#bbzone_CH16_JUN = nonzero((lon16_JUN >= 120.23) & (lon16_JUN <= 126.23) & (lat16_JUN >= -78.06) & (lat16_JUN <= -72.06))
+
+
+## FEBUARY ##
+emis17_FEB_moy = np.zeros([28],float)
+tb17_FEB_moy = np.zeros([28],float)
+ts17_FEB_moy = np.zeros([28],float)
+orog17_FEB_moy = np.zeros([28],float)
+for ijr in range (0,28):
+     print 'jour=', ijr+1
+     ## ch17 ##
+     ind_jr17_FEB = np.where(jjr17_FEB[bbzone_CH17_FEB]==ijr+1)[0]
+     a = emis17_FEB[bbzone_CH17_FEB][ind_jr17_FEB]
+     b = tb17_FEB[bbzone_CH17_FEB][ind_jr17_FEB]
+     c = ts17_FEB[bbzone_CH17_FEB][ind_jr17_FEB]
+     d = orog17_FEB[bbzone_CH17_FEB][ind_jr17_FEB]
+     emis17_FEB_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))])
+     tb17_FEB_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))])
+     ts17_FEB_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))])
+     orog17_FEB_moy[ijr] = mean(d[nonzero((d!=-500.))])
+
+
+## APRIL ##
+emis17_APR_moy = np.zeros([30],float)
+tb17_APR_moy = np.zeros([30],float)
+ts17_APR_moy = np.zeros([30],float)
+orog17_APR_moy = np.zeros([30],float)
+for ijr in range (0,30):
+     print 'jour=', ijr+1
+     ind_jr17_APR = np.where(jjr17_APR[bbzone_CH17_APR]==ijr+1)[0]
+     a = emis17_APR[bbzone_CH17_APR][ind_jr17_APR]
+     b = tb17_APR[bbzone_CH17_APR][ind_jr17_APR]
+     c = ts17_APR[bbzone_CH17_APR][ind_jr17_APR]
+     d = orog17_APR[bbzone_CH17_APR][ind_jr17_APR]
+     emis17_APR_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))])
+     tb17_APR_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))])
+     ts17_APR_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))])
+     orog17_APR_moy[ijr] = mean(d[nonzero((d!=-500.))])
+
+
+## MAY ##
+emis17_MAY_moy = np.zeros([31],float)
+tb17_MAY_moy = np.zeros([31],float)
+ts17_MAY_moy = np.zeros([31],float)
+orog17_MAY_moy = np.zeros([31],float)
+for ijr in range (0,31):
+     print 'jour=', ijr+1
+     ind_jr17_MAY = np.where(jjr17_MAY[bbzone_CH17_MAY]==ijr+1)[0]
+     a = emis17_MAY[bbzone_CH17_MAY][ind_jr17_MAY]
+     b = tb17_MAY[bbzone_CH17_MAY][ind_jr17_MAY]
+     c = ts17_MAY[bbzone_CH17_MAY][ind_jr17_MAY]
+     d = orog17_MAY[bbzone_CH17_MAY][ind_jr17_MAY]
+     emis17_MAY_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))])
+     tb17_MAY_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))])
+     ts17_MAY_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))])
+     orog17_MAY_moy[ijr] = mean(d[nonzero((d!=-500.))])
+
+
+## JUNE ##
 emis17_JUN_moy = np.zeros([30],float)
 tb17_JUN_moy = np.zeros([30],float)
 ts17_JUN_moy = np.zeros([30],float)
-
-## ch18 ##
-#bbzone_CH18_JUN = nonzero((lon18_JUN >= 120.23) & (lon18_JUN <= 126.23) & (lat18_JUN >= -78.06) & (lat18_JUN <= -72.06))
-bbzone_CH18_JUN = nonzero((lon18_JUN >= -60.) & (lon18_JUN <= -40.) & (lat18_JUN >= -85.) & (lat18_JUN <= -75.))
-emis18_JUN_moy = np.zeros([30],float)
-tb18_JUN_moy = np.zeros([30],float)
-ts18_JUN_moy = np.zeros([30],float)
-
-## ch 16 ##
-#bbzone_CH16_JUN = nonzero((lon16_JUN >= 120.23) & (lon16_JUN <= 126.23) & (lat16_JUN >= -78.06) & (lat16_JUN <= -72.06))
+#orog17_JUN_moy = np.zeros([30],float)
+#emis18_JUN_moy = np.zeros([30],float)
+#tb18_JUN_moy = np.zeros([30],float)
+#ts18_JUN_moy = np.zeros([30],float)
 #emis16_JUN_moy = np.zeros([30],float)
 #tb16_JUN_moy = np.zeros([30],float)
 #ts16_JUN_moy = np.zeros([30],float)
-
 for ijr in range (0,30):
      print 'jour=', ijr+1
-     ## ch17 ##
-     ind_jr17 = np.where(jjr17_JUN[bbzone_CH17_JUN]==ijr+1)[0]
-     a = emis17_JUN[bbzone_CH17_JUN][ind_jr17]
-     b = tb17_JUN[bbzone_CH17_JUN][ind_jr17]
-     c = ts17_JUN[bbzone_CH17_JUN][ind_jr17]
+     ind_jr17_JUN = np.where(jjr17_JUN[bbzone_CH17_JUN]==ijr+1)[0]
+     a = emis17_JUN[bbzone_CH17_JUN][ind_jr17_JUN]
+     b = tb17_JUN[bbzone_CH17_JUN][ind_jr17_JUN]
+     c = ts17_JUN[bbzone_CH17_JUN][ind_jr17_JUN]
+     d = orog17_JUN[bbzone_CH17_JUN][ind_jr17_JUN]
      emis17_JUN_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))])
      tb17_JUN_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))])
      ts17_JUN_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))])
+#     orog17_JUN_moy[ijr] = mean(d[nonzero((d!=-500.))])
      ## ch18 ##
-     ind_jr18 = np.where(jjr18_JUN[bbzone_CH18_JUN]==ijr+1)[0]
-     d = emis18_JUN[bbzone_CH18_JUN][ind_jr18]
-     e = tb18_JUN[bbzone_CH18_JUN][ind_jr18]
-     f = ts18_JUN[bbzone_CH18_JUN][ind_jr18]
-     emis18_JUN_moy[ijr] = mean(d[nonzero((d != -500.) & (d <= 1.))])
-     tb18_JUN_moy[ijr] = mean(e[nonzero((e != -500.) & (e != 0.))])
-     ts18_JUN_moy[ijr] = mean(f[nonzero((f != -500.) & (f != 0.))])
+#     ind_jr18 = np.where(jjr18_JUN[bbzone_CH18_JUN]==ijr+1)[0]
+#     d = emis18_JUN[bbzone_CH18_JUN][ind_jr18]
+#     e = tb18_JUN[bbzone_CH18_JUN][ind_jr18]
+#     f = ts18_JUN[bbzone_CH18_JUN][ind_jr18]
+#     emis18_JUN_moy[ijr] = mean(d[nonzero((d != -500.) & (d <= 1.))])
+#     tb18_JUN_moy[ijr] = mean(e[nonzero((e != -500.) & (e != 0.))])
+#     ts18_JUN_moy[ijr] = mean(f[nonzero((f != -500.) & (f != 0.))])
      ## ch16 ##
 #     ind_jr16 = np.where(jjr16_JUN[bbzone_CH16_JUN]==ijr+1)[0]
@@ -237 +445 @@
 
 
+## JULY ##
+emis17_JUL_moy = np.zeros([31],float)
+tb17_JUL_moy = np.zeros([31],float)
+ts17_JUL_moy = np.zeros([31],float)
+orog17_JUL_moy = np.zeros([31],float)
+for ijr in range (0,31):
+     print 'jour=', ijr+1
+     ind_jr17_JUL = np.where(jjr17_JUL[bbzone_CH17_JUL]==ijr+1)[0]
+     a = emis17_JUL[bbzone_CH17_JUL][ind_jr17_JUL]
+     b = tb17_JUL[bbzone_CH17_JUL][ind_jr17_JUL]
+     c = ts17_JUL[bbzone_CH17_JUL][ind_jr17_JUL]
+     d = orog17_JUL[bbzone_CH17_JUL][ind_jr17_JUL]
+     emis17_JUL_moy[ijr] = mean(a[nonzero((a!=-500.)&(a<=1.))])
+     tb17_JUL_moy[ijr] = mean(b[nonzero((b!=-500.)&(b!=0.))])
+     ts17_JUL_moy[ijr] = mean(c[nonzero((c!=-500.)&(c!=0.))])
+     orog17_JUL_moy[ijr] = mean(d[nonzero((d!=-500.))])
+
+
+## plot evolution temporelle ##
+## FEBUARY ##
+fig = plt.figure()
+plt.subplot(3, 1, 1)
+plt.title('FEBUARY')
+plt.plot(arange(0, 28, 1), emis17_FEB_moy, c = 'b')
+plt.ylabel('Emissivity')
+xticks(arange(0, 28, 1), arange(1, 29, 1))
+xlim (0, 28)
+grid(True)
+plt.subplot(3, 1, 2)
+plt.plot(arange(0, 28, 1), tb17_FEB_moy, c = 'g', label = 'Tb')
+ylabel('Brightness temperature')
+xticks(arange(0, 28, 1), arange(1, 29, 1))
+xlim (0, 28)
+grid(True)
+plt.subplot(3, 1, 3)
+plt.plot(arange(0, 28, 1), ts17_FEB_moy, c = 'r', label = 'Ts')
+plt.ylabel('Surface temperature')
+xticks(arange(0, 28, 1), arange(1, 29, 1))
+xlim (0, 28)
+grid(True)
+fig.show()
+
+## from APRIL to JULY ##
+emis1 = np.append(emis17_APR_moy, emis17_MAY_moy)
+emis2 = np.append(emis1, emis17_JUN_moy)
+emis17_moy = np.append(emis2, emis17_JUL_moy)
+tb1 = np.append(tb17_APR_moy, tb17_MAY_moy)
+tb2 = np.append(tb1, tb17_JUN_moy)
+tb17_moy = np.append(tb2, tb17_JUL_moy)
+ts1 = np.append(ts17_APR_moy, ts17_MAY_moy)
+ts2 = np.append(ts1, ts17_JUN_moy)
+ts17_moy = np.append(ts2, ts17_JUL_moy)
+fig = plt.figure()
+plt.subplot(3, 1, 1)
+plt.plot(arange(0, 122, 1), emis17_moy, c = 'b')
+plt.ylabel('Emissivity')
+xticks(array([0, 30, 61, 91, 122]), date[1:])
+xlim (0, 122)
+grid(True)
+plt.subplot(3, 1, 2)
+plt.plot(arange(0, 122, 1), tb17_moy, c = 'g', label = 'Tb')
+ylabel('Brightness temperature')
+xticks(array([0, 30, 61, 91, 122]), date[1:])
+xlim (0, 122)
+grid(True)
+plt.subplot(3, 1, 3)
+plt.plot(arange(0, 122, 1), ts17_moy, c = 'r', label = 'Ts')
+plt.ylabel('Surface temperature')
+xticks(array([0, 30, 61, 91, 122]), date[1:])
+xlim (0, 122)
+grid(True)
+fig.show()
+
+## JUNE ##
 ## ch17 - ch18 ##
 fig=plt.figure()
@@ -262 +544 @@
 legend(loc=7)
 fig.show()
-
-## ch16 ##
+## ch16 - ch17 ##
 fig=plt.figure()
 plt.subplot(3,1,1)
@@ -290 +571 @@
 
 
-
-tb17_JUN_anom = np.zeros([30],float)
-tb18_JUN_anom = np.zeros([30],float)
+## calcul anomalie de Tb
+## FEBUARY ##
+tb17_FEB_anom = np.zeros([28], float)
+for ijr in range (0,28):
+     print 'jour=', ijr + 1
+     tb17_FEB_anom[ijr]=tb17_FEB_moy[ijr]-mean(tb17_FEB_moy)
+
+## APRIL ##
+tb17_APR_anom = np.zeros([30], float)
 for ijr in range (0,30):
      print 'jour=', ijr + 1
-     tb17_JUN_anom[ijr]=tb17_JUN_moy[ijr]-tb17_JUN_moy.mean()
-     tb18_JUN_anom[ijr]=tb18_JUN_moy[ijr]-tb18_JUN_moy.mean()
-
-
+     tb17_APR_anom[ijr]=tb17_APR_moy[ijr]-mean(tb17_APR_moy)
+
+## MAY ##
+bbnan = nonzero(isnan(tb17_MAY_moy) == False)
+tb17_MAY_anom = np.zeros([31], float)
+for ijr in range (0,31):
+     print 'jour=', ijr + 1
+     tb17_MAY_anom[ijr]=tb17_MAY_moy[ijr]-mean(tb17_MAY_moy[bbnan])
+
+## JUNE ##
+tb17_JUN_anom = np.zeros([30],float)
+#tb18_JUN_anom = np.zeros([30],float)
+for ijr in range (0,30):
+     print 'jour=', ijr + 1
+     tb17_JUN_anom[ijr]=tb17_JUN_moy[ijr]-mean(tb17_JUN_moy)
+#     tb18_JUN_anom[ijr]=tb18_JUN_moy[ijr]-mean(tb18_JUN_moy)
+
+## JULY ##
+tb17_JUL_anom = np.zeros([31],float)
+for ijr in range (0,30):
+     print 'jour=', ijr + 1
+     tb17_JUL_anom[ijr]=tb17_JUL_moy[ijr]-mean(tb17_JUL_moy)
+
+
+
+a1 = np.append(tb17_APR_anom, tb17_MAY_anom)
+a2 = np.append(a1, tb17_JUN_anom)
+tb17_month_anom = np.append(a2, tb17_JUL_anom)
+
+
+## for APRIL to JULY ##
+fig = plt.figure()
+plt.plot(arange(0, 122, 1), tb17_month_anom, c = 'k')
+plt.plot(arange(0, 122, 1), np.zeros([122]), '--', c = 'r')
+ylabel('Tb anomaly (K)')
+xlim(0, 122)
+xticks(array([0, 30, 61, 91 122]), date[1:])
+grid(True)
+plt.title('SSMIS - zone2')
+fig.show()
+
+## FEBRUARY ##
+fig = plt.figure()
+plt.plot(arange(0, 28, 1), tb17_FEB_anom, c = 'k')
+plt.plot(arange(0, 28, 1), np.zeros([28]), '--', c = 'r')
+ylabel('Tb anomaly (K)')
+xlim(0, 27)
+xticks(arange(0, 28, 1), arange(1, 29, 1))
+grid(True)
+plt.xlabel('FABRUARY')
+plt.title('SSMIS - zone2')
+fig.show()
+
+## JUNE ##
 fig=plt.figure()
 plt.plot(arange(0,30,1),tb17_JUN_anom,label='91.66Ghz-V',c='r')
@@ -324 +661 @@
 y0, y1 = -90, -30
 
+
+## FEBRUARY ##
+bbtb_ch17_FEB = nonzero((tb17_FEB != -500.) & (tb17_FEB != 0.))
+OUTZCH17_FEB = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),28], float)
+outzch17_FEB = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float)
+lonch17_FEB = np.zeros([len(np.arange(x0, x1+1, dx))], float)
+latch17_FEB = np.zeros([len(np.arange(y0, y1+1, dy))], float)
+for ijr in range(0,28):
+    print 'jour=', ijr+1
+    ## ch17 ##
+    ind_jr17_FEB = np.where(jjr17_FEB[bbtb_ch17_FEB] == ijr+1)[0]
+    xx = lon17_FEB[bbtb_ch17_FEB][ind_jr17_FEB]
+    yy = lat17_FEB[bbtb_ch17_FEB][ind_jr17_FEB]
+    zz = tb17_FEB[bbtb_ch17_FEB][ind_jr17_FEB]
+    zz0 = min(zz)
+    zz1 = max(zz)
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch17_FEB = outz
+    lonch17_FEB = outx
+    latch17_FEB = outy
+    OUTZCH17_FEB[:,:,ijr] = outzch17_FEB[:,:]
+
+## APRIL ##
+bbtb_ch17_APR = nonzero((tb17_APR != -500.) & (tb17_APR != 0.))
+OUTZCH17_APR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float)
+outzch17_APR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float)
+lonch17_APR = np.zeros([len(np.arange(x0, x1+1, dx))], float)
+latch17_APR = np.zeros([len(np.arange(y0, y1+1, dy))], float)
+for ijr in range(0,30):
+    print 'jour=', ijr+1
+    ## ch17 ##
+    ind_jr17_APR = np.where(jjr17_APR[bbtb_ch17_APR] == ijr+1)[0]
+    xx = lon17_APR[bbtb_ch17_APR][ind_jr17_APR]
+    yy = lat17_APR[bbtb_ch17_APR][ind_jr17_APR]
+    zz = tb17_APR[bbtb_ch17_APR][ind_jr17_APR]
+    zz0 = min(zz)
+    zz1 = max(zz)
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch17_APR = outz
+    lonch17_APR = outx
+    latch17_APR = outy
+    OUTZCH17_APR[:,:,ijr] = outzch17_APR[:,:]
+
+## MAY ##
+bbtb_ch17_MAY = nonzero((tb17_MAY != -500.) & (tb17_MAY != 0.))
+OUTZCH17_MAY = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),31], float)
+outzch17_MAY = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float)
+lonch17_MAY = np.zeros([len(np.arange(x0, x1+1, dx))], float)
+latch17_MAY = np.zeros([len(np.arange(y0, y1+1, dy))], float)
+for ijr in range(0,31):
+    print 'jour=', ijr+1
+    ## ch17 ##
+    ind_jr17_MAY = np.where(jjr17_MAY[bbtb_ch17_MAY] == ijr+1)[0]
+    xx = lon17_MAY[bbtb_ch17_MAY][ind_jr17_MAY]
+    yy = lat17_MAY[bbtb_ch17_MAY][ind_jr17_MAY]
+    zz = tb17_MAY[bbtb_ch17_MAY][ind_jr17_MAY]
+    zz0 = min(zz)
+    zz1 = max(zz)
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch17_MAY = outz
+    lonch17_MAY = outx
+    latch17_MAY = outy
+    OUTZCH17_MAY[:,:,ijr] = outzch17_MAY[:,:]
+
+## JUNE ##
 ## ch17 ##
-bbemis_ch17_JUN=nonzero((emis17_JUN!=-500.)&(emis17_JUN<1.)&(emis17_JUN>0.))
-bbtb_ch17_JUN=nonzero((tb17_JUN!=-500.)&(tb17_JUN!=0.))
+bbtb_ch17_JUN = nonzero((tb17_JUN != -500.) & (tb17_JUN != 0.))
 OUTZCH17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float)
 outzch17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float)
-lonch17_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float)
-latch17_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float)
-
+lonch17_JUN = np.zeros([len(np.arange(x0, x1+1, dx))], float)
+latch17_JUN = np.zeros([len(np.arange(y0, y1+1, dy))], float)
 ## ch18 ##
 bbemis_ch18_JUN=nonzero((emis18_JUN!=-500.)&(emis18_JUN<1.)&(emis18_JUN>0.))
@@ -339 +739 @@
 lonch18_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float)
 latch18_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float)
-
 for ijr in range(0,30):
     print 'jour=', ijr+1
@@ -367 +766 @@
     OUTZCH18_JUN[:,:,ijr] = outzch18_JUN[:,:]
 
-    
+## JULY ##
+bbtb_ch17_JUL = nonzero((tb17_JUL != -500.) & (tb17_JUL != 0.))
+OUTZCH17_JUL = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),31], float)
+outzch17_JUL = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float)
+lonch17_JUL = np.zeros([len(np.arange(x0, x1+1, dx))], float)
+latch17_JUL = np.zeros([len(np.arange(y0, y1+1, dy))], float)
+for ijr in range(0,31):
+    print 'jour=', ijr+1
+    ## ch17 ##
+    ind_jr17_JUL = np.where(jjr17_JUL[bbtb_ch17_JUL] == ijr+1)[0]
+    xx = lon17_JUL[bbtb_ch17_JUL][ind_jr17_JUL]
+    yy = lat17_JUL[bbtb_ch17_JUL][ind_jr17_JUL]
+    zz = tb17_JUL[bbtb_ch17_JUL][ind_jr17_JUL]
+    zz0 = min(zz)
+    zz1 = max(zz)
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch17_JUL = outz
+    lonch17_JUL = outx
+    latch17_JUL = outy
+    OUTZCH17_JUL[:,:,ijr] = outzch17_JUL[:,:]
+
+
+## calcul de la climatologie moyenne sur le mois et anomalie en chaque lon/lat ##
+## FEBRUARY ##
+mean_outzch17_FEB = np.zeros([len(latch17_FEB), len(lonch17_FEB)], float)
+for ilon in range (0,len(lonch17_FEB)):
+    for ilat in range (0,len(latch17_FEB)):
+        mean_outzch17_FEB[ilat,ilon] = mean(OUTZCH17_FEB[ilat,ilon,:])
+
+tbch17_anom_FEB = np.zeros([len(latch17_FEB),len(lonch17_FEB),28], float)
+for ijr in range (0,28):
+    tbch17_anom_FEB[:,:,ijr] = OUTZCH17_FEB[:,:,ijr] - mean_outzch17_FEB[:,:]
+
+## APRIL ##
+mean_outzch17_APR = np.zeros([len(latch17_APR), len(lonch17_APR)], float)
+for ilon in range (0,len(lonch17_APR)):
+    for ilat in range (0,len(latch17_APR)):
+        mean_outzch17_APR[ilat,ilon] = mean(OUTZCH17_APR[ilat,ilon,:])
+
+tbch17_anom_APR = np.zeros([len(latch17_APR),len(lonch17_APR),30], float)
+for ijr in range (0,30):
+    tbch17_anom_APR[:,:,ijr] = OUTZCH17_APR[:,:,ijr] - mean_outzch17_APR[:,:]
+
+## MAY ##
+mean_outzch17_MAY = np.zeros([len(latch17_MAY), len(lonch17_MAY)], float)
+for ilon in range (0,len(lonch17_MAY)):
+    for ilat in range (0,len(latch17_MAY)):
+        mean_outzch17_MAY[ilat,ilon] = mean(OUTZCH17_MAY[ilat,ilon,:])
+
+tbch17_anom_MAY = np.zeros([len(latch17_MAY),len(lonch17_MAY),31], float)
+for ijr in range (0,31):
+    tbch17_anom_MAY[:,:,ijr] = OUTZCH17_MAY[:,:,ijr] - mean_outzch17_MAY[:,:]
+
+
+## JUNE ##
 ## ch17 ##
 mean_outzch17_JUN = np.zeros([len(latch17_JUN), len(lonch17_JUN)], float)
@@ -374 +827 @@
         mean_outzch17_JUN[ilat,ilon] = OUTZCH17_JUN[ilat,ilon,:].mean()
 
+tbch17_anom_JUN = np.zeros([len(latch17_JUN),len(lonch17_JUN),30], float)
+for ijr in range (0,30):
+    tbch17_anom_JUN[:,:,ijr] = OUTZCH17_JUN[:,:,ijr] - mean_outzch17_JUN[:,:]
 
 ## ch18 ##
@@ -381 +837 @@
         mean_outzch18_JUN[ilat,ilon] = OUTZCH18_JUN[ilat,ilon,:].mean()
 
-
-## ch17 ##
-tbch17_anom_JUN = np.zeros([len(latch17_JUN),len(lonch17_JUN),30], float)
-for ijr in range (0,30):
-    tbch17_anom_JUN[:,:,ijr] = OUTZCH17_JUN[:,:,ijr] - mean_outzch17_JUN[:,:]
-
-
-## ch18 ##
 tbch18_anom_JUN = np.zeros([len(latch18_JUN),len(lonch18_JUN),30], float)
 for ijr in range (0,30):
     tbch18_anom_JUN[:,:,ijr] = OUTZCH18_JUN[:,:,ijr] - mean_outzch18_JUN[:,:]
+
 
 
@@ -416 +865 @@
 m = Basemap(llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-90, urcrnrlat=-30, projection='cyl', resolution='c', fix_aspect=True)
 m.drawcoastlines(linewidth = 1)
-m.drawparallels(np.arange(-90., 90., 20))
+m.drawparallels(np.arange(-90., -30., 20))
 m.drawmeridians(np.arange(-180., 180., 20))
 #m.fillcontinents()
-xii,yii = m(*np.meshgrid(lonch17_JUN, latch17_JUN))
-
+xii,yii = m(*np.meshgrid(lonch17_FEB, latch17_FEB))
+plt.xticks(arange(-180, 200, 20))
+plt.yticks(arange(-90, -30, 20))
+## ch17 ##
+clevs = arange(150., 270., 1.)
+cs = m.contourf(xii, yii, mean_outzch17_FEB, clevs, cmap=cm.s3pcpn_l_r)
+cbar = colorbar(cs)
+cbar.set_label('Mean Tb FEB [CH17] - SSMIS')
 ## BIAIS ch17-ch18 ## 
 biais_JUN = mean_outzch17_JUN - mean_outzch18_JUN
@@ -447 +902 @@
 # DIAGRAMME DE HOVMOLLER #
 ##########################
-
 # shape(tbch17_anom_JUN) = [ilat, ilon, ijr]
-
-## tranche de latitude étudiée ##
-#bbtranche17_JUN = nonzero((latch17_JUN == -75.))
+## FEBRUARY ##
+bbtranche17_FEB = nonzero((latch17_FEB >= -85.) & (latch17_FEB <= -75))
+mean_tbch17_anom_FEB = np.zeros([len(lonch17_FEB), 28], float)
+for ilon in range (0, len(lonch17_FEB)):
+    for ijr in range (0,28):
+        mean_tbch17_anom_FEB[ilon,ijr] = mean(tbch17_anom_FEB[bbtranche17_FEB][:,ilon,ijr])
+
+y_time, x_space = np.meshgrid(arange(0,28,1), lonch17_FEB)
+fig = plt.figure()
+plt.pcolor(x_space, y_time, mean_tbch17_anom_FEB, cmap=cm.s3pcpn_l_r, vmin = -10., vmax = 15.)
+plt.axis([-180., 180., 0, 28])
+cb = plt.colorbar()
+cb.set_label('Tb anomaly - SSMIS CH17')
+plt.xticks(arange(-180.,200.,40))
+plt.yticks(arange(0, 28, 1), arange (1, 29, 1))
+plt.xlabel('longitude')
+plt.ylabel('FEBRUARY 2010')
+
+## MAY ##
+bbtranche17_MAY = nonzero((latch17_MAY >= -85.) & (latch17_MAY <= -75))
+mean_tbch17_anom_MAY = np.zeros([len(lonch17_MAY), 31], float)
+for ilon in range (0, len(lonch17_MAY)):
+    for ijr in range (0,31):
+        mean_tbch17_anom_MAY[ilon,ijr] = mean(tbch17_anom_MAY[bbtranche17_MAY][:,ilon,ijr])
+
+y_time, x_space = np.meshgrid(arange(0,31,1), lonch17_MAY)
+fig = plt.figure()
+plt.pcolor(x_space, y_time, mean_tbch17_anom_MAY, cmap=cm.s3pcpn_l_r, vmin = -12., vmax = 25.)
+plt.axis([-180., 180., 0, 30])
+cb = plt.colorbar()
+cb.set_label('Tb anomaly - SSMIS CH17')
+plt.xticks(arange(-180.,200.,40))
+plt.yticks(arange(0, 30, 1), arange(1,31,1))
+plt.xlabel('longitude')
+plt.ylabel('MAY 2010')
+
+## JUNE ##
 bbtranche17_JUN = nonzero((latch17_JUN >= -85.) & (latch17_JUN <= -75))
-#bbtranche17_JUN = nonzero((latch17_JUN >= -90.) & (latch17_JUN <= -85))
 mean_tbch17_anom_JUN = np.zeros([len(lonch17_JUN), 30], float)
 for ilon in range (0, len(lonch17_JUN)):
@@ -460 +947 @@
 
 
-y_time, x_space = m(*np.meshgrid(arange(0,30,1), lonch17_JUN))
+y_time, x_space = np.meshgrid(arange(0,30,1), lonch17_APR)
 fig = plt.figure()
-plt.pcolor(x_space, y_time, mean_tbch17_anom_JUN, cmap=cm.s3pcpn_l_r, vmin = -30., vmax = 25.)
+plt.pcolor(x_space, y_time, mean_tbch17_anom_APR, cmap=cm.s3pcpn_l_r, vmin = -30., vmax = 25.)
 plt.axis([-180., 180., 0, 29])
 cb = plt.colorbar()