00001 subroutine save_solution_BNS_mpt(icycle)
00002 use phys_constant, only : long
00003 use def_bh_parameter
00004 use def_quantities
00005 use grid_parameter, only : rgin
00006 implicit none
00007 integer :: icycle
00008 character(40) :: char1, char2, char3, char4, char5
00009 character(100) :: dircommand
00010
00011 write(char1, '(i5)') icycle
00012 char2 = adjustl(char1)
00013 if (icycle < 10) then
00014 char5 = 'cycle_0' // trim(char2)
00015 else
00016 char5 = 'cycle_' // trim(char2)
00017 endif
00018
00019 dircommand = 'mkdir ' // char5
00020 call system(dircommand)
00021
00022
00023
00024
00025 call chdir(char5)
00026
00027
00028 call write_last_physq_BNS_mpt
00029 char3 = 'main_bnsphys_all_mpt.txt'
00030 call printout_physq_BNS_all_mpt(char3)
00031
00032 call chdir('../')
00033
00034 end subroutine save_solution_BNS_mpt