RICHARDSON
heading goes here.
RICHARDSON
is a FORTRAN77 (program/library) which
(description goes here).
Usage:
Licensing:
The computer code and data files described and made available on this web page
are distributed under
the GNU LGPL license.
Related Data and Programs:
Reference:
-
Peter Lynch,
The Emergence of Numerical Weather Prediction: Richardson's Dream,
Cambridge University Press, 2006,
ISBN: 978-0521857291,
LC: QC995.L93.
Source Code:
Examples and Tests:
List of Routines:
-
DIMENSION VORT(IM,JM,LM)
-
COMMON / PGRAD / DPDX(IM,JM2,LM), DPDY(IM,JM2,LM)
-
real ZPLOT(IM,2*JM)
-
DATA WHOUR / 00 , 0.25, 06, 12, 24, 48, 4*0 /
-
OPEN(UNIT=NDISC,FILE='RICHARDSON.DAT',FORM='FORMATTED')
-
print*, ' READ in the INITIAL fields '
-
READ(NDISC,777) P,MU,MV,TTOP,OROG
-
PMID2(I,J,1) = P2(I,J,1) + ( RGAS*TTOP2(I,J)/GRAV )
-
DO 1002 I=3,IMM2
-
DO 1002 J=3,JM2M3
-
print9455, L, RMID
-
SURFAC = 0.2 !!! Now read in on a card.
-
cc WLO = W2(I,J-1,L)*RHO2I(I,J-1,L)
-
cc WHI = W2(I,J-1,L-1)*RHO2I(I,J-1,L)
-
cc WHI = W2(I,J-1,L-1)*RHO2I(I,J,L-1)
-
cc WLO = W2(I,J-1,L)*RHO2I(I,J-1,L)
-
YMIN = 0.
-
YMAX = 0.
-
CALL PLOTLN (TIME,w,0,M,YMIN,YMAX,' w ')
-
print*, (TIME2(nT), nT=0,2*M)
-
print*, (w2(nT), nT=0,2*M)
-
YMIN = 0.
-
YMAX = 0.
-
CALL PLOTLN (TIME2(0),w2(0),0,2*M,YMIN,YMAX,' w2 ')
-
CALL RESPONSE(w2,0, N+1,2*NMAX,HH,PI)
-
print98, ( ARRAY(I,J) , I=2,IMM1,2 )
-
print99, ( ARRAY(I,J) , I=1,IM ,2 )
-
print*,'P CHECK ', P(I,J,LM) , PS(I,J)
-
print*,' PTOZ:: ZDYNAM, PBAR, TBAR : ', ZL,PL,TL
-
print*,' PTOZ:: ZDYNAM, PBAR, TBAR : ', ZL,PL,TL
-
print*,'P CHECK ', P2(I,J,LM) , PSEA2(I,J)
-
print*,' Z1000:: PL,TL,HS,DZ:', PL,TL,HS,DZ
-
CCCCC these scaling factors make printed values
-
CCCCC "look the same size" as in LFR's table.
-
print*,'P CHECK ', PS2(I,J) , PSEA2(I,J)
-
From NUMERICAL RECIPES (Press, et al.)
-
DETERMINE CUTOFF FREQUENCY/WAVENUMBER.
-
PI = 2*ASIN(1.)
-
THETAC = 2*ATAN(EPSILON**(-1./6.))
-
XLENC = 2*PI/THETAC
-
TYPE *,' XLENC, THETAC, EPSILON ',XLENC,THETAC,EPSILON
You can go up one level to
the FORTRAN77 source codes.
Last revised on Thu Jan 14 11:08:34 2010