AtHeartEngineer/kaiju
1
0
Fork 0
mirror of https://github.com/JHUAPL/kaiju.git synced 2026-01-09 15:17:56 -05:00
Code Issues Packages Projects Releases Wiki Activity

Correcting dipole test code in ebsquish, and fixing non-mpi dipole squish test

Browse Source
This commit is contained in:
Jeffrey Garretson
2024-09-09 22:00:19 -06:00
parent 4d576122f2
commit 4ace1f469a
2 changed files with 20 additions and 19 deletions
Download Patch File Download Diff File Expand all files Collapse all files

24
tests/voltron/testebsquish.pf
View File

@@ -153,18 +153,20 @@ contains
call Squish(voltronApp)
call SquishEnd(voltronApp)
@assertEqual(COUNT(voltronApp%chmp2mhd%isGood),112344,'Squish Fake Projection Good Count is wrong. Check Squish Processing and Output.')
@assertEqual(COUNT(norm2(ebGr%xyz(ebGr%is:voltronApp%iDeep+1,ebGr%js:ebGr%je+1,ebGr%ks:ebGr%ke+1,:),4) <= voltronApp%rTrc),COUNT(voltronApp%chmp2mhd%isGood),'Squish Fake Projection Good Count is wrong. Check Squish Processing and Output.')
do k=ebGr%ks,ebGr%ke
do j=ebGr%js,ebGr%je
do i=ebGr%is,ebGr%ie
xyz = ebGr%xyz(i,j,k,XDIR:ZDIR)
invlat = InvLatitude(xyz)
invlon = katan2(xyz(YDIR),xyz(XDIR))
write (checkMessage,'(A,I0,A,I0,A,I0,A)'), 'Squish Fake Projection Latitude is wrong at (',i,',',j,',',k,'). Check Squish Processing and Output.'
@assertEqual(voltronApp%chmp2mhd%xyzSquish(i,j,k,1),invlat,1e-17_rp,trim(checkMessage))
write (checkMessage,'(A,I0,A,I0,A,I0,A)'), 'Squish Fake Projection Longitude is wrong at (',i,',',j,',',k,'). Check Squish Processing and Output.'
@assertEqual(voltronApp%chmp2mhd%xyzSquish(i,j,k,2),invlon,1e-17_rp,trim(checkMessage))
do k=ebGr%ks,ebGr%ke+1
do j=ebGr%js+1,ebGr%je ! js and je+1 are smoothed singularities, not checked
do i=ebGr%is,voltronApp%iDeep+1
if(voltronApp%chmp2mhd%isGood(i,j,k)) then
xyz = ebGr%xyz(i,j,k,XDIR:ZDIR)
invlat = InvLatitude(xyz)
invlon = katan2(xyz(YDIR),xyz(XDIR))
write (checkMessage,'(A,I0,A,I0,A,I0,A)'), 'Squish Fake Projection Latitude is wrong at (',i,',',j,',',k,'). Check Squish Processing and Output.'
@assertEqual(invlat,voltronApp%chmp2mhd%xyzSquish(i,j,k,1),1e-17_rp,trim(checkMessage))
write (checkMessage,'(A,I0,A,I0,A,I0,A)'), 'Squish Fake Projection Longitude is wrong at (',i,',',j,',',k,'). Check Squish Processing and Output.'
@assertEqual(invlon,voltronApp%chmp2mhd%xyzSquish(i,j,k,2),1e-17_rp,trim(checkMessage))
endif
enddo
enddo
enddo