debug/goprof.m

   1 % debugging canopy interpolation
   2  step=2
   3  file='wrfout_d05_2012-11-11_12:00:00'
   4  long=-86.73
   5  lat=30.53
   6  p=wrfatm2struct(file,step);
   7  pp=nc2struct(file,{'CAN_TOP','CUF','CVF','FWH','FZ0','UF','VF','UAH','VAH'},{},step);
   8  p.can_top=pp.can_top;
   9  p.cuf=pp.cuf;
  10  p.cvf=pp.cvf;
  11  p.fwh=pp.fwh;
  12  p.fz0=pp.fz0;
  13  p.uf=pp.uf;
  14  p.vf=pp.vf;
  15  p.uah=pp.uah;
  16  p.vah=pp.vah;
  17 prof(p,long,lat)