wind rotated ncl

From: Diego Saśl Carrió Carrió <carrio.carrio_at_nyahnyahspammersnyahnyah>
Date: Wed Jan 29 2014 - 08:27:49 MST

Dear All,

I have an important doubt. I would like to get wind's field referenced at
earth coordinates. I have found a way to achieve this:

Uearth = U*cosalpha + V*sinalpha
Vearth = V*cosalpha - U*sinalpha.

, where:
 sinalpha = wrf_user_getvar(a,"SINALPHA",0)
 cosalpha = wrf_user_getvar(a,"COSALPHA",0)
 u = wrf_user_getvar(a,"ua",-1) ; u at mass point
 v = wrf_user_getvar(a,"va",-1) ; v at mass point

But I see, there are another way to do this. Using this function:

  uvm = wrf_user_getvar(a,"uvmet",-1)

The problem is uvm(0, , , , ) {component u} isn't the same than Uearth.

I don't know what is the problem.

Could you help me, please?

Best Regards,

Diego

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jan 29 08:28:02 2014

This archive was generated by hypermail 2.1.8 : Fri Feb 07 2014 - 16:39:11 MST