Changeset 27
- Timestamp:
- 06/10/14 18:13:05 (10 years ago)
- Location:
- trunk/src/scripts_Laura
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/scripts_Laura/diff_frequ_AMSUA_test.py
r25 r27 27 27 imo = 0 28 28 ## ch1 ## 29 fovZen_ch1_JAN = np.where(fov1_JAN == 15.)29 fovZen_ch1_JAN = np.where(fov1_JAN == 21.) 30 30 bbemis_ch1_JAN = nonzero((emis1_JAN[fovZen_ch1_JAN] != -500.) & (emis1_JAN[fovZen_ch1_JAN] <= 1.)) 31 31 OUTZCH1_JAN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 34 34 latch1_JAN = np.zeros([len(np.arange(y0, y1+1, dy))], float) 35 35 ## ch15 ## 36 fovZen_ch15_JAN = np.where(fov15_JAN == 15.)36 fovZen_ch15_JAN = np.where(fov15_JAN == 21.) 37 37 bbemis_ch15_JAN = nonzero((emis15_JAN[fovZen_ch15_JAN] != -500.) & (emis15_JAN[fovZen_ch15_JAN] <= 1.)) 38 38 OUTZCH15_JAN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 73 73 imo = 1 74 74 ## ch1 ## 75 fovZen_ch1_FEB = np.where(fov1_FEB == 15.)75 fovZen_ch1_FEB = np.where(fov1_FEB == 21.) 76 76 bbemis_ch1_FEB = nonzero((emis1_FEB[fovZen_ch1_FEB] != -500.) & (emis1_FEB[fovZen_ch1_FEB] <= 1.)) 77 77 OUTZCH1_FEB = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 80 80 latch1_FEB = np.zeros([len(np.arange(y0, y1+1, dy))], float) 81 81 ## ch15 ## 82 fovZen_ch15_FEB = np.where(fov15_FEB == 15.)82 fovZen_ch15_FEB = np.where(fov15_FEB == 21.) 83 83 bbemis_ch15_FEB = nonzero((emis15_FEB[fovZen_ch15_FEB] != -500.) & (emis15_FEB[fovZen_ch15_FEB] <= 1.)) 84 84 OUTZCH15_FEB = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 120 120 imo = 2 121 121 ## ch1 ## 122 fovZen_ch1_MAR = np.where(fov1_MAR == 15.)122 fovZen_ch1_MAR = np.where(fov1_MAR == 21.) 123 123 bbemis_ch1_MAR = nonzero((emis1_MAR[fovZen_ch1_MAR] != -500.) & (emis1_MAR[fovZen_ch1_MAR] <= 1.)) 124 124 OUTZCH1_MAR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 127 127 latch1_MAR = np.zeros([len(np.arange(y0, y1+1, dy))], float) 128 128 ## ch15 ## 129 fovZen_ch15_MAR = np.where(fov15_MAR == 15.)129 fovZen_ch15_MAR = np.where(fov15_MAR == 21.) 130 130 bbemis_ch15_MAR = nonzero((emis15_MAR[fovZen_ch15_MAR] != -500.) & (emis15_MAR[fovZen_ch15_MAR] <= 1.)) 131 131 OUTZCH15_MAR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 166 166 imo = 3 167 167 ## ch1 ## 168 fovZen_ch1_APR = np.where(fov1_APR == 15.)168 fovZen_ch1_APR = np.where(fov1_APR == 21.) 169 169 bbemis_ch1_APR = nonzero((emis1_APR[fovZen_ch1_APR] != -500.) & (emis1_APR[fovZen_ch1_APR] <= 1.)) 170 170 OUTZCH1_APR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 173 173 latch1_APR = np.zeros([len(np.arange(y0, y1+1, dy))], float) 174 174 ## ch15 ## 175 fovZen_ch15_APR = np.where(fov15_APR == 15.)175 fovZen_ch15_APR = np.where(fov15_APR == 21.) 176 176 bbemis_ch15_APR = nonzero((emis15_APR[fovZen_ch15_APR] != -500.) & (emis15_APR[fovZen_ch15_APR] <= 1.)) 177 177 OUTZCH15_APR = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 212 212 imo = 4 213 213 ## ch1 ## 214 fovZen_ch1_MAY = np.where(fov1_MAY == 15.)214 fovZen_ch1_MAY = np.where(fov1_MAY == 21.) 215 215 bbemis_ch1_MAY = nonzero((emis1_MAY[fovZen_ch1_MAY] != -500.) & (emis1_MAY[fovZen_ch1_MAY] <= 1.)) 216 216 OUTZCH1_MAY = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 219 219 latch1_MAY = np.zeros([len(np.arange(y0, y1+1, dy))], float) 220 220 ## ch15 ## 221 fovZen_ch15_MAY = np.where(fov15_MAY == 15.)221 fovZen_ch15_MAY = np.where(fov15_MAY == 21.) 222 222 bbemis_ch15_MAY = nonzero((emis15_MAY[fovZen_ch15_MAY] != -500.) & (emis15_MAY[fovZen_ch15_MAY] <= 1.)) 223 223 OUTZCH15_MAY = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 258 258 imo = 5 259 259 ## ch1 ## 260 fovZen_ch1_JUN = np.where(fov1_JUN == 15.)260 fovZen_ch1_JUN = np.where(fov1_JUN == 21.) 261 261 bbemis_ch1_JUN = nonzero((emis1_JUN[fovZen_ch1_JUN] != -500.) & (emis1_JUN[fovZen_ch1_JUN] <= 1.)) 262 262 OUTZCH1_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 265 265 latch1_JUN = np.zeros([len(np.arange(y0, y1+1, dy))], float) 266 266 ## ch15 ## 267 fovZen_ch15_JUN = np.where(fov15_JUN == 15.)267 fovZen_ch15_JUN = np.where(fov15_JUN == 21.) 268 268 bbemis_ch15_JUN = nonzero((emis15_JUN[fovZen_ch15_JUN] != -500.) & (emis15_JUN[fovZen_ch15_JUN] <= 1.)) 269 269 OUTZCH15_JUN = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 304 304 imo = 6 305 305 ## ch1 ## 306 fovZen_ch1_JUL = np.where(fov1_JUL == 15.)306 fovZen_ch1_JUL = np.where(fov1_JUL == 21.) 307 307 bbemis_ch1_JUL = nonzero((emis1_JUL[fovZen_ch1_JUL] != -500.) & (emis1_JUL[fovZen_ch1_JUL] <= 1.)) 308 308 OUTZCH1_JUL = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 311 311 latch1_JUL = np.zeros([len(np.arange(y0, y1+1, dy))], float) 312 312 ## ch15 ## 313 fovZen_ch15_JUL = np.where(fov15_JUL == 15.)313 fovZen_ch15_JUL = np.where(fov15_JUL == 21.) 314 314 bbemis_ch15_JUL = nonzero((emis15_JUL[fovZen_ch15_JUL] != -500.) & (emis15_JUL[fovZen_ch15_JUL] <= 1.)) 315 315 OUTZCH15_JUL = np.zeros([len(np.arange(y0, y1+1, dy)),len(np.arange(x0, x1+1, dx)),len_month[imo]], float) … … 473 473 m.drawmeridians(np.arange(-180., 180., 20)) 474 474 xii,yii = m(*np.meshgrid(lon, lat)) 475 #clevs = (arange(0., 0.006, 0.00001))476 cs = m.contourf(xii, yii, monthly_emis_JUL[0, :, :], cmap=cm.s3pcpn_l_r)475 clevs = (arange(0.5, 0.9, 0.001)) 476 cs = m.contourf(xii, yii, OUTZ_JAN[0, :, :,20], clevs, cmap=cm.s3pcpn_l_r) 477 477 cbar = colorbar(cs) 478 478 cbar.set_label('CH1 - 23.8GHz') -
trunk/src/scripts_Laura/write_netcdf_AMSUA_test.py
r25 r27 46 46 ## JANUARY ## 47 47 imo = 0 48 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')48 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 49 49 rootgrp.createDimension('longitude', len(lon)) 50 50 rootgrp.createDimension('latitude', len(lat)) … … 67 67 ## FEBRUARY ## 68 68 imo = 1 69 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')69 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 70 70 rootgrp.createDimension('longitude', len(lon)) 71 71 rootgrp.createDimension('latitude', len(lat)) … … 88 88 ## MARCH ## 89 89 imo = 2 90 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')90 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 91 91 rootgrp.createDimension('longitude', len(lon)) 92 92 rootgrp.createDimension('latitude', len(lat)) … … 109 109 ## APRIL ## 110 110 imo = 3 111 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')111 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 112 112 rootgrp.createDimension('longitude', len(lon)) 113 113 rootgrp.createDimension('latitude', len(lat)) … … 130 130 ## MAY ## 131 131 imo = 4 132 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')132 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 133 133 rootgrp.createDimension('longitude', len(lon)) 134 134 rootgrp.createDimension('latitude', len(lat)) … … 151 151 ## JUNE ## 152 152 imo = 5 153 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')153 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 154 154 rootgrp.createDimension('longitude', len(lon)) 155 155 rootgrp.createDimension('latitude', len(lat)) … … 172 172 ## JUlY## 173 173 imo = 6 174 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_ ANTARC.nc', 'w', format = 'NETCDF4')174 rootgrp = Dataset ('/net/argos/data/parvati/lahlod/AMUSUA_CH1-CH15_' + month[imo] +'_fov21_ANTARC.nc', 'w', format = 'NETCDF4') 175 175 rootgrp.createDimension('longitude', len(lon)) 176 176 rootgrp.createDimension('latitude', len(lat)) … … 193 193 194 194 195 196 197
Note: See TracChangeset
for help on using the changeset viewer.