dl On surface DA (reEer to Ficiure 3.24) HL = H& Now sinœ j, = ii x a' The wmnt density can be computeâ at any point on the cylinder.
3.4.4 Code and Resub Based on the formulations in sections 3.4.1, 3.4.2 and 3.4.3 code MGTDBX.FOR was created. The code covers angles of incidence 6 behneen O and 90'. It uses the generic UTD formulation for the di-on coefficient. Results from the moment method solution describeci in d o n 3.5 are plotted (in small cirdes) also for cornparison. The results frorn MGTDBX.FOR chosen to be reproduœâ hem, are foi a square cylinder illuminatd by a T G plane wave. The radius of the cylinder = a. The plots of the cuvent density am against pla. whem p is a linear distance on the surface of the cylinder from edge A moving dockwise. This means thus at pla=l. the position is halfiiiwy between A and B, if #a* the position is at C etc ... The idea behind this method of plotting is to compare the results obtained to plots already existent in the open literature. Figures 3.25, 3.26, 3.27 and 3.28 show the current density distribution on the surface of a square cylinder (of radius a) Mai ka= 10. 5.2 an 1 resmvely and the T G plane wave incident at 89.9'. These plots are in agreement mth pl]  m ub (Note:  corrects some of the plots found in ). Differenœs are notiœable at edges for the UTD solution but this is expcted so close to the edge. For a = 0. i5Qlk, the differenœ h mm noticeable on the shadow side with the UTD method. The differenœ be-n UT0 and moment method resuk (as a pemntage of the momnt method rwwlt) is still hmver below 5% at 0.05k from dge C. Figures 3.29 and 3.30 show the current density distribution with 60' angle of incidenœ for k-2 and 1 mspsctively.