Encyclopaedia Index

```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
* 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

AZYV=GRND
NAMSAT=MOSG
* PLANT blocks begin
REGION(1,1,1,1) /IZ.LE.20
REGION(1,1,1,1) /IZ.GT.20.AND.IZ.LE.40
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
```