Hi, I use ncl-6.0.0
Sorry to make a trouble
fl =new(7, file)
fl(0) = addfile ("$HOME/lu_modif_wrf/base-run/wrfout_d01_2012-07-01_00:00:00.nc","r")
......
fl(6) = addfile ("$HOME/lu_modif_wrf/allurban-run/wrfout_d01_2012-07-01_00:00:00.nc","r")
do i=0, 6
U10(i,:,:,:)=wrf_user_getvar(fl(i),"U10",-1)
.....
end do
But when I plot the time series of U10, the outcome is about ~100m/s, the units seemed to be exaggerated.And there is no error info.
I then used fl0, fl1, .... (it is inconvenient since cannot use do loop )to replace fl(i) in wrf_user_getvar, the outcome is right(compared with ncview).
I doubt whether fl(i) is the source of problem. Is there any good way to be more efficient in my case?
Regards.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Apr 27 22:37:55 2013
This archive was generated by hypermail 2.1.8 : Thu May 02 2013 - 11:16:45 MDT