wiki:u/erica/ProblemAfterStepRoutine

Version 1 (modified by Erica Kaminski, 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

Note: See TracWiki for help on using the wiki.