00001 subroutine initializePhysics 00002 00003 00004 00005 use global 00006 00007 00008 implicit none 00009 00010 00011 integer i,j 00012 00013 00014 allocate(myModel%physics(numberOfPhysics)) 00015 00016 00017 00018 allocate(myModel%pVec(numberOfPhysics)) 00019 00020 pVec=>myModel%pVec 00021 00022 do i=1,numberOfPhysics 00023 00024 pVec(i)=i 00025 00026 end do 00027 00028 00029 00030 do j=1,numberOfPhysics 00031 00032 00033 currentPhysics=j 00034 00035 call physicsInfo 00036 00037 myModel%physics(currentPhysics)%totalDof=0 00038 00039 allocate(myModel%physics(currentPhysics)%domain(numberOfDomains)) 00040 00041 dFX=>myModel%physics(currentPhysics)%dFX 00042 00043 do i=1,numberOfDomains 00044 00045 currentDomain=i 00046 00047 !call domainInfo 00048 call initializeDomains 00049 00050 00051 end do 00052 00053 00054 00055 00056 00057 00058 end do 00059 00060 00061 00062 00063 00064 00065 00066 end subroutine initializePhysics
1.6.1