Opened 9 years ago

Closed 7 years ago

Last modified 3 years ago

#856 closed Enhancement (fixed)

Enhanced initialization for LIM3

Reported by: vancop Owned by: vancop
Priority: normal Milestone:
Component: LIM3 Version: release-3.3
Severity: Keywords: LIM*
Cc: Review:
MP ready?:
Progress:

Description

Current initialization for LIM3 generates negative concentrations, which are filtered out after one time step; but generate alerts in the code that can be misinterpreted.

A cleaner routine has now been implemented for v3_1.

Based on namelist values of initial ice concentration and thickness, the new initialization computes values for the thickness and area for each category. A gaussian distribution that respects the following principles is assumed:
1) volume and area are conserved
2) thickness of categories has to be within the prescribed bounds (hi_max's)
3) concentration of ice in each category has to be positive

The documentation is attached.
Any question to be addressed to Martin!

namelist_ice has to be changed because some of the initialization-related variables in the former code are useless (deformed ice thickness and area).

Best wishes to everyone,
Martin

Commit History (0)

(No commits)

Attachments (3)

limistate.F90 (23.3 KB) - added by vancop 9 years ago.
routine
namelist_ice (13.6 KB) - added by vancop 9 years ago.
namelist_ice
doc_newinit_lim3.pdf (114.5 KB) - added by vancop 9 years ago.
doc

Download all attachments as: .zip

Change History (7)

Changed 9 years ago by vancop

routine

Changed 9 years ago by vancop

namelist_ice

Changed 9 years ago by vancop

doc

comment:1 Changed 8 years ago by clevy

  • Owner changed from NEMO team to vancop

comment:2 Changed 7 years ago by barthele

The new version of limistate.F90 proposed here by Martin is not perfect yet, because :

1) the initial ice thickness provided seems to be (initial thickness in namelist) * (initial concentration in namelist) ;
2) it should be cleaned to avoid disturbing warnings that appear even if the routine does its job as expected.

But, it is already way better than the current version, in which the actual initial thickness is totally different from what it should be.

Hence, in my opinion, it's a step in the right direction towards a reliable initialization of ice in LIM3.

comment:3 Changed 7 years ago by clem

  • Resolution set to fixed
  • Status changed from new to closed

Fixed in branch: dev_r4028_CNRS_LIM3 to be included in the next merge party (Nov. 2013)

comment:4 Changed 3 years ago by nemo

  • Keywords LIM* added
Note: See TracTickets for help on using tickets.