← Back to branch summary

~fluidity-core/fluidity/sediment_restructure

« back to all changes in this revision

Viewing changes to parameterisation/Gls_vertical_turbulence_model.F90

merge with branch including remap_surface_to_field but then removed this change and approached it a different way which seems to work

Show diffs side-by-side

added added

removed removed

Lines of Context:
parameterisation/Gls_vertical_turbulence_model.F90
1504
1504
    real, dimension(2)                   :: temp_vector_2D
1505
1505
    real, dimension(3)                   :: temp_vector_3D
1506
1506
    logical                              :: surface_allocated
 
1507
    integer                              :: stat
1507
1508
 
1508
1509
    MaxIter = 10
1509
1510
    z0s_min = 0.003
1524
1525
            call create_surface_mesh(ocean_mesh, top_surface_nodes, tke%mesh, &
1525
1526
                                     top_surface_element_list, 'OceanTop')
1526
1527
            call allocate(surface_forcing, wind_surface_field%dim, ocean_mesh, name="surface_velocity")
1527
 
            call allocate(surface_pos, positions%dim, ocean_mesh, name="surface_positions")
 
1528
            surface_pos = get_coordinates_remapped_to_surface(positions, ocean_mesh, top_surface_element_list) 
1528
1529
            call deallocate(ocean_mesh)
1529
 
            call remap_field_to_surface(positions, surface_pos, top_surface_element_list)
1530
1530
 
1531
1531
            if (tke%mesh%continuity == velocity%mesh%continuity) then
1532
1532
                call set(surface_forcing,wind_surface_field)