Changes in / [20:30]


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

Legend:

Unmodified
Added
Removed

  • /trunk/src/scripts_Laura/read_AMSUA_test.py

    r20 r30  
    99import netCDF4 
    1010 
    11 ################ fichiers par mois - CH1 ################################################### 
    12  
    13 f1='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH1_ANTARC_' 
    14 f2='2010.DAT' 
    15 date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    16  
    17 numlines_month=np.zeros([7],float) 
    18  
    19 for 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 
     11################ fichiers par mois  ################################################### 
     12 
     13######### 
     14## CH1 ## 
     15######### 
     16 
     17f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH1_ANTARC_' 
     18f2 = '2010.DAT' 
     19date = np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     20 
     21numlines_month = np.zeros([7],float) 
     22 
     23for 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 
    2429 
    2530    fichier.close() 
    2631 
    2732 
    28 i=0 #JANUARY 
    29 f=f1+str(date[i])+f2 
    30 fichier=open(f,'r') 
     33i = 0 #JANUARY 
     34f = f1 + str(date[i]) + f2 
     35fichier = open(f,'r') 
     36amch = np.zeros([20, numlines_month[i]], float) 
     37for 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 
     43fichier.close() 
     44amch1_JAN = amch 
     45lon1_JAN = amch1_JAN[0,:] 
     46lat1_JAN = amch1_JAN[1,:] 
     47jjr1_JAN = amch1_JAN[4,:] 
     48fov1_JAN = amch1_JAN[7,:] 
     49ts1_JAN = amch1_JAN[10,:] 
     50emis1_JAN = amch1_JAN[16,:] 
     51tb1_JAN = amch1_JAN[15,:] 
     52tup1_JAN = amch1_JAN[18,:] 
     53tdn1_JAN = amch1_JAN[17,:] 
     54trans1_JAN = amch1_JAN[19,:] 
     55orog1_JAN = amch1_JAN[13,:] 
     56 
     57 
     58i = 1 #FEBRUARY 
     59f = f1 + str(date[i]) + f2 
     60fichier = open(f,'r') 
     61amch = np.zeros([20, numlines_month[i]], float) 
     62for 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 
     68fichier.close() 
     69amch1_FEB = amch 
     70lon1_FEB = amch1_FEB[0,:] 
     71lat1_FEB = amch1_FEB[1,:] 
     72jjr1_FEB = amch1_FEB[4,:] 
     73fov1_FEB = amch1_FEB[7,:] 
     74ts1_FEB = amch1_FEB[10,:] 
     75emis1_FEB = amch1_FEB[16,:] 
     76tb1_FEB = amch1_FEB[15,:] 
     77tup1_FEB = amch1_FEB[18,:] 
     78tdn1_FEB = amch1_FEB[17,:] 
     79trans1_FEB = amch1_FEB[19,:] 
     80orog1_FEB = amch1_FEB[13,:] 
     81 
     82 
     83i = 2 #MARCH 
     84f = f1 + str(date[i]) + f2 
     85fichier = open(f, 'r') 
     86amch = np.zeros([20, numlines_month[i]], float) 
     87for 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 
     93fichier.close() 
     94amch1_MAR = amch 
     95lon1_MAR = amch1_MAR[0,:] 
     96lat1_MAR = amch1_MAR[1,:] 
     97jjr1_MAR = amch1_MAR[4,:] 
     98fov1_MAR = amch1_MAR[7,:] 
     99ts1_MAR = amch1_MAR[10,:] 
     100emis1_MAR = amch1_MAR[16,:] 
     101tb1_MAR = amch1_MAR[15,:] 
     102tup1_MAR = amch1_MAR[18,:] 
     103tdn1_MAR = amch1_MAR[17,:] 
     104trans1_MAR = amch1_MAR[19,:] 
     105orog1_MAR = amch1_MAR[13,:] 
     106 
     107 
     108i = 3 #APRIL 
     109f = f1 + str(date[i]) + f2 
     110fichier = open(f, 'r') 
     111amch = np.zeros([20, numlines_month[i]], float) 
     112for 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 
     118fichier.close() 
     119amch1_APR = amch 
     120lon1_APR = amch1_APR[0,:] 
     121lat1_APR = amch1_APR[1,:] 
     122jjr1_APR = amch1_APR[4,:] 
     123fov1_APR = amch1_APR[7,:] 
     124ts1_APR = amch1_APR[10,:] 
     125emis1_APR = amch1_APR[16,:] 
     126tb1_APR = amch1_APR[15,:] 
     127tup1_APR = amch1_APR[18,:] 
     128tdn1_APR = amch1_APR[17,:] 
     129trans1_APR = amch1_APR[19,:] 
     130orog1_APR = amch1_APR[13,:] 
     131 
     132 
     133i = 4 #MAY 
     134f = f1 + str(date[i]) + f2 
     135fichier = open(f, 'r') 
    31136amch=np.zeros([20,numlines_month[i]],float) 
    32137for iligne in range (0,numlines_month[i]): 
    … …  
    36141        amch[j,iligne] = float(liste[j]) 
    37142 
    38  
    39     fichier.close() 
    40  
    41  
    42 amch1_JAN=amch 
    43 #### def des variables  
    44 lon1_JAN=amch1_JAN[0,:] 
    45 lat1_JAN=amch1_JAN[1,:] 
    46 jjr1_JAN=amch1_JAN[4,:] 
    47 ts1_JAN=amch1_JAN[10,:] 
    48 emis1_JAN=amch1_JAN[16,:] 
    49 tb1_JAN=amch1_JAN[15,:] 
    50 tup1_JAN=amch1_JAN[18,:] 
    51 tdn1_JAN=amch1_JAN[17,:] 
    52 trans1_JAN=amch1_JAN[19,:] 
    53 orog1_JAN=amch1_JAN[13,:] 
    54  
    55  
    56 i=1 #FEBRUARY 
    57 f=f1+str(date[i])+f2 
    58 fichier=open(f,'r') 
    59 amch=np.zeros([20,numlines_month[i]],float) 
    60 for 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  
    70 amch1_FEB=amch 
    71 #### def des variables CH1 corrigees 
    72 lon1_FEB=amch1_FEB[0,:] 
    73 lat1_FEB=amch1_FEB[1,:] 
    74 jjr1_FEB=amch1_FEB[4,:] 
    75 ts1_FEB=amch1_FEB[10,:] 
    76 emis1_FEB=amch1_FEB[16,:] 
    77 tb1_FEB=amch1_FEB[15,:] 
    78 tup1_FEB=amch1_FEB[18,:] 
    79 tdn1_FEB=amch1_FEB[17,:] 
    80 trans1_FEB=amch1_FEB[19,:] 
    81 orog1_FEB=amch1_FEB[13,:] 
    82  
    83  
    84 i=2 #MARCH 
    85 f=f1+str(date[i])+f2 
    86 fichier=open(f,'r') 
    87 amch=np.zeros([20,numlines_month[i]],float) 
    88 for 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  
    98 amch1_MAR=amch 
    99 #### def des variables 
    100 lon1_MAR=amch1_MAR[0,:] 
    101 lat1_MAR=amch1_MAR[1,:] 
    102 jjr1_MAR=amch1_MAR[4,:] 
    103 ts1_MAR=amch1_MAR[10,:] 
    104 emis1_MAR=amch1_MAR[16,:] 
    105 tb1_MAR=amch1_MAR[15,:] 
    106 tup1_MAR=amch1_MAR[18,:] 
    107 tdn1_MAR=amch1_MAR[17,:] 
    108 trans1_MAR=amch1_MAR[19,:] 
    109 orog1_MAR=amch1_MAR[13,:] 
    110  
    111  
    112 i=3 #APRIL 
    113 f=f1+str(date[i])+f2 
    114 fichier=open(f,'r') 
    115 amch=np.zeros([20,numlines_month[i]],float) 
    116 for 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  
    126 amch1_APR=amch 
    127 #### def des variables  
    128 lon1_APR=amch1_APR[0,:] 
    129 lat1_APR=amch1_APR[1,:] 
    130 jjr1_APR=amch1_APR[4,:] 
    131 ts1_APR=amch1_APR[10,:] 
    132 emis1_APR=amch1_APR[16,:] 
    133 tb1_APR=amch1_APR[15,:] 
    134 tup1_APR=amch1_APR[18,:] 
    135 tdn1_APR=amch1_APR[17,:] 
    136 trans1_APR=amch1_APR[19,:] 
    137 orog1_APR=amch1_APR[13,:] 
    138  
    139  
    140 i=4 #MAY 
    141 f=f1+str(date[i])+f2 
    142 fichier=open(f,'r') 
    143 amch=np.zeros([20,numlines_month[i]],float) 
    144 for 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  
     143fichier.close() 
    154144amch1_MAY=amch 
    155 #### def des variables 
    156145lon1_MAY=amch1_MAY[0,:] 
    157146lat1_MAY=amch1_MAY[1,:] 
    158147jjr1_MAY=amch1_MAY[4,:] 
     148fov1_MAY = amch1_MAY[7,:] 
    159149ts1_MAY=amch1_MAY[10,:] 
    160150emis1_MAY=amch1_MAY[16,:] 
    … …  
    165155orog1_MAY=amch1_MAY[13,:] 
    166156 
     157 
    167158i=5 #JUNE 
    168159f=f1+str(date[i])+f2 
    … …  
    175166        amch[j,iligne] = float(liste[j]) 
    176167 
    177  
    178     fichier.close() 
    179  
    180  
     168fichier.close() 
    181169amch1_JUN=amch 
    182 #### def des variables 
    183170lon1_JUN=amch1_JUN[0,:] 
    184171lat1_JUN=amch1_JUN[1,:] 
    185172jjr1_JUN=amch1_JUN[4,:] 
     173fov1_JUN = amch1_JUN[7,:] 
    186174ts1_JUN=amch1_JUN[10,:] 
    187175emis1_JUN=amch1_JUN[16,:] 
    … …  
    192180orog1_JUN=amch1_JUN[13,:] 
    193181 
    194 i=6 #JULLY 
     182 
     183i=6 #JULY 
    195184f=f1+str(date[i])+f2 
    196185fichier=open(f,'r') 
    … …  
    202191        amch[j,iligne] = float(liste[j]) 
    203192 
    204  
    205     fichier.close() 
    206  
    207  
     193fichier.close() 
    208194amch1_JUL=amch 
    209 #### def des variables  
    210195lon1_JUL=amch1_JUL[0,:] 
    211196lat1_JUL=amch1_JUL[1,:] 
    212197jjr1_JUL=amch1_JUL[4,:] 
     198fov1_JUL = amch1_JUL[7,:] 
    213199ts1_JUL=amch1_JUL[10,:] 
    214200emis1_JUL=amch1_JUL[16,:] 
    … …  
    220206 
    221207 
     208########## 
     209## CH15 ## 
     210########## 
     211 
     212f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH15_ANTARC_' 
     213f2 = '2010.DAT' 
     214date = np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
     215 
     216numlines_month = np.zeros([7],float) 
     217 
     218for 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 
     228i = 0 #JANUARY 
     229f = f1 + str(date[i]) + f2 
     230fichier = open(f,'r') 
     231amch = np.zeros([20, numlines_month[i]], float) 
     232for 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 
     238fichier.close() 
     239amch15_JAN = amch 
     240lon15_JAN = amch15_JAN[0,:] 
     241lat15_JAN = amch15_JAN[1,:] 
     242jjr15_JAN = amch15_JAN[4,:] 
     243fov15_JAN = amch15_JAN[7,:] 
     244ts15_JAN = amch15_JAN[10,:] 
     245emis15_JAN = amch15_JAN[16,:] 
     246tb15_JAN = amch15_JAN[15,:] 
     247tup15_JAN = amch15_JAN[18,:] 
     248tdn15_JAN = amch15_JAN[17,:] 
     249trans15_JAN = amch15_JAN[19,:] 
     250orog15_JAN = amch15_JAN[13,:] 
     251 
     252 
     253i = 1 #FEBRUARY 
     254f = f1 + str(date[i]) + f2 
     255fichier = open(f,'r') 
     256amch = np.zeros([20, numlines_month[i]], float) 
     257for 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 
     263fichier.close() 
     264amch15_FEB = amch 
     265lon15_FEB = amch15_FEB[0,:] 
     266lat15_FEB = amch15_FEB[1,:] 
     267jjr15_FEB = amch15_FEB[4,:] 
     268fov15_FEB = amch15_FEB[7,:] 
     269ts15_FEB = amch15_FEB[10,:] 
     270emis15_FEB = amch15_FEB[16,:] 
     271tb15_FEB = amch15_FEB[15,:] 
     272tup15_FEB = amch15_FEB[18,:] 
     273tdn15_FEB = amch15_FEB[17,:] 
     274trans15_FEB = amch15_FEB[19,:] 
     275orog15_FEB = amch15_FEB[13,:] 
     276 
     277 
     278i = 2 #MARCH 
     279f = f1 + str(date[i]) + f2 
     280fichier = open(f, 'r') 
     281amch = np.zeros([20, numlines_month[i]], float) 
     282for 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 
     288fichier.close() 
     289amch15_MAR = amch 
     290lon15_MAR = amch15_MAR[0,:] 
     291lat15_MAR = amch15_MAR[1,:] 
     292jjr15_MAR = amch15_MAR[4,:] 
     293fov15_MAR = amch15_MAR[7,:] 
     294ts15_MAR = amch15_MAR[10,:] 
     295emis15_MAR = amch15_MAR[16,:] 
     296tb15_MAR = amch15_MAR[15,:] 
     297tup15_MAR = amch15_MAR[18,:] 
     298tdn15_MAR = amch15_MAR[17,:] 
     299trans15_MAR = amch15_MAR[19,:] 
     300orog15_MAR = amch15_MAR[13,:] 
     301 
     302 
     303i = 3 #APRIL 
     304f = f1 + str(date[i]) + f2 
     305fichier = open(f, 'r') 
     306amch = np.zeros([20, numlines_month[i]], float) 
     307for 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 
     313fichier.close() 
     314amch15_APR = amch 
     315lon15_APR = amch15_APR[0,:] 
     316lat15_APR = amch15_APR[1,:] 
     317jjr15_APR = amch15_APR[4,:] 
     318fov15_APR = amch15_APR[7,:] 
     319ts15_APR = amch15_APR[10,:] 
     320emis15_APR = amch15_APR[16,:] 
     321tb15_APR = amch15_APR[15,:] 
     322tup15_APR = amch15_APR[18,:] 
     323tdn15_APR = amch15_APR[17,:] 
     324trans15_APR = amch15_APR[19,:] 
     325orog15_APR = amch15_APR[13,:] 
     326 
     327 
     328i = 4 #MAY 
     329f = f1 + str(date[i]) + f2 
     330fichier = open(f, 'r') 
     331amch=np.zeros([20,numlines_month[i]],float) 
     332for 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 
     338fichier.close() 
     339amch15_MAY=amch 
     340lon15_MAY=amch15_MAY[0,:] 
     341lat15_MAY=amch15_MAY[1,:] 
     342jjr15_MAY=amch15_MAY[4,:] 
     343fov15_MAY = amch15_MAY[7,:] 
     344ts15_MAY=amch15_MAY[10,:] 
     345emis15_MAY=amch15_MAY[16,:] 
     346tb15_MAY=amch15_MAY[15,:] 
     347tup15_MAY=amch15_MAY[18,:] 
     348tdn15_MAY=amch15_MAY[17,:] 
     349trans15_MAY=amch15_MAY[19,:] 
     350orog15_MAY=amch15_MAY[13,:] 
     351 
     352 
     353i=5 #JUNE 
     354f=f1+str(date[i])+f2 
     355fichier=open(f,'r') 
     356amch=np.zeros([20,numlines_month[i]],float) 
     357for 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 
     363fichier.close() 
     364amch15_JUN=amch 
     365lon15_JUN=amch15_JUN[0,:] 
     366lat15_JUN=amch15_JUN[1,:] 
     367jjr15_JUN=amch15_JUN[4,:] 
     368fov15_JUN = amch15_JUN[7,:] 
     369ts15_JUN=amch15_JUN[10,:] 
     370emis15_JUN=amch15_JUN[16,:] 
     371tb15_JUN=amch15_JUN[15,:] 
     372tup15_JUN=amch15_JUN[18,:] 
     373tdn15_JUN=amch15_JUN[17,:] 
     374trans15_JUN=amch15_JUN[19,:] 
     375orog15_JUN=amch15_JUN[13,:] 
     376 
     377 
     378i=6 #JULY 
     379f=f1+str(date[i])+f2 
     380fichier=open(f,'r') 
     381amch=np.zeros([20,numlines_month[i]],float) 
     382for 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 
     388fichier.close() 
     389amch15_JUL = amch 
     390lon15_JUL = amch15_JUL[0,:] 
     391lat15_JUL = amch15_JUL[1,:] 
     392jjr15_JUL = amch15_JUL[4,:] 
     393fov15_JUL = amch15_JUL[7,:] 
     394ts15_JUL = amch15_JUL[10,:] 
     395emis15_JUL = amch15_JUL[16,:] 
     396tb15_JUL = amch15_JUL[15,:] 
     397tup15_JUL = amch15_JUL[18,:] 
     398tdn15_JUL = amch15_JUL[17,:] 
     399trans15_JUL = amch15_JUL[19,:] 
     400orog15_JUL = amch15_JUL[13,:] 
    222401 
    223402 

  • /trunk/src/scripts_Laura/read_SSMIS_test.py

    r20 r30  
    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.