C C User subroutine vufield C subroutine VUField( C Write only - * rUserField, C Read only - * nBlock, nField, kField, nComp, * kStep, kInc, jNodeUid, time, * coords, U, V, A ) * include 'vaba_param.inc' dimension rUserField(nBlock,nComp,nField) dimension jNodeUid(nBlock), time(4), coords(3,nBlock) dimension U(8,nBlock), V(8,nBlock), A(8,nBlock) * parameter ( i_ufld_Current = 1, * i_ufld_Increment = 2, * i_ufld_Period = 3, * i_ufld_Total = 4 ) * parameter ( i_ufld_CoordX = 1, * i_ufld_CoordY = 2, * i_ufld_CoordZ = 3 ) * parameter ( i_ufld_SpaDisplX = 1, * i_ufld_SpaDisplY = 2, * i_ufld_SpaDisplZ = 3, * i_ufld_RotDisplX = 4, * i_ufld_RotDisplY = 5, * i_ufld_RotDisplZ = 6, * i_ufld_AcoPress = 7, * i_ufld_Temp = 8 ) * parameter ( i_ufld_SpaVelX = 1, * i_ufld_SpaVelY = 2, * i_ufld_SpaVelZ = 3, * i_ufld_RotVelX = 4, * i_ufld_RotVelY = 5, * i_ufld_RotVelZ = 6, * i_ufld_DAcoPress = 7, * i_ufld_DTemp = 8 ) * parameter ( i_ufld_SpaAccelX = 1, * i_ufld_SpaAccelY = 2, * i_ufld_SpaAccelZ = 3, * i_ufld_RotAccelX = 4, * i_ufld_RotAccelY = 5, * i_ufld_RotAccelZ = 6, * i_ufld_DDAcoPress = 7, * i_ufld_DDTemp = 8 ) parameter (ptFive = 0.5d0) * if (kField .eq. 2) then do kComp = 1, nComp do kNod = 1, nBlock rUserField(kNod,kComp,1) = ptFive end do end do end if * return end