Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...

10.STIMA DEI PARAMETRI DELLA TURBOLENZA ATMOSFERICA—————————————————————⎯⎯⎯⎯⎯⎯⎯⎯————————APPENDICE 10.2La Subroutine FORTRAN SPE_MESA produce la determinazione dello spettro a Massima Entropia diun segnale x di Ndata. La routine princiaple e le routines secondarie sono scritte ispirandosi a Press eal., 1992.! ==================================================================Subroutine SPE_MESA(dt,Ndata,x,mpoli,Nf,F,Sx)! ==================================================================! Stima dello Spettro col Metodo di Massima Entropia!! Input:! dt = intervallo di tempo tra due misure consecutive (sec.)! Ndata = numero di dati della serie di misure! x = vettore delle misure! Mpoli = numero di poli p! Nf = numero di frequenze per cui si vuole lo spettro! F = vettore delle frequenze! Output:! Sx = vettore delle densità spettrali! da !Press e al., (1992)! ------------------------------------------------------------------REAL x(Ndata)REAL F(Nf),Sx(Nf)REAL dxx(mpoli)! Calcolo dei coefficienti del processo AR(Mpoli)Call MEMCOF(x,Ndata,mpoli,xmsx,dxx)!! Calcolo delle densità spettrali.DO i=1,Nffdt = F(i)*dtSx(i) = 2.*dt*EVLMEM(fdt,dxx,mpoli,xmsx)EndDORETURNEND!! ==================================================================SUBROUTINE MEMCOF(data,n,m,xms,d)! ==================================================================INTEGERm,n,MMAX,NMAXREALxms,d(m),data(n)PARAMETER (MMAX=100,NMAX=40000)REALwk1(NMAX),wk2(NMAX),wkm(MMAX)! ------------------------------------------------------------------IF (m.gt.MMAX.or.n.gt.NMAX) PAUSE 'workspace too small in memcof'p = 0.DO j=1,np = p + data(j)**2EndDOxms = p/nwk1(1) = data(1)wk2(n-1) = data(n)————————————————————————————————————————- 503 -