00001 subroutine read_parameter_drot
00002 use def_matter_parameter
00003 implicit none
00004 open(1,file='rnspar_drot.dat',status='old')
00005 read(1,'(2x,a2,10x,1p,e14.6)') ROT_LAW, DRAT_A2DR
00006 if (ROT_LAW.eq."DR") then
00007 read(1,'(1p,2e14.6)') index_DR, A2DR
00008 end if
00009 if (ROT_LAW.eq."OJ") then
00010 read(1,'(1p,2e14.6)') index_DRq, A2DR
00011 read(1,'(1p,2e14.6)') index_DRp, B2DR
00012 read(1,'(1p,2e14.6)') DRAT_A2DR, DRAT_B2DR
00013 if (index_DRq.ne.3.0) stop 'NOT KEPLER'
00014 end if
00015 close(1)
00016 end subroutine read_parameter_drot
