<div dir="ltr"><div><div><div>Hi,<br><br></div>I have a problem with hpts() function to interpolating data. I was using hpts() function without any problem in userchk and by input it was giving me the interpolation, but after some days I don't what change I did, whatever I interpolate now it just gives me multiplication of certain points(like in one case is 107 points and in other it is 90), I think there is some problems with setting parameters on SIZE, below is my SIZE file for one of my cases which I can just get interpolation with multiplication of 107 points (if I have less than 107 points in <a href="http://hpts.in">hpts.in</a>, it just gives me nothing in htps.out!!):<br>
<br>C Dimension file to be included <br>C<br>C HCUBE array dimensions<br>C<br> parameter (ldim=3)<br> parameter (lx1=6,ly1=lx1,lz1=lx1,lelt=410,lelv=lelt)<br> parameter (lxd=10,lyd=lxd,lzd=lxd)<br>
parameter (lelx=1,lely=1,lelz=1)<br> <br> parameter (lzl=3 + 2*(ldim-3))<br> <br> parameter (lx2=lx1-2)<br> parameter (ly2=ly1-2)<br> parameter (lz2=lz1-2)<br> parameter (lx3=lx1)<br> parameter (ly3=ly1)<br>
parameter (lz3=lz1)<br><br> parameter (lp = 512)<br> parameter (lelg = 16000)<br>c<br>c parameter (lpelv=lelv,lpelt=lelt,lpert=3) ! perturbation<br>c parameter (lpx1=lx1,lpy1=ly1,lpz1=lz1) ! array sizes<br>
c parameter (lpx2=lx2,lpy2=ly2,lpz2=lz2)<br>c<br> parameter (lpelv=1,lpelt=1,lpert=1) ! perturbation<br> parameter (lpx1=1,lpy1=1,lpz1=1) ! array sizes<br> parameter (lpx2=1,lpy2=1,lpz2=1)<br>
c<br>c parameter (lbelv=lelv,lbelt=lelt) ! MHD<br>c parameter (lbx1=lx1,lby1=ly1,lbz1=lz1) ! array sizes<br>c parameter (lbx2=lx2,lby2=ly2,lbz2=lz2)<br>c<br> parameter (lbelv=1,lbelt=1) ! MHD<br>
parameter (lbx1=1,lby1=1,lbz1=1) ! array sizes<br> parameter (lbx2=1,lby2=1,lbz2=1)<br> <br>C LX1M=LX1 when there are moving meshes; =1 otherwise<br> parameter (lx1m=1,ly1m=1,lz1m=1)<br> parameter (ldimt= 2) ! 3 passive scalars + T<br>
parameter (ldimt1=ldimt+1)<br> parameter (ldimt3=ldimt+3)<br>c<br>c Note: In the new code, LELGEC should be about sqrt(LELG)<br>c<br> PARAMETER (LELGEC = 1)<br> PARAMETER (LXYZ2 = 1)<br> PARAMETER (LXZ21 = 1)<br>
<br> PARAMETER (LMAXV=LX1*LY1*LZ1*LELV)<br> PARAMETER (LMAXT=LX1*LY1*LZ1*LELT)<br> PARAMETER (LMAXP=LX2*LY2*LZ2*LELV)<br> PARAMETER (LXZ=LX1*LZ1)<br> PARAMETER (LORDER=3)<br> PARAMETER (MAXOBJ=4,MAXMBR=LELT*6)<br>
PARAMETER (lhis=200) ! # of pts a proc reads from <a href="http://hpts.in">hpts.in</a><br> ! Note: lhis*np > npoints in <a href="http://hpts.in">hpts.in</a><br>C<br>C Common Block Dimensions<br>
C<br> PARAMETER (LCTMP0 =2*LX1*LY1*LZ1*LELT)<br> PARAMETER (LCTMP1 =4*LX1*LY1*LZ1*LELT)<br>C<br>C The parameter LVEC controls whether an additional 42 field arrays<br>C are required for Steady State Solutions. If you are not using<br>
C Steady State, it is recommended that LVEC=1.<br>C<br> PARAMETER (LVEC=1)<br>C<br>C Uzawa projection array dimensions<br>C<br> parameter (mxprev = 20)<br> parameter (lgmres = 30)<br>C<br>C Split projection array dimensions<br>
C<br> parameter(lmvec = 1)<br> parameter(lsvec = 1)<br> parameter(lstore=lmvec*lsvec)<br>c<br>c NONCONFORMING STUFF<br>c<br> parameter (maxmor = lelt)<br>C<br>C Array dimensions<br>C<br> COMMON/DIMN/NELV,NELT,NX1,NY1,NZ1,NX2,NY2,NZ2<br>
$,NX3,NY3,NZ3,NDIM,NFIELD,NPERT,NID<br> $,NXD,NYD,NZD<br><br>c automatically added by makenek<br> parameter(lxo = lx1) ! max output grid size (lxo>=lx1)<br><br>c automatically added by makenek<br> parameter(lpart = 1 ) ! max number of particles<br>
<br>c automatically added by makenek<br> integer ax1,ay1,az1,ax2,ay2,az2<br> parameter (ax1=lx1,ay1=ly1,az1=lz1,ax2=lx2,ay2=ly2,az2=lz2) ! running averages<br><br>c automatically added by makenek<br> parameter (lxs=1,lys=lxs,lzs=(lxs-1)*(ldim-2)+1) !New Pressure Preconditioner<br>
<br>c automatically added by makenek<br> parameter (lfdm=0) ! == 1 for fast diagonalization method<br><br></div>Is anyone has idea what is the problem in this case?!<br><br></div>Thanks,<br>Ami<br><div><div><div><br>
</div></div></div></div>