Changeset 20

Timestamp:

05/26/14 18:28:33 (10 years ago)

Author:

lahlod

Message:

modif_script_SSMIS

Location:

trunk/src/scripts_Laura

Files:

• read_AMSUB_test.py (modified) (4 diffs)
• read_SSMIS_test.py (modified) (10 diffs)

trunk/src/scripts_Laura/read_AMSUB_test.py

@@ -11 +11 @@
 ################ fichiers par mois - CH1 ###################################################
 
-f1='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUB_CH1_ANTARC_'
-f2='2010.DAT'
-date=np.array(['JANUARY', 'FEBRUARY', 'MARCH'])
-
-numlines = np.zeros([3],float)
+f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUB_CH1_ANTARC_'
+f2 = '2010.DAT'
+#date = np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
+date = np.array(['APRIL', 'MAY', 'JUNE', 'JULY'])
+
+numlines = np.zeros(len(date),float)
 
 for i in range (0,len(numlines)):
@@ -27 +28 @@
 
 
-i=0 #JANUARY
+i=0 # JANUARY
 f=f1+str(date[i])+f2
 fichier=open(f,'r')
@@ -111 +112 @@
 
 
+i=0 #APRIL
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch1_APR=amch
+#### def des variables
+lon1_APR=amch1_APR[0,:]
+lat1_APR=amch1_APR[1,:]
+jjr1_APR=amch1_APR[4,:]
+ts1_APR=amch1_APR[10,:]
+emis1_APR=amch1_APR[16,:]
+tb1_APR=amch1_APR[15,:]
+tup1_APR=amch1_APR[18,:]
+tdn1_APR=amch1_APR[17,:]
+trans1_APR=amch1_APR[19,:]
+orog1_APR=amch1_APR[13,:]
+
+
+i=1 #MAY
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch1_MAY=amch
+#### def des variables
+lon1_MAY=amch1_MAY[0,:]
+lat1_MAY=amch1_MAY[1,:]
+jjr1_MAY=amch1_MAY[4,:]
+ts1_MAY=amch1_MAY[10,:]
+emis1_MAY=amch1_MAY[16,:]
+tb1_MAY=amch1_MAY[15,:]
+tup1_MAY=amch1_MAY[18,:]
+tdn1_MAY=amch1_MAY[17,:]
+trans1_MAY=amch1_MAY[19,:]
+orog1_MAY=amch1_MAY[13,:]
+
+
+i=2 #JUNE
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch1_JUN=amch
+#### def des variables
+lon1_JUN=amch1_JUN[0,:]
+lat1_JUN=amch1_JUN[1,:]
+jjr1_JUN=amch1_JUN[4,:]
+ts1_JUN=amch1_JUN[10,:]
+emis1_JUN=amch1_JUN[16,:]
+tb1_JUN=amch1_JUN[15,:]
+tup1_JUN=amch1_JUN[18,:]
+tdn1_JUN=amch1_JUN[17,:]
+trans1_JUN=amch1_JUN[19,:]
+orog1_JUN=amch1_JUN[13,:]
+
+
+i=0 #JULY
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch1_JUL=amch
+#### def des variables
+lon1_JUL=amch1_JUL[0,:]
+lat1_JUL=amch1_JUL[1,:]
+jjr1_JUL=amch1_JUL[4,:]
+ts1_JUL=amch1_JUL[10,:]
+emis1_JUL=amch1_JUL[16,:]
+tb1_JUL=amch1_JUL[15,:]
+tup1_JUL=amch1_JUL[18,:]
+tdn1_JUL=amch1_JUL[17,:]
+trans1_JUL=amch1_JUL[19,:]
+orog1_JUL=amch1_JUL[13,:]
+
+
+
 ################ fichiers par mois - CH2  ###################################################
 
 f1='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUB_CH2_ANTARC_'
 f2='2010.DAT'
