
import load_dump

def test_wind_rotation():
    """
    Test that winds are being rotated correctly.

    Look particular geographic locations near the North pole in both the
    atmosphere, ice and ocean. Get the wind directly at these points and
    normalise. Make sure they are the same.
    """

    load_dump.load_coupling_fields('om_360x300-unit_test', ['uwnd', 'vwnd'])
    routine = 'atmo_boundar_layer'
    inputs = ['uatm', 'vatm']
    inputs, _ = load_dump.load('om_360x300-unit_test', routine, inputs, [])
    
    # Choose our point, should be a place where grid cells are going to have 
    # different rotations. 

    lat = 1.22317  
    lon = -1.77761

    # Check that there is a non-zero angle at this point. 

    # Get the atm and ice grid boxes that corresponds to this location. 

    # U and V winds at this location in both atmos and ice. 

    # The should not be the same, the ice should have been rotated. 

@@ -299,16 +299,27 @@
 %22%22%0A%0A    
+pass%0A%0A    #
 load_dum
@@ -387,16 +387,17 @@
 '%5D)%0A    
+#
 routine 
@@ -419,24 +419,25 @@
 _layer'%0A    
+#
 inputs = %5B'u
@@ -454,16 +454,17 @@
 m'%5D%0A    
+#
 inputs, 
@@ -527,20 +527,16 @@
 ts, %5B%5D)%0A
-    
 %0A    # C
@@ -604,17 +604,16 @@
  to have
- 
 %0A    # d
@@ -631,23 +631,23 @@
 tations.
- 
 %0A%0A    
+#
 lat = 1.
@@ -655,15 +655,14 @@
 2317
-  
 %0A    
+#
 lon 
@@ -730,17 +730,16 @@
 s point.
- 
 %0A%0A    # 
@@ -803,17 +803,16 @@
 ocation.
- 
 %0A%0A    # 
@@ -864,17 +864,16 @@
 and ice.
- 
 %0A%0A    # 
@@ -933,10 +933,9 @@
 rotated.
- 
 %0A

Comment out incomplete test.