source: trunk/src/read_AMSUA_test.py @ 54

#!/usr/bin/env python
# -*- coding: utf-8 -*-
import string
import numpy as np
import matplotlib.pyplot as plt
from pylab import *
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.basemap import shiftgrid, cm
import netCDF4

################ fichiers par mois - CH1 ###################################################

f1='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH1_ANTARC_'
f2='2010.DAT'
date=np.array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])

numlines_month=np.zeros([7],float)

for i in range (0,7):
    f=f1+str(date[i])+f2
    fichier=open(f,'r')
    numlines_month[i]=0
    for line in fichier: numlines_month[i]+=1

    fichier.close()


i=0 #JANUARY
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_JAN=amch

#### enlever les valeurs a -500.000
#bb_ts_JAN=nonzero(amch_JAN[10,:]!=-500.)
#bb_emis_JAN=nonzero((amch_JAN[16,:]!=-500.)&(amch_JAN[16,:]<=1.))
#bb_tb_JAN=nonzero(amch_JAN[15,:]!=-500.)
#bb_tup_JAN=nonzero(amch_JAN[18,:]!=-500.)
#bb_tdn_JAN=nonzero(amch_JAN[17,:]!=-500.)
#bb_trans_JAN=nonzero(amch_JAN[19,:]!=-500.)
#bb_orog_JAN=nonzero(amch_JAN[13,:]>=0.)

#### def des variables 
lon1_JAN=amch1_JAN[0,:]
lat1_JAN=amch1_JAN[1,:]
jjr1_JAN=amch1_JAN[4,:]
ts1_JAN=amch1_JAN[10,:]
emis1_JAN=amch1_JAN[16,:]
tb1_JAN=amch1_JAN[15,:]
tup1_JAN=amch1_JAN[18,:]
tdn1_JAN=amch1_JAN[17,:]
trans1_JAN=amch1_JAN[19,:]
orog1_JAN=amch1_JAN[13,:]


i=1 #FEBRUARY
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_FEB=amch

#### enlever les valeurs a -500.000
#bb_ts_FEB=nonzero(amch_FEB[10,:]!=-500)
#bb_emis_FEB=nonzero((amch_FEB[16,:]!=-500)&(amch_FEB[16,:]<=1.))
#bb_tb_FEB=nonzero(amch_FEB[15,:]!=-500)
#bb_tup_FEB=nonzero(amch_FEB[18,:]!=-500)
#bb_tdn_FEB=nonzero(amch_FEB[17,:]!=-500)
#bb_trans_FEB=nonzero(amch_FEB[19,:]!=-500)
#bb_orog_FEB=nonzero(amch_FEB[13,:]>=0.)

#### def des variables CH1 corrigees
lon1_FEB=amch1_FEB[0,:]
lat1_FEB=amch1_FEB[1,:]
jjr1_FEB=amch1_FEB[4,:]
ts1_FEB=amch1_FEB[10,:]
emis1_FEB=amch1_FEB[16,:]
tb1_FEB=amch1_FEB[15,:]
tup1_FEB=amch1_FEB[18,:]
tdn1_FEB=amch1_FEB[17,:]
trans1_FEB=amch1_FEB[19,:]
orog1_FEB=amch1_FEB[13,:]


i=2 #MARCH
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_MAR=amch

#### enlever les valeurs a -500.000
#bb_ts_MAR=nonzero(amch_MAR[10,:]!=-500)
#bb_emis_MAR=nonzero((amch_MAR[16,:]!=-500)&(amch_MAR[16,:]<=1.))
#bb_tb_MAR=nonzero(amch_MAR[15,:]!=-500)
#bb_tup_MAR=nonzero(amch_MAR[18,:]!=-500)
#bb_tdn_MER=nonzero(amch_MAR[17,:]!=-500)
#bb_trans_MAR=nonzero(amch_MAR[19,:]!=-500)
#bb_orog_MAR=nonzero(amch_MAR[13,:]>=0.)

#### def des variables
lon1_MAR=amch1_MAR[0,:]
lat1_MAR=amch1_MAR[1,:]
jjr1_MAR=amch1_MAR[4,:]
ts1_MAR=amch1_MAR[10,:]
emis1_MAR=amch1_MAR[16,:]
tb1_MAR=amch1_MAR[15,:]
tup1_MAR=amch1_MAR[18,:]
tdn1_MAR=amch1_MAR[17,:]
trans1_MAR=amch1_MAR[19,:]
orog1_MAR=amch1_MAR[13,:]


i=3 #APRIL
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_APR=amch

#### def des variables 
lon1_APR=amch1_APR[0,:]
lat1_APR=amch1_APR[1,:]
jjr1_APR=amch1_APR[4,:]
ts1_APR=amch1_APR[10,:]
emis1_APR=amch1_APR[16,:]
tb1_APR=amch1_APR[15,:]
tup1_APR=amch1_APR[18,:]
tdn1_APR=amch1_APR[17,:]
trans1_APR=amch1_APR[19,:]
orog1_APR=amch1_APR[13,:]


i=4 #MAY
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_MAY=amch

#### def des variables
lon1_MAY=amch1_MAY[0,:]
lat1_MAY=amch1_MAY[1,:]
jjr1_MAY=amch1_MAY[4,:]
ts1_MAY=amch1_MAY[10,:]
emis1_MAY=amch1_MAY[16,:]
tb1_MAY=amch1_MAY[15,:]
tup1_MAY=amch1_MAY[18,:]
tdn1_MAY=amch1_MAY[17,:]
trans1_MAY=amch1_MAY[19,:]
orog1_MAY=amch1_MAY[13,:]

i=5 #JUNE
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_JUN=amch

