Changeset 19


Ignore:
Timestamp:
05/23/14 18:34:24 (10 years ago)
Author:
lahlod
Message:

nouveaux scripts

Location:
trunk/src/scripts_Laura
Files:
3 added
2 deleted
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/scripts_Laura/read_AMSUA_test.py

    r18 r19  
    3737 
    3838 
    39 fichier.close 
     39    fichier.close() 
     40 
    4041 
    4142amch1_JAN=amch 
    42  
    43 #### enlever les valeurs a -500.000 
    44 #bb_ts_JAN=nonzero(amch_JAN[10,:]!=-500.) 
    45 #bb_emis_JAN=nonzero((amch_JAN[16,:]!=-500.)&(amch_JAN[16,:]<=1.)) 
    46 #bb_tb_JAN=nonzero(amch_JAN[15,:]!=-500.) 
    47 #bb_tup_JAN=nonzero(amch_JAN[18,:]!=-500.) 
    48 #bb_tdn_JAN=nonzero(amch_JAN[17,:]!=-500.) 
    49 #bb_trans_JAN=nonzero(amch_JAN[19,:]!=-500.) 
    50 #bb_orog_JAN=nonzero(amch_JAN[13,:]>=0.) 
    51  
    5243#### def des variables  
    5344lon1_JAN=amch1_JAN[0,:] 
    … …  
    7465 
    7566 
    76 fichier.close 
     67    fichier.close() 
     68 
    7769 
    7870amch1_FEB=amch 
    79  
    80 #### enlever les valeurs a -500.000 
    81 #bb_ts_FEB=nonzero(amch_FEB[10,:]!=-500) 
    82 #bb_emis_FEB=nonzero((amch_FEB[16,:]!=-500)&(amch_FEB[16,:]<=1.)) 
    83 #bb_tb_FEB=nonzero(amch_FEB[15,:]!=-500) 
    84 #bb_tup_FEB=nonzero(amch_FEB[18,:]!=-500) 
    85 #bb_tdn_FEB=nonzero(amch_FEB[17,:]!=-500) 
    86 #bb_trans_FEB=nonzero(amch_FEB[19,:]!=-500) 
    87 #bb_orog_FEB=nonzero(amch_FEB[13,:]>=0.) 
    88  
    8971#### def des variables CH1 corrigees 
    9072lon1_FEB=amch1_FEB[0,:] 
    … …  
    11193 
    11294 
    113 fichier.close 
     95    fichier.close() 
     96 
    11497 
    11598amch1_MAR=amch 
    116  
    117 #### enlever les valeurs a -500.000 
    118 #bb_ts_MAR=nonzero(amch_MAR[10,:]!=-500) 
    119 #bb_emis_MAR=nonzero((amch_MAR[16,:]!=-500)&(amch_MAR[16,:]<=1.)) 
    120 #bb_tb_MAR=nonzero(amch_MAR[15,:]!=-500) 
    121 #bb_tup_MAR=nonzero(amch_MAR[18,:]!=-500) 
    122 #bb_tdn_MER=nonzero(amch_MAR[17,:]!=-500) 
    123 #bb_trans_MAR=nonzero(amch_MAR[19,:]!=-500) 
    124 #bb_orog_MAR=nonzero(amch_MAR[13,:]>=0.) 
    125  
    12699#### def des variables 
    127100lon1_MAR=amch1_MAR[0,:] 
    … …  
    148121 
    149122 
    150 fichier.close 
     123    fichier.close() 
     124 
    151125 
    152126amch1_APR=amch 
    153  
    154127#### def des variables  
    155128lon1_APR=amch1_APR[0,:] 
    … …  
    176149 
    177150 
    178 fichier.close 
     151    fichier.close() 
     152 
    179153 
    180154amch1_MAY=amch 
    181  
    182155#### def des variables 
    183156lon1_MAY=amch1_MAY[0,:] 
    … …  
    203176 
    204177 
    205 fichier.close 
     178    fichier.close() 
     179 
    206180 
    207181amch1_JUN=amch 
    208  
    209182#### def des variables 
    210183lon1_JUN=amch1_JUN[0,:] 
    … …  
    230203 
    231204 
    232 fichier.close 
     205    fichier.close() 
     206 
    233207 
    234208amch1_JUL=amch 
    235  
    236209#### def des variables  
    237210lon1_JUL=amch1_JUL[0,:] 
    … …  
    247220 
    248221 
     222 
     223 
     224 
    249225################ fichiers par frequence ################# 
    250226 

  • trunk/src/scripts_Laura/read_SSMIS_test.py

    r18 r19  
    1616f3 = '_ANTARC_JUNE2010.DAT' 
    1717#date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY']) 
    18 channel = np.array([17, 18, 15, 16]) 
     18channel = np.array([12, 13, 15, 16]) 
    1919numlines = np.zeros([len(channel)],int) 
    2020 
    … …  
    3232fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    3333ssmis = np.zeros([18, numlines[ich]], float) 
    34 for iligne in range (0,numlines[ich]-1): 
     34for iligne in range (0,numlines[ich]): 
    3535    line = fichier.readline() 
    3636    liste = line.split() 
    … …  
    3939 
    4040 
    41 fichier.close 
     41    fichier.close 
     42 
     43 
    4244ssch15_JUN=ssmis 
    4345lon15_JUN=ssch15_JUN[0,:] 
    … …  
    5658fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    5759ssmis = np.zeros([18, numlines[ich]], float) 
    58 for iligne in range (0,numlines[ich]-1): 
     60for iligne in range (0,numlines[ich]): 
    5961    line = fichier.readline() 
    6062    liste = line.split() 
    … …  
    6365 
    6466 
    65 fichier.close 
     67    fichier.close 
     68 
     69 
    6670ssch16_JUN=ssmis 
    6771lon16_JUN=ssch16_JUN[0,:] 
    … …  
    7781 
    7882 
    79 ich = 0 # 91.66GHz, V polar 
     83ich = 0 # 19.35GHz, H polar 
    8084fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    8185ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    8791 
    8892 
    89 fichier.close 
    90 ssch17_JUN=ssmis 
    91 lon17_JUN=ssch17_JUN[0,:] 
    92 lat17_JUN=ssch17_JUN[1,:] 
    93 jjr17_JUN=ssch17_JUN[4,:] 
    94 ts17_JUN=ssch17_JUN[8,:] 
    95 emis17_JUN=ssch17_JUN[14,:] 
    96 tb17_JUN=ssch17_JUN[13,:] 
    97 tup17_JUN=ssch17_JUN[16,:] 
    98 tdn17_JUN=ssch17_JUN[15,:] 
    99 trans17_JUN=ssch17_JUN[17,:] 
    100 orog17_JUN=ssch17_JUN[11,:] 
    101  
    102  
    103 ich = 1 # 91.66GHz, H polar 
     93    fichier.close 
     94 
     95 
     96ssch12_JUN=ssmis 
     97lon12_JUN=ssch12_JUN[0,:] 
     98lat12_JUN=ssch12_JUN[1,:] 
     99jjr12_JUN=ssch12_JUN[4,:] 
     100ts12_JUN=ssch12_JUN[8,:] 
     101emis12_JUN=ssch12_JUN[14,:] 
     102tb12_JUN=ssch12_JUN[13,:] 
     103tup12_JUN=ssch12_JUN[16,:] 
     104tdn12_JUN=ssch12_JUN[15,:] 
     105trans12_JUN=ssch12_JUN[17,:] 
     106orog12_JUN=ssch12_JUN[11,:] 
     107 
     108 
     109ich = 1 # 19.35GHz, V polar 
    104110fichier = open(f1 + str(channel[ich]) + f3, 'r') 
    105111ssmis = np.zeros([18, numlines[ich]], float) 
    … …  
    111117 
    112118 
    113 fichier.close 
    114 ssch18_JUN=ssmis 
    115 lon18_JUN=ssch18_JUN[0,:] 
    116 lat18_JUN=ssch18_JUN[1,:] 
    117 jjr18_JUN=ssch18_JUN[4,:] 
    118 ts18_JUN=ssch18_JUN[8,:] 
    119 emis18_JUN=ssch18_JUN[14,:] 
    120 tb18_JUN=ssch18_JUN[13,:] 
    121 tup18_JUN=ssch18_JUN[16,:] 
    122 tdn18_JUN=ssch18_JUN[15,:] 
    123 trans18_JUN=ssch18_JUN[17,:] 
    124 orog18_JUN=ssch18_JUN[11,:] 
     119    fichier.close 
     120 
     121 
     122ssch13_JUN=ssmis 
     123lon13_JUN=ssch13_JUN[0,:] 
     124lat13_JUN=ssch13_JUN[1,:] 
     125jjr13_JUN=ssch13_JUN[4,:] 
     126ts13_JUN=ssch13_JUN[8,:] 
     127emis13_JUN=ssch13_JUN[14,:] 
     128tb13_JUN=ssch13_JUN[13,:] 
     129tup13_JUN=ssch13_JUN[16,:] 
     130tdn13_JUN=ssch13_JUN[15,:] 
     131trans13_JUN=ssch13_JUN[17,:] 
     132orog13_JUN=ssch13_JUN[11,:] 
    125133 
    126134 
Note: See TracChangeset for help on using the changeset viewer.