arpi3d pc/t3d type input file.....*******FOR F90 CODE********
doz45t  ! runname 
0       ! restart: =0 no restart, =1 restart run
0       ! detrend: 0= no detrending, 1=detrending(pressure) NOT ACTIVE
1       ! cordopt: 0=ARPS trans, 1=SIMPLE transformation SIMPLE ONLY
3       ! dimopt  dimension opt...
3       ! mpiopt  Option for mpi code, =0 no mpi, =1 x only, =2, y, =3 x,y
1       ! passopt MPI passing option, =0 single pass, =1 opt. passing model
4       ! nprocs0 Total number of processors nprocs = nprocx*nprocy
2       ! nprocx  Number of processors in the x direction.
2       ! nprocy  Number of processors in the y direction.
0       ! ngridopt Nested grid option, =0 no nesting, =1 nesting...
1       ! nlevel  Number of nest levels (including base grid).
1       ! ngrid   Number of grid nests in each level. (one per nlevel)
1       ! xratio nested grid dx ratio...
1       ! yratio nested grid dy ratio...
1       ! zratio nested grid dz ratio...
105     ! nx for the first (coarse grid) NOTE: add additional nx,ny,nz's   
55     ! ny for the first (coarse grid)       for more than 1 grid...   
83      ! nz for the first (coarse grid)  
0       ! interopt Interpol. option, =0 none, = 1 inter. from outer grid
0.0     ! x0 x coordinate of the southwest corner of model grid in meters...
0.0     ! y0 y coordinate of the southwest corner of model grid in meters...
2       ! inibasopt:  = 2 isent, =3 isothermal, = 4, constant static stab.
1       ! pt0opt: 1=bubble,2=horiz flds of random pert,3=symm random pert
3.0     ! ptpert0: Magnitude of initial pot temp perturbation (K)
2000.   ! pt0radx: Radius of initial bubble; infinity if pt0radx < 0
2000.   ! pt0rady: Radius of initial bubble; infinity if pt0rady < 0
1000.0  ! pt0radz: Radius of initial bubble; infinity if pt0radz < 0
25125.  ! pt0ctrx: x-coordinate of bubble center (m)
13125.  ! pt0ctry: y-coordinate of bubble center (m)
750.0   ! pt0ctrz: z-coordinate of bubble center (m)
0.0     ! ubarc base sate u wind.
0.0     ! vbarc base sate v wind.
0.0     ! wbarc base sate w wind.
250.0   ! t0 base state temperature.
1.096e-4 ! nsqrd static stability  (N).
300.0   ! pt0 base state potential temperature.
0       ! teropt: 0=flat, 1=bell mountain, 2=terrain data,6= mpi ter
0       ! ridgeopt: 0=no 2-d ridge in 3-d run, 1=2-d ridge y=0,2=centered 
500.0   ! hmount:  height of bell-shaped mountain (m)
10000.  ! mtnctrx: x-coordinate of terrain center (m)
10000.  ! mtnctry: y-coordinate of terrain center (m)
10000.0 ! mtnwidx: half-width of mount in x dir.; infinity if mntwidx < 0
10000.0 ! mtnwidy: half-width of mount in y dir.; infinity if mntwidy < 0
1.0     ! Ridgf ridge length is ridgf*mtnwidy
1       ! soundopt init.option,=0 analytical,=1 sounding,=2 restart file.
0       ! sndtype  = 0 arps, = 1 nws-gempack...
'doz4.snd' ! sndfile Name of file containing the base state sounding.
'doz701.rst007200g01' ! rstfile Name of file containing the restart data dump.
250.0   ! dx: Grid spacing in the x-direction (m)
250.0   ! dy: Grid spacing in the x-direction (m)
250.0   ! dz: Average grid spacing in z-dir in stretched or non-str. case 
0.0     ! zrefsfc Height of the surface ....
13000.  ! zflat   height of the level comp surface...
0       ! stretch option = 0, no stretch, = 1 stretch use pgradopt=1!
10.0    ! dzmin
600.0   ! dlayer1
100000. ! dlayer2
1.0     ! strhtune
0       ! initstrt slow startup: =0 none, =1 ubarc, =2 gravity.
0       ! bousopt =0 solve full equation set, =1 bous approx. applied
0       ! plinopt =0 pressure grad term full strength, =1 linearized
1       ! solvwp  =0 explicit solver, =1 vertical implicit solver...
0       ! pgradopt=0 coord trans pgrad, =1 const z inter. for stretch!
0.6     ! beta    =0.6 im, =1.0 EXP.  time average of p for wp solver.
1       ! totpdiv =0 coef. div linearized, =1 full strength.
1       ! psopt   =0 heating source term off, =1 included in p-eq.
0       ! vmixopt =0 explicit vertical tmix, =1 implicit tmixing.
0.5     ! crank1  = value 0<crank1<1.0 recommended =0.5 for 2nd order.
2.5     ! dtbig  Big timestep for non-acoustic modes (seconds)
0.5     ! dtsml   Small timestep for acoustic modes (seconds)
00.0    !  tstart  start time of the model run
7200.0  ! tstop:  Termination time of model integration (seconds)
1       ! hadvopt = 0 2nd order horz., = 1 fourth order adv. (ONLY)
1       ! vadvopt = 0 2nd order vert., = 1 fourth order adv.
0       ! coropt  = 0 no coriolis, =1 unbalanced cor, =2 balan cor.
43.2    ! flat    Latitude used to set the coriolis force f.
1       ! moist  water option  = 0 dry run, = 1 moisture on (Kessler).
0       ! rhopt  set the rh at a specified value everywhere.
70.0    ! rhval  specified rh in percent used for testing.
0       ! qvsfc     =0 no surface qv flux, = 1 constant sfc qvflux
0.0000  ! qvsfcflx  value for the constant qv sfc flux option..
0       ! rainopt   option for including/generating rain =0 no, =1 yes.
1       ! dozhell  option for const. heat at the surface =0 no, =1 yes.
5       ! ebc  lat bc 1=kw model,3=zero grad,4=rad(ORL),5=env,6=wall
5       ! wbc  lat bc 1=kw model,3=zero grad,4=rad(ORL),5=env,6=wall
5       ! nbc  lat bc 1=kw model,3=zero grad,4=rad(ORL),5=env,6=wall
5       ! sbc  lat bc 1=kw model,3=zero grad,4=rad(ORL),5=env,6=wall
1       ! lbctime =0 big time step, =1 small time step (for lbc=2,4)
1       ! ortype  =0 arps orlanski form., =1 durran form (for lbc=2,4)
1       ! tbc    Top bc, = 1 rigid lid , 4=linear w-p rad. condition
3       ! fftopt  = 1 periodic fft, 3= even fft for upper rad condition
1       ! bbc    Bottom bc: 1=rigid, 3=zero grad (not functional)
0       ! hgradopt =0 lat. bc not comp. in adv. part =1 assume zero grd
0       ! rlxopt Relaxation option =0 none, =1 nsew, =2 ws, =3 en.
0.0     ! rlxlbc Relaxation coefficient  0.0<rlxlbc<1.0
4.0    ! cphase  phase speed for constant option...used in lbc=1,2,4
1.5     ! efactin  factor mult. the inflow environmental phase speed
1.5     ! efactout  factor mult. the outflow environmental phase speed
1.0     ! cflimit  time step limitation =cdt/dx...
0       ! fricopt = 0 no friction, = 1 surface fric. from Miller/Durran
0       ! basebal  =0 base state wind not frictionally balanced, =1 yes
2       ! turbopt  =0 no turb. applied, =2 base Sullivan et.al. (1994)
2       ! ktop    initial guess to the top of the boundary layer.
650.0   ! lmax  mixing length max for sun and chang/fiedler method...
50.0    ! lmin  CAREFULL!! mix len min for sun/chang/fiedler method.
2000.0  ! estpblmax  = estimated max level for top of pbl...(meters) 
1.0e-05 ! tkeprimer  minumum value of tke to start the turbulence mixing.
0.0010  ! tkemin minumum value used to set the minimum tke.
0.0     ! stablim minumum value of N used to find the pbl top. USED!!!
0       ! heatopt surf heat option. 0 = no, 1 = flux, 2= heating type
0       ! lesopt =0 no LES s/c ,=1 deard+ptprt rand once,=2 heat once
1       ! lesperi  index used to set periodic cond. in the PBL (not used)
1       ! heatvar  =0 constant with time, =1 sine wave in time.
000.0   ! heatstart  starting time in seconds the heating is applied.
86400.0 ! heatper   FULL period of the time dep. heating function.
216.0   ! heaval    value for option heatopt = 1 flux, =2 (watts/m**2).
-50.0   ! heatmin   minimum value of the heating function.
2       ! heatlev   k-index at which the heat(heatopt=2) is applied.
3       ! cmixopt   =0 no mixing, =1 horiz., =2 vert. =3 horiz and vert.
1       ! cmono  =0 original mixing, =1 monotonic limited (Xue, 2000)
5.2e-4  ! cfcmh4  ND Coeff. of 4th-order HORIZONTAL comp mixing (m**4/s)
5.2e-4  ! cfcmv4  ND Coeff. of 4th-order VERTICAL comp mixing (m**4/s)
0       ! divdmp  Acous.-wave div. damping (0=off, 1=on) (removed)
0.050   ! divdmpnd Non-dimensional divergence damping coefficient
1       ! raydmp  0=off,=1 vertical arps,=2 horz,3 h+v? arps,=4 vert durran
0.0033  ! cfrdmpv Vertical Rayleigh sponge coefficient (inverse seconds)
12000.0 ! zbrdmp  Altitude (meters) of the base of the Rayleigh sponge
0.0010  ! cfrdmph  Horz. rayleigh sponge coefficient (inverse seconds)
1000.   ! raydstrt Starting point in meters for horizontal damper
5000.   ! raydend  Ending point in meters for horizontal damper
0.20    ! flteps   Asselin time filter coef. Non-dimensional.
1200.0   ! trstdmp  Time interval (seconds) between restart data dumps
60.0   ! thisdmp  Time interval (seconds) between history data dumps
100.0   ! tmaxmin  Time interval (seconds) between max/min statistics
100.    ! tdiagnum Time interval (seconds) between diagnostics comps.
1       ! fopen    Number of reading/writing files at once(for mpi)
1       ! sdopt    Option for writing the surf wave drag,=0 no, =1 yes
1       ! sdindex  j index for computing the surface drag.
1       ! massopt  Option for computing/printing the model mass.
1       ! umaxopt  Option for computing/writing the max wind.
2       ! umaxk    Level at which the max x-y slice wind is computed.
1       ! dwater   Option for computing/writing the total water.