! NB - to print field at x(degN) and y(km) with all longitudes use code
! at GPRINT.48 and just zonally symm number use code at XSECT.95
&INPB
RNTAPE=0.0	!run identifier on history records
TSPD=96!192!96		!time steps per day (in run)
KITS=3		!number of initial short timesteps (in run)
TMEAN=41*250.0!41	!NL*basic state temperature (in run)
P0=1.E5		!surface pressure
TDISS=0.25	!diffusion timescale (in run)
RESTIM=0. !4.62963290414011,4.62964739451667,4.62966230909385,4.62968974521928,4.62974021527448,4.62983305700905,4.63000384207439,4.63031800240416,4.63089588983697,4.63195885147918,4.63391390941225,4.63750926835309,4.64411948388083,4.65626701811761,4.67857147270056,4.71946153071006,4.79421009814248,4.93014223828935,5.17500889927798,5.60867764707162,6.35409999369992,7.57186265964173,9.40504213040493,11.8508785451700,14.6372099692440,17.2934891221783,19.4288463053302,20.9222702033051,21.8679094007900,22.4295242646568,22.7504538756229 !newtonian cooling (for NL=31 and ZHI,ZLOW as in run)
DRAG=41*0.!41 !3.92151776545858,4.43657941852848,4.68147896395137,4.97090497023231,5.31767250872592,5.74002927350022,6.26487629805510,6.93358819324983,7.81326119631138,9.02038131689040,10.7765211990690,13.5613581043730,18.6429242879512,30.8267888278709,99.2299491552724,16*0.0 !TFRC (in run, for NL=31 and ZHI,ZLOW as in run)     
BEGDAY=0.0	!starting day (in run)
ENDDAY=100.!2920.!600.!365.!2455.!3370.!365.!600.	!final day for which data is to be processed
KOUNTD=96	!fields analysed every KOUNTD/TSPD days
		!(corresp to KOUNTH in run I think)
KOUNTP=96	!frequency of printed output (multiple of KOUNTD)
KOUNTF=0	!frequency of writing to film (multiple of KOUNTD)
KOUNTE=0	!frequency of output of KE spectra (multiple of KOUNTD)
KOUNTA=0	!frequency of writing to data transfer file (multiple of KOUNTD)
IZGTYP=1	!projection for lat-long plots. 4 = polar stereographic, 
		!1 = regular latitude-longitude.
LROG=.T.	!include orography?
		!CHECK THIS IS WHAT WE WANT!!
NCOEFF=0	!highest total wavenumber of spectral coeffs. printed.
LPC2=.F.	!plot fields on isobaric levels?
NPC2=61		!channel on which fields are printed
LSGEDY=.F.	!removes zonal mean from the fields
LSGGRZ=.F.	!prints zonal average of the fields
LGPO=17*.F.,.T.,23*.F. !12*.F.,.T.,18*.F.,.T.,7*.F.,.T.,.F.		!may need this to print more than just the surface geopotential...
LSGGR=4*.F.,.T.,7*.F.!7*.F.,.T.,4*.F. !variables to print
!use *with zonally symm part removed* to print GSHAPE over longitudes and 
!height.  Also use (see above) for printing zonal mean on selected levels.
!LSGPL		!variables to plot
SGFAC=7*0.1,100.!10.	!dimensional factor
SGINC=0.0!200.	!dimensional contour interval
LPC3=.T.	!plot zonally averaged fields?
NPC3=60		!channel on which zonal averages are printed
LXDGR=.T.,27*.F.! .T.,.F.,2*.T.,10*.F.,.T.,6*.F.,.T.,6*.F. !variables to be printed
LXDPL=28*.F. !variables to be plotted
XDFAC=0.1!1.E19!0.1!1.E19	!dimensional factor (from documentation)
XDINC=5.0!0.!5.0!-0.1!0.	!dimensional contour interval (from documentation)
LPC5=.F.	!plot global energetics?
LPC6=.F.	!plot fields on isentropic levels?
NPC6=60
THSURD=900.,750.,550.,300.,200. !1650.
CPHASE=0.
LTHGR=2*.F.,.T.,.F.
LTHPL=4*.F.
THFAC=2*.1,1.,1.
THINC=0.5
LPC8=.F.	!plot tendencies of physical processes?
&END

&INPHYS
LDRY=.T.
BEGDP=0.0
&END