More advanced example

From CrewWiki

(Difference between revisions)
Jump to: navigation, search
(Calculate multi algorithm ensemble average)
(Calculate multi algorithm ensemble average)
Line 36: Line 36:
==Calculate multi algorithm ensemble average==
==Calculate multi algorithm ensemble average==
 +
Note, if you want to have
<code>
<code>
o=obj_new('cws_stats', scale=1) <br>
o=obj_new('cws_stats', scale=1) <br>
Line 42: Line 43:
avg = o->algorithm_average(/calc, /save, mask='com')<br>
avg = o->algorithm_average(/calc, /save, mask='com')<br>
o->set_group,'AVG'<br>
o->set_group,'AVG'<br>
-
o->make_image, /single, /no_track<br>
+
o->make_image, /single<br>
 +
o->set_group,'N_OBS'<br>
 +
o->make_image, /single<br>
obj_destroy, o<br>
obj_destroy, o<br>
</code><br>
</code><br>

Revision as of 10:14, 21 February 2014

Go back to CREW software package

Contents

Display a cloud property of one algorithm on the MSG disk

o=obj_new('cws_stats')
o->set_product,'cth'
o->set_group,'CMS'
o->set_date,2008,6,13,12,15
o->make_image, /single
obj_destroy, o

CMS cloud top height

Difference of one data set to another (one picture per algorithm pair)

o=obj_new('cws_stats')
o->set_product,'ctt'
o->set_date,2008,6,13,12,00
o->diff_maps, group='AWG'
obj_destroy, o

cloud top temperature difference

Display a cloud property of all algorithms on the MSG disk

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->make_image, /no_track
obj_destroy, o

overview cloud top pressure

Calculate multi algorithm ensemble average

Note, if you want to have o=obj_new('cws_stats', scale=1)
o->set_product,'ctt'
o->set_date,2008,6,13,12,00
avg = o->algorithm_average(/calc, /save, mask='com')
o->set_group,'AVG'
o->make_image, /single
o->set_group,'N_OBS'
o->make_image, /single
obj_destroy, o

multi algorithm average

1d histogram

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->hist1d
obj_destroy, o

cloud top pressure histogram

dependency on latitude

o=obj_new('cws_stats')
o->set_product,'cth'
o->set_date,2008,6,13,12,00
o->dep
obj_destroy, o

latitudinal mean of cloud top temperature

dependency on viewing zenith angle

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->dep, /vza
obj_destroy, o

cloud top pressure per viewing zenith angle

Simple statistics as a table

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->make_tables
obj_destroy, o

cloud top pressure statistics

Two algorithm scatter plot

o=obj_new('cws_stats')
o->set_product,'ctt'
o->set_date,2008,6,13,12,00
o->make_corSingle
obj_destroy, o

cloud top temperature scatter plots

Thumbnails of all scatter plots

o=obj_new('cws_stats')
o->set_product,'ctt'
o->set_date,2008,6,13,12,00
o->make_corThumb
obj_destroy, o

overview cloud top temperature correlation plots

Make ISCCP like plots

o=obj_new('cws_stats')
o->cod_ctp
obj_destroy, o

ISCCP like CTP COD histogram