** radflux **
if(nx.gt.1) then
 store(qrx)
endif
if(ny.gt.1) then
 store(qry)
endif
if(nz.gt.1) then
 store(qrz)
endif