source: trunk/src/scripts_Laura/read_SSMIS_test.py @ 23

Last change on this file since 23 was 23, checked in by lahlod, 10 years ago

modifs
File size: 13.8 KB
Line 
1#!/usr/bin/env python
2# -*- coding: utf-8 -*-
3import string
4import numpy as np
5import matplotlib.pyplot as plt
6from pylab import *
7from mpl_toolkits.basemap import Basemap
8from mpl_toolkits.basemap import shiftgrid, cm
9from netCDF4 import Dataset
10import ffgrid2
11
12
13
14################ fichiers par canal - mois de juin ###################################################
15
16
17f1 = '/net/dedale/usr/dedale/surf/lelod/ANTARC/SSMIS_CH'
18f3 = '_ANTARC_JUNE2010.DAT'
19channel = np.array([12, 13, 15, 16, 17, 18])
20numlines = np.zeros([len(channel)],int)
21
22for ich in range (0, len(channel)):
23     print channel[ich] 
24     f = f1 + str(channel[ich]) + f3
25     fichier = open(f, 'r')
26     numlines[ich] = 0
27     for line in fichier: numlines[ich] += 1
28
29      fichier.close()
30
31
32
33
34
35ich = 0 # 37GHz, H
36fichier = open(f1 + str(channel[ich]) + f3, 'r')
37ssmis = np.zeros([18, numlines[ich]], float)
38for iligne in range (0,numlines[ich]):
39    line = fichier.readline()
40    liste = line.split()
41    for j in range(0,18):
42        ssmis[j,iligne] = float(liste[j])
43
44
45    fichier.close
46
47
48ssch15_JUN=ssmis
49lon15_JUN=ssch15_JUN[0,:]
50lat15_JUN=ssch15_JUN[1,:]
51jjr15_JUN=ssch15_JUN[4,:]
52ts15_JUN=ssch15_JUN[8,:]
53emis15_JUN=ssch15_JUN[14,:]
54tb15_JUN=ssch15_JUN[13,:]
55tup15_JUN=ssch15_JUN[16,:]
56tdn15_JUN=ssch15_JUN[15,:]
57trans15_JUN=ssch15_JUN[17,:]
58orog15_JUN=ssch15_JUN[11,:]
59
60
61ich = 3 # 37GHz, V polar
62fichier = open(f1 + str(channel[ich]) + f3, 'r')
63ssmis = np.zeros([18, numlines[ich]], float)
64for iligne in range (0,numlines[ich]):
65    line = fichier.readline()
66    liste = line.split()
67    for j in range(0,18):
68        ssmis[j,iligne] = float(liste[j])
69
70
71    fichier.close
72
73
74ssch16_JUN=ssmis
75lon16_JUN=ssch16_JUN[0,:]
76lat16_JUN=ssch16_JUN[1,:]
77jjr16_JUN=ssch16_JUN[4,:]
78ts16_JUN=ssch16_JUN[8,:]
79emis16_JUN=ssch16_JUN[14,:]
80tb16_JUN=ssch16_JUN[13,:]
81tup16_JUN=ssch16_JUN[16,:]
82tdn16_JUN=ssch16_JUN[15,:]
83trans16_JUN=ssch16_JUN[17,:]
84orog16_JUN=ssch16_JUN[11,:]
85
86
87ich = 0 # 19.35GHz, H polar
88fichier = open(f1 + str(channel[ich]) + f3, 'r')
89ssmis = np.zeros([18, numlines[ich]], float)
90for iligne in range (0,numlines[ich]-1):
91    line = fichier.readline()
92    liste = line.split()
93    for j in range(0,18):
94        ssmis[j,iligne] = float(liste[j])
95
96
97    fichier.close
98
99
100ssch12_JUN=ssmis
101lon12_JUN=ssch12_JUN[0,:]
102lat12_JUN=ssch12_JUN[1,:]
103jjr12_JUN=ssch12_JUN[4,:]
104ts12_JUN=ssch12_JUN[8,:]
105emis12_JUN=ssch12_JUN[14,:]
106tb12_JUN=ssch12_JUN[13,:]
107tup12_JUN=ssch12_JUN[16,:]
108tdn12_JUN=ssch12_JUN[15,:]
109trans12_JUN=ssch12_JUN[17,:]
110orog12_JUN=ssch12_JUN[11,:]
111
112
113ich = 1 # 19.35GHz, V polar
114fichier = open(f1 + str(channel[ich]) + f3, 'r')
115ssmis = np.zeros([18, numlines[ich]], float)
116for iligne in range (0,numlines[ich]-1):
117    line = fichier.readline()
118    liste = line.split()
119    for j in range(0,18):
120        ssmis[j,iligne] = float(liste[j])
121
122
123    fichier.close
124
125
126ssch13_JUN=ssmis
127lon13_JUN=ssch13_JUN[0,:]
128lat13_JUN=ssch13_JUN[1,:]
129jjr13_JUN=ssch13_JUN[4,:]
130ts13_JUN=ssch13_JUN[8,:]
131emis13_JUN=ssch13_JUN[14,:]
132tb13_JUN=ssch13_JUN[13,:]
133tup13_JUN=ssch13_JUN[16,:]
134tdn13_JUN=ssch13_JUN[15,:]
135trans13_JUN=ssch13_JUN[17,:]
136orog13_JUN=ssch13_JUN[11,:]
137
138
139ich = 4 # 91.66GHz, V polar
140fichier = open(f1 + str(channel[ich]) + f3, 'r')
141ssmis = np.zeros([18, numlines[ich]], float)
142for iligne in range (0,numlines[ich]-1):
143    line = fichier.readline()
144    liste = line.split()
145    for j in range(0,18):
146        ssmis[j,iligne] = float(liste[j])
147
148
149    fichier.close
150
151
152ssch17_JUN=ssmis
153lon17_JUN=ssch17_JUN[0,:]
154lat17_JUN=ssch17_JUN[1,:]
155jjr17_JUN=ssch17_JUN[4,:]
156ts17_JUN=ssch17_JUN[8,:]
157emis17_JUN=ssch17_JUN[14,:]
158tb17_JUN=ssch17_JUN[13,:]
159tup17_JUN=ssch17_JUN[16,:]
160tdn17_JUN=ssch17_JUN[15,:]
161trans17_JUN=ssch17_JUN[17,:]
162orog17_JUN=ssch17_JUN[11,:]
163
164
165ich = 5 # 91.66GHz, H polar
166fichier = open(f1 + str(channel[ich]) + f3, 'r')
167ssmis = np.zeros([18, numlines[ich]], float)
168for iligne in range (0,numlines[ich]-1):
169    line = fichier.readline()
170    liste = line.split()
171    for j in range(0,18):
172        ssmis[j,iligne] = float(liste[j])
173
174
175    fichier.close
176
177
178ssch18_JUN=ssmis
179lon18_JUN=ssch18_JUN[0,:]
180lat18_JUN=ssch18_JUN[1,:]
181jjr18_JUN=ssch18_JUN[4,:]
182ts18_JUN=ssch18_JUN[8,:]
183emis18_JUN=ssch18_JUN[14,:]
184tb18_JUN=ssch18_JUN[13,:]
185tup18_JUN=ssch18_JUN[16,:]
186tdn18_JUN=ssch18_JUN[15,:]
187trans18_JUN=ssch18_JUN[17,:]
188orog18_JUN=ssch18_JUN[11,:]
189
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 TracBrowser for help on using the repository browser.