PARAMETER (mbig=1000000000,mz=0,fac=1.d0/mbig)

      real*8 mbig,mseed,mz

      PARAMETER (mbig=4000000.,mz=0.,fac=1.d0/mbig)

      IF (firstCall .OR. seed.GE.mz) THEN

        IF (seed.GE.mz) mseed = seed

          IF (mk.LT.mz) mk=mk+mbig

            IF (ma(i).LT.mz) ma(i)=ma(i)+mbig

      IF (seed.GE.mz) THEN

        IF (mj.LT.mz) mj=mj+mbig

     & + Uinflow*sin(mz*PI*(float(k)-0.5 _d 0)/float(Nr))

     & *sin(mz*PI*(float(k)-0.5 _d 0)/float(Nr))

      _RL mz,strat,kx

      mz=1.0 _d 0

       vertStructWst(K)=cos(mz*PI* (rC(K)/rF(Nr+1)) )

       kx=mz*2. _d 0*pi/400.0 _d 0