1 function w=read_wrfout_tign(f)
2 % w=read_wrfout_tign(f)
3 % read variables from wrfout f for detect_fit_level2
4 t=nc2struct(f,{'Times'},{}); nframes=size(t.times,2)
5 w=nc2struct(f,{'Times','TIGN_G','FXLONG','FXLAT','UNIT_FXLAT','UNIT_FXLONG',...
6 'XLONG','XLAT','NFUEL_CAT','ITIMESTEP'},{'DX','DY','DT'},nframes);