Changes in / [30:20]


Ignore:
Location:
/trunk/src
Files:
10 deleted
2 edited

Legend:

Unmodified
Added
Removed

  • /trunk/src/scripts_Laura/read_AMSUA_test.py

    r30 r20  
    99import netCDF4 
    1010 
    11 ################ fichiers par mois  ################################################### 
    12  
    13 ######### 
    14 ## CH1 ## 
    15 ######### 
    16  
    17 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH1_ANTARC_' 
    18 f2 = '2010.DAT' 
    19 date = np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    20  
    21 numlines_month = np.zeros([7],float) 
    22  
    23 for i in range (0, 7): 
    24     print 'mois =', date[i] 
    25     f = f1 + str(date[i]) + f2 
    26     fichier = open(f,'r') 
    27     numlines_month[i] = 0 
    28     for line in fichier : numlines_month[i] += 1 
    29  
    30     fichier.close() 
    31  
    32  
    33 i = 0 #JANUARY 
    34 f = f1 + str(date[i]) + f2 
    35 fichier = open(f,'r') 
    36 amch = np.zeros([20, numlines_month[i]], float) 
    37 for iligne in range (0, numlines_month[i]): 
    38     line = fichier.readline() 
    39     liste = line.split() 
    40     for j in range (0, 20): 
    41         amch[j,iligne] = float(liste[j]) 
    42  
    43 fichier.close() 
    44 amch1_JAN = amch 
    45 lon1_JAN = amch1_JAN[0,:] 
    46 lat1_JAN = amch1_JAN[1,:] 
    47 jjr1_JAN = amch1_JAN[4,:] 
    48 fov1_JAN = amch1_JAN[7,:] 
    49 ts1_JAN = amch1_JAN[10,:] 
    50 emis1_JAN = amch1_JAN[16,:] 
    51 tb1_JAN = amch1_JAN[15,:] 
    52 tup1_JAN = amch1_JAN[18,:] 
    53 tdn1_JAN = amch1_JAN[17,:] 
    54 trans1_JAN = amch1_JAN[19,:] 
    55 orog1_JAN = amch1_JAN[13,:] 
    56  
    57  
    58 i = 1 #FEBRUARY 
    59 f = f1 + str(date[i]) + f2 
    60 fichier = open(f,'r') 
    61 amch = np.zeros([20, numlines_month[i]], float) 
    62 for iligne in range (0, numlines_month[i]): 
    63     line = fichier.readline() 
    64     liste = line.split() 
    65     for j in range(0,20): 
    66         amch[j,iligne] = float(liste[j]) 
    67  
    68 fichier.close() 
    69 amch1_FEB = amch 
    70 lon1_FEB = amch1_FEB[0,:] 
    71 lat1_FEB = amch1_FEB[1,:] 
    72 jjr1_FEB = amch1_FEB[4,:] 
    73 fov1_FEB = amch1_FEB[7,:] 
    74 ts1_FEB = amch1_FEB[10,:] 
    75 emis1_FEB = amch1_FEB[16,:] 
    76 tb1_FEB = amch1_FEB[15,:] 
    77 tup1_FEB = amch1_FEB[18,:] 
    78 tdn1_FEB = amch1_FEB[17,:] 
    79 trans1_FEB = amch1_FEB[19,:] 
    80 orog1_FEB = amch1_FEB[13,:] 
    81  
    82  
    83 i = 2 #MARCH 
    84 f = f1 + str(date[i]) + f2 
    85 fichier = open(f, 'r') 
    86 amch = np.zeros([20, numlines_month[i]], float) 
    87 for iligne in range (0, numlines_month[i]): 
    88     line = fichier.readline() 
    89     liste = line.split() 
    90     for j in range (0, 20): 
    91         amch[j,iligne] = float(liste[j]) 
    92  
    93 fichier.close() 
    94 amch1_MAR = amch 
    95 lon1_MAR = amch1_MAR[0,:] 
    96 lat1_MAR = amch1_MAR[1,:] 
    97 jjr1_MAR = amch1_MAR[4,:] 
    98 fov1_MAR = amch1_MAR[7,:] 
    99 ts1_MAR = amch1_MAR[10,:] 
    100 emis1_MAR = amch1_MAR[16,:] 
    101 tb1_MAR = amch1_MAR[15,:] 
    102 tup1_MAR = amch1_MAR[18,:] 
    103 tdn1_MAR = amch1_MAR[17,:] 
    104 trans1_MAR = amch1_MAR[19,:] 
    105 orog1_MAR = amch1_MAR[13,:] 
    106  
    107  
    108 i = 3 #APRIL 
    109 f = f1 + str(date[i]) + f2 
    110 fichier = open(f, 'r') 
    111 amch = np.zeros([20, numlines_month[i]], float) 
    112 for iligne in range (0, numlines_month[i]): 
    113     line = fichier.readline() 
    114     liste = line.split() 
    115     for j in range (0, 20): 
    116         amch[j,iligne] = float(liste[j]) 
    117  
    118 fichier.close() 
    119 amch1_APR = amch 
    120 lon1_APR = amch1_APR[0,:] 
    121 lat1_APR = amch1_APR[1,:] 
    122 jjr1_APR = amch1_APR[4,:] 
    123 fov1_APR = amch1_APR[7,:] 
    124 ts1_APR = amch1_APR[10,:] 
    125 emis1_APR = amch1_APR[16,:] 
    126 tb1_APR = amch1_APR[15,:] 
    127 tup1_APR = amch1_APR[18,:] 
    128 tdn1_APR = amch1_APR[17,:] 
    129 trans1_APR = amch1_APR[19,:] 
    130 orog1_APR = amch1_APR[13,:] 
    131  
    132  
    133 i = 4 #MAY 
    134 f = f1 + str(date[i]) + f2 
    135 fichier = open(f, 'r') 
    136 amch=np.zeros([20,numlines_month[i]],float) 
    137 for iligne in range (0,numlines_month[i]): 
    138     line=fichier.readline() 
    139     liste = line.split() 
    140     for j in range(0,20): 
    141         amch[j,iligne] = float(liste[j]) 
    142  
    143 fichier.close() 
     11################ fichiers par mois - CH1 ################################################### 
     12 
     13f1='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH1_ANTARC_' 
     14f2='2010.DAT' 
     15date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     16 
     17numlines_month=np.zeros([7],float) 
     18 
     19for i in range (0,7): 
     20    f=f1+str(date[i])+f2 
     21    fichier=open(f,'r') 
     22    numlines_month[i]=0 
     23    for line in fichier: numlines_month[i]+=1 
     24 
     25    fichier.close() 
     26 
     27 
     28i=0 #JANUARY 
     29f=f1+str(date[i])+f2 
     30fichier=open(f,'r') 
     31amch=np.zeros([20,numlines_month[i]],float) 
     32for iligne in range (0,numlines_month[i]): 
     33    line=fichier.readline() 
     34    liste = line.split() 
     35    for j in range(0,20): 
     36        amch[j,iligne] = float(liste[j]) 
     37 
     38 
     39    fichier.close() 
     40 
     41 
     42amch1_JAN=amch 
     43#### def des variables  
     44lon1_JAN=amch1_JAN[0,:] 
     45lat1_JAN=amch1_JAN[1,:] 
     46jjr1_JAN=amch1_JAN[4,:] 
     47ts1_JAN=amch1_JAN[10,:] 
     48emis1_JAN=amch1_JAN[16,:] 
     49tb1_JAN=amch1_JAN[15,:] 
     50tup1_JAN=amch1_JAN[18,:] 
     51tdn1_JAN=amch1_JAN[17,:] 
     52trans1_JAN=amch1_JAN[19,:] 
     53orog1_JAN=amch1_JAN[13,:] 
     54 
     55 
     56i=1 #FEBRUARY 
     57f=f1+str(date[i])+f2 
     58fichier=open(f,'r') 
     59amch=np.zeros([20,numlines_month[i]],float) 
     60for iligne in range (0,numlines_month[i]): 
     61    line=fichier.readline() 
     62    liste = line.split() 
     63    for j in range(0,20): 
     64        amch[j,iligne] = float(liste[j]) 
     65 
     66 
     67    fichier.close() 
     68 
     69 
     70amch1_FEB=amch 
     71#### def des variables CH1 corrigees 
     72lon1_FEB=amch1_FEB[0,:] 
     73lat1_FEB=amch1_FEB[1,:] 
     74jjr1_FEB=amch1_FEB[4,:] 
     75ts1_FEB=amch1_FEB[10,:] 
     76emis1_FEB=amch1_FEB[16,:] 
     77tb1_FEB=amch1_FEB[15,:] 
     78tup1_FEB=amch1_FEB[18,:] 
     79tdn1_FEB=amch1_FEB[17,:] 
     80trans1_FEB=amch1_FEB[19,:] 
     81orog1_FEB=amch1_FEB[13,:] 
     82 
     83 
     84i=2 #MARCH 
     85f=f1+str(date[i])+f2 
     86fichier=open(f,'r') 
     87amch=np.zeros([20,numlines_month[i]],float) 
     88for iligne in range (0,numlines_month[i]): 
     89    line=fichier.readline() 
     90    liste = line.split() 
     91    for j in range(0,20): 
     92        amch[j,iligne] = float(liste[j]) 
     93 
     94 
     95    fichier.close() 
     96 
     97 
     98amch1_MAR=amch 
     99#### def des variables 
     100lon1_MAR=amch1_MAR[0,:] 
     101lat1_MAR=amch1_MAR[1,:] 
     102jjr1_MAR=amch1_MAR[4,:] 
     103ts1_MAR=amch1_MAR[10,:] 
     104emis1_MAR=amch1_MAR[16,:] 
     105tb1_MAR=amch1_MAR[15,:] 
     106tup1_MAR=amch1_MAR[18,:] 
     107tdn1_MAR=amch1_MAR[17,:] 
     108trans1_MAR=amch1_MAR[19,:] 
     109orog1_MAR=amch1_MAR[13,:] 
     110 
     111 
     112i=3 #APRIL 
     113f=f1+str(date[i])+f2 
     114fichier=open(f,'r') 
     115amch=np.zeros([20,numlines_month[i]],float) 
     116for iligne in range (0,numlines_month[i]): 
     117    line=fichier.readline() 
     118    liste = line.split() 
     119    for j in range(0,20): 
     120        amch[j,iligne] = float(liste[j]) 
     121 
     122 
     123    fichier.close() 
     124 
     125 
     126amch1_APR=amch 
     127#### def des variables  
     128lon1_APR=amch1_APR[0,:] 
     129lat1_APR=amch1_APR[1,:] 
     130jjr1_APR=amch1_APR[4,:] 
     131ts1_APR=amch1_APR[10,:] 
     132emis1_APR=amch1_APR[16,:] 
     133tb1_APR=amch1_APR[15,:] 
     134tup1_APR=amch1_APR[18,:] 
     135tdn1_APR=amch1_APR[17,:] 
     136trans1_APR=amch1_APR[19,:] 
     137orog1_APR=amch1_APR[13,:] 
     138 
     139 
     140i=4 #MAY 
     141f=f1+str(date[i])+f2 
     142fichier=open(f,'r') 
     143amch=np.zeros([20,numlines_month[i]],float) 
     144for iligne in range (0,numlines_month[i]): 
     145    line=fichier.readline() 
     146    liste = line.split() 
     147    for j in range(0,20): 
     148        amch[j,iligne] = float(liste[j]) 
     149 
     150 
     151    fichier.close() 
     152 
     153 
    144154amch1_MAY=amch 
     155#### def des variables 
    145156lon1_MAY=amch1_MAY[0,:] 
    146157lat1_MAY=amch1_MAY[1,:] 
    147158jjr1_MAY=amch1_MAY[4,:] 
    148 fov1_MAY = amch1_MAY[7,:] 
    149159ts1_MAY=amch1_MAY[10,:] 
    150160emis1_MAY=amch1_MAY[16,:] 
    … …  
    155165orog1_MAY=amch1_MAY[13,:] 
    156166 
    157  
    158167i=5 #JUNE 
    159168f=f1+str(date[i])+f2 
    … …  
    166175        amch[j,iligne] = float(liste[j]) 
    167176 
    168 fichier.close() 
     177 
     178    fichier.close() 
     179 
     180 
    169181amch1_JUN=amch 
     182#### def des variables 
    170183lon1_JUN=amch1_JUN[0,:] 
    171184lat1_JUN=amch1_JUN[1,:] 
    172185jjr1_JUN=amch1_JUN[4,:] 
    173 fov1_JUN = amch1_JUN[7,:] 
    174186ts1_JUN=amch1_JUN[10,:] 
    175187emis1_JUN=amch1_JUN[16,:] 
    … …  
    180192orog1_JUN=amch1_JUN[13,:] 
    181193 
    182  
    183 i=6 #JULY 
    184 f=f1+str(date[i])+f2 
    185 fichier=open(f,'r') 
    186 amch=np.zeros([20,numlines_month[i]],float) 
    187 for iligne in range (0,numlines_month[i]): 
    188     line=fichier.readline() 
    189     liste = line.split() 
    190     for j in range(0,20): 
    191         amch[j,iligne] = float(liste[j]) 
    192  
    193 fichier.close() 
     194i=6 #JULLY 
     195f=f1+str(date[i])+f2 
     196fichier=open(f,'r') 
     197amch=np.zeros([20,numlines_month[i]],float) 
     198for iligne in range (0,numlines_month[i]): 
     199    line=fichier.readline() 
     200    liste = line.split() 
     201    for j in range(0,20): 
     202        amch[j,iligne] = float(liste[j]) 
     203 
     204 
     205    fichier.close() 
     206 
     207 
    194208amch1_JUL=amch 
     209#### def des variables  
    195210lon1_JUL=amch1_JUL[0,:] 
    196211lat1_JUL=amch1_JUL[1,:] 
    197212jjr1_JUL=amch1_JUL[4,:] 
    198 fov1_JUL = amch1_JUL[7,:] 
    199213ts1_JUL=amch1_JUL[10,:] 
    200214emis1_JUL=amch1_JUL[16,:] 
    … …  
    206220 
    207221 
    208 ########## 
    209 ## CH15 ## 
    210 ########## 
    211  
    212 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH15_ANTARC_' 
    213 f2 = '2010.DAT' 
    214 date = np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    215  
    216 numlines_month = np.zeros([7],float) 
    217  
    218 for i in range (0, 7): 
    219     print 'mois =', date[i] 
    220     f = f1 + str(date[i]) + f2 
    221     fichier = open(f,'r') 
    222     numlines_month[i] = 0 
    223     for line in fichier : numlines_month[i] += 1 
    224  
    225     fichier.close() 
    226  
    227  
    228 i = 0 #JANUARY 
    229 f = f1 + str(date[i]) + f2 
    230 fichier = open(f,'r') 
    231 amch = np.zeros([20, numlines_month[i]], float) 
    232 for iligne in range (0, numlines_month[i]): 
    233     line = fichier.readline() 
    234     liste = line.split() 
    235     for j in range (0, 20): 
    236         amch[j,iligne] = float(liste[j]) 
    237  
    238 fichier.close() 
    239 amch15_JAN = amch 
    240 lon15_JAN = amch15_JAN[0,:] 
    241 lat15_JAN = amch15_JAN[1,:] 
    242 jjr15_JAN = amch15_JAN[4,:] 
    243 fov15_JAN = amch15_JAN[7,:] 
    244 ts15_JAN = amch15_JAN[10,:] 
    245 emis15_JAN = amch15_JAN[16,:] 
    246 tb15_JAN = amch15_JAN[15,:] 
    247 tup15_JAN = amch15_JAN[18,:] 
    248 tdn15_JAN = amch15_JAN[17,:] 
    249 trans15_JAN = amch15_JAN[19,:] 
    250 orog15_JAN = amch15_JAN[13,:] 
    251  
    252  
    253 i = 1 #FEBRUARY 
    254 f = f1 + str(date[i]) + f2 
    255 fichier = open(f,'r') 
    256 amch = np.zeros([20, numlines_month[i]], float) 
    257 for iligne in range (0, numlines_month[i]): 
    258     line = fichier.readline() 
    259     liste = line.split() 
    260     for j in range(0,20): 
    261         amch[j,iligne] = float(liste[j]) 
    262  
    263 fichier.close() 
    264 amch15_FEB = amch 
    265 lon15_FEB = amch15_FEB[0,:] 
    266 lat15_FEB = amch15_FEB[1,:] 
    267 jjr15_FEB = amch15_FEB[4,:] 
    268 fov15_FEB = amch15_FEB[7,:] 
    269 ts15_FEB = amch15_FEB[10,:] 
    270 emis15_FEB = amch15_FEB[16,:] 
    271 tb15_FEB = amch15_FEB[15,:] 
    272 tup15_FEB = amch15_FEB[18,:] 
    273 tdn15_FEB = amch15_FEB[17,:] 
    274 trans15_FEB = amch15_FEB[19,:] 
    275 orog15_FEB = amch15_FEB[13,:] 
    276  
    277  
    278 i = 2 #MARCH 
    279 f = f1 + str(date[i]) + f2 
    280 fichier = open(f, 'r') 
    281 amch = np.zeros([20, numlines_month[i]], float) 
    282 for iligne in range (0, numlines_month[i]): 
    283     line = fichier.readline() 
    284     liste = line.split() 
    285     for j in range (0, 20): 
    286         amch[j,iligne] = float(liste[j]) 
    287  
    288 fichier.close() 
    289 amch15_MAR = amch 
    290 lon15_MAR = amch15_MAR[0,:] 
    291 lat15_MAR = amch15_MAR[1,:] 
    292 jjr15_MAR = amch15_MAR[4,:] 
    293 fov15_MAR = amch15_MAR[7,:] 
    294 ts15_MAR = amch15_MAR[10,:] 
    295 emis15_MAR = amch15_MAR[16,:] 
    296 tb15_MAR = amch15_MAR[15,:] 
    297 tup15_MAR = amch15_MAR[18,:] 
    298 tdn15_MAR = amch15_MAR[17,:] 
    299 trans15_MAR = amch15_MAR[19,:] 
    300 orog15_MAR = amch15_MAR[13,:] 
    301  
    302  
    303 i = 3 #APRIL 
    304 f = f1 + str(date[i]) + f2 
    305 fichier = open(f, 'r') 
    306 amch = np.zeros([20, numlines_month[i]], float) 
    307 for iligne in range (0, numlines_month[i]): 
    308     line = fichier.readline() 
    309     liste = line.split() 
    310     for j in range (0, 20): 
    311         amch[j,iligne] = float(liste[j]) 
    312  
    313 fichier.close() 
    314 amch15_APR = amch 
    315 lon15_APR = amch15_APR[0,:] 
    316 lat15_APR = amch15_APR[1,:] 
    317 jjr15_APR = amch15_APR[4,:] 
    318 fov15_APR = amch15_APR[7,:] 
    319 ts15_APR = amch15_APR[10,:] 
    320 emis15_APR = amch15_APR[16,:] 
    321 tb15_APR = amch15_APR[15,:] 
    322 tup15_APR = amch15_APR[18,:] 
    323 tdn15_APR = amch15_APR[17,:] 
    324 trans15_APR = amch15_APR[19,:] 
    325 orog15_APR = amch15_APR[13,:] 
    326  
    327  
    328 i = 4 #MAY 
    329 f = f1 + str(date[i]) + f2 
    330 fichier = open(f, 'r') 
    331 amch=np.zeros([20,numlines_month[i]],float) 
    332 for iligne in range (0,numlines_month[i]): 
    333     line=fichier.readline() 
    334     liste = line.split() 
    335     for j in range(0,20): 
    336         amch[j,iligne] = float(liste[j]) 
    337  
    338 fichier.close() 
    339 amch15_MAY=amch 
    340 lon15_MAY=amch15_MAY[0,:] 
    341 lat15_MAY=amch15_MAY[1,:] 
    342 jjr15_MAY=amch15_MAY[4,:] 
    343 fov15_MAY = amch15_MAY[7,:] 
    344 ts15_MAY=amch15_MAY[10,:] 
    345 emis15_MAY=amch15_MAY[16,:] 
    346 tb15_MAY=amch15_MAY[15,:] 
    347 tup15_MAY=amch15_MAY[18,:] 
    348 tdn15_MAY=amch15_MAY[17,:] 
    349 trans15_MAY=amch15_MAY[19,:] 
    350 orog15_MAY=amch15_MAY[13,:] 
    351  
    352  
    353 i=5 #JUNE 
    354 f=f1+str(date[i])+f2 
    355 fichier=open(f,'r') 
    356 amch=np.zeros([20,numlines_month[i]],float) 
    357 for iligne in range (0,numlines_month[i]): 
    358     line=fichier.readline() 
    359     liste = line.split() 
    360     for j in range(0,20): 
    361         amch[j,iligne] = float(liste[j]) 
    362  
    363 fichier.close() 
    364 amch15_JUN=amch 
    365 lon15_JUN=amch15_JUN[0,:] 
    366 lat15_JUN=amch15_JUN[1,:] 
    367 jjr15_JUN=amch15_JUN[4,:] 
    368 fov15_JUN = amch15_JUN[7,:] 
    369 ts15_JUN=amch15_JUN[10,:] 
    370 emis15_JUN=amch15_JUN[16,:] 
    371 tb15_JUN=amch15_JUN[15,:] 
    372 tup15_JUN=amch15_JUN[18,:] 
    373 tdn15_JUN=amch15_JUN[17,:] 
    374 trans15_JUN=amch15_JUN[19,:] 
    375 orog15_JUN=amch15_JUN[13,:] 
    376  
    377  
    378 i=6 #JULY 
    379 f=f1+str(date[i])+f2 
    380 fichier=open(f,'r') 
    381 amch=np.zeros([20,numlines_month[i]],float) 
    382 for iligne in range (0,numlines_month[i]): 
    383     line=fichier.readline() 
    384     liste = line.split() 
    385     for j in range(0,20): 
    386         amch[j,iligne] = float(liste[j]) 
    387  
    388 fichier.close() 
    389 amch15_JUL = amch 
    390 lon15_JUL = amch15_JUL[0,:] 
    391 lat15_JUL = amch15_JUL[1,:] 
    392 jjr15_JUL = amch15_JUL[4,:] 
    393 fov15_JUL = amch15_JUL[7,:] 
    394 ts15_JUL = amch15_JUL[10,:] 
    395 emis15_JUL = amch15_JUL[16,:] 
    396 tb15_JUL = amch15_JUL[15,:] 
    397 tup15_JUL = amch15_JUL[18,:] 
    398 tdn15_JUL = amch15_JUL[17,:] 
    399 trans15_JUL = amch15_JUL[19,:] 
    400 orog15_JUL = amch15_JUL[13,:] 
    401222 
    402223 

  • /trunk/src/scripts_Laura/read_SSMIS_test.py

    r30 r20  
    77from mpl_toolkits.basemap import Basemap 
    88from mpl_toolkits.basemap import shiftgrid, cm 
    9 from netCDF4 import Dataset 
     9import netCDF4 
    1010import ffgrid2 
    1111 
    1212 
    13  
    14 ################ fichiers par canal - mois de juin ################################################### 
    15  
     13################ fichiers par mois - CH1 ################################################### 
    1614 
    1715f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH' 
    1816f3 = '_ANTARC_JUNE2010.DAT' 
     17#date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    1918channel = np.array([12, 13, 15, 16, 17, 18]) 
    2019numlines = np.zeros([len(channel)],int) 
    … …  
    2726     for line in fichier: numlines[ich] += 1 
    2827 
    29       fichier.close() 
    30  
    31  
    32  
    33  
    34  
    35 ich = 0 # 37GHz, H  
     28 
     29     fichier.close() 
     30 
     31ich = 2 # 37GHz, H polar 
    3632fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    3733ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    163159 
    164160 
    165 ich = 5 # 91.66GHz, H polar 
     161ich = 5 # 91.66GHz, V polar 
    166162fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    167163ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    188184orog18_JUN=ssch18_JUN[11,:] 
    189185 
    190  
    191  
    192  
    193 ############### fichiers par mois pour deux canaux (polars H et V) ################################## 
    194  
    195 ######## 
    196 # ch17 # 
    197 ######## 
    198  
    199 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH17_ANTARC_' 
    200 f3 = '2010.DAT' 
    201 date=np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    202 numlines = np.zeros([len(date)],int) 
    203  
    204 for 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  
    214 imo = 0 # FEBUARY 
    215 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    216 ssmis = np.zeros([18, numlines[imo]], float) 
    217 for 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  
    225 ssch17_FEB=ssmis 
    226 lon17_FEB=ssch17_FEB[0,:] 
    227 lat17_FEB=ssch17_FEB[1,:] 
    228 jjr17_FEB=ssch17_FEB[4,:] 
    229 ts17_FEB=ssch17_FEB[8,:] 
    230 emis17_FEB=ssch17_FEB[14,:] 
    231 tb17_FEB=ssch17_FEB[13,:] 
    232 tup17_FEB=ssch17_FEB[16,:] 
    233 tdn17_FEB=ssch17_FEB[15,:] 
    234 trans17_FEB=ssch17_FEB[17,:] 
    235 orog17_FEB=ssch17_FEB[11,:] 
    236  
    237  
    238 imo = 1 # APRIL 
    239 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    240 ssmis = np.zeros([18, numlines[imo]], float) 
    241 for 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  
    249 ssch17_APR=ssmis 
    250 lon17_APR=ssch17_APR[0,:] 
    251 lat17_APR=ssch17_APR[1,:] 
    252 jjr17_APR=ssch17_APR[4,:] 
    253 ts17_APR=ssch17_APR[8,:] 
    254 emis17_APR=ssch17_APR[14,:] 
    255 tb17_APR=ssch17_APR[13,:] 
    256 tup17_APR=ssch17_APR[16,:] 
    257 tdn17_APR=ssch17_APR[15,:] 
    258 trans17_APR=ssch17_APR[17,:] 
    259 orog17_APR=ssch17_APR[11,:] 
    260  
    261  
    262 imo = 2 # MAY 
    263 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    264 ssmis = np.zeros([18, numlines[imo]], float) 
    265 for 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  
    273 ssch17_MAY=ssmis 
    274 lon17_MAY=ssch17_MAY[0,:] 
    275 lat17_MAY=ssch17_MAY[1,:] 
    276 jjr17_MAY=ssch17_MAY[4,:] 
    277 ts17_MAY=ssch17_MAY[8,:] 
    278 emis17_MAY=ssch17_MAY[14,:] 
    279 tb17_MAY=ssch17_MAY[13,:] 
    280 tup17_MAY=ssch17_MAY[16,:] 
    281 tdn17_MAY=ssch17_MAY[15,:] 
    282 trans17_MAY=ssch17_MAY[17,:] 
    283 orog17_MAY=ssch17_MAY[11,:] 
    284  
    285  
    286 imo = 3 # JUNE 
    287 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    288 ssmis = np.zeros([18, numlines[imo]], float) 
    289 for 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  
    297 ssch17_JUN=ssmis 
    298 lon17_JUN=ssch17_JUN[0,:] 
    299 lat17_JUN=ssch17_JUN[1,:] 
    300 jjr17_JUN=ssch17_JUN[4,:] 
    301 ts17_JUN=ssch17_JUN[8,:] 
    302 emis17_JUN=ssch17_JUN[14,:] 
    303 tb17_JUN=ssch17_JUN[13,:] 
    304 tup17_JUN=ssch17_JUN[16,:] 
    305 tdn17_JUN=ssch17_JUN[15,:] 
    306 trans17_JUN=ssch17_JUN[17,:] 
    307 orog17_JUN=ssch17_JUN[11,:] 
    308  
    309  
    310 imo = 4 # JULY 
    311 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    312 ssmis = np.zeros([18, numlines[imo]], float) 
    313 for 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  
    321 ssch17_JUL=ssmis 
    322 lon17_JUL=ssch17_JUL[0,:] 
    323 lat17_JUL=ssch17_JUL[1,:] 
    324 jjr17_JUL=ssch17_JUL[4,:] 
    325 ts17_JUL=ssch17_JUL[8,:] 
    326 emis17_JUL=ssch17_JUL[14,:] 
    327 tb17_JUL=ssch17_JUL[13,:] 
    328 tup17_JUL=ssch17_JUL[16,:] 
    329 tdn17_JUL=ssch17_JUL[15,:] 
    330 trans17_JUL=ssch17_JUL[17,:] 
    331 orog17_JUL=ssch17_JUL[11,:] 
    332  
    333  
    334 ######## 
    335 # ch18 # 
    336 ######## 
    337  
    338 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH18_ANTARC_' 
    339 f3 = '2010.DAT' 
    340 date=np.array(['FEBRUARY', 'APRIL', 'JULY']) 
    341 numlines = np.zeros([len(date)],int) 
    342  
    343 for 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  
    353 imo = 0 # FEBUARY 
    354 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    355 ssmis = np.zeros([18, numlines[imo]], float) 
    356 for 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  
    364 ssch18_FEB=ssmis 
    365 lon18_FEB=ssch18_FEB[0,:] 
    366 lat18_FEB=ssch18_FEB[1,:] 
    367 jjr18_FEB=ssch18_FEB[4,:] 
    368 ts18_FEB=ssch18_FEB[8,:] 
    369 emis18_FEB=ssch18_FEB[14,:] 
    370 tb18_FEB=ssch18_FEB[13,:] 
    371 tup18_FEB=ssch18_FEB[16,:] 
    372 tdn18_FEB=ssch18_FEB[15,:] 
    373 trans18_FEB=ssch18_FEB[17,:] 
    374 orog18_FEB=ssch18_FEB[11,:] 
    375  
    376  
    377 imo = 1 # APRIL 
    378 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    379 ssmis = np.zeros([18, numlines[imo]], float) 
    380 for 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  
    388 ssch18_APR=ssmis 
    389 lon18_APR=ssch18_APR[0,:] 
    390 lat18_APR=ssch18_APR[1,:] 
    391 jjr18_APR=ssch18_APR[4,:] 
    392 ts18_APR=ssch18_APR[8,:] 
    393 emis18_APR=ssch18_APR[14,:] 
    394 tb18_APR=ssch18_APR[13,:] 
    395 tup18_APR=ssch18_APR[16,:] 
    396 tdn18_APR=ssch18_APR[15,:] 
    397 trans18_APR=ssch18_APR[17,:] 
    398 orog18_APR=ssch18_APR[11,:] 
    399  
    400  
    401  
    402 imo = 2 # JULY 
    403 fichier = open(f1 + str(date[imo]) + f3, 'r') 
    404 ssmis = np.zeros([18, numlines[imo]], float) 
    405 for 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  
    413 ssch18_JUL=ssmis 
    414 lon18_JUL=ssch18_JUL[0,:] 
    415 lat18_JUL=ssch18_JUL[1,:] 
    416 jjr18_JUL=ssch18_JUL[4,:] 
    417 ts18_JUL=ssch18_JUL[8,:] 
    418 emis18_JUL=ssch18_JUL[14,:] 
    419 tb18_JUL=ssch18_JUL[13,:] 
    420 tup18_JUL=ssch18_JUL[16,:] 
    421 tdn18_JUL=ssch18_JUL[15,:] 
    422 trans18_JUL=ssch18_JUL[17,:] 
    423 orog18_JUL=ssch18_JUL[11,:] 
    424  
    425  
    426 ########## 
    427 ## ch12 ## 
    428 ########## 
    429  
    430 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH12_ANTARC_' 
    431 f3 = '2010.DAT' 
    432 month = np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    433 numlines = np.zeros([len(month)],int) 
    434  
    435 for 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  
    445 imo = 0 # FEB  
    446 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    447 ssmis = np.zeros([18, numlines[imo]], float) 
    448 for 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  
    458 ssch12_FEB=ssmis 
    459 lon12_FEB=ssch12_FEB[0,:] 
    460 lat12_FEB=ssch12_FEB[1,:] 
    461 jjr12_FEB=ssch12_FEB[4,:] 
    462 ts12_FEB=ssch12_FEB[8,:] 
    463 emis12_FEB=ssch12_FEB[14,:] 
    464 tb12_FEB=ssch12_FEB[13,:] 
    465 tup12_FEB=ssch12_FEB[16,:] 
    466 tdn12_FEB=ssch12_FEB[15,:] 
    467 trans12_FEB=ssch12_FEB[17,:] 
    468 orog12_FEB=ssch12_FEB[11,:] 
    469  
    470 imo = 1 # APR  
    471 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    472 ssmis = np.zeros([18, numlines[imo]], float) 
    473 for 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  
    482 ssch12_APR=ssmis 
    483 lon12_APR=ssch12_APR[0,:] 
    484 lat12_APR=ssch12_APR[1,:] 
    485 jjr12_APR=ssch12_APR[4,:] 
    486 ts12_APR=ssch12_APR[8,:] 
    487 emis12_APR=ssch12_APR[14,:] 
    488 tb12_APR=ssch12_APR[13,:] 
    489 tup12_APR=ssch12_APR[16,:] 
    490 tdn12_APR=ssch12_APR[15,:] 
    491 trans12_APR=ssch12_APR[17,:] 
    492 orog12_APR=ssch12_APR[11,:] 
    493  
    494  
    495 imo = 4 # JUL  
    496 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    497 ssmis = np.zeros([18, numlines[imo]], float) 
    498 for 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  
    507 ssch12_JUL=ssmis 
    508 lon12_JUL=ssch12_JUL[0,:] 
    509 lat12_JUL=ssch12_JUL[1,:] 
    510 jjr12_JUL=ssch12_JUL[4,:] 
    511 ts12_JUL=ssch12_JUL[8,:] 
    512 emis12_JUL=ssch12_JUL[14,:] 
    513 tb12_JUL=ssch12_JUL[13,:] 
    514 tup12_JUL=ssch12_JUL[16,:] 
    515 tdn12_JUL=ssch12_JUL[15,:] 
    516 trans12_JUL=ssch12_JUL[17,:] 
    517 orog12_JUL=ssch12_JUL[11,:] 
    518  
    519  
    520 ########## 
    521 ## ch13 ## 
    522 ########## 
    523  
    524 f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH13_ANTARC_' 
    525 f3 = '2010.DAT' 
    526 month = np.array(['FEBRUARY', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    527 numlines = np.zeros([len(month)],int) 
    528  
    529 for 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  
    539 imo = 0 # FEB  
    540 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    541 ssmis = np.zeros([18, numlines[imo]], float) 
    542 for 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  
    552 ssch13_FEB=ssmis 
    553 lon13_FEB=ssch13_FEB[0,:] 
    554 lat13_FEB=ssch13_FEB[1,:] 
    555 jjr13_FEB=ssch13_FEB[4,:] 
    556 ts13_FEB=ssch13_FEB[8,:] 
    557 emis13_FEB=ssch13_FEB[14,:] 
    558 tb13_FEB=ssch13_FEB[13,:] 
    559 tup13_FEB=ssch13_FEB[16,:] 
    560 tdn13_FEB=ssch13_FEB[15,:] 
    561 trans13_FEB=ssch13_FEB[17,:] 
    562 orog13_FEB=ssch13_FEB[11,:] 
    563  
    564 imo = 1 # APR  
    565 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    566 ssmis = np.zeros([18, numlines[imo]], float) 
    567 for 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  
    576 ssch13_APR=ssmis 
    577 lon13_APR=ssch13_APR[0,:] 
    578 lat13_APR=ssch13_APR[1,:] 
    579 jjr13_APR=ssch13_APR[4,:] 
    580 ts13_APR=ssch13_APR[8,:] 
    581 emis13_APR=ssch13_APR[14,:] 
    582 tb13_APR=ssch13_APR[13,:] 
    583 tup13_APR=ssch13_APR[16,:] 
    584 tdn13_APR=ssch13_APR[15,:] 
    585 trans13_APR=ssch13_APR[17,:] 
    586 orog13_APR=ssch13_APR[11,:] 
    587  
    588  
    589 imo = 4 # JUL  
    590 fichier = open(f1 + str(month[imo]) + f3, 'r') 
    591 ssmis = np.zeros([18, numlines[imo]], float) 
    592 for 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  
    601 ssch13_JUL=ssmis 
    602 lon13_JUL=ssch13_JUL[0,:] 
    603 lat13_JUL=ssch13_JUL[1,:] 
    604 jjr13_JUL=ssch13_JUL[4,:] 
    605 ts13_JUL=ssch13_JUL[8,:] 
    606 emis13_JUL=ssch13_JUL[14,:] 
    607 tb13_JUL=ssch13_JUL[13,:] 
    608 tup13_JUL=ssch13_JUL[16,:] 
    609 tdn13_JUL=ssch13_JUL[15,:] 
    610 trans13_JUL=ssch13_JUL[17,:] 
    611 orog13_JUL=ssch13_JUL[11,:] 
     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)) 
     193bbzone_CH17_JUN = nonzero((lon17_JUN >= -60.) & (lon17_JUN <= -40.) & (lat17_JUN >= -85.) & (lat17_JUN <= -75.)) 
     194emis17_JUN_moy = np.zeros([30],float) 
     195tb17_JUN_moy = np.zeros([30],float) 
     196ts17_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)) 
     200bbzone_CH18_JUN = nonzero((lon18_JUN >= -60.) & (lon18_JUN <= -40.) & (lat18_JUN >= -85.) & (lat18_JUN <= -75.)) 
     201emis18_JUN_moy = np.zeros([30],float) 
     202tb18_JUN_moy = np.zeros([30],float) 
     203ts18_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 
     211for 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 ## 
     240fig=plt.figure() 
     241plt.subplot(3,1,1) 
     242plt.plot(arange(0,30,1),emis17_JUN_moy,label='91.66GHz-V',c='r') 
     243plt.plot(arange(0,30,1),emis18_JUN_moy,label='91.66GHz-H',c='b') 
     244plt.xticks(arange(1,31,1)) 
     245grid(True) 
     246ylabel('emissivity') 
     247legend(loc=7) 
     248plt.title('SSMIS - JUNE 2010') 
     249plt.subplot(3,1,2) 
     250plt.plot(arange(0,30,1),tb17_JUN_moy,label='91.66GHz-V',c='r') 
     251plt.plot(arange(0,30,1),tb18_JUN_moy,label='91.66GHz-H',c='b') 
     252plt.xticks(arange(1,31,1)) 
     253grid(True) 
     254ylabel('brightness temperature') 
     255legend(loc=7) 
     256plt.subplot(3,1,3) 
     257plt.plot(arange(0,30,1),ts17_JUN_moy,label='91.66GHz-V',c='r') 
     258plt.plot(arange(0,30,1),ts18_JUN_moy,label='91.66GHz-H',c='b') 
     259plt.xticks(arange(1,31,1)) 
     260grid(True) 
     261ylabel('surface temperature') 
     262legend(loc=7) 
     263fig.show() 
     264 
     265## ch16 ## 
     266fig=plt.figure() 
     267plt.subplot(3,1,1) 
     268plt.plot(arange(0,30,1),emis16_JUN_moy,label='37GHz-V',c='r') 
     269plt.plot(arange(0,30,1),emis17_JUN_moy,label='91.66GHz-V',c='g') 
     270plt.xticks(arange(1,31,1)) 
     271grid(True) 
     272ylabel('emissivity') 
     273legend(loc=7) 
     274plt.title('SSMIS - JUNE 2010') 
     275plt.subplot(3,1,2) 
     276plt.plot(arange(0,30,1),tb16_JUN_moy,label='37GHz-V',c='r') 
     277plt.plot(arange(0,30,1),tb17_JUN_moy,label='91.66GHz-V',c='g') 
     278plt.xticks(arange(1,31,1)) 
     279grid(True) 
     280ylabel('brightness temperature') 
     281legend(loc=7) 
     282plt.subplot(3,1,3) 
     283plt.plot(arange(0,30,1),ts16_JUN_moy,label='37GHz-V',c='r') 
     284plt.plot(arange(0,30,1),ts17_JUN_moy,label='91.66GHz-V',c='g') 
     285plt.xticks(arange(1,31,1)) 
     286grid(True) 
     287ylabel('surface temperature') 
     288legend(loc=7) 
     289fig.show() 
     290 
     291 
     292 
     293tb17_JUN_anom = np.zeros([30],float) 
     294tb18_JUN_anom = np.zeros([30],float) 
     295for 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 
     301fig=plt.figure() 
     302plt.plot(arange(0,30,1),tb17_JUN_anom,label='91.66Ghz-V',c='r') 
     303plt.plot(arange(0,30,1),tb18_JUN_anom,label='91.66GHz-H',c='b') 
     304plt.plot(arange(0,31,1),np.zeros([31]),'--',c='k') 
     305ylabel('Tb anomaly') 
     306legend() 
     307twinx().plot(arange(0,30,1),ts17_JUN_moy,label='surface temperature',c='g') 
     308ylabel('surface temperature') 
     309plt.xticks(arange(0, 30, 1), arange(1, 31, 1)) 
     310grid(True) 
     311plt.title('SSMIS - JUNE 2010') 
     312fig.show() 
     313 
     314 
     315 
     316 
     317################ 
     318# CARTOGRAPHIE # 
     319################ 
     320## Etude sur l'Antarctique ## 
     321dx=5. 
     322dy=5. 
     323x0, x1 = -180, 180 
     324y0, y1 = -90, -30 
     325 
     326## ch17 ## 
     327bbemis_ch17_JUN=nonzero((emis17_JUN!=-500.)&(emis17_JUN<1.)&(emis17_JUN>0.)) 
     328bbtb_ch17_JUN=nonzero((tb17_JUN!=-500.)&(tb17_JUN!=0.)) 
     329OUTZCH17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float) 
     330outzch17_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float) 
     331lonch17_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float) 
     332latch17_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float) 
     333 
     334## ch18 ## 
     335bbemis_ch18_JUN=nonzero((emis18_JUN!=-500.)&(emis18_JUN<1.)&(emis18_JUN>0.)) 
     336bbtb_ch18_JUN=nonzero((tb18_JUN!=-500.)&(tb18_JUN!=0.)) 
     337OUTZCH18_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),30], float) 
     338outzch18_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx))], float) 
     339lonch18_JUN=np.zeros([len(np.arange(x0, x1+1, dx))], float) 
     340latch18_JUN=np.zeros([len(np.arange(y0, y1+1, dy))], float) 
     341 
     342for 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 ## 
     371mean_outzch17_JUN = np.zeros([len(latch17_JUN), len(lonch17_JUN)], float) 
     372for 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 ## 
     378mean_outzch18_JUN = np.zeros([len(latch18_JUN), len(lonch18_JUN)], float) 
     379for 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 ## 
     385tbch17_anom_JUN = np.zeros([len(latch17_JUN),len(lonch17_JUN),30], float) 
     386for ijr in range (0,30): 
     387    tbch17_anom_JUN[:,:,ijr] = OUTZCH17_JUN[:,:,ijr] - mean_outzch17_JUN[:,:] 
     388 
     389 
     390## ch18 ## 
     391tbch18_anom_JUN = np.zeros([len(latch18_JUN),len(lonch18_JUN),30], float) 
     392for ijr in range (0,30): 
     393    tbch18_anom_JUN[:,:,ijr] = OUTZCH18_JUN[:,:,ijr] - mean_outzch18_JUN[:,:] 
     394 
     395 
     396for 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 
     414figure() 
     415plt.ion() 
     416m = Basemap(llcrnrlon=-180, urcrnrlon=180, llcrnrlat=-90, urcrnrlat=-30, projection='cyl', resolution='c', fix_aspect=True) 
     417m.drawcoastlines(linewidth = 1) 
     418m.drawparallels(np.arange(-90., 90., 20)) 
     419m.drawmeridians(np.arange(-180., 180., 20)) 
     420#m.fillcontinents() 
     421xii,yii = m(*np.meshgrid(lonch17_JUN, latch17_JUN)) 
     422 
     423## BIAIS ch17-ch18 ##  
     424biais_JUN = mean_outzch17_JUN - mean_outzch18_JUN 
     425clevs = arange(0., 45., 1.) 
     426cs = m.contourf(xii, yii, biais_JUN, clevs, cmap=cm.s3pcpn_l_r) 
     427cbar = colorbar(cs) 
     428cbar.set_label('Bias [CH17-CH18] of Mean Tb JUNE - SSMIS') 
     429 
     430## VARIANCE ch17 - ch18 ## 
     431std1 = np.zeros([len(latch17_JUN), len(lonch17_JUN)], float) 
     432for 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 
     437N = len(lonch17_JUN)*len(latch17_JUN) 
     438std_JUN = std1/N 
     439clevs = arange(0., 25., 0.1) 
     440cs = m.contourf(xii, yii, std_JUN, clevs, cmap=cm.s3pcpn_l_r) 
     441cbar = colorbar(cs) 
     442cbar.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.)) 
     454bbtranche17_JUN = nonzero((latch17_JUN >= -85.) & (latch17_JUN <= -75)) 
     455#bbtranche17_JUN = nonzero((latch17_JUN >= -90.) & (latch17_JUN <= -85)) 
     456mean_tbch17_anom_JUN = np.zeros([len(lonch17_JUN), 30], float) 
     457for 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 
     462y_time, x_space = m(*np.meshgrid(arange(0,30,1), lonch17_JUN)) 
     463fig = plt.figure() 
     464plt.pcolor(x_space, y_time, mean_tbch17_anom_JUN, cmap=cm.s3pcpn_l_r, vmin = -30., vmax = 25.) 
     465plt.axis([-180., 180., 0, 29]) 
     466cb = plt.colorbar() 
     467cb.set_label('Tb anomaly - SSMIS CH17') 
     468plt.xticks(arange(-180.,220.,40)) 
     469plt.yticks(arange(0, 30, 1), arange(1,31,1)) 
     470plt.xlabel('longitude') 
     471plt.ylabel('JUNE 2010') 
     472 
Note: See TracChangeset for help on using the changeset viewer.