<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p>Hi Stefan,</p>
<p><br>
</p>
<p>In case the attachments didn't go through, please see below for the .usr file and log file.<br>
</p>
<p>I only used 'userchk' and 'usrdat' in the .usr file for interpolation. The other subroutines are empty.<br>
</p>
<p>(I added several print lines in 'gfldr.f' to see where the interpolation halted.)</p>
<p><br>
</p>
<p>Best regards,</p>
<p>Tony<br>
</p>
<br>
<p>.USR FILE:</p>
<p></p>
<div>c-----------------------------------------------------------------------<br>
      subroutine userchk()<br>
      include 'SIZE'<br>
      include 'TOTAL'<br>
      include 'ZPER'<br>
      include 'CTIMER'<br>
      include 'RESTART'<br>
<br>
      character*132 sourcefld<br>
<br>
      sourcefld='source.f00001'<br>
      call gfldr(sourcefld)<br>
      call outpost(vx,vy,vz,pr,t,'rst')<br>
      call exitt()<br>
<br>
      return<br>
      end<br>
c-----------------------------------------------------------------------<br>
      subroutine usrdat()   ! This routine to modify element vertices<br>
      include 'SIZE'<br>
      include 'TOTAL'<br>
<br>
      parameter (ielx=150)<br>
      parameter (iely=22)<br>
      parameter (ielz=142)<br>
      integer i,e,eg,ex,ey,ez<br>
      real cony<br>
<br>
      cony = 2.36d0<br>
<br>
      do e=1,nelv<br>
        eg = lglel(e)<br>
       if(eg<=ielx*iely*ielz) then<br>
        call get_exyz(ex,ey,ez,eg,ielx,iely,ielz)<br>
