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:

1. Peter Lynch,
   The Emergence of Numerical Weather Prediction: Richardson's Dream,
   Cambridge University Press, 2006,
   ISBN: 978-0521857291,
   LC: QC995.L93.

Source Code:

• richardson.f, the source code.
• richardson.csh, commands to compile the source code.

Examples and Tests:

• richardson_prb.f, a sample calling program.
• richardson_prb.csh, commands to compile and run the sample program.
• richardson_prb_output.txt, the output file.

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.