Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/write_parameter

00001 subroutine write_parameter_drot
00002   use def_matter_parameter
00003   implicit none
00004   character(len=28) dum
00005   character(len=50) char_50(4)
00006 !
00007   open(1,file='rnspar_drot.dat',status='old')
00008   read(1,'(a28,a50)') dum, char_50(1)  ! For ROT_LAW, DROT_ratio
00009   if (ROT_LAW.eq."DR") then
00010     read(1,'(a28,a50)') dum, char_50(2)  ! For index_DR, A2DR
00011   end if
00012   if (ROT_LAW.eq."OJ") then
00013     read(1,'(a28,a50)') dum, char_50(2)  ! For index_DR, A2DR
00014     read(1,'(a28,a50)') dum, char_50(3)  ! For index_DRp, B2DR
00015     read(1,'(a28,a50)') dum, char_50(4)  ! For DRAT_A2DR, DRAT_B2DR
00016   end if
00017   close(1)
00018 !
00019   open(1,file='rnspar_drot.las',status='unknown')
00020   write(1,'(2x,a2,10x,1p,e14.6,a50)') ROT_LAW, DRAT_A2DR, char_50(1)
00021   if (ROT_LAW.eq."DR") then
00022     write(1,'(1p,2e14.6,a50)')      index_DR, A2DR,       char_50(2)
00023   end if
00024   if (ROT_LAW.eq."OJ") then
00025     write(1,'(1p,2e14.6,a50)')      index_DRq, A2DR,      char_50(2)
00026     write(1,'(1p,2e14.6,a50)')      index_DRp, B2DR,      char_50(3)
00027     write(1,'(1p,2e14.6,a50)')      DRAT_A2DR, DRAT_B2DR, char_50(4)
00028   end if
00029   close(1)
00030 end subroutine write_parameter_drot

/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/write_parameter_drot.f90