New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
2010WP/2010Stream3 (diff) – NEMO

Changes between Version 1 and Version 2 of 2010WP/2010Stream3


Ignore:
Timestamp:
2009-11-04T18:53:56+01:00 (14 years ago)
Author:
acc
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • 2010WP/2010Stream3

    v1 v2  
    2222 
    2323=== Description and participants === 
     24---- 
     25 
     26== '''S3.2 : from NOCS-NEMO team (short-term)''' == 
     27[[BR]] 
     28 
     29=== S3.2-a — Revisit of mass and salt fluxes between ice, ocean and air === 
     30         '''Motivation:''' As a consequence of the embedding of sea-ice into the ocean (see S3.2-b), the total mass of salt of the combined sea-ice/ocean system is  
     31         conserved by removing a volume of water from the ocean when ice forms or adding a volume of water when it melts.  When ice forms, salt contained 
     32         in the water is returned to the ocean, less the salt retained in the ice itself.  When ice melts, the salt contained is added back to the ocean.[[BR]] 
     33         '''Status :''' Carried over from 2009 but work is well advanced. Full testing and documentation will be complete for 3.3 release [[BR]] 
     34         '''Principal Investigator :''' Yevgeny Aksenov (yka@noc.soton.ac.uk) [[BR]] 
     35=== S3.2-b — Embedding sea-ice inside the ocean === 
     36         '''Motivation:''' To improve the coupling between the sea-ice and ocean models, we embed the sea-ice into the ocean by changing the volume of  
     37         the ocean and depressing the ocean interface by a depth equal to the mean ice draft.   
     38         The latter is calculated as the water equivalent of the total ice and snow mass in the model grid box using in-situ local densities of seawater.   
     39         The pressure gradient that takes into account the total weight of sea ice and snow in the top model grid box is added to the ocean momentum budget.  
     40         The vertical coordinate is redefined for the whole water column.[[BR]] 
     41         '''Status :''' During 2009 the above scheme was implemented and tested at NOCS in a 100-year long integration of an ORCA2-LIM2 configuration using both  
     42         linear and non-linear free surface options.   
     43         It was also tested for short runs in the ORCA2-LIM3 configurations.   
     44         A longer run in ORCA2-LIM3 will take place early in 2010 and the code committed for the NEMO 3.3 release. 
     45         '''Principal Investigator :''' Yevgeny Aksenov (yka@noc.soton.ac.uk) [[BR]] 
     46=== S3.2-c — Embedding sea-ice inside the ocean for CICE === 
     47        '''Motivation:'''  Similar changes will have to be made in order to fully embed CICE-modelled sea-ice. Subtle differences and problems may arise because of the B-grid basis of CICE. This work is listed separately to ensure that possible delays will not impact on the delivery of an embedded solution for LIM in 3.3 
     48         '''Principal Investigator :''' Yevgeny Aksenov (yka@noc.soton.ac.uk) [[BR]] 
     49=== S3.2-d — Lateral ocean physics === 
     50 
     51        '''Motivation:'''  Griffies's implementation of the Gent and  McWilliams eddy transport as a skew flux is being implemented. This has the advantage that  its numerical discretization can be written in terms of contributions  from quarter cells – ‘triads’. This gives a tighter stencil,  disallowing 2-gridpoint numerical noise that is permitted by the  advective discretization. Isopycnal diffusion is also being  discretized in terms of these triads, obviating the need to smooth  isopycnal slopes horizontally with a Shapiro filter (as currently  implemented in NEMO), or apply a background horizontal diffusivity  that would mix diapycnally. Visbeck et al.'s formulation of spatially varying  diffusivities is also being implemented, as an alternative to the current formulation based on  the Held and Larichev timescale and the Rossby radius as a  lengthscale. Mixed-layer depth is sensitive to the slope limiting that  is employed, and a slope limiting algorithm will be added that  behaves satisfactorily within and immediately below the mixed layer. [[BR]] 
     52 
     53        '''Status: ''' A working prototype code for the iso-neutral and skew-flux operator was carefully reviewed in 2009. An final version will be ready for the 3.3 release. Considerable care has been taken in the formulation to ensure the tensorial representation is consistent with the variable volume layers (S-coordinate/z* representation). This has been documented and will be provided with the release.   [[BR]] 
     54        '''Principal Investigator :''' George Nurser (agn@noc.soton.ac.uk) [[BR]] 
    2455 
    2556===  Advancement  ===