demona: D:/demona/2.0/manual/doc/2-31set/preprocessor/initializeDomains.f90 Source File

00001 subroutine initializeDomains
00002 
00003 
00004 
00005 use global
00006 
00007 
00008 implicit none
00009 
00010 
00011 integer,pointer :: minI,minJ,maxI,maxJ
00012 
00013 integer,pointer :: nx,ny
00014 
00015 !allocate(myModel%physics(numberOfPhysics))
00016 
00017 
00018 
00019 
00020 call domainInfo
00021 
00022 minI=>myModel%physics(currentPhysics)%domain(currentDomain)%minI
00023 maxI=>myModel%physics(currentPhysics)%domain(currentDomain)%maxI
00024 minJ=>myModel%physics(currentPhysics)%domain(currentDomain)%minJ
00025 maxJ=>myModel%physics(currentPhysics)%domain(currentDomain)%maxJ
00026 
00027 dof=>myModel%physics(currentPhysics)%dof
00028 
00029 
00030 nx=>myModel%physics(currentPhysics)%domain(currentDomain)%nx
00031 ny=>myModel%physics(currentPhysics)%domain(currentDomain)%ny
00032 
00033 
00034 allocate(myModel%physics(currentPhysics)%domain(currentDomain)%dofMatrix(minI:maxI,minJ:maxJ,1:dof,4))
00035 allocate(myModel%physics(currentPhysics)%domain(currentDomain)%nodeMatrix(nx,ny,2))
00036 
00037 
00038 myModel%physics(currentPhysics)%domain(currentDomain)%dofMatrix=0
00039 myModel%physics(currentPhysics)%domain(currentDomain)%nodeMatrix=0
00040 
00041 
00042 
00043 
00044 ! if skip function is to be applied, it should be called here.
00045 
00046 
00047 
00048 
00049 
00050 
00051 end subroutine initializeDomains