UoN / SACAMOS

Blame view

SRC/WRITE_FH2_IPFILE/create_grid.F90 1.06 KB
Edit Raw Normal View History
197629dd   Chris Smartt   Update proximity_... 
1
2
3
4
5
 

!    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
ca8ab9e9   Chris Smartt   Update to proximi... 
6
7
8
9
10
11
12
13
14
15
16
 

    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
197629dd   Chris Smartt   Update proximity_... 
17
18
 
    
    write(*,*)'Create grid: nxmax=',nxmax,' nymax=',nymax
ca8ab9e9   Chris Smartt   Update to proximi... 
19
20
21
22
23
24
25
 
        
    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