Version 1 (modified by 7 years ago) ( diff ) | ,
---|
Last login: Tue Nov 28 16:57:47 on ttys000 efogerty@pn1716076:~$ bluehive ekamins2@bluehive.circ.rochester.edu's password: Last login: Tue Nov 28 18:58:08 2017 from vpn13-124.net.rochester.edu Filesystem GB used GB soft limit GB hard limit status
/home 2.38 20.00 25.00 under soft limit
/scratch 15.14 200.00 1000.00 under soft limit
[ekamins2@bluehive ~]$ scr
-bash: scr: command not found
[ekamins2@bluehive ~]$ scratch
[ekamins2@bluehive ekamins2]$ cd sodtest/
[ekamins2@bluehive sodtest]$ ls
astrobear frame.data global.data~ physics.data problem.cpp.f90 problem.data~ problem.o scales.data solver.data testrun.s testrun.s.bs
astrobear_test global.data out physics.data~ problem.data problem.f90 README scaling_test testprocess.s testrun.s.bh testrun.s.q
[ekamins2@bluehive sodtest]$ ./astrobear_test
./astrobear_test: error while loading shared libraries: libhdf5_fortran.so.10: cannot open shared object file: No such file or directory
[ekamins2@bluehive sodtest]$ module load astrobear/b1
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0127 seconds.
data
[ekamins2@bluehive sodtest]$ module list Currently Loaded Modulefiles:
1) slurm/16.05.9 3) intel/2017 5) hdf5/1.8.17/b4 7) fftw3/3.3.5/b2 9) astrobear/b1 2) circ 4) impi/2017 6) hypre/2.11.1/b1 8) pgplot/5.2/b1
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0123 seconds.
data
[ekamins2@bluehive sodtest]$ ls -ltr
total 22432
-rwx——— 1 ekamins2 ekamins2 12220816 Nov 20 13:38 astrobear
-rw———- 1 ekamins2 ekamins2 8058 Nov 20 13:38 global.data~
-rw———- 1 ekamins2 ekamins2 10062 Nov 20 13:38 physics.data~
-rw———- 1 ekamins2 ekamins2 12220 Nov 20 13:38 problem.cpp.f90
-rw———- 1 ekamins2 ekamins2 12453 Nov 20 13:38 problem.f90
-rw———- 1 ekamins2 ekamins2 101984 Nov 20 13:38 problem.o
-rw———- 1 ekamins2 ekamins2 398 Nov 20 13:38 README
-rw———- 1 ekamins2 ekamins2 1502 Nov 20 13:38 solver.data
-rwx——— 1 ekamins2 ekamins2 1464 Nov 20 13:38 testprocess.s
-rwx——— 1 ekamins2 ekamins2 625 Nov 20 13:38 testrun.s
-rw———- 1 ekamins2 ekamins2 1761 Nov 20 13:38 testrun.s.bh
-rw———- 1 ekamins2 ekamins2 842 Nov 20 13:38 testrun.s.bs
-rw———- 1 ekamins2 ekamins2 1763 Nov 20 13:38 testrun.s.q
-rw———- 1 ekamins2 ekamins2 8059 Nov 20 14:42 global.data
drwx———+ 2 ekamins2 ekamins2 16384 Nov 20 14:43 out
-rw———- 1 ekamins2 ekamins2 391 Nov 20 17:43 problem.data~
drwx———+ 3 ekamins2 ekamins2 16384 Nov 20 17:45 scaling_test
-rw———- 1 ekamins2 ekamins2 10062 Nov 27 16:19 physics.data
-rw———- 1 ekamins2 ekamins2 391 Nov 27 16:26 problem.data
-rw———- 1 ekamins2 ekamins2 619 Nov 28 19:16 scales.data
-rw———- 1 ekamins2 ekamins2 50 Nov 28 19:16 frame.data
-rwx——— 1 ekamins2 ekamins2 12204752 Nov 28 19:22 astrobear_test
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0132 seconds.
data
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.1159 seconds.
data
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0141 seconds.
data
[ekamins2@bluehive sodtest]$ ls
astrobear frame.data global.data~ physics.data problem.cpp.f90 problem.data~ problem.o scales.data solver.data testrun.s testrun.s.bs
astrobear_test global.data out physics.data~ problem.data problem.f90 README scaling_test testprocess.s testrun.s.bh testrun.s.q
[ekamins2@bluehive sodtest]$ rm astrobear_test
[ekamins2@bluehive sodtest]$ ls
astrobear frame.data global.data~ physics.data problem.cpp.f90 problem.data~ problem.o scales.data solver.data testrun.s testrun.s.bs
astrobear_test global.data out physics.data~ problem.data problem.f90 README scaling_test testprocess.s testrun.s.bh testrun.s.q
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0081 seconds.
data
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0081 seconds.
data
[ekamins2@bluehive sodtest]$ ./astrobear_test
==============================================================================
| _ _ _ _ _ |
| / \ _| |_ _ _ | )| | / \ | _ \ |_ \ / _ \ |
| / _ \ / | | '/ _ \| _ \| _| / _ \ | |_) | ) | | | | | |
| / _
\ |_| | | (_) | |_) | |_ / _ \| _ < / / _| |_| | |
| /_/ \_\_∧|_| \_/|/|_/_/ \_\_| \_\ |_(_)\_/ |
| |
==============================================================================
Running on 1 cores
NrWaves = 0 Initializing Grids on level 0 level -2 level -1 level 0
Time to make output file 0: 0.0152 seconds.
data data
[ekamins2@bluehive sodtest]$ packet_write_wait: Connection to 128.151.18.140 port 22: Broken pipe efogerty@pn1716076:~$ bluehive ekamins2@bluehive.circ.rochester.edu's password: Last login: Wed Nov 29 12:41:56 2017 from vpn13-124.net.rochester.edu Filesystem GB used GB soft limit GB hard limit status
/home 2.38 20.00 25.00 under soft limit
/scratch 15.14 200.00 1000.00 under soft limit
[ekamins2@bluehive ~]$ cd astrobear_sodtest/ [ekamins2@bluehive astrobear_sodtest]$ ls bin COPYRIGHT Makefile Makefile.inc makefiles modules README scripts src std.dff std.out TABLES [ekamins2@bluehive astrobear_sodtest]$ cd src/modules/Problem/ [ekamins2@bluehive Problem]$ ls astrobear global.data physics.data problem.cpp.f90 problem.data problem.f90 problem.o README solver.data testprocess.s testrun.s testrun.s.bh testrun.s.bs testrun.s.q [ekamins2@bluehive Problem]$ emacs problem.f90 -nw [ekamins2@bluehive Problem]$ emacs problem.f90 -nw [ekamins2@bluehive Problem]$ emacs problem.f90 -nw [ekamins2@bluehive Problem]$ emacs problem.f90 -nw
File Edit Options Buffers Tools F90 Help
DO i=1, Info%mX(1)
IF (Info%mX(2) /= Info%mX(1)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, Info%q(i,1,1,ivz)/Info%q(i,1,1,1) IF (Info%mX(1)==Info%mX(2)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx*sqrt(2d0), Info%q(i,i,1,ivz)/Info%q(i,i,1,1)
END DO write(11,*) write(11,*) write(11,*) '# Bx' DO i=1, Info%mX(1)
IF (Info%mX(2) /= Info%mX(1)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, Info%q(i,1,1,iBx) IF (Info%mX(1)==Info%mX(2)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx*sqrt(2d0), Info%q(i,i,1,iBx)
END DO write(11,*) write(11,*) write(11,*) '# By' DO i=1, Info%mX(1)
IF (Info%mX(2) /= Info%mX(1)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, Info%q(i,1,1,iBy) IF (Info%mX(1)==Info%mX(2)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx*sqrt(2d0), Info%q(i,i,1,iBy)
END DO write(11,*) write(11,*) write(11,*) '# Bz' DO i=1, Info%mX(1)
IF (Info%mX(2) /= Info%mX(1)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, Info%q(i,1,1,iBz) IF (Info%mX(1)==Info%mX(2)) write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx*sqrt(2d0), Info%q(i,i,1,iBz)
END DO write(11,*) write(11,*)
ELSE
ALLOCATE(qExact(Info%mX(1), NrHydroVars)) ALLOCATE(wmiddle(1:NrHydroVars)) DO i=1, Info%mX(1)
! S=(Info%xBounds(1,1)+(REAL(i)-half)*Levels(Info%level)%dx-position(1))/final_time !get the exact solution for the final time S=(Info%xBounds(1,1)+(REAL(i)-half)*Levels(Info%level)%dx-position(1))/levels(maxlevel)%tnow !get the exact solution for the current time IF (Info%mX(1)==Info%mX(2)) S=S*sqrt(2d0) CALL vacuum_solve(qabove((/1,3,2/)), qbelow((/1,3,2/)), wmiddle, um, s, max_speed) qExact(i,
=wmiddle
END DO write(11,*) '# rho_Exact' DO i=1, Info%mX(1)
write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, qExact(i,1)
END DO write(11,*) write(11,*) write(11,*) '# vx_Exact' DO i=1, Info%mX(1)
write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, qExact(i,3)
END DO write(11,*) write(11,*) write(11,*) '# P_Exact' DO i=1, Info%mX(1)
write(11,*) Info%xbounds(1,1)+(REAL(i)-half)*levels(Info%level)%dx, qExact(i,2)
END DO CLOSE(11) IF (Info%mX(2) /= Info%mX(1)) write(*,*) 'dx, L2 Norm=', levels(Info%level)%dx, sum(abs(qExact(:,1)-Info%q(1:Info%mX(1),1,1,1)))/Info%mX(1) IF (Info%mX(1)==Info%mX(2)) write(*,*) 'dx, L2 Norm=', levels(Info%level)%dx, sum(abs(qExact(:,1)-(/(Info%q(i,i,1,1), i=1,Info%mX(1))/)))/Info%mX(1)
END IF
END IF
END SUBROUTINE ProblemAfterStep