Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/scripts_Laura/read_AMSUA_test.py

    r19 r24  
    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 
Note: See TracChangeset for help on using the changeset viewer.