TALK=T;RUN(1,1)
  PHOTON USE
  p



  GR OU Y     2 X     1     1 Z     3     3 COL     2
  GR OU Y     1 X     1     1 Z     3     3 COL     2
  GR OU Y     3 X     1     1 Z     4     4 COL     3
  GR OU Y     1 X     1     1 Z     4     4 COL     3
  GR OU Y     4 X     1     1 Z     5     5 COL     4
  GR OU Y     1 X     1     1 Z     5     5 COL     4
  GR OU Y     5 X     1     1 Z     6    13 COL     5
  GR OU Y     1 X     1     1 Z     6    13 COL     5
  GR OU Y     7 X     1     1 Z    14    14 COL     6
  GR OU Y     1 X     1     1 Z    14    14 COL     6
  GR OU Y     8 X     1     1 Z    15    17 COL     7
  GR OU Y     1 X     1     1 Z    15    17 COL     7
  GR OU Y     7 X     1     1 Z    18    19 COL     8
  GR OU Y     1 X     1     1 Z    18    19 COL     8
  GR OU Y     6 X     1     1 Z    20    30 COL     9
  GR OU Y     1 X     1     1 Z    20    30 COL     9
  GR OU Y     5 X     1     1 Z    31    35 COL    10
  GR OU Y     1 X     1     1 Z    31    35 COL    10
  GR OU Y     4 X     1     1 Z    36    36 COL    11
  GR OU Y     1 X     1     1 Z    36    36 COL    11
  GR OU X     2 Y    10    16 Z    16    16 COL     5
  GR OU X     1 Y    10    16 Z    16    16 COL     5
  GR OU Y    17 X     1     1 Z    16    16 COL     5
  GR OU Y    10 X     1     1 Z    16    16 COL     5
  GR OU X     2 Y     9    14 Z    29    29 COL     6
  GR OU X     1 Y     9    14 Z    29    29 COL     6
  GR OU Y    15 X     1     1 Z    29    29 COL     6
  GR OU Y     9 X     1     1 Z    29    29 COL     6
  GR OU X 1; GR OU X M; GR OU Y 1; GR OU Y M
  red
  msg                 ON-LINE WIND TURBINES
  msg                      Geometry
  msg Press Enter to continue
  pause
  gr OFF;red
  msg                 ON-LINE WIND TURBINES
  MSG                         Grid
  GR X 1
  msg Press Enter to continue
  pause
  gr OFF;red
  gr ou x 1
  msg                 ON-LINE WIND TURBINES
  msg                    Velocity field
  VE X 1 SH
  msg Press Enter to continue
  pause
  ve OFF;red
  msg                 ON-LINE WIND TURBINES
  msg           Streamwise velocity contours
  con w1 x 1 sh;int 10
  msg Press Enter to continue
  pause
  con OFF;red
  msg                 ON-LINE WIND TURBINES
  msg               Turbulence energy contours
  con ke x 1 fil;0.01
  msg Press E to end
  ENDUSE
  DISPLAY
   This run  introduces the influence of the stand of widely-spaced
   trees on the wake fields generated by two on-line wind  turbines
   located  on  the complex terrain of Valaam island of Ladoga lake
   in Nothern Russia.  The standard K-E model is modified  to  take
   into   account  the  plant  airflow  interactions  resulting  in
   momentum sinks. The corresponding GROUND codings are made by way
   of  PLANT  option.  The  main  details  can be found in paper by
   S.Zhubrin  et  al."Practical  evalution  of  turbulence   models
   applied   to  on-line  wind  turbines  on  a  complex  terrain",
   presented  at  International  Seminar  of  AirFlows  on  Complex
   Terrain, University of Karlsrue, February 1995
  ENDDIS
  PLANTBEGIN
   * Plant airflow interactions
     - form darg
PATCH(INTER,PHASEM,1,1,5,9,7,12,1,1)
    CO=0.5*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,V1,GRND,0.0)
    CO=0.5*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,W1,GRND,0.0)
     - wake production
    VAL=0.5*0.6*10.5*(V1*V1+W1*W1)**1.5
COVAL(INTER,KE,fixflu,GRND)
     - turbulence energy loss
    CO=2.*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,KE,GRND,0.0)
     - dissipation production
    VAL=0.75*EP/KE*0.6*10.5*(V1*V1+W1*W1)**1.5
