TALK=f;RUN( 1, 1) TEXT( Library case Y602: EXPANDING AND CONTRACTING THE GRID >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>> PLANT information : * Data input groups used: 4 * Ground groups planted : 4 * Headings used : SCYS?? * Functions used : None * Commands used : REGION <<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<< GROUP 1. Run title and other preliminaries REAL(POWER, HSLOT) HSLOT=0.1 GROUP 4. Y-direction grid specification NY=20;YVLAST=HSLOT;YFRAC(1)=-20.;YFRAC(2)=1.0/20. ** Power law grid expansion and contraction POWER=0.50;ZWADD=HSLOT AZYV=GRND NAMSAT=MOSG * PLANT blocks beginYRAT=1.+:POWER:*DZ/(ZWADD+ZW) REGION(1,1,1,1) /IZ.LE.20 YRAT=1.+:POWER:*(-1.)*DZ/(ZWADD+ZW) REGION(1,1,1,1) /IZ.GT.20.AND.IZ.LE.40 YRAT=1.+:POWER:*1.*DZ/(ZWADD+ZW) REGION(1,1,1,1) /IZ.GT.40 * PLANT blocks end >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>>>>>>>> The above three statements make the domain width a power law function of the downstream distance of the current slab (ZW). The power is changed with Z-direction in step-wise manner. It is positive over first 20 slabs, negative over next 20 slabs and positive again for the remainder of the domain. <<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<<<<<<<< GROUP 5. Z-direction grid specification PARAB=T ** The z-direction step size is made proportional to YVLAST by setting AZDZ=PROPY and DZW1 here NZ=100; AZDZ=PROPY; DZW1=0.05 GROUP 7. Variables stored, solved & named STORE(V1,W1) LITHYD=1 ** grid will be dumped in XYZ format file, PARXYZ. TSTSWP=-1;IDISPA=5;IDISPB=1;IDISPC=NZ GROUP 23. Field print-out and plot control PHOTON USE x; parxyz rot y ang 180 gr x 1 msg EXPANDING - CONTRACTING GRID, Y-Z msg Press E to end enduse STOP