create_grid.F90
1019 Bytes
if (conductor_data(conductor)%auto_grid_density) then
required_dl=conductor_data(conductor)%skin_depth ! could add reduction factor here e.g. /2.0
conductor_data(conductor)%dl=required_dl
end if
conductor_data(conductor)%nxmax=NINT(conductor_data(conductor)%grid_dim/conductor_data(conductor)%dl)+1
conductor_data(conductor)%nxmin=-conductor_data(conductor)%nxmax
conductor_data(conductor)%nymax=conductor_data(conductor)%nxmax
conductor_data(conductor)%nymin=conductor_data(conductor)%nxmin
nxmin=conductor_data(conductor)%nxmin
nxmax=conductor_data(conductor)%nxmax
nymin=conductor_data(conductor)%nymin
nymax=conductor_data(conductor)%nymax
ALLOCATE( conductor_data(conductor)%grid(nxmin:nxmax,nymin:nymax) )
ALLOCATE( conductor_data(conductor)%depth(nxmin:nxmax,nymin:nymax) )
! Reset the grid
conductor_data(conductor)%grid(nxmin:nxmax,nymin:nymax)=0
conductor_data(conductor)%depth(nxmin:nxmax,nymin:nymax)=0.0