More advanced example

From CrewWiki

(Difference between revisions)
Jump to: navigation, search
Line 33: Line 33:
</code><br>
</code><br>
[[File:group_image2.jpg|200px|overview cloud top pressure]]
[[File:group_image2.jpg|200px|overview cloud top pressure]]
 +
<br>
 +
 +
==Calculate multi algorithm ensemble average==
 +
<code>
 +
o=obj_new('cws_stats') <br>
 +
o->set_product,'ctt'<br>
 +
o->set_date,2008,6,13,12,00<br>
 +
avg = o->algorithm_average(/calc, /save, mask='com')
 +
o->set_group,'AVG'<br>
 +
o->make_image, /no_track<br>
 +
obj_destroy, o<br>
 +
</code><br>
 +
[[File:multi_algorithm_average.jpg|200px|multi algorithm average]]
<br>
<br>

Revision as of 09:52, 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

o=obj_new('cws_stats')
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, /no_track
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