def __init__(self, emissivity_sfc=<NUM_LIT:0.>, **kwargs):
super(FourBandSW, self).__init__(emissivity_sfc=emissivity_sfc, **kwargs)<EOL>self.band_fraction = np.array([<NUM_LIT>, <NUM_LIT>, <NUM_LIT>, <NUM_LIT>])<EOL>if '<STR_LIT>' not in self.absorber_vmr:<EOL><INDENT>self.absorber_vmr['<STR_LIT>'] = <NUM_LIT> * np.ones_like(self.Tatm)<EOL><DEDENT>if '<STR_LIT>' not in self.absorber_vmr:<EOL><INDENT>self.absorber_vmr['<STR_LIT>'] = np.zeros_like(self.Tatm)<EOL><DEDENT>if '<STR_LIT>' not in self.absorber_vmr:<EOL><INDENT>self.absorber_vmr['<STR_LIT>'] = self.q<EOL><DEDENT>O3 = np.array([<NUM_LIT>, <NUM_LIT:0.>, <NUM_LIT>, <NUM_LIT:0.>]) * const.Rd / const.kBoltzmann<EOL>H2O = np.array([<NUM_LIT:0.>, <NUM_LIT:0.>, <NUM_LIT:0.>, <NUM_LIT>])<EOL>for n in range(self.Tatm.domain.numdims):<EOL><INDENT>H2O = H2O[:, np.newaxis]<EOL>O3 = O3[:, np.newaxis]<EOL><DEDENT>self.absorption_cross_section['<STR_LIT>'] = O3<EOL>self.absorption_cross_section['<STR_LIT>'] = H2O<EOL>self.absorption_cross_section['<STR_LIT>'] =np.zeros_like(self.absorption_cross_section['<STR_LIT>'])<EOL>self.cosZen = <NUM_LIT:0.5><EOL>
A four-band mdoel for shortwave radiation.

        The spectral decomposition used here is largely based on the
        "Moist Radiative-Convective Model" by Aarnout van Delden, Utrecht University
        a.j.vandelden@uu.nl
        http://www.staff.science.uu.nl/~delde102/RCM.htm

        Four SW channels:
            channel 0 is Hartley and Huggins band (UV, 6%, <340 nm)
            channel 1 is part of visible with no O3 absorption (14%, 340 - 500 nm)
            channel 2 is Chappuis band (27%, 500 - 700 nm)
            channel 3 is near-infrared (53%, > 700 nm)