cgxomeg.htm

c
c
c

file-name gxomeg.htm 011102

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 '/phoenics/d_includ/farray' INCLUDE '/phoenics/d_includ/satear' INCLUDE '/phoenics/d_includ/grdloc' INCLUDE '/phoenics/d_includ/satgrd' INCLUDE '/phoenics/d_includ/grdear' INCLUDE '/phoenics/d_includ/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 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