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)