demona: D:/demona/2.0/manual/doc/2-31set/processor/sets/globalSet.f90 Source File

00001 subroutine globalSet
00002 
00003 use global
00004 
00005 implicit none
00006 
00007 integer i,j,t
00008 
00009 integer myP,myD
00010 
00011 integer,pointer :: IM(:,:,:,:)
00012 
00013 
00014 integer,pointer :: minI,minJ,maxI,maxJ
00015 
00016 
00017 
00018 
00019 do myP=1,size(pVec)
00020 
00021         currentPhysics=pVEc(myP)
00022 
00023         dVec=>myModel%physics(currentPhysics)%dVec
00024 
00025         do myD=1,size(dVec)
00026 
00027                 currentDomain=dVec(myD)
00028 
00029                 IM=>myModel%physics(currentPhysics)%domain(currentDomain)%dofMatrix
00030 
00031                 minI=>myModel%physics(currentPhysics)%domain(currentDomain)%minI
00032                 maxI=>myModel%physics(currentPhysics)%domain(currentDomain)%maxI
00033                 minJ=>myModel%physics(currentPhysics)%domain(currentDomain)%minJ
00034                 maxJ=>myModel%physics(currentPhysics)%domain(currentDomain)%maxJ
00035 
00036                                                 
00037                 !IM(:,:,:,flag)=1
00038                 FORALL(i = minI:maxI, j = minJ:maxJ, t=1:myModel%physics(currentPhysics)%dof, IM(i,j,t,absolute).ne.0) IM(i,j,t,flag)=1
00039 
00040         end do
00041 
00042 end do
00043 
00044 call currentSet
00045 
00046 
00047 
00048 
00049 
00050 end subroutine globalSet