January 2009, mirin Created using an updated version of CAM's interpic code (interpic_new): On bangkok: cd models/atm/cam/tools/interpic_new gmake set dir=/fs/cgd/csm/inputdata/atm/cam interpic -v -t $dir/coords/fv_0.23x0.31.nc $dir/dst/dst_1x1.25_c060901.nc dst_0.23x0.31_c090117.nc define dim from template file: lat len= 768 define dim from template file: lon len= 1152 define dim from template file: slat len= 767 define dim from template file: slon len= 1152 copy CVS_Id = $Id$ copy creation_date = 2006/08/22 11:35:36.159 copy prg_ID = Source file unknown Version unknown Date unknown copy cmd_ln = bds --bsn_fct_hrz=F --flg_cst_tpr=T --flg_rgn_xcl=T --bkf_itr_nbr=1 -x 288 -y 181 -m 33 -t 19801980 -o /data/zender/data/dst_1x1.25.nc copy history = Fri Sep 1 12:52:09 2006: ncrename -v rdb_fct_gmr,mbl_bsn_fct_geo dst_1x1.25.nc Fri Sep 1 12:52:09 2006: ncks -v rdb_fct_gmr from_charlie/dst_1x1.25.nc dst_1x1.25.nc Tue Aug 22 11:39:12 2006: ncap -O -s mbl_bsn_fct=rdb_fct_gmr /data/zender/data/dst_1x1.25.nc /data/zender/data/dst_1x1.25.nc Tue Aug 22 11:39:12 2006: ncks -O -v area,bsn_enm,bsn_sz,lat_grd,lon_grd,lat,lat_fv,lon,lon_sz,lat_sz,oro,lnd_frc,lnd_msk,rdb_fct_gmr,rdb_fct_unity,rdb_fct_tpg,rdb_fct_hyd,rdb_fct_rfl_mds_lnr,rdb_fct_rfl_mds_sqr,slat,slon /data/zender/data/dst_1x1.25.nc /data/zender/data/dst_1x1.25.nc Tue Aug 22 11:39:12 2006: ncap -O -s rdb_fct_unity=0.0*mbl_bsn_fct+0.1614 -s rdb_fct_tpg=mbl_bsn_fct -s rdb_fct_gmr=5.707*sfc_acm_fct -s rdb_fct_hyd=22.429*flw_acm_fct -s rdb_fct_rfl_mds_lnr=rfl_sfc_lnd_nrm_mds_lnr -s rdb_fct_rfl_mds_sqr=rfl_sfc_lnd_nrm_mds_sqr /data/zender/data/dst_1x1.25.nc /data/zender/data/dst_1x1.25.nc handle_special_cases from template file: Copying special var gw Copying special var lat Copying special var lon Copying special var slat Copying special var slon Copying special var w_stag skipping special case var lat skipping special case var lon create var mbl_bsn_fct_geo ndims= 2 dids= 2 1 vid= 7 fillvar: mbl_bsn_fct_geo(lon= 288 lat= 181 ) fillvar: mbl_bsn_fct_geo(lon= 1152 lat= 768 ) Copy and interpolate from input to output. ntime= 1 Starting time sample 1 Interpolating mbl_bsn_fct_geo input: shape= xy nx= 288 ny= 181 nz= 1 output: shape= xy nx= 1152 ny= 768 nz= 1 input min, avg, max = 0.000000000000000E+00 1.865315806208183E-02 5.707000000000000 doing horizontal interpolation output min, avg, max= 0.000000000000000E+00 1.864974213400789E-02 5.677375118411894 Finished.