L(223)
    GROUP 1. Run title
TEXT(Polar Wind 360 Degree             
TITLE
  The 180 degree domain is extended over 360 degrees in this
  example.
INTEGER(XFLOW1,XFLOW2);XFLOW1=6;XFLOW2=15
 
    GROUP 3. X-direction grid specification
GRDPWR(X,20,2.0*PI,1.0)
 
    GROUP 4. Y-direction grid specification
GRDPWR(Y,20,2.5,1.)
 
    GROUP 7. Variables stored, solved & named
SOLVE(P1,U1,V1);SOLUTN( P1,Y,Y,Y,N,N,N )
 
    GROUP 9. Properties of the medium (or media)
RHO1=1.205;ENUL=1.8E-5/RHO1
 
    GROUP 13. Boundary conditions and special sources
XCYCLE=T
PATCH(PINF,CELL,1,NX,NY,NY,1,NZ,1,1)
COVAL(PINF,P1,1.0E6,0.0 )
COVAL(PINF,U1,ONLYMS, 0.0);COVAL(PINF,V1,ONLYMS, 0.0)
PATCH(UPOL,CELL,XFLOW1-1,XFLOW2,NY,NY,1,NZ,1,1)
COVAL(UPOL,U1,FIXVAL,SETSPEED)
PATCH(VPOL,CELL,XFLOW1,XFLOW2,NY-1,NY-1,1,NZ,1,1)
COVAL(VPOL,V1,FIXVAL,SETSPEED)
 
    GROUP 16. Termination of iterations
LITER(P1) = 12;LITER(U1)=1;LITER(V1)=1
 
    GROUP 17. Under-relaxation devices
RELAX(U1,FALSDT,0.15);RELAX(V1,FALSDT,0.15)
 
    GROUP 19. Special calls from EARTH to GROUND
   **Inlet region extends from pi/2 to 3pi/2 and V1=VIN*CosX ,which
     gives a negative value (-Y direction) in this region
POLRA=VIN
 
    GROUP 21. Print-out of variables
 
 
    GROUP 22. Spot-value print-out
IXMON=2;IYMON=12;IZMON=1
 
    GROUP 23. Field print-out and plot control
NYPRIN=NY/5