c
C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
c
C SUBROUTINE GXOMEG provides for the calculation of the vorticity
c
C parameter in the MMK K-E turbulence Model:
C
C F = OMEG/S = 2*SQRT[ (O12**2 + O13**2 + O23**2)/GEN1 ]
C
C where GEN1 = 2.Sij*Sij with Sij=0.5*(Ui,j+Uj,i)
C and Oij=0.5*(Ui,j-Uj,i).
C
C F is used as a multiplier in the expression for ENUT, i.e.
C ENUT=CMUCD*F*KE**2/EP.
C
C The Subroutine is called from Group 19 Section 4 of the EARTH
C subroutine GR.
C-------------------------------------------------------------------
SUBROUTINE GXOMEG
INCLUDE 'farray'
INCLUDE 'satear'
INCLUDE 'grdloc'
INCLUDE 'satgrd'
INCLUDE 'grdear'
INCLUDE 'grdbfc'
COMMON/GENI/NXNY,NXM1NY,IGFIL1(7),NFM,IGFIL2(50)
COMMON/GENL/LGFIL1(14),DEBGTZ,LGFIL2(45)
COMMON/MMKKE/LBFOMG,L0FOMG
COMMON/GENFL/ LDUDX,LDUDY,LDUDZ,LDVDX,LDVDY,LDVDZ,LDWDX,LDWDY,
1 LDWDZ,LGNFIL
1 /LBDFDL/IDUDX,IDUDY,IDUDZ,IDVDX,IDVDY,IDVDZ,IDWDX,IDWDY,
1 IDWDZ,IDSDX,IDSDY,IDSDZ,IDU2X,IDU2Y,IDU2Z,IDV2X,
1 IDV2Y,IDV2Z,IDW2X,IDW2Y,IDW2Z
COMMON/NAMFN/NAMFUN,NAMSUB
CHARACTER*6 NAMFUN,NAMSUB
LOGICAL LGFIL1,DEBGTZ,LGFIL2
LOGICAL LDUDX,LDUDY,LDUDZ,LDVDX,LDVDY,LDVDZ,LDWDX,LDWDY,LDWDZ,
1 LGNFIL
NAMSUB='GXOMEG'
C IF(IGR.EQ.1.AND.ISC.EQ.1) THEN
C allocate storage
C ELSEIF(IGR.EQ.19.AND.ISC.EQ.4) THEN
CALL FN1(LBFOMG,0.0)
C.... omega,xy = 0.5*(DU/DR - DV/DX)
IF(LDUDY.OR.LDVDX) THEN
CALL FN1(EASP1,0.0)
IF(LDUDY) CALL FN34(EASP1,IDUDY,0.5)
IF(LDVDX) CALL FN34(EASP1,IDVDX,-0.5)
CALL FN51(EASP1,2.0)
CALL FN34(LBFOMG,EASP1,1.0)
ENDIF
C.... omega,xz = 0.5*(DU/DZ - DW/DX)
IF(LDUDZ.OR.LDWDX) THEN
CALL FN1(EASP1,0.0)
c call writ1i('isweep ',isweep)
IF(LDUDZ) CALL FN34(EASP1,IDUDZ,0.5)
IF(LDWDX) CALL FN34(EASP1,IDWDX,-0.5)
CALL FN51(EASP1,2.0)
CALL FN34(LBFOMG,EASP1,1.0)
ENDIF
C.... omega,yz = 0.5*(DV/DZ - DW/DY)
IF(LDVDZ.OR.LDWDY) THEN
CALL FN1(EASP1,0.0)
IF(LDVDZ) CALL FN34(EASP1,IDVDZ,0.5)
IF(LDWDY) CALL FN34(EASP1,IDWDY,-0.5)
CALL FN51(EASP1,2.0)
CALL FN34(LBFOMG,EASP1,1.0)
ENDIF
CALL FN25(LBFOMG,4.0)
CALL FN27(LBFOMG,GEN1)
CALL FN30(LBFOMG)
IF(IENUTA.EQ.12) CALL FN23(LBFOMG,1.0)
NAMSUB = 'gxomeg'
END
c