<br>
        do i=1,8<br>
         if(i<=2 .or. i==5 .or. i==6) then<br>
          if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR.<br>
     &      (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then<br>
           yc(i,e) = 2.0e0*Real(ey-1)/Real(iely)<br>
          end if<br>
         else if(i==3 .or. i==4 .or. i>=7) then<br>
          if((yc(i,e)>0.0e0.and.yc(i,e)<1.0e0).OR.<br>
     &      (yc(i,e)>1.0e0.and.yc(i,e)<2.0e0)) then<br>
           yc(i,e) = 2.0e0*Real(ey)/Real(iely)<br>
          end if<br>
         end if<br>
        end do<br>
<br>
       end if<br>
      end do<br>
<br>
      do e=1,nelv<br>
        eg = lglel(e)<br>
        if(eg==1) print*,yc(1,e),yc(2,e),yc(5,e),yc(6,e)<br>
        if(eg==1) print*,yc(3,e),yc(4,e),yc(7,e),yc(8,e)<br>
        do i=1,8<br>
         if(yc(i,e)>0.0e0 .AND. yc(i,e)<1.0e0) then<br>
           yc(i,e) = yc(i,e)-1.0e0<br>
           yc(i,e) = tanh(cony*yc(i,e))/tanh(cony)<br>
           yc(i,e) = yc(i,e)+1.0e0<br>
         else if(yc(i,e)>1.0e0 .AND. yc(i,e)<2.0e0) then<br>
           yc(i,e) = yc(i,e)-1.0e0<br>
           yc(i,e) = tanh(cony*yc(i,e))/tanh(cony)<br>
           yc(i,e) = yc(i,e)+1.0e0<br>
         end if<br>
        end do<br>
      end do<br>
<br>
      return<br>
      end<br>
c-----------------------------------------------------------------------<br>
</div>
<p></p>
<br>
<br>
<p>LOG FILE:<br>
</p>
<br>
<div>/----------------------------------------------------------\\<br>
|      _   __ ______ __ __  ______  ____   ____   ____     |<br>
|     / | / // ____// //_/ / ____/ / __ \\/ __ \\/ __ \\   |<br>
|    /  |/ // __/  / ,<   /___ \\ / / / // / / // / / /    |<br>
|   / /|  // /___ / /| | ____/ / / /_/ // /_/ // /_/ /     |<br>
|  /_/ |_//_____//_/ |_|/_____/  \\___/ \\___/ \\___/      |<br>
|                                                          |<br>
|----------------------------------------------------------|<br>
|                                                          |<br>
| NEK5000:  Open Source Spectral Element Solver            |<br>
| COPYRIGHT (c) 2008-2017 UCHICAGO ARGONNE, LLC            |<br>
| Version:  17.0.0                                         |<br>
| Web:      http://nek5000.mcs.anl.gov                     |<br>
|                                                          |<br>
\\----------------------------------------------------------/<br>
<br>
<br>
 Number of processors:         720<br>
 REAL    wdsize      :           8<br>
 INTEGER wdsize      :           4<br>
 Timer accuracy      : 9.54E-08<br>
 Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.rea<br>
<br>
 mapping elements to processors<br>
 Reading /work/e01/e01/tony/geometry_nonosci_cavity/interp.map<br>
 element load imbalance:            1         652         653<br>
 done :: mapping elements to processors<br>
<br>
           0  objects found<br>
 118   Parameters from file:<br>
   1      1.00000     p001 DENSITY<br>
   2     -3150.00     p002 VISCOS<br>
  12     -0.11750E-01 p012 DT<br>
  15      500.000     p015 IOSTEP<br>
  21     0.100000E-07 p021 DIVERGENCE<br>
  22     0.100000E-07 p022 HELMHOLTZ<br>
  26      2.00000     p026 COURANT/NTAU<br>
  27      3.00000     p027 TORDER<br>
  59      0.00000     p059 !=0 --> full Jac. eval. for each el.<br>
  63      8.00000     p063 =8 --> force 8-byte output<br>
  65      1.00000     p065 #iofiles (eg, 0 or 64); <0 --> sep. dirs<br>
  66      6.00000     p066 output : <0=ascii, else binary<br>
  67      6.00000     p067 restart: <0=ascii, else binary<br>
  93      20.0000     p093 Number of previous pressure solns saved<br>
  94      5.00000     p094 start projecting velocity after p94 step<br>
  95      5.00000     p095 start projecting pressure after p95 step<br>
  99      3.00000     p099 dealiasing: <0--> off/3--> old/4--> new<br>
 103      0.10000E-01 p103 weight of stabilizing filter (.01)<br>
<br>
 nelgt/nelgv/lelt:      469500      469500         653<br>
 lx1  /lx2  /lx3 :           8           8           8<br>
 done :: read .rea file    45.222     sec<br>
<br>
 setup mesh topology<br>
   Right-handed check complete for  469500 elements. OK.<br>
   setvert3d:   8    60855000   162267000    60855000    60855000<br>
 call usrsetvert<br>
 done :: usrsetvert<br>
<br>
gs_setup: 10325502 unique labels shared<br>
   pairwise times (avg, min, max): 0.00044986 0.00033021 0.000500798<br>
   crystal router                : 0.00145524 0.00143301 0.0014781<br>
   used all_to_all method: pairwise<br>
   handle bytes (avg, min, max): 2.64212e+06 2463020 2796252<br>
   buffer bytes (avg, min, max): 491593 352288 632640<br>
   setupds time 5.3504E-01 seconds   0  8    60855000      469500<br>
           8  max multiplicity<br>
 done :: setup mesh topology<br>
<br>
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000<br>
  0.000000000000000E+000<br>
  9.090909090909091E-002  9.090909090909091E-002  9.090909090909091E-002<br>
  9.090909090909091E-002<br>
 call usrdat<br>
 done :: usrdat<br>
<br>
 generate geometry data<br>
 NOTE: All elements deformed , param(59) ^=0<br>
 done :: generate geometry data<br>
<br>
 call usrdat2<br>
 done :: usrdat2<br>
<br>
 regenerate geometry data           1<br>
 NOTE: All elements deformed , param(59) ^=0<br>
 done :: regenerate geometry data           1<br>
<br>
 verify mesh topology<br>
  0.000000000000000E+000   6.00000000000000       Xrange<br>
 -0.120000000000000        2.00000000000000       Yrange<br>
  0.000000000000000E+000   49.7000000000000       Zrange<br>
 done :: verify mesh topology<br>
<br>
 IFTRAN   = T<br>
 IFFLOW   = T<br>
 IFHEAT   = F<br>
 IFSPLIT  = T<br>
 IFLOMACH = F<br>
 IFUSERVP = F<br>
 IFUSERMV = F<br>
 IFSTRS   = F<br>
 IFCHAR   = T<br>
 IFCYCLIC = F<br>
 IFAXIS   = F<br>
 IFMVBD   = F<br>
 IFMELT   = F<br>
 IFMODEL  = F<br>
 IFKEPS   = F<br>
 IFMOAB   = F<br>
 IFNEKNEK = F<br>
 IFSYNC   = T<br>
<br>
 IFVCOR   = F<br>
 IFINTQ   = F<br>
 IFCWUZ   = F<br>
 IFSWALL  = F<br>
 IFGEOM   = F<br>
 IFSURT   = F<br>
 IFWCNO   = F<br>
<br>
 IFTMSH for field           1    =  F<br>
 IFADVC for field           1    =  T<br>
 IFNONL for field           1    =  F<br>
<br>
 Dealiasing enabled, lxd=          12<br>
 verify mesh topology<br>
  0.000000000000000E+000   6.00000000000000       Xrange<br>
 -0.120000000000000        2.00000000000000       Yrange<br>
  0.000000000000000E+000   49.7000000000000       Zrange<br>
 done :: verify mesh topology<br>
<br>
 call usrdat3<br>
 done :: usrdat3<br>
</div>
<br>
<div> set initial conditions<br>
 nekuic (1) for ifld            1<br>
 call nekuic for vel<br>
 xyz min     0.0000     -0.12000       0.0000<br>
 uvwpt min   0.0000       0.0000       0.0000       0.0000       0.0000<br>
 xyz max     6.0000       2.0000       49.700<br>
 uvwpt max   0.0000       0.0000       0.0000       0.0000       0.0000<br>
 done :: set initial conditions<br>
<br>
gridpoints unique/tot:     162267000    240384000<br>
  dofs:             159972750             162104250<br>
<br>
 Initialization successfully completed    46.588     sec<br>
<br>
 nsteps=0 -> skip time loop<br>
 running solver in post processing mode<br>
<br>
 call userchk<br>
 call gfldr<br>
 byte swap: F   6.543210     -2.9312772E+35<br>
 sourcefld X  0.000000000000000E+000   6.00000000000000<br>
 sourcefld Y -0.120000000000000        2.00000000000000<br>
 sourcefld Z  0.000000000000000E+000   49.7000000000000<br>
 CHECK interp prepare<br>
 CHECK interp read<br>
<br>
</div>
<div style="color: rgb(19, 19, 18);">
<div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" color="#000000" face="Calibri, sans-serif"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users-request@lists.mcs.anl.gov <nek5000-users-request@lists.mcs.anl.gov><br>
<b>Sent:</b> 28 March 2017 20:51<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Nek5000-users Digest, Vol 97, Issue 33</font> </div>
</div>
<font size="2"><span style="font-size:10pt;"></span></font><font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
Message: 3<br>
Date: Tue, 28 Mar 2017 19:51:04 +0000<br>
From: nek5000-users@lists.mcs.anl.gov<br>
To: "nek5000-users@lists.mcs.anl.gov"<br>
        <nek5000-users@lists.mcs.anl.gov><br>
Subject: Re: [Nek5000-users] Interpolation between two meshes<br>
Message-ID:<br>
        <mailman.3915.1490730701.2967.nek5000-users@lists.mcs.anl.gov><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi Stefan,<br>
<br>
<br>
Many thanks for your reply. I'm running my simulations on ARCHER (Cray XC30 supercomputer) and I'm using the Intel compiler.<br>
<br>
<br>
The current simulation has 469500 elements with polynomial order of 7. The source field I used for interpolation has 383700 elements with the same polynomial order. The size of the source field is about 10GB. I'm not sure if it can be easily transferred. But
 I have attached my .usr and log files herewith (also the gfldr.f with few added checking lines). I ran the interpolation using 720 processors.<br>
<br>
<br>
Please have a look and let me know if you need further information about the case. Thank you very much for you help.<br>
<br>
<br>
Best regards,<br>
<br>
Tony<br>
</div>
</span></font></div>
</div>
</body>
</html>