-date=np.array(['JANUARY', 'FEBRUARY', 'MARCH'])
-
-numlines = np.zeros([3],float)
+#date=np.array(['JANUARY', 'FEBRUARY', 'MARCH'])
+date = np.array(['APRIL', 'MAY', 'JUNE', 'JULY'])
+
+numlines = np.zeros([len(date)],float)
 
 for i in range (0,len(numlines)):
@@ -212 +327 @@
 orog2_MAR=amch2_MAR[13,:]
 
+
+i=0 #APRIL
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch2_APR=amch
+#### def des variables
+lon2_APR=amch2_APR[0,:]
+lat2_APR=amch2_APR[1,:]
+jjr2_APR=amch2_APR[4,:]
+ts2_APR=amch2_APR[10,:]
+emis2_APR=amch2_APR[16,:]
+tb2_APR=amch2_APR[15,:]
+tup2_APR=amch2_APR[18,:]
+tdn2_APR=amch2_APR[17,:]
+trans2_APR=amch2_APR[19,:]
+orog2_APR=amch2_APR[13,:]
+
+
+i=1 #MAY
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch2_MAY=amch
+#### def des variables
+lon2_MAY=amch2_MAY[0,:]
+lat2_MAY=amch2_MAY[1,:]
+jjr2_MAY=amch2_MAY[4,:]
+ts2_MAY=amch2_MAY[10,:]
+emis2_MAY=amch2_MAY[16,:]
+tb2_MAY=amch2_MAY[15,:]
+tup2_MAY=amch2_MAY[18,:]
+tdn2_MAY=amch2_MAY[17,:]
+trans2_MAY=amch2_MAY[19,:]
+orog2_MAY=amch2_MAY[13,:]
+
+
+i=2 #JUNE
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch2_JUN=amch
+#### def des variables
+lon2_JUN=amch2_JUN[0,:]
+lat2_JUN=amch2_JUN[1,:]
+jjr2_JUN=amch2_JUN[4,:]
+ts2_JUN=amch2_JUN[10,:]
+emis2_JUN=amch2_JUN[16,:]
+tb2_JUN=amch2_JUN[15,:]
+tup2_JUN=amch2_JUN[18,:]
+tdn2_JUN=amch2_JUN[17,:]
+trans2_JUN=amch2_JUN[19,:]
+orog2_JUN=amch2_JUN[13,:]
+
+
+i=0 #JULY
+f=f1+str(date[i])+f2
+fichier=open(f,'r')
+amch=np.zeros([20,numlines[i]],float)
+for iligne in range (0,numlines[i]):
+    line=fichier.readline()
+    liste = line.split()
+    for j in range(0,20):
+        amch[j,iligne] = float(liste[j])
+
+
+fichier.close()
+
+
+amch1_JUL=amch
+#### def des variables
+lon1_JUL=amch1_JUL[0,:]
+lat1_JUL=amch1_JUL[1,:]
+jjr1_JUL=amch1_JUL[4,:]
+ts1_JUL=amch1_JUL[10,:]
+emis1_JUL=amch1_JUL[16,:]
+tb1_JUL=amch1_JUL[15,:]
+tup1_JUL=amch1_JUL[18,:]
+tdn1_JUL=amch1_JUL[17,:]
+trans1_JUL=amch1_JUL[19,:]
+orog1_JUL=amch1_JUL[13,:]
+

trunk/src/scripts_Laura/read_SSMIS_test.py

@@ -16 +16 @@
 f3 = '_ANTARC_JUNE2010.DAT'
 #date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
-channel = np.array([12, 13, 15, 16])
+channel = np.array([12, 13, 15, 16, 17, 18])
 numlines = np.zeros([len(channel)],int)
 
@@ -27 +27 @@
 
 
-fichier.close()
+     fichier.close()
 
 ich = 2 # 37GHz, H polar
