2 function ros=read_data_from_wrfout(wrfout,time)
4 % for witch its 3100 2600
8 p=nc2struct(wrfout,{'F_ROS11','F_ROS12','F_ROS13','F_ROS11','F_ROS21','F_ROS23','F_ROS31','F_ROS32','F_ROS33'},{},time);
10 ros=zeros(size(p.f_ros11,1),size(p.f_ros11,2),3,3);
12 ros(:,:,1,1)=p.f_ros11;
14 ros(:,:,1,2)=p.f_ros12;
16 ros(:,:,1,3)=p.f_ros13;
18 ros(:,:,3,1)=p.f_ros31;
20 ros(:,:,3,2)=p.f_ros32;
22 ros(:,:,3,3)=p.f_ros33;
24 ros(:,:,2,1)=p.f_ros21;
26 ros(:,:,2,3)=p.f_ros23;