[Fortran] 纯文本查看 复制代码
module hydromodelUtil
use iso_c_binding
implicit none
double precision, parameter :: pi=3.141593d0
!这部分为输入值
type ,Bind(C)::unitinfo
integer(c_int) :: controlno
character(c_char) :: unitname
integer(c_int) :: unittype
integer(c_int) :: solid
real (c_double):: coearea
real(c_double) ,dimension(3) :: shapepara
integer(c_int):: pn=0
double precision, dimension(:,:), allocatable :: plist
real(c_double), dimension(3) :: rotationpara, tranpara
end type unitinfo
!这部分为输出值
type ,Bind(C):: singleinfo
integer(c_int) :: controlno
integer(c_int) :: unittype
integer(c_int) :: solid
real(c_double) :: coearea
real(c_double):: volume
real(c_double), dimension(3) :: center
integer (c_int):: pointcount, trianglecount
double precision, dimension(:,:), allocatable :: pointchain
double precision, dimension(:,:,:), allocatable :: trianglechain
end type singleinfo
contains
subroutine cylindermeshing(r, h, cylinderinfo) Bind(C,Name="cylindermeshing")
!DEC$ ATTRIBUTES DLLEXPORT::cylindermeshing
implicit none
real(c_double) :: r, h
type(singleinfo) :: cylinderinfo
integer(c_int) :: polyedge=60
double precision, dimension(:,:), allocatable :: point
real(c_double) :: dtheta
integer(c_int) :: meshi
dtheta=2*pi/real(polyedge)
cylinderinfo%volume=(0.5d0*polyedge*sin(dtheta)*r**2)*h
cylinderinfo%center=(/ 0.0d0, 0.0d0, 0.5d0*h /)
cylinderinfo%pointcount=2*polyedge
cylinderinfo%trianglecount=4*polyedge
allocate(point(2,polyedge+1) )
allocate(cylinderinfo%pointchain(3,cylinderinfo%pointcount) )
do meshi=1, polyedge, 1
point(1,meshi)=r*cos( (real(meshi)-1.0d0)*dtheta)
point(2,meshi)=r*sin( (real(meshi)-1.0d0)*dtheta)
end do
point(:, polyedge+1)=point(:, 1)
do meshi=1, polyedge, 1
cylinderinfo%pointchain(1:2,meshi)=point(1:2,meshi)
cylinderinfo%pointchain(3,meshi)=0.0d0
cylinderinfo%pointchain(1:2,polyedge+meshi)=point(1:2,meshi)
cylinderinfo%pointchain(3,polyedge+meshi)=h
end do
end subroutine cylindermeshing
end module hydromodelUtil