@@ -101,6 +101,11 @@ def collapse_dimensions_for_plotting(longitude, latitude, vname, vx, vd, dims):
101
101
xticks_labels = None
102
102
xticks = np .arange (- 90 ,91 ,30 )
103
103
xlabel = 'Latitude (deg)'
104
+ if vd ['xaxis_type' ] == 'CFODD_NDBZE' :
105
+ x = np .arange (- 30 ,21 ,2 )
106
+ xticks = x
107
+ xticks_labels = None
108
+ xlabel = 'CloudSat equivalent reflectivity factor (dBZ)'
104
109
if vd ['yaxis_type' ] == 'pres7' :
105
110
yticks_labels = ('1000' , '800' , '680' , '560' , '440' , '310' , '180' ,'' )
106
111
yticks = y
@@ -119,6 +124,11 @@ def collapse_dimensions_for_plotting(longitude, latitude, vname, vx, vd, dims):
119
124
yticks_labels = ('0' ,'4' , '8' , '10' , '12.5' , '15' , '20' , '30' , '10000' )
120
125
yticks = y
121
126
ylabel = 'Liquid particle size (micron)'
127
+ if vd ['yaxis_type' ] == 'CFODD_NICOD' :
128
+ y = np .arange (0 ,61 ,2 )
129
+ yticks_labels = None
130
+ yticks = y
131
+ ylabel = 'MODIS in-cloud optical depth'
122
132
if vd ['yaxis_type' ] == 'levStat' :
123
133
if not dims ['levStat' ].any ():
124
134
# For diagnostics on model levels, all elements in levStat
@@ -298,7 +308,8 @@ def variable2D_metadata(var_list, fname):
298
308
('hgt16' , 'tau7' , 'loc' ),
299
309
('REICE_MODIS' , 'tau7' , 'loc' ),
300
310
('RELIQ_MODIS' , 'tau7' , 'loc' ),
301
- ('levStat' , 'SR_BINS' , 'loc' ))
311
+ ('levStat' , 'SR_BINS' , 'loc' ),
312
+ ('CFODD_NICOD' , 'CFODD_NDBZE' , 'loc' ))
302
313
zcs_dims = (('levStat' ,'loc' ), ('lev' ,'loc' ))
303
314
f_id = netCDF4 .Dataset (fname , 'r' )
304
315
vmeta = {}
@@ -352,11 +363,15 @@ def variable2D_metadata(var_list, fname):
352
363
'albisccp' , 'misr_meanztop' , 'misr_cldarea' , 'cltmodis' , 'clwmodis' , 'climodis' , 'clhmodis' ,
353
364
'clmmodis' , 'cllmodis' , 'tautmodis' , 'tauwmodis' , 'tauimodis' , 'tautlogmodis' , 'tauwlogmodis' ,
354
365
'tauilogmodis' , 'reffclwmodis' , 'reffclimodis' , 'pctmodis' , 'lwpmodis' , 'iwpmodis' ,
355
- 'cltlidarradar' , 'cloudsat_tcc' , 'cloudsat_tcc2' ,'parasolGrid_refl' ]
366
+ 'cltlidarradar' , 'cloudsat_tcc' , 'cloudsat_tcc2' ,'parasolGrid_refl' ,
367
+ 'ptcloudsatflag0' , 'ptcloudsatflag1' , 'ptcloudsatflag2' , 'ptcloudsatflag3' , 'ptcloudsatflag4' ,
368
+ 'ptcloudsatflag5' , 'ptcloudsatflag6' , 'ptcloudsatflag7' , 'ptcloudsatflag8' , 'ptcloudsatflag9' ,
369
+ 'cloudsatpia' , 'npdfcld' , 'npdfdrz' , 'npdfrain' ]
356
370
v2D_hists_names = ['clisccp' , 'clmodis' , 'cfadDbze94' , 'clMISR' ,
357
371
'modis_Optical_Thickness_vs_ReffICE' ,
358
372
'modis_Optical_Thickness_vs_ReffLIQ' ,
359
- 'cfadLidarsr532' , 'cfadLidarsr532gr' , 'cfadLidarsr355' ]
373
+ 'cfadLidarsr532' , 'cfadLidarsr532gr' , 'cfadLidarsr355' ,
374
+ 'ncfodd1' , 'ncfodd2' , 'ncfodd3' ]
360
375
v2D_zcs_names = ['clcalipsoice' ,'clcalipsoliq' ,'clcalipsoun' ,'clcalipsotmp' ,'clcalipsotmpice' ,'clcalipsotmpliq' ,
361
376
'clcalipsotmpun' ,'clcalipso' ,'clcalipsoopaque' ,'clcalipsothin' ,'clcalipsozopaque' ,
362
377
'clcalipsoopacity' ,'clgrLidar532' ,'clatlid' ,'clcalipso2' ,
0 commit comments