#### def des variables
lon1_JUN=amch1_JUN[0,:]
lat1_JUN=amch1_JUN[1,:]
jjr1_JUN=amch1_JUN[4,:]
ts1_JUN=amch1_JUN[10,:]
emis1_JUN=amch1_JUN[16,:]
tb1_JUN=amch1_JUN[15,:]
tup1_JUN=amch1_JUN[18,:]
tdn1_JUN=amch1_JUN[17,:]
trans1_JUN=amch1_JUN[19,:]
orog1_JUN=amch1_JUN[13,:]

i=6 #JULLY
f=f1+str(date[i])+f2
fichier=open(f,'r')
amch=np.zeros([20,numlines_month[i]],float)
for iligne in range (0,numlines_month[i]):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


fichier.close

amch1_JUL=amch

#### def des variables 
lon1_JUL=amch1_JUL[0,:]
lat1_JUL=amch1_JUL[1,:]
jjr1_JUL=amch1_JUL[4,:]
ts1_JUL=amch1_JUL[10,:]
emis1_JUL=amch1_JUL[16,:]
tb1_JUL=amch1_JUL[15,:]
tup1_JUL=amch1_JUL[18,:]
tdn1_JUL=amch1_JUL[17,:]
trans1_JUL=amch1_JUL[19,:]
orog1_JUL=amch1_JUL[13,:]


################ fichiers par frequence #################


f11='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH'
f22='_ANTARC_JANUARY2010.DAT'
chan=np.array([1, 2, 3, 4, 5, 6, 7, 15])

numlines_channel=np.zeros([8],float)

for i in range (0,len(chan)):
    f=f11+str(chan[i])+f22
    fichier=open(f,'r')
    numlines_channel[i]=0
    for line in fichier: numlines_channel[i]+=1

    fichier.close()

i=1 #CHANNEL2
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch2_JAN=amch

#### def des variables 
lon2_JAN=amch2_JAN[0,:]
lat2_JAN=amch2_JAN[1,:]
jjr2_JAN=amch2_JAN[4,:]
ts2_JAN=amch2_JAN[10,:]
emis2_JAN=amch2_JAN[16,:]
tb2_JAN=amch2_JAN[15,:]
tup2_JAN=amch2_JAN[18,:]
tdn2_JAN=amch2_JAN[17,:]
trans2_JAN=amch2_JAN[19,:]
orog2_JAN=amch2_JAN[13,:]


i=2 #CHANNEL3
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch3_JAN=amch

#### def des variables 
lon3_JAN=amch3_JAN[0,:]
lat3_JAN=amch3_JAN[1,:]
jjr3_JAN=amch3_JAN[4,:]
ts3_JAN=amch3_JAN[10,:]
emis3_JAN=amch3_JAN[16,:]
tb3_JAN=amch3_JAN[15,:]
tup3_JAN=amch3_JAN[18,:]
tdn3_JAN=amch3_JAN[17,:]
trans3_JAN=amch3_JAN[19,:]
orog3_JAN=amch3_JAN[13,:]


i=3 #CHANNEL4
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch4_JAN=amch

#### def des variables 
lon4_JAN=amch4_JAN[0,:]
lat4_JAN=amch4_JAN[1,:]
jjr4_JAN=amch4_JAN[4,:]
ts4_JAN=amch4_JAN[10,:]
emis4_JAN=amch4_JAN[16,:]
tb4_JAN=amch4_JAN[15,:]
tup4_JAN=amch4_JAN[18,:]
tdn4_JAN=amch4_JAN[17,:]
trans4_JAN=amch4_JAN[19,:]
orog4_JAN=amch4_JAN[13,:]


################ fichiers par frequence - JULY #################


f11='/net/dedale/usr/dedale/surf/lelod/ANTARC/AMSUA_CH'
f22='_ANTARC_JULY2010.DAT'
chan=np.array([1, 2, 3, 4, 5, 6, 7, 15])

numlines_channel=np.zeros([8],float)

for i in range (0,len(chan)):
    f=f11+str(chan[i])+f22
    fichier=open(f,'r')
    numlines_channel[i]=0
    for line in fichier: numlines_channel[i]+=1

    fichier.close()

i=1 #CHANNEL2
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch2_JUL=amch

#### def des variables 
lon2_JUL=amch2_JUL[0,:]
lat2_JUL=amch2_JUL[1,:]
jjr2_JUL=amch2_JUL[4,:]
ts2_JUL=amch2_JUL[10,:]
emis2_JUL=amch2_JUL[16,:]
tb2_JUL=amch2_JUL[15,:]
tup2_JUL=amch2_JUL[18,:]
tdn2_JUL=amch2_JUL[17,:]
trans2_JUL=amch2_JUL[19,:]
orog2_JUL=amch2_JUL[13,:]


i=2 #CHANNEL3
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch3_JUL=amch

#### def des variables 
lon3_JUL=amch3_JUL[0,:]
lat3_JUL=amch3_JUL[1,:]
jjr3_JUL=amch3_JUL[4,:]
ts3_JUL=amch3_JUL[10,:]
emis3_JUL=amch3_JUL[16,:]
tb3_JUL=amch3_JUL[15,:]
tup3_JUL=amch3_JUL[18,:]
tdn3_JUL=amch3_JUL[17,:]
trans3_JUL=amch3_JUL[19,:]
orog3_JUL=amch3_JUL[13,:]


i=3 #CHANNEL4
f=f11+str(chan[i])+f22
fichier=open(f,'r')
amch=np.zeros([20,numlines_channel[i]],float)
for iligne in range (0,numlines_channel[i]-1):
    line=fichier.readline()
    liste = line.split()
    for j in range(0,20):
        amch[j,iligne] = float(liste[j])


    fichier.close

amch4_JUL=amch