&mksrc_input
;
; This code interpolates in horizontal and/or vertical, but is not
; fully generic. It will work for either or both of the following resolutions:
;   horizontal (lat x lon): from 5.0 to 2.5 degrees
;   vertical (ln(p0/p)): from 0.5 to 0.25
;
; Namelist inputs:
;
; histfile_read: Local input disk file to read. If this file does not exist, 
;   msspath_read will be read from mss to this disk file). If histfile_read
;   is not provided, msspath_read must be provided (in this case, msspath_read 
;   will be read from mss to the cwd).
;
; msspath_read: Mss path of file to read. If histfile_read is not
;   found, msrcp will be called to read this file from mss to disk file
;   histfile_read. If msspath_read is not provided, file will not be
;   read from mss (in this case, histfile_read must exist).
;
; histfile_write: Local disk file to write output (if it pre-exists, it 
;   will be overwritten). If histfile_write is not provided, msspath_write
;   must be provided (histfile_write in cwd will be constructed from 
;   msspath_write). Be sure to provide a path with enough disk space (or 
;   run the program from a directory with plenty of disk space).
;
; msspath_write: Mss path to which output file is to be disposed.
;   (will dispose histfile_write to this path). If msspath_write is not
;   provided, the disk file will not be disposed to the mss.
;
; horizontal_in: Horizontal resolution (delta degrees lat and lon) of
;   input file.
; 
; horizontal_out: Horizontal resolution (delta degrees lat and lon) of
;   output file. If horizontal_in == horizontal_out, then no interpolation
;   will be performed in the horizontal dimension.
;
; vertical_in: Vertical resolution (ln(p/p0)) of input file, e.g., 0.5.
;
; vertical_out: Vertical resolution (ln(p/p0)) of output file, e.g., 0.25.
;   If vertical_in == vertical_out, then no inpterpolation will be performed
;   in the vertical dimension.
;
; mtimes: model time(s) (day,hour,min) of input and output histories.
;
;------------------------------------------------------------------------
;
 histfile_read  = '/hao/aim1/tgcm/data/ABRECHT.V10.pequmin035nr2.ubc18.nc'
 histfile_write = '/hao/aim1/tgcm/data/ABRECHT.V10.pequmin035nr2.ubc18.dres.nc'
;
 dlon_in        = 5.00 ; longitude resolution of histfile_read  (delta degrees)
 dlon_out       = 2.50 ; longitude resolution of histfile_write (delta degrees)
 dlat_in        = 5.00 ; latitude resolution of histfile_read  (delta degrees)
 dlat_out       = 2.50 ; latitude resolution of histfile_write (delta degrees)
 vertical_in    = 0.50
 vertical_out   = 0.25
 mtimes = 211,0,0
/