D:/demona/2.0/manual/doc/2-31set/problems/bratu/1/info/physicsInfo.f90
Go to the documentation of this file.00001 subroutine physicsInfo 00002 00003 00004 use global 00005 00006 00007 implicit none 00008 00009 00010 00011 numberOfDomains=>myModel%physics(currentPhysics)%numberOfDomains 00012 dof=>myModel%physics(currentPhysics)%dof 00013 dFX=>myModel%physics(currentPhysics)%dFX 00014 00015 00016 refXcoor=>myModel%physics(currentPhysics)%refXcoor 00017 refYcoor=>myModel%physics(currentPhysics)%refYcoor 00018 00019 Lx=>myModel%physics(currentPhysics)%Lx 00020 Ly=>myModel%physics(currentPhysics)%Ly 00021 00022 numberOfPhysicalParameters=>myModel%physics(currentPhysics)%numberOfPhysicalParameters 00023 00024 00025 !##################################################################### 00026 !user supplied variables 00027 !##################################################################### 00028 select case(currentPhysics) 00029 00030 case(1) 00031 00032 numberOfDomains=1 00033 dof=1 00034 !dof=1 00035 dFX=1 00036 TecplotVariables='VARIABLES="x","y","t"' 00037 refXcoor=0.d0 00038 refYcoor=0.d0 00039 Lx=1.d0 00040 Ly=1.d0 00041 00042 allocate(myModel%physics(currentPhysics)%damping(dof)) 00043 damping=>myModel%physics(currentPhysics)%damping 00044 damping(1)=1.0d0 00045 !damping(2)=0.01d0 00046 !damping(2)=1.d0 00047 00048 numberOfPhysicalParameters=1 00049 allocate(myModel%physics(currentPhysics)%physicalParameters(numberOfPhysicalParameters)) 00050 allocate(myModel%physics(currentPhysics)%physicalParametersInitial(numberOfPhysicalParameters)) 00051 allocate(myModel%physics(currentPhysics)%physicalParametersFinal(numberOfPhysicalParameters)) 00052 myModel%physics(currentPhysics)%physicalParameters=0.d0 00053 myModel%physics(currentPhysics)%physicalParametersInitial=200.d0 00054 myModel%physics(currentPhysics)%physicalParametersFinal=6.8d0 00055 00056 00057 00058 case default 00059 00060 print*,'invalid physics definition' 00061 print*,'demona stop' 00062 00063 end select 00064 !##################################################################### 00065 00066 !##################################################################### 00067 !user supplied variables 00068 !##################################################################### 00069 ! if different pVec are to be used, something similar to the code lines 00070 ! below should be implemented. 00071 !deallocate(myModel%pVec) 00072 !allocate(myModel%pVec(*)) 00073 !pVec(1)=2 00074 !pVec(2)=5 00075 !##################################################################### 00076 00077 00078 myModel%physics(currentPhysics)%tecVariables=TecplotVariables 00079 00080 allocate(myModel%physics(currentPhysics)%endIndex(numberOfDomains)) 00081 myModel%physics(currentPhysics)%endIndex=1 00082 00083 00084 00085 00086 end subroutine physicsInfo
