D:/demona/2.0/manual/doc/2-31set/problems/test/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 !##################################################################### 00023 !user supplied variables 00024 !##################################################################### 00025 select case(currentPhysics) 00026 00027 case(1) 00028 00029 numberOfDomains=1 00030 !dof=2 00031 dof=1 00032 dFX=1 00033 TecplotVariables='VARIABLES="x","y","T"' 00034 !TecplotVariables='VARIABLES="x","y","T1","U1"' 00035 !TecplotVariables='VARIABLES="x","y","str","vor"' 00036 refXcoor=0.d0 00037 refYcoor=0.d0 00038 Lx=1.d0 00039 Ly=1.d0 00040 00041 allocate(myModel%physics(currentPhysics)%damping(dof)) 00042 damping=>myModel%physics(currentPhysics)%damping 00043 damping(1)=1.0 00044 !damping(2)=0.001 00045 00046 00047 00048 case(2) 00049 numberOfDomains=1 00050 dof=1 00051 dFX=1 00052 TecplotVariables='VARIABLES="x","y","T"' 00053 !TecplotVariables='VARIABLES="x","y","T2","U2"' 00054 refXcoor=0.875d0 00055 refYcoor=-0.125d0 00056 Lx=0.5d0 00057 Ly=1.25d0 00058 00059 00060 allocate(myModel%physics(currentPhysics)%damping(dof)) 00061 damping=>myModel%physics(currentPhysics)%damping 00062 damping(1)=1.0 00063 !damping(2)=0.01 00064 00065 case(3) 00066 numberOfDomains=1 00067 dof=2 00068 dFX=1 00069 TecplotVariables='VARIABLES="x","y","str","vor"' 00070 refXcoor=0.d0 00071 refYcoor=0.d0 00072 00073 allocate(myModel%physics(currentPhysics)%damping(dof)) 00074 damping=>myModel%physics(currentPhysics)%damping 00075 damping(1)=1.0 00076 damping(2)=0.01 00077 00078 !< the format below could also been used 00079 !case(1:2) 00080 ! numberOfDomains=1 00081 ! dof=1 00082 ! dFX=1 00083 00084 case default 00085 00086 print*,'invalid physics definition' 00087 print*,'demona stop' 00088 00089 end select 00090 !##################################################################### 00091 00092 !##################################################################### 00093 !user supplied variables 00094 !##################################################################### 00095 ! if different pVec are to be used, something similar to the code lines 00096 ! below should be implemented. 00097 !deallocate(myModel%pVec) 00098 !allocate(myModel%pVec(*)) 00099 !pVec(1)=2 00100 !pVec(2)=5 00101 !##################################################################### 00102 00103 00104 myModel%physics(currentPhysics)%tecVariables=TecplotVariables 00105 00106 allocate(myModel%physics(currentPhysics)%endIndex(numberOfDomains)) 00107 myModel%physics(currentPhysics)%endIndex=1 00108 00109 00110 00111 00112 end subroutine physicsInfo
