SUBROUTINE UPDATE CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C In this revised version (February, 1998), the Newtonian C C damping parameter for the tangential wind (toward a C C reference value) at each gridpoint is computed with C C empirical formulas. The damping is not made above SSTDX C C sigma level, beyond the radius RSTDX and in the planetary C C boundary layer. C C Time integration period for the vortex generation is C C 60 hours, including one hour free-up period. C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC c ccc USE setparms c include 'RESOLUTION.h' include 'PARAMETERS.h' PARAMETER (ima=250) PARAMETER (NP=26, nmx=24) PARAMETER (IBLK95=ima*20 + kmax*ima*4) COMMON /PSTOWE/ P1,P2 COMMON /FACT/ HURR,DEEP COMMON /NING/ TNGAS,NWRA COMMON /FLEXK/ B(KMAX),BI(KMAX),SIGT,KTOP,KTOPM,KTOPP COMMON /SURFC/ TSMEAN,WTMEAN,WSIN,PC,PMIN,IOFF COMMON /TANG/ VT(KMAX,IMX),TBAR(KMAX),RBAR(KMAX),PN,F * ,PMXMN,IENV COMMON /NEWF/ typ1f(ima,20),typ2f(kmax,ima,4) COMMON /POSIT/ XOLD,YOLD,XCORN,YCORN,RO,XV,YV,ROMM,TSIZE COMMON /SM/ CRW(NP),RC(NP),VC(NP),RR(IMX),VCR(IMX),WC(IMX) COMMON /INGCUT/ TNGCUT,VMXCUT COMMON /MAXVAL/ RMAX, VMAX, SGMST, SCAPCT COMMON /STRMSZ/ TCHGHT, TCRDUS COMMON /TIMET/ AFACT(KMAX,IMX),AA(KMAX,IMX) common /storu/ stpsf(imx),stdlf(imx) common /storet/ FK(KMAX),XK(IMX) C C include 'BKINFO.h' include 'COMMUNICATE.h' include 'GDINFO.h' include 'LIMIT.h' include 'CONMLEV.h' include 'CONSLEV.h' include 'FILEIF.h' include 'FLAGS.h' include 'QLOGS.h' include 'TIME.h' CC DIMENSION ro(nmx),typ1c(imx,lgi),typ2c(kmax,imx,nwat3) DIMENSION TM(KMAX,IMX),TN(KMAX,IMX),RLPS(IMX),RPS(IMX) DIMENSION RT(KMAX,IMX),FILC(IBLK95) DIMENSION TABL(8),TABS(8),TAXL(9),TAXS(9) CCC DIMENSION SGMS(KMAX) CCC EQUIVALENCE (typ1f(1,1),FILC(1)),(stpf(1),typ1c(1,1)) equivalence (uf(1,1),typ2c(1,1,1) ) CC DATA TABL /1.E30,1500.,800.,300.,150.,40.,0.,0./ DATA TABS /.15,.2,.3,.4,.5,.7,.85,1./ DATA TAXL /1.E-20,5.,25.,50.,100.,200.,400.,800.,1000./ DATA TAXS /8.,15.,20.,25.,30.,35.,40.,56.,200./ C do i = 1 , imx stpsf(i) = 0.0 stdlf(i) = 0.0 enddo do i = i1 , i2 stpsf(i) = psf(i) stdlf(i) = dlsf(i) enddo C call MPP_SUM_LOCAL(stpsf,imx) call MPP_SUM_LOCAL(stdlf,imx) C do i = 1 , imax psf(i) = stpsf(i) dlsf(i) = stdlf(i) enddo C CC REMOVE THE FREEING UP PERIOD FOR A TROPICAL STORM INTENSITY CC CC TMFT = TNGAS - 1.0 CC IF(HURR.EQ.0.0)TMFT = TNGAS CC TOW = TMFT*3600 C CCC Period of freeing up in the time integration (following C forced damping) is one hour regardless of intensity. C HTOW = TNGAS HTOW0 = TNGCUT TOW = HTOW*3600. TOW0 = HTOW0*3600. GGCNST = TOW0/TOW CCC PSUM = 0.0 IF(NSTEP.LE.1)THEN P2 = stpsf(IMAX) P1 = stpsf(IMAX-1) ENDIF CC CC********************* HOUR = TTIME(1) + DELTAT CC IF(HOUR.LE.TOW)THEN CC GG = EXP(GGCNST - (TOW0/HOUR)) CCC ELSE GG = 1.0 ENDIF CC PSUM =(P2-psf(IMAX))*dlsf(IMAX)+dlsf(IMAX-1)*(P1-psf(IMAX-1)) ASUM = dlsf(IMAX)+dlsf(IMAX-1) PCORR = PSUM/ASUM DO 83 I = i1 , i2 cccc psf(I) = psf(I) + PCORR 83 CONTINUE CCC ====================================================== CCC Specify a functional form for Newtonian damping IF(NSTEP-1)5005,5005,5100 5005 CONTINUE DO 5010 K = 1 , KMAX SGMS(K) = Q(K) 5010 CONTINUE CCC First, sigma-dependent part FK(K). TC height TCHGHT(=SSTD) CCC and two constants SAMP and SDEL are used. No damping CCC for SGMS