D:/demona/2.0/manual/doc/2-31set/problems/ldc_fv/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=3 00034 !dof=1 00035 dFX=2 00036 TecplotVariables='VARIABLES="x","y","p","u","v"' 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)=0.01d0 00045 damping(2)=0.01d0 00046 damping(3)=0.01d0 00047 !damping(2)=1.d0 00048 00049 numberOfPhysicalParameters=1 00050 allocate(myModel%physics(currentPhysics)%physicalParameters(numberOfPhysicalParameters)) 00051 allocate(myModel%physics(currentPhysics)%physicalParametersInitial(numberOfPhysicalParameters)) 00052 allocate(myModel%physics(currentPhysics)%physicalParametersFinal(numberOfPhysicalParameters)) 00053 myModel%physics(currentPhysics)%physicalParameters=0.d0 00054 myModel%physics(currentPhysics)%physicalParametersInitial=100.d0 00055 myModel%physics(currentPhysics)%physicalParametersFinal=1000.d0 00056 00057 myModel%physics(currentPhysics)%discOrder=highOrder 00058 myModel%physics(currentPhysics)%evalOrder=>myModel%physics(currentPhysics)%discOrder 00059 00060 relativeModel=relativeBC 00061 00062 case default 00063 00064 print*,'invalid physics definition' 00065 print*,'demona stop' 00066 00067 end select 00068 !##################################################################### 00069 00070 !##################################################################### 00071 !user supplied variables 00072 !##################################################################### 00073 ! if different pVec are to be used, something similar to the code lines 00074 ! below should be implemented. 00075 !deallocate(myModel%pVec) 00076 !allocate(myModel%pVec(*)) 00077 !pVec(1)=2 00078 !pVec(2)=5 00079 !##################################################################### 00080 00081 00082 myModel%physics(currentPhysics)%tecVariables=TecplotVariables 00083 00084 allocate(myModel%physics(currentPhysics)%endIndex(numberOfDomains)) 00085 myModel%physics(currentPhysics)%endIndex=1 00086 00087 00088 00089 00090 end subroutine physicsInfo
