&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 =  'ROBLE.tgcm24r.peqmb02.nc'
;msspath_read  =  '/ROBLE/tgcm24r/peqmb02.nc'
;histfile_write = 'FOSTER.timegcm-dhvres.peqmb02.nc'
;msspath_write  = '/FOSTER/timegcm_dhvres/peqmb02.nc'
;
;histfile_read =  '/hao/tgcm/timegcm1.2_yr02/py2j/ROBLE.timegcm1.2.py2j008.nc'
;msspath_read  =  '/ROBLE/tgcm24r/peqmb02.nc'
;histfile_write = 'FOSTER.timegcmdr.py2j008.nc'
;msspath_write  = '/FOSTER/timegcm_dhvres/peqmb02.nc'
;
;histfile_read =  'ROBLE.timegcm1.2.peqmb015.nc'
;histfile_write = 'TGCM.timegcmdr.peqmb015.nc'
;
; 11/2/06: make double-res solstice file for timegcm1.3:
;histfile_read =  'ROBLE.timegcm1.2.pdsnc021.nc'
;histfile_write = 'TGCM.timegcmdr.pdsnc021.nc'
;
; 1/31/07: make another double-res solstice file for timegcm1.3:
;histfile_read =  'ROBLE.timegcm1.2.py2j007.nc'
;histfile_write = 'TGCM.timegcmdr.py2j007.nc'
;
;histfile_read =  'ROBLE.timegcm1.2.py5j034.nc'
;histfile_write = 'TGCM.timegcmdr.py5j034.nc'
;
; 9/5/07: make day 121 (May 1), 1997 history for Hanli:
;histfile_read =  'ROBLE.timegcm1.2.py97ur13.nc'
;histfile_write = 'TGCM.timegcmdr.py97ur13.nc'
;
;histfile_read =  'ROBLE.timegcm1.2dp7.peqnc005.nc'
;histfile_write = 'ROBLE.timegcm1.2dres.peqnc005.nc'
;
; 11/4/08: make 3d double-res files from timegcm1.2 files:
;
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:21 peqdc004.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:22 peqgc004.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:21 peqhc004.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:21 peqnc004.nc
;
;histfile_read =  'peqdc004.nc'
;histfile_write = 'peqdc004_dres.nc'
;histfile_read =  'peqgc004.nc'
;histfile_write = 'peqgc004_dres.nc'
;histfile_read =  'peqhc004.nc'
;histfile_write = 'peqhc004_dres.nc'
;histfile_read =  'peqnc004.nc'
;histfile_write = 'peqnc004_dres.nc'
;
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:42 pdsdc004.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:42 pdsgc004.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 15:07 pdshc005.nc
;-rw-r--r--    1 foster   ncar       52980440 Nov 04 14:41 pdsnc004.nc
;
;histfile_read =  'pdsdc004.nc'
;histfile_write = 'pdsdc004_dres.nc'
;histfile_read =  'pdsgc004.nc'
;histfile_write = 'pdsgc004_dres.nc'
;histfile_read =  'pdshc005.nc'
;histfile_write = 'pdshc005_dres.nc'
;histfile_read =  'pdsnc004.nc'
;histfile_write = 'pdsnc004_dres.nc'
;histfile_read =  '/ptmp/wbwang/Airglow/pyr2007031.nc'
;histfile_write = 'dpyr2007031.nc'
;histfile_read =  'ROBLE.timegcm1.41r.p96sept002.nc'
;histfile_write = 'ROBLE.timegcm1.41d.p96sept002.nc'

 histfile_read =  '/hao/tgcm/data/tiegcm1.93/TGCM.tiegcm1.93.pcntr_jsol_smin.nc'
 histfile_write = 'TGCM.tiegcm1.93_dres.pcntr_jsol_smin.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 = 172,0,0
/