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=3 00031 !dof=1 00032 dFX=1 00033 TecplotVariables='VARIABLES="x","y","str","vor","temp"' 00034 refXcoor=0.d0 00035 refYcoor=0.d0 00036 Lx=1.d0 00037 Ly=1.d0 00038 00039 allocate(myModel%physics(currentPhysics)%damping(dof)) 00040 damping=>myModel%physics(currentPhysics)%damping 00041 damping(1)=0.01d0 00042 damping(2)=0.01d0 00043 damping(3)=0.01d0 00044 !damping(2)=1.d0 00045 00046 00047 00048 00049 case default 00050 00051 print*,'invalid physics definition' 00052 print*,'demona stop' 00053 00054 end select 00055 !##################################################################### 00056 00057 !##################################################################### 00058 !user supplied variables 00059 !##################################################################### 00060 ! if different pVec are to be used, something similar to the code lines 00061 ! below should be implemented. 00062 !deallocate(myModel%pVec) 00063 !allocate(myModel%pVec(*)) 00064 !pVec(1)=2 00065 !pVec(2)=5 00066 !##################################################################### 00067 00068 00069 myModel%physics(currentPhysics)%tecVariables=TecplotVariables 00070 00071 allocate(myModel%physics(currentPhysics)%endIndex(numberOfDomains)) 00072 myModel%physics(currentPhysics)%endIndex=1 00073 00074 00075 00076 00077 end subroutine physicsInfo
1.6.1