no error info., but there is an error in the units of the output variable of the function wrf_user_getvar

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