00001 subroutine calc_3velocity_corot
00002 use phys_constant, only : long
00003 use grid_parameter, only : nrf, ntf, npf
00004 use def_matter_parameter, only : ome
00005 use def_matter_velocity, only : vxu, vyu, vzu
00006 use def_vector_phi, only : vec_phif
00007 implicit none
00008 integer :: ir, it, ip
00009
00010 do ip = 0, npf
00011 do it = 0, ntf
00012 do ir = 0, nrf
00013
00014 vxu(ir,it,ip) = ome*vec_phif(ir,it,ip,1)
00015 vyu(ir,it,ip) = ome*vec_phif(ir,it,ip,2)
00016 vzu(ir,it,ip) = 0.0d0
00017
00018 end do
00019 end do
00020 end do
00021 end subroutine calc_3velocity_corot