COVAL(INTER,EP,fixflu,GRND)
     - dissipation loss
    CO=3.*0.6*10.5*(V1*V1+W1*W1)**0.5
COVAL(INTER,EP,GRND,0.0)
   PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(ON-LINE WIND TURBINES:144               )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 102
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 25
 YVLAST =470.
 YFRAC(1)=0.04255 ;YFRAC(2)=0.06383
 YFRAC(3)=0.08511 ;YFRAC(4)=0.1064
 YFRAC(5)=0.1277 ;YFRAC(6)=0.1383
 YFRAC(7)=0.1596 ;YFRAC(8)=0.1809
 YFRAC(9)=0.2234 ;YFRAC(10)=0.2319
 YFRAC(11)=0.2404 ;YFRAC(12)=0.2489
 YFRAC(13)=0.2574 ;YFRAC(14)=0.266
 YFRAC(15)=0.2819 ;YFRAC(16)=0.2979
 YFRAC(17)=0.3617 ;YFRAC(18)=0.4415
 YFRAC(19)=0.5213 ;YFRAC(20)=0.6011
 YFRAC(21)=0.6809 ;YFRAC(22)=0.7606
 YFRAC(23)=0.8404 ;YFRAC(24)=0.9202
 YFRAC(25)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 36
 ZWLAST =1300.
 ZFRAC(1)=0.03077 ;ZFRAC(2)=0.06154
 ZFRAC(3)=0.08462 ;ZFRAC(4)=0.1154
 ZFRAC(5)=0.1538 ;ZFRAC(6)=0.1923
 ZFRAC(7)=0.2308 ;ZFRAC(8)=0.2692
 ZFRAC(9)=0.3077 ;ZFRAC(10)=0.3462
 ZFRAC(11)=0.3846 ;ZFRAC(12)=0.4231
 ZFRAC(13)=0.4615 ;ZFRAC(14)=0.4846
 ZFRAC(15)=0.5077 ;ZFRAC(16)=0.5154
 ZFRAC(17)=0.5231 ;ZFRAC(18)=0.5385
 ZFRAC(19)=0.5538 ;ZFRAC(20)=0.5692
 ZFRAC(21)=0.5846 ;ZFRAC(22)=0.6
 ZFRAC(23)=0.6154 ;ZFRAC(24)=0.6308
 ZFRAC(25)=0.6462 ;ZFRAC(26)=0.6615
 ZFRAC(27)=0.6769 ;ZFRAC(28)=0.6923
 ZFRAC(29)=0.7 ;ZFRAC(30)=0.7308
 ZFRAC(31)=0.7692 ;ZFRAC(32)=0.8077
 ZFRAC(33)=0.8462 ;ZFRAC(34)=0.8846
 ZFRAC(35)=0.9231 ;ZFRAC(36)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(149)=PRPS
 NAME(150)=ENUT
    * Y in SOLUTN argument list denotes:
    * 1-stored 2-solved 3-whole-field
    * 4-point-by-point 5-explicit 6-harmonic averaging 
 SOLUTN(P1,Y,Y,Y,N,N,N)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 VIST = 150
 PRPS = 149
 ************************************************************
  Group 8. Terms & Devices
    * Y in TERMS argument list denotes:
    * 1-built-in source 2-convection 3-diffusion 4-transient
    * 5-first phase variable 6-interphase transport         
 TERMS(P1,Y,Y,Y,N,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENT = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1.22 ;TMP1 =0.
 EL1 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 ENUL =1.465E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(W1)=4. ;FIINIT(KE)=1.
 FIINIT(EP)=3.0E-04 ;FIINIT(PRPS)=-1.
 FIINIT(ENUT)=1.0E-10
 
 PATCH(CMP1 ,INIVAL, 1, 1, 1, 1, 3, 3, 1, 1)
 INIT(CMP1 ,PRPS,0. ,198. )
 
 PATCH(CMP2 ,INIVAL, 1, 1, 1, 2, 4, 4, 1, 1)
 INIT(CMP2 ,PRPS,0. ,198. )
 
 PATCH(CMP3 ,INIVAL, 1, 1, 1, 3, 5, 5, 1, 1)
 INIT(CMP3 ,PRPS,0. ,198. )
 
 PATCH(CMP4 ,INIVAL, 1, 1, 1, 4, 6, 13, 1, 1)
 INIT(CMP4 ,PRPS,0. ,198. )
 
 PATCH(CMP5 ,INIVAL, 1, 1, 1, 6, 14, 14, 1, 1)
 INIT(CMP5 ,PRPS,0. ,198. )
 
 PATCH(CMP6 ,INIVAL, 1, 1, 1, 7, 15, 17, 1, 1)
 INIT(CMP6 ,PRPS,0. ,198. )
 
 PATCH(CMP7 ,INIVAL, 1, 1, 1, 6, 18, 19, 1, 1)
 INIT(CMP7 ,PRPS,0. ,198. )
 
 PATCH(CMP8 ,INIVAL, 1, 1, 1, 5, 20, 30, 1, 1)
 INIT(CMP8 ,PRPS,0. ,198. )
 
 PATCH(CMP9 ,INIVAL, 1, 1, 1, 4, 31, 35, 1, 1)
 INIT(CMP9 ,PRPS,0. ,198. )
 
 PATCH(CMP10 ,INIVAL, 1, 1, 1, 3, 36, 36, 1, 1)
 INIT(CMP10 ,PRPS,0. ,198. )
 INIADD = F
 FSWEEP = 1
 NAMFI =CHAM
 ************************************************************
  Group 12. Patchwise adjustment of terms
  Patches for this group are printed with those
  for Group 13.
  Their names begin either with GP12 or &
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(INL ,LOW , 1, 1, 1, 25, 1, 1, 1, 1)
 COVAL(INL ,P1 , FIXFLU ,4.88 )
 COVAL(INL ,W1 ,0. ,4. )
 COVAL(INL ,KE ,0. ,1. )
 COVAL(INL ,EP ,0. ,3.0E-04 )
 
 PATCH(WAKE1 ,PHASEM, 1, 1, 10, 16, 16, 16, 1, 1)
 COVAL(WAKE1 ,W1 ,-0.6 ,0. )
 
 PATCH(WAKE2 ,PHASEM, 1, 1, 9, 14, 29, 29, 1, 1)
 COVAL(WAKE2 ,W1 ,-0.6 ,0. )
 
 PATCH(FREE ,NORTH , 1, 1, 25, 25, 1, 36, 1, 1)
 COVAL(FREE ,P1 ,1. ,0. )
 COVAL(FREE ,V1 ,0. , SAME )
 COVAL(FREE ,W1 ,0. , SAME )
 COVAL(FREE ,KE ,0. , SAME )
 COVAL(FREE ,EP ,0. , SAME )
 
 PATCH(OUT ,HIGH , 1, 1, 4, 25, 36, 36, 1, 1)
 COVAL(OUT ,P1 ,1000. ,0. )
 COVAL(OUT ,V1 ,0. , SAME )
 COVAL(OUT ,W1 ,0. , SAME )
 COVAL(OUT ,KE ,0. , SAME )
 COVAL(OUT ,EP ,0. , SAME )
 
 PATCH(INTER ,PHASEM, 1, 1, 5, 9, 7, 12, 1, 1)
 COVAL(INTER ,V1 , GRND ,0. )
 COVAL(INTER ,W1 , GRND ,0. )
 COVAL(INTER ,KE , GRND ,0. )
 COVAL(INTER ,EP , GRND ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 250 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(EP)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.4)
 RELAX(V1,FALSDT,3.333333)
 RELAX(W1,FALSDT,3.333333)
 RELAX(KE,FALSDT,3.333333)
 RELAX(EP,FALSDT,3.333333)
 RELAX(PRPS,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 KELIN = 0
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=1.96 ;EX(V1)=0.3132
 EX(W1)=2.452 ;EX(KE)=1.561
 EX(EP)=8.638E-03 ;EX(PRPS)=34.93
 EX(ENUT)=47.310001
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = F
    * Y in OUTPUT argument list denotes:
    * 1-field 2-correction-eq. monitor 3-selective dumping      
    * 4-whole-field residual 5-spot-value table 6-residual table
 OUTPUT(P1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 12 ;IZMON = 23
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP