vis3d/hillframe.m

   1 disp('1. start matlab in wrf/WRFV3/test/em_fire to set up the paths properly')
   2 disp('   or cd there and run startup while in the directory')
   3 disp('2. read a wrfrst file into the workspace as in')
   4 disp('   >> ncload wrfrst_d01_0001-01-01_00:01:00')
   5
   6 % note the variables are read here without the WRF permutation of dimensions
   7 ideal=1;
   8 rst=1;
   9
  10 if rst,
  11     u=u_2;
  12     v=v_2;
  13     w=w_2;
  14     ph = ph_2;
  15 end
  16
  17 clf
  18 amin=[15,15,1];  % the atm grid part to show
  19 amax=[30,30,2];
  20 amin=[1,1,1];  % the atm grid part to show
  21 amax=[41,41,2];
  22
  23 qstep=[5,5];        % quiver step for wind on the surface fire grid
  24 astep=[1,1,1];      % quiver step for wind on the atmosphere grid
  25 qs=1.5;             % scaling for quiver arrows
  26 r=[10,10];          % refinement factor
  27 swind=1;            % plot surface wind
  28
  29 %-------------------------------------------------
  30
  31 frame3d(swind,amin,amax,astep,qstep,qs,...
  32     fxlong,fxlat,xlong,xlat,zsf,fgrnhfx,uf,vf,u,v,w,ph,phb,hgt)