@@ -133 +133 @@
 
 
+ich = 4 # 91.66GHz, V polar
+fichier = open(f1 + str(channel[ich]) + f3, 'r')
+ssmis = np.zeros([18, numlines[ich]], float)
+for iligne in range (0,numlines[ich]-1):
+    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,:]
+
+
+ich = 5 # 91.66GHz, V polar
+fichier = open(f1 + str(channel[ich]) + f3, 'r')
+ssmis = np.zeros([18, numlines[ich]], float)
+for iligne in range (0,numlines[ich]-1):
+    line = fichier.readline()
+    liste = line.split()
+    for j in range(0,18):
+        ssmis[j,iligne] = float(liste[j])
+
+
+    fichier.close
+
+
+ssch18_JUN=ssmis
+lon18_JUN=ssch18_JUN[0,:]
+lat18_JUN=ssch18_JUN[1,:]
+jjr18_JUN=ssch18_JUN[4,:]
+ts18_JUN=ssch18_JUN[8,:]
+emis18_JUN=ssch18_JUN[14,:]
+tb18_JUN=ssch18_JUN[13,:]
+tup18_JUN=ssch18_JUN[16,:]
+tdn18_JUN=ssch18_JUN[15,:]
+trans18_JUN=ssch18_JUN[17,:]
+orog18_JUN=ssch18_JUN[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)## 
 ## 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))
+bbzone_CH17_JUN = nonzero((lon17_JUN >= -60.) & (lon17_JUN <= -40.) & (lat17_JUN >= -85.) & (lat17_JUN <= -75.))
 emis17_JUN_moy = np.zeros([30],float)
 tb17_JUN_moy = np.zeros([30],float)
@@ -144 +197 @@
 
 ## 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 >= 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)
@@ -150 +204 @@
 
 ## ch 16 ##
-bbzone_CH16_JUN = nonzero((lon16_JUN >= 120.23) & (lon16_JUN <= 126.23) & (lat16_JUN >= -78.06) & (lat16_JUN <= -72.06))
-emis16_JUN_moy = np.zeros([30],float)
-tb16_JUN_moy = np.zeros([30],float)
-ts16_JUN_moy = np.zeros([30],float)
+#bbzone_CH16_JUN = nonzero((lon16_JUN >= 120.23) & (lon16_JUN <= 126.23) & (lat16_JUN >= -78.06) & (lat16_JUN <= -72.06))
+#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(jjr16_JUN[bbzone_CH16_JUN]==ijr+1)[0]
-#     a = emis17_JUN[bbzone_CH16_JUN][ind_jr16]
-#     b = tb17_JUN[bbzone_CH16_JUN][ind_jr16]
-#     c = ts17_JUN[bbzone_CH16_JUN][ind_jr16]
-#     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.))])
+     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]
+     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.))])
      ## 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]
-     g = emis16_JUN[bbzone_CH16_JUN][ind_jr16]
-     h = tb16_JUN[bbzone_CH16_JUN][ind_jr16]
-     l = ts16_JUN[bbzone_CH16_JUN][ind_jr16]
-     emis16_JUN_moy[ijr] = mean(g[nonzero((g != -500.) & (g <= 1.))])
-     tb16_JUN_moy[ijr] = mean(h[nonzero((h != -500.) & (h != 0.))])
-     ts16_JUN_moy[ijr] = mean(l[nonzero((l != -500.) & (l != 0.))])
+#     ind_jr16 = np.where(jjr16_JUN[bbzone_CH16_JUN]==ijr+1)[0]
+#     g = emis16_JUN[bbzone_CH16_JUN][ind_jr16]
+#     h = tb16_JUN[bbzone_CH16_JUN][ind_jr16]
+#     l = ts16_JUN[bbzone_CH16_JUN][ind_jr16]
+#     emis16_JUN_moy[ijr] = mean(g[nonzero((g != -500.) & (g <= 1.))])
+#     tb16_JUN_moy[ijr] = mean(h[nonzero((h != -500.) & (h != 0.))])
+#     ts16_JUN_moy[ijr] = mean(l[nonzero((l != -500.) & (l != 0.))])
 
 
@@ -253 +307 @@
 twinx().plot(arange(0,30,1),ts17_JUN_moy,label='surface temperature',c='g')
 ylabel('surface temperature')
-plt.xticks(arange(1,31,1))
+plt.xticks(arange(0, 30, 1), arange(1, 31, 1))
 grid(True)
 plt.title('SSMIS - JUNE 2010')
+fig.show()
 
 
@@ -292 +347 @@
     yy = lat17_JUN[bbtb_ch17_JUN][ind_jr17_JUN]
     zz = tb17_JUN[bbtb_ch17_JUN][ind_jr17_JUN]
-#    zz0 = min(zz)
-#    zz1 = max(zz)
-#    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,zz0, zz1)
-#    outzch17_JUN=outz
-#    lonch17_JUN=outx
-#    latch17_JUN=outy
-#    OUTZCH17_JUN[:,:,ijr] = outzch17_JUN[:,:]
+    zz0 = min(zz)
+    zz1 = max(zz)
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch17_JUN = outz
+    lonch17_JUN = outx
+    latch17_JUN = outy
+    OUTZCH17_JUN[:,:,ijr] = outzch17_JUN[:,:]
     ## ch18 ##
     ind_jr18_JUN = np.where(jjr18_JUN[bbtb_ch18_JUN] == ijr+1)[0]
@@ -306 +361 @@
     zz0 = min(zz)
     zz1 = max(zz)
-    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0,x1,y0,y1,zz0, zz1)
-    outzch18_JUN=outz
-    lonch18_JUN=outx
-    latch18_JUN=outy
+    outz, outx, outy = ffgrid2.ffgrid(xx, yy, zz, dx, dy, x0, x1, y0, y1, zz0, zz1)
+    outzch18_JUN = outz
+    lonch18_JUN = outx
+    latch18_JUN = outy
     OUTZCH18_JUN[:,:,ijr] = outzch18_JUN[:,:]
 
@@ -339 +394 @@
 
 
-for ijr in range (12,23):
+for ijr in range (23, 30):
      figure()
      plt.ion()
-     m = Basemap(llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-90, urcrnrlat=-30, projection='cyl', resolution='c', fix_aspect=True)
+     m = Basemap(llcrnrlon=-60, urcrnrlon=-40, llcrnrlat=-85, urcrnrlat=-75, projection='cyl', resolution='c', fix_aspect=True)
      m.drawcoastlines(linewidth = 1)
-     m.drawparallels(np.arange(-90., 90., 20))
-     m.drawmeridians(np.arange(-180., 180., 20))
+     m.drawparallels(np.arange(-85., 75., 2))
+     m.drawmeridians(np.arange(-60., -40., 2))
      #m.fillcontinents()
-     clevs = arange(-30., 30., 1.)
+     clevs = arange(-25., 5., 0.1)
      xii,yii = m(*np.meshgrid(lonch17_JUN, latch17_JUN))
      cs = m.contourf(xii, yii, tbch17_anom_JUN[:,:,ijr], clevs, cmap=cm.s3pcpn_l_r)
      cbar = colorbar(cs)
      cbar.set_label('Tb anomaly SSMIS CH17 - JUNE')
-     plt.savefig('/usr/home/lahlod/twice_d/figures_output_ANTARC/mean_tb_anomaly_map_'+str(ijr+1)+'JUN_ch17_SSMIS')
+     xticks(np.arange(-60., -40., 2))
+     yticks(np.arange(-85., 75., 2))
+     plt.savefig('/usr/home/lahlod/twice_d/figures_output_ANTARC/SSMIS/mean_tb_anomaly_map_'+str(ijr+1)+'JUN_ch17_SSMIS_Zoom_zone2')
 
 
@@ -395 +452 @@
 ## tranche de latitude étudiée ##
 #bbtranche17_JUN = nonzero((latch17_JUN == -75.))
-#bbtranche17_JUN = nonzero((latch17_JUN >= -80.) & (latch17_JUN <= -75))
+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)