tr= getParam( 'timerange', '2012-11-15', 'timerange to load' ) mag= getDataSet('vap+cdf:http://emfisis.physics.uiowa.edu/Flight/RBSP-B/L3/$Y/$m/$d/rbsp-b_magnetometer_4sec-gei_emfisis-L3_$Y$m$d_v$(v,sep).cdf?timerange=%s&Magnitude' % tr,monitor ) typ= getParam( 'typ', 'fce', 'the parameter to use', [ 'bmag', 'fce', '2fce', '0.5fce', 'fci', 'flhr' ] ) if ( typ=='fce' ): result= mag*28 result.putProperty( QDataSet.UNITS, Units.hertz ) result.putProperty( QDataSet.LABEL, 'f_ce' ) elif ( typ=='2fce'): result= mag*28*2 result.putProperty( QDataSet.UNITS, Units.hertz ) result.putProperty( QDataSet.LABEL, '2*f_ce' ) elif ( typ=='0.5fce'): result= mag*14 result.putProperty( QDataSet.UNITS, Units.hertz ) result.putProperty( QDataSet.LABEL, '0.5*f_ce' ) elif ( typ=='fci' ): result = (mag*28)/1837. result.putProperty( QDataSec.UNITS, Units.hertz ) result.putProperty( QDataSet.LABEL, 'f_ci' ) elif ( typ=='flhr' ): result = (mag*28)/42 result.putProperty( QDataSec.UNITS, Units.hertz ) result.putProperty( QDataSet.LABEL, 'f_lhr' ) elif ( typ=='bmag' ): result= mag elif ( typ=='fuh' ): print else: print 'huh'