thanks
I defined CURVGRID and tried rotating the winds in ana_winds:
I took a look at set_data.f and found how the clockwise rotation matrix is applied. Looks something like this:
Code: Select all
cff1=Uwind*CosAngler + Vwind*SinAngler
cff2=Vwind*CosAngler - Uwind*SinAngler
Uwind(i,j)=cff1
Vwind(i,j)=cff2
----
So I modified ana_winds.h to call for 'angler' and wrote:
Code: Select all
u_wind=0.0_r8
v_wind=10.0_r8
DO j=JstrT,JendT
DO i=IstrP,IendT
Uwind(i,j)=u_wind*COS(angler(i,j))+v_wind*SIN(angler(i,j))
Vwind(i,j)=v_wind*COS(angler(i,j))-v_wind*SIN(angler(i,j))
END DO
END DO
Figure from history file drawn as roms_quivergrd((Uwind(:,:)),(Vwind(:,:)),grd,20,0.01);
Magnitude increasing landward?
----
For reference, this is a figure of the grid angle from the history output (grid 'grd' obtained with roms_get_grid)
Colorbar shows angles in degrees.
[han,data]=roms_quivergrd(cos(grd.angle),sin(grd.angle),grd,20,0.1)
![Image](https://i.imgur.com/J8daC5n.jpg?1)