source: trunk/src/scripts_Laura/write_netcdf_AMSUA_test.py @ 54

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

modifs
File size: 7.1 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
10
11
12############################
13## monthly mean emisivity ##
14############################
15month = array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
16#rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_FEB-APR-JUL_ANTARC.nc', 'w', format = 'NETCDF4')
17rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_JAN-MAR-MAY-JUN_ANTARC.nc', 'w', format = 'NETCDF4')
18rootgrp.createDimension('longitude', len(lon))
19rootgrp.createDimension('latitude', len(lat))
20rootgrp.createDimension('channels', 2)
21rootgrp.createDimension('month', 4)
22nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
23nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
24ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
25ncmon = rootgrp.createVariable('month', 'f8', ('month',))
26ncemis = rootgrp.createVariable('emissivity', 'f8', ('month', 'channels', 'latitude', 'longitude'))
27nclon[:] =lon
28nclat[:] = lat
29ncchan[:] = [23.8, 89]
30ncmon[:] = [1., 2., 3., 4.]
31ncemis[:] = mean_outz2
32
33rootgrp.close()
34
35
36
37#####################
38## daily emisivity ##
39#####################
40month = array(['JANUARY', 'FEBRUARY', 'MARCH', 'APRIL', 'MAY', 'JUNE', 'JULY'])
41len_month = np.array([31, 28, 31, 30, 31, 30, 31])
42chan = array(['CH1', 'CH2'])
43lon = lonch1_JAN
44lat = latch1_JAN
45
46## JANUARY ##
47imo = 0
48rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
49rootgrp.createDimension('longitude', len(lon))
50rootgrp.createDimension('latitude', len(lat))
51rootgrp.createDimension('channels', len(chan))
52rootgrp.createDimension('day', len_month[imo])
53nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
54nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
55ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
56ncjours = rootgrp.createVariable('day', 'f8', ('day',))
57ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
58nclon[:] =lon
59nclat[:] = lat
60ncchan[:] = [23.8, 89]
61ncjours[:] = np.arange(0, len_month[imo], 1)
62ncemis[:] = OUTZ_JAN
63
64rootgrp.close()
65
66
67## FEBRUARY ##
68imo = 1
69rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
70rootgrp.createDimension('longitude', len(lon))
71rootgrp.createDimension('latitude', len(lat))
72rootgrp.createDimension('channels', len(chan))
73rootgrp.createDimension('day', len_month[imo])
74nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
75nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
76ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
77ncjours = rootgrp.createVariable('day', 'f8', ('day',))
78ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
79nclon[:] =lon
80nclat[:] = lat
81ncchan[:] = [23.8, 89]
82ncjours[:] = np.arange(0, len_month[imo], 1)
83ncemis[:] = OUTZ_FEB
84
85rootgrp.close()
86
87
88## MARCH ##
89imo = 2
90rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
91rootgrp.createDimension('longitude', len(lon))
92rootgrp.createDimension('latitude', len(lat))
93rootgrp.createDimension('channels', len(chan))
94rootgrp.createDimension('day', len_month[imo])
95nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
96nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
97ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
98ncjours = rootgrp.createVariable('day', 'f8', ('day',))
99ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
100nclon[:] =lon
101nclat[:] = lat
102ncchan[:] = [23.8, 89]
103ncjours[:] = np.arange(0, len_month[imo], 1)
104ncemis[:] = OUTZ_MAR
105
106rootgrp.close()
107
108
109## APRIL ##
110imo = 3
111rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
112rootgrp.createDimension('longitude', len(lon))
113rootgrp.createDimension('latitude', len(lat))
114rootgrp.createDimension('channels', len(chan))
115rootgrp.createDimension('day', len_month[imo])
116nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
117nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
118ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
119ncjours = rootgrp.createVariable('day', 'f8', ('day',))
120ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
121nclon[:] =lon
122nclat[:] = lat
123ncchan[:] = [23.8, 89]
124ncjours[:] = np.arange(0, len_month[imo], 1)
125ncemis[:] = OUTZ_APR
126
127rootgrp.close()
128
129
130## MAY ##
131imo = 4
132rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
133rootgrp.createDimension('longitude', len(lon))
134rootgrp.createDimension('latitude', len(lat))
135rootgrp.createDimension('channels', len(chan))
136rootgrp.createDimension('day', len_month[imo])
137nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
138nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
139ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
140ncjours = rootgrp.createVariable('day', 'f8', ('day',))
141ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
142nclon[:] =lon
143nclat[:] = lat
144ncchan[:] = [23.8, 89]
145ncjours[:] = np.arange(0, len_month[imo], 1)
146ncemis[:] = OUTZ_MAY
147
148rootgrp.close()
149
150
151## JUNE ##
152imo = 5
153rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
154rootgrp.createDimension('longitude', len(lon))
155rootgrp.createDimension('latitude', len(lat))
156rootgrp.createDimension('channels', len(chan))
157rootgrp.createDimension('day', len_month[imo])
158nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
159nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
160ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
161ncjours = rootgrp.createVariable('day', 'f8', ('day',))
162ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
163nclon[:] =lon
164nclat[:] = lat
165ncchan[:] = [23.8, 89]
166ncjours[:] = np.arange(0, len_month[imo], 1)
167ncemis[:] = OUTZ_JUN
168
169rootgrp.close()
170
171
172## JUlY##
173imo = 6
174rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4')
175rootgrp.createDimension('longitude', len(lon))
176rootgrp.createDimension('latitude', len(lat))
177rootgrp.createDimension('channels', len(chan))
178rootgrp.createDimension('day', len_month[imo])
179nclon = rootgrp.createVariable('longitude', 'f8', ('longitude',))
180nclat = rootgrp.createVariable('latitude', 'f8', ('latitude',))
181ncchan = rootgrp.createVariable('channels', 'f8', ('channels',))
182ncjours = rootgrp.createVariable('day', 'f8', ('day',))
183ncemis = rootgrp.createVariable('emissivity', 'f8', ('channels', 'latitude', 'longitude', 'day'))
184nclon[:] =lon
185nclat[:] = lat
186ncchan[:] = [23.8, 89]
187ncjours[:] = np.arange(0, len_month[imo], 1)
188ncemis[:] = OUTZ_JUL
189
190rootgrp.close()
191
192
193
194
Note: See TracBrowser for help on using the repository browser.