!---------------------------------------------------------------------- ! Copyright (c) 2010--2011 Shin-ichi Takehiro. All rights reserved. !---------------------------------------------------------------------- ! !表題 wtq_mpi_module テストプログラム ! ! 微分のテスト ! !履歴 2010/04/18 竹広真一 wt_module_sjpack_deriv_test1.f90 より改造 ! 2011/09/13 竹広真一 MPI 並列化 ! 2011/09/14 竹広真一 wtq_mpi_module 用に改造 ! program wtq_mpi_module_derivative_wt1 use dc_message, only : MessageNotify use dc_test, only : AssertEqual use wtq_mpi_module use mpi implicit none integer,parameter :: im=32, jm=16 ! 格子点の設定(経度, 緯度, 動径) integer,parameter :: kmo=16, kmi=10 ! 格子点の設定(球殻動径, 球動径) integer,parameter :: nm=10 ! 切断波数の設定(水平) integer,parameter :: lmo=10, lmi=16 ! 切断波数の設定(球殻動径, 球動径) real(8),parameter :: ri=0.5, ro=1.5 ! 内外半径 real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_Data real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_DData real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_VLon real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_VLat real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_VRad real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_Div real(8), dimension(0:im-1,1:jm,0:kmo) :: xyz_Psi real(8), allocatable :: xvz_Data(:,:,:) real(8), allocatable :: xvz_DData(:,:,:) real(8), allocatable :: xvz_VLon(:,:,:) real(8), allocatable :: xvz_VLat(:,:,:) real(8), allocatable :: xvz_VRad(:,:,:) real(8), allocatable :: xvz_Div(:,:,:) real(8), allocatable :: xvz_Psi(:,:,:) ! 判定誤差設定 integer, parameter :: check_digits = 10 integer, parameter :: ignore = -11 integer, parameter :: nmin=1, nmax=9 integer :: n integer :: iproc, np, ierr !---------------- MPI スタート --------------------- call MPI_INIT(IERR) call MPI_COMM_RANK(MPI_COMM_WORLD,IPROC,IERR) call MPI_COMM_SIZE(MPI_COMM_WORLD,NP,IERR) call MessageNotify('M','wtq_mpi_module_derivative_wt_test1', & 'wtq_mpi_module wt-derivative function test #1') call wtq_mpi_Initial(im,jm,kmi,kmo,nm,lmi,lmo,ri,ro) allocate(xvz_Data(0:im-1,1:jc,0:kmo)) allocate(xvz_DData(0:im-1,1:jc,0:kmo)) allocate(xvz_VLon(0:im-1,1:jc,0:kmo)) allocate(xvz_VLat(0:im-1,1:jc,0:kmo)) allocate(xvz_VRad(0:im-1,1:jc,0:kmo)) allocate(xvz_Div(0:im-1,1:jc,0:kmo)) allocate(xvz_Psi(0:im-1,1:jc,0:kmo)) !==================== wt_DivRad_wt ==================== do n=nmin,nmax write(6,*) 'n=',n xyz_Data = xyz_Rad**n xyz_DData = (n+2)*xyz_Rad**(n-1) xvz_Data = xvz_Rad**n xvz_DData = (n+2)*xvz_Rad**(n-1) call AssertEqual(& message='wt_Divrat_wt with r^n', & answer = xyz_DData, & check = xyz_wt(wt_DivRad_wt(wt_xyz(xyz_Data))), & significant_digits = check_digits, ignore_digits = ignore & ) call AssertEqual(& message='wt_Divrat_wt with r^n', & answer = xvz_DData, & check = xvz_wt(wt_DivRad_wt(wt_xvz(xvz_Data))), & significant_digits = check_digits, ignore_digits = ignore & ) enddo !========== xyz_GradLon_wt, xyz_GradLat_wt, wt_Div_xyz_xyz_xyz ========== ! ------ r^2 Y_1^1 Toridal field ----- n = 2 write(6,*) 'n=',n xyz_Psi = xyz_Rad**n * cos(xyz_Lat)*sin(xyz_Lon) ! r^2 Y_1^1 xyz_VRad = 0 xvz_Psi = xvz_Rad**n * cos(xvz_Lat)*sin(xvz_Lon) ! r^2 Y_1^1 xvz_VRad = 0 xyz_VLon = xyz_Rad**(n-1)*cos(xyz_Lon) call AssertEqual(& message='xyz_GradLon_wt with r^2 Y_1^-1', & answer = xyz_VLon, & check = xyz_GradLon_wt(wt_xyz(xyz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLon = xvz_Rad**(n-1)*cos(xvz_Lon) call AssertEqual(& message='xvz_GradLon_wt with r^2 Y_1^-1', & answer = xvz_VLon, & check = xvz_GradLon_wt(wt_xvz(xvz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xyz_VLat = -xyz_Rad**(n-1)*sin(xyz_Lat)*sin(xyz_Lon) call AssertEqual(& message='xyz_GradLat_wt with r^2 Y_1^-1', & answer = xyz_VLat, & check = xyz_GradLat_wt(wt_xyz(xyz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLat = -xvz_Rad**(n-1)*sin(xvz_Lat)*sin(xvz_Lon) call AssertEqual(& message='xvz_GradLat_wt with r^2 Y_1^-1', & answer = xvz_VLat, & check = xvz_GradLat_wt(wt_xvz(xvz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xyz_Div = - 2 * xyz_Psi/xyz_Rad**2 call AssertEqual(& message='wt_Div_xyz_xyz_xyz_wt with Toroiral r^2 Y_1^-1', & answer = xyz_Div, & check = xyz_wt(wt_Div_xyz_xyz_xyz(xyz_VLon,xyz_VLat,xyz_VRad)),& significant_digits = check_digits, ignore_digits = ignore & ) xvz_Div = - 2 * xvz_Psi/xvz_Rad**2 call AssertEqual(& message='wt_Div_xvz_xvz_xvz_wt with Toroiral r^2 Y_1^-1', & answer = xvz_Div, & check = xvz_wt(wt_Div_xvz_xvz_xvz(xvz_VLon,xvz_VLat,xvz_VRad)),& significant_digits = check_digits, ignore_digits = ignore & ) ! ------ r^2Y_2^-1 Toridal field ----- n = 2 xyz_VRad = 0 xyz_Psi = xyz_Rad**n * cos(xyz_Lat)*sin(xyz_Lat) * sin(xyz_Lon) ! r^2Y_2^-1 xvz_VRad = 0 xvz_Psi = xvz_Rad**n * cos(xvz_Lat)*sin(xvz_Lat) * sin(xvz_Lon) ! r^2Y_2^-1 xyz_VLon = xyz_Rad**(n-1)*sin(xyz_Lat)*cos(xyz_Lon) call AssertEqual(& message='xyz_GradLon_wt with r^2 Y_2^-1', & answer = xyz_VLon, & check = xyz_GradLon_wt(wt_xyz(xyz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLon = xvz_Rad**(n-1)*sin(xvz_Lat)*cos(xvz_Lon) call AssertEqual(& message='xvz_GradLon_wt with r^2 Y_2^-1', & answer = xvz_VLon, & check = xvz_GradLon_wt(wt_xvz(xvz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xyz_VLat = xyz_Rad**(n-1)*cos(2*xyz_Lat)*sin(xyz_Lon) call AssertEqual(& message='xyz_GradLat_wt with r^2 Y_2^-1', & answer = xyz_VLat, & check = xyz_GradLat_wt(wt_xyz(xyz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLat = xvz_Rad**(n-1)*cos(2*xvz_Lat)*sin(xvz_Lon) call AssertEqual(& message='xvz_GradLat_wt with r^2 Y_2^-1', & answer = xvz_VLat, & check = xvz_GradLat_wt(wt_xvz(xvz_Psi)), & significant_digits = check_digits, ignore_digits = ignore & ) xyz_Div = - 6* xyz_Psi/xyz_Rad**2 call AssertEqual(& message='wt_Div_xyz_xyz_xyz_wt with Toroiral r2 Y_1^-1', & answer = xyz_Div, & check = xyz_wt(wt_Div_xyz_xyz_xyz(xyz_VLon,xyz_VLat,xyz_VRad)),& significant_digits = check_digits, ignore_digits = ignore & ) xvz_Div = - 6* xvz_Psi/xvz_Rad**2 call AssertEqual(& message='wt_Div_xvz_xvz_xvz_wt with Toroiral r2 Y_1^-1', & answer = xvz_Div, & check = xvz_wt(wt_Div_xvz_xvz_xvz(xvz_VLon,xvz_VLat,xvz_VRad)),& significant_digits = check_digits, ignore_digits = ignore & ) !==================== wt_DivLon_xyz, wt_DivLat_xyz ==================== n=2 write(6,*)'n =',n xyz_VLon = xyz_rad**n * cos(xyz_Lat)**2*sin(xyz_Lon) xyz_DData = xyz_rad**(n-1)*cos(xyz_Lat)*cos(xyz_Lon) call AssertEqual(& message='wt_DivLon_xyz with r^n cos^2 φ sinλ', & answer = xyz_DData, & check = xyz_wt(wt_DivLon_xyz(xyz_VLon)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLon = xvz_rad**n * cos(xvz_Lat)**2*sin(xvz_Lon) xvz_DData = xvz_rad**(n-1)*cos(xvz_Lat)*cos(xvz_Lon) call AssertEqual(& message='wt_DivLon_xvz with r^n cos^2 φ sinλ', & answer = xvz_DData, & check = xvz_wt(wt_DivLon_xvz(xvz_VLon)), & significant_digits = check_digits, ignore_digits = ignore & ) xyz_VLat = xyz_rad**n * cos(xyz_Lat)**2*sin(xyz_Lon) xyz_DData = -3*xyz_rad**(n-1)*cos(xyz_Lat)*sin(xyz_Lat)*sin(xyz_Lon) call AssertEqual(& message='wt_DivLat_xyz with r^n cos^2 φ sinλ', & answer = xyz_DData, & check = xyz_wt(wt_DivLat_xyz(xyz_VLat)), & significant_digits = check_digits, ignore_digits = ignore & ) xvz_VLat = xvz_rad**n * cos(xvz_Lat)**2*sin(xvz_Lon) xvz_DData = -3*xvz_rad**(n-1)*cos(xvz_Lat)*sin(xvz_Lat)*sin(xvz_Lon) call AssertEqual(& message='wt_DivLat_xvz with r^n cos^2 φ sinλ', & answer = xvz_DData, & check = xvz_wt(wt_DivLat_xvz(xvz_VLat)), & significant_digits = check_digits, ignore_digits = ignore & ) call MessageNotify('M','wtq_mpi_module_derivative_wt_test1', & 'wtq_mpi_module wt-derivative function test #1 succeeded!') !------ MPIの終了 ------ call MPI_FINALIZE(IERR) end program wtq_mpi_module_derivative_wt1