diff options
Diffstat (limited to 'benchmarks/CUDA/WP/compare_snaps.F')
| -rw-r--r-- | benchmarks/CUDA/WP/compare_snaps.F | 83 |
1 files changed, 0 insertions, 83 deletions
diff --git a/benchmarks/CUDA/WP/compare_snaps.F b/benchmarks/CUDA/WP/compare_snaps.F deleted file mode 100644 index cedbd39..0000000 --- a/benchmarks/CUDA/WP/compare_snaps.F +++ /dev/null @@ -1,83 +0,0 @@ -! file format -!write(46) ids,ide,jds,jde,kds,kde,ims,ime,jms,jme,kms,kme,its,ite,jts,jte,kts,kte -!write(46) th(its:ite,kts:kte,jts:jte) -!write(46) q(its:ite,kts:kte,jts:jte) -!write(46) qc(its:ite,kts:kte,jts:jte) -!write(46) qr(its:ite,kts:kte,jts:jte) -!write(46) qi(its:ite,kts:kte,jts:jte) -!write(46) qs(its:ite,kts:kte,jts:jte) -!!INOUT 2D -!write(46) rain(its:ite,jts:jte) -!write(46) rainncv(its:ite,jts:jte) -!write(46) sr(its:ite,jts:jte) -!write(46) snow(its:ite,jts:jte) -!write(46) snowncv(its:ite,jts:jte) - - program compare_snaps - implicit none - integer, parameter :: unit1 = 31 - integer, parameter :: unit2 = 32 - integer ids,ide,jds,jde,kds,kde,ims,ime,jms,jme,kms,kme,its,ite,jts,jte,kts,kte - - open(unit=unit1,file="snap_out_010",status="old",form="unformatted",err=22) - open(unit=unit2,file="snap_gpu_010",status="old",form="unformatted",err=23) - - read(unit=unit1) ids,ide,jds,jde,kds,kde,ims,ime,jms,jme,kms,kme,its,ite,jts,jte,kts,kte - write(*,*)ids,ide,jds,jde,kds,kde - write(*,*)ims,ime,jms,jme,kms,kme - write(*,*)its,ite,jts,jte,kts,kte - read(unit=unit2) ids,ide,jds,jde,kds,kde,ims,ime,jms,jme,kms,kme,its,ite,jts,jte,kts,kte - write(*,*)ids,ide,jds,jde,kds,kde - write(*,*)ims,ime,jms,jme,kms,kme - write(*,*)its,ite,jts,jte,kts,kte - - call read_and_output3( 'th', unit1, 88 ) - call read_and_output3( 'th', unit2, 98 ) - call read_and_output3( 'q', unit1, 88 ) - call read_and_output3( 'q', unit2, 98 ) - call read_and_output3( 'qc', unit1, 88 ) - call read_and_output3( 'qc', unit2, 98 ) - call read_and_output3( 'qr', unit1, 88 ) - call read_and_output3( 'qr', unit2, 98 ) - call read_and_output3( 'qi', unit1, 88 ) - call read_and_output3( 'qi', unit2, 98 ) - call read_and_output3( 'qs', unit1, 88 ) - call read_and_output3( 'qs', unit2, 98 ) - - close(unit1) - close(unit2) - stop -22 write(*,*)'error opening ',unit1 - stop -23 write(*,*)'error opening ',unit2 - stop - end - - subroutine read_and_output3 ( varname,iu, ou ) - implicit none - integer iu,ou - character*(*) varname -! 2 72 2 59 1 27 - integer, parameter :: ix = 71 - integer, parameter :: jx = 58 - integer, parameter :: kx = 27 - integer i,j,k - real s3(ix,kx,jx) - - read(iu)s3 - write(ou,*)ix,jx,' ',trim(varname) - do j = 1,jx - do i = 1,ix - write(ou,*)s3(i,1,j) - enddo - enddo - return - end - - - - - - - - |