00001 subroutine generateCoordinates(dx,dy)
00002
00003 use global
00004
00005 integer myDomaim
00006
00007
00008 integer,pointer :: myfullNx,myfullNy
00009
00010 integer,pointer :: myNumberOfDomains
00011
00012 integer,pointer :: myLocalDof
00013
00014 real(8) dx,dy
00015
00016 real(8) myX,myY
00017
00018
00019
00020 myNumberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains
00021
00022
00023
00024 myfullNx=>myModel%physics(currentPhysics)%domain(currentDomain)%Nx
00025 myfullNy=>myModel%physics(currentPhysics)%domain(currentDomain)%Ny
00026
00027
00028 do j=1,myfullNy
00029
00030
00031 myY=(j-1)*dy+myModel%physics(currentPhysics)%refYcoor+myModel%physics(currentPhysics)%domain(currentDomain)%relYcoor
00032
00033 do i=1,myfullNx
00034
00035
00036 myX=(i-1)*dx+myModel%physics(currentPhysics)%refXcoor+myModel%physics(currentPhysics)%domain(currentDomain)%relXcoor
00037
00038
00039 if (myModel%physics(currentPhysics)%domain(currentDomain)%nodeMatrix(i,j,2).ne.0) then
00040
00041
00042
00043 myModel%physics(currentPhysics)%domain(currentDomain)%xCoor(myModel%physics(currentPhysics)%domain(currentDomain)%nodeMatrix(i,j,2))=myX
00044 myModel%physics(currentPhysics)%domain(currentDomain)%yCoor(myModel%physics(currentPhysics)%domain(currentDomain)%nodeMatrix(i,j,2))=myY
00045
00046
00047 end if
00048
00049 end do
00050
00051 end do
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064 end subroutine generateCoordinates