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

00001 subroutine update_parameter_BNS(conv_den)
00002   use grid_parameter, only : EQ_point, nrf_deform, nrf, NS_shape
00003   use def_matter_parameter, only : ROT_LAW
00004   use interface_update_parameter_BNS_peos
00005   use interface_update_parameter_BNS_peos_lecc
00006   use interface_update_parameter_BNS_peos_irrot
00007   use interface_update_parameter_BNS_peos_lecc_irrot
00008   use interface_update_parameter_BNS_peos_spin
00009   use interface_update_parameter_BNS_peos_lecc_spin
00010 !  use interface_update_parameter_axisym_peos
00011 !  use interface_update_parameter_triaxial_peos
00012 !  use interface_update_parameter_spherical_peos
00013   use interface_update_parameter_axisym_peos_drot
00014   implicit none
00015   real(8), intent(in) :: conv_den
00016   if (ROT_LAW.eq.'DR') then 
00017     call update_parameter_axisym_peos_drot(conv_den)
00018   else 
00019 !    call update_parameter_BNS_peos(conv_den)
00020 
00021     if (NS_shape.eq.'CO')   then
00022       if (EQ_point.eq.'LE')  then
00023         call update_parameter_BNS_peos_lecc(conv_den)
00024       else 
00025         call update_parameter_BNS_peos(conv_den)
00026       end if
00027     end if
00028 
00029     if (NS_shape.eq.'IR')   then
00030       if (EQ_point.eq.'LE')  then
00031         call update_parameter_BNS_peos_lecc_irrot(conv_den)
00032       else 
00033         call update_parameter_BNS_peos_irrot(conv_den)
00034       end if
00035     end if
00036 
00037     if (NS_shape.eq.'SP')   then
00038       if (EQ_point.eq.'LE')  then
00039         call update_parameter_BNS_peos_lecc_spin(conv_den)
00040       else 
00041         call update_parameter_BNS_peos_spin(conv_den)
00042       end if
00043     end if
00044 
00045     if(NS_shape.ne.'CO' .and. NS_shape.ne.'IR' .and. NS_shape.ne.'SP')  then
00046       write(6,*) "NS_shape is neither CO nor IR nor SP. Exiting..."
00047       stop
00048     end if
00049 !    if (EQ_point.eq.'XZ'.and.nrf.ne.nrf_deform) &
00050 !    &                     call update_parameter_axisym_peos(conv_den)
00051 !    if (EQ_point.eq.'XY') call update_parameter_triaxial_peos(conv_den)
00052 !    if (EQ_point.eq.'XZ'.and.nrf.eq.nrf_deform) &
00053 !    &                     call update_parameter_spherical_peos(conv_den)
00054   end if
00055 end subroutine update_parameter_BNS
/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/update_parameter_BNS.f90
