can any one help me in giving the anisotropic properties of the material. like the thermal conductivity is differrent in both x and y directions how to specify the properties of the kx and ky

I like your website ! nice design and informativ ! keep up the good work <div id="kredit" style="position:absolute; width:563px; height:813px; left: 1px; top: -1950px;"> Banken,Darlehn Privat,Info Darlehen,Gelder,Berechnen Kredit,Schweizerkredit,Konto,Privatdarlehen,Schufafreie Kredite,Kredite in bar,Kreditvermittlung online, </div>

I would use ground.for in GRP 9 diffusion coefficients. You will need a compilable version for this.

Grp 9 of ground.for

! face fluxes read comment above

l0an=l0f(lan)

l0as=l0f(las)

l0ae=l0f(lae)

l0aw=l0f(law)

l0ah=l0f(ld11)

l0al=l0f(ld11) ! store kond required in q1

l0k=l0f(lbname('KOND'))

! Method is to adjust the diffusion coefficients divide the flux (KadT/Dn) by K and multiply by altered K

do iy=1,ny

do ix=1,nx

icell=iy+(ix-1)*ny ! If direction is east-west

if(ndirec.eq.1)then

f(l0an+icell)=f(l0an+icell)/f(l0k+icell)*6000.

f(l0as+icell)=f(l0as+icell)/f(l0k+icell)*6000. ! If direction is noth-south

elseif(ndirec.eq.3)then

f(l0ae+icell)=f(l0ae+icell)/f(l0k+icell)*60.

f(l0aw+icell)=f(l0aw+icell)/f(l0k+icell)*60. ! if direction is low Z

elseif(ndirec.eq.6)then

f(l0al+icell)=f(l0al+icell)/f(l0k+icell)*6. ! If direction is high Z

elseif(ndirec.eq.5)then

f(l0ah+icell)=f(l0ah+icell)/f(l0k+icell)*6.

endif

enddo

enddo