Version 19 (modified by ctlod, 16 years ago) (diff) |
---|
ORCA2_LIM reference runs
A - Purposes:
Document the changes from ORCA2-LIM based on OPA8.2 and ORCA2-LIM base on last NEMO version
Key changes to be documented:
# Existing one :
- partial step representation of bottom topography
- updated TKE physics (surface/internal wave breaking)
- new sea-ice component (LIM-3)
# Existing but not included in NEMO trunk :
- Kz tides including ITF specific formulation (Laurent & Ariane PhD)
- mesoscale eddy parametrisation : added in momentum equations (Aiki) (Claude/Gurvan? Japan)
- 3 band light penetration (RGB) (Sebastian Japan)
- Ocean color used to specified time-space varying light penetration (Sebastian Japan)
- Diurnal cycle on qsr (from daily mean value) (Sebastian Japan)
# Forth coming one :
- mesoscale eddy parametrisation : Griffies scheme + Visbeck coefficient (MetOffice?, NOCS)
- npc reactivation (Claude/Gurvan?)
- conservative leap-frog (Matthieu)
- vvl + conservative leap-frog (Matthieu)
B - To do list/ Remaining Questions:
- check climatological CORE bulk forcing fields computed during 1st and 2nd year
- do moficiation to perfom annual mean outputs except the last year: OK
- check with Martin that all ice parameters for LIM 3.0 are the best one for ORCA2_LIM.
- ensure the upload of the STRAIT.dat file (transport through straits) : OK
- rebuild the meshmask since the scale factors have been changed on some straits: OK
Questions to solve when setting the control simulation:
- What are we doing with runoffs ? A: We used the standard climatology
- freswater budget, annual correction ? A: yes, in using the previous year ssh drift.
- Do we take into account the ice melting of Greenland and Antarctic into runoffs directly ?
- advection scheme to use ? A: the TVD for tracers
- define the way to manage the wind speed .vs. the sea-ice/ocean velocity, i.e do the difference of fields
- define the relaxation time scale for the SSS restoring term and define the area i.e. (global even under sea-ice ?) Done everywhere even under sea-ice, same restorint time scale (12 days) as for SST which require a value of deds= - 864 mm/day
- use the corrected interpolated temp/ salinity fields in Arctic from PHC for initialization suggested by Sylvain Bouillon see ticket: #227; Conkright et al. (2002) everywhere except in Arctic and Steel et al. (2001) for Arctic region data for the initialization ? (See Core paper)
- bathymetry review ( special treatment of straits Gibraltar, Bab El ..ec ) A: yes around Gibraltar and Bab El Mandeb straits with new depth, and a reduing of scale factors.
C - Few characteristics of runs as a reminder:
- perform internal damping (key_tradmp) ndmp = -1 around Gibraltat & Bab El Mandeb straits
- SSS damping ln_sssr = true with a restoring time scale of 12 days (same as for SST) deds = -864 mm/day
D - Last results:
Plot of the differences between: wind module climatology and wind module computed with climatology wind speed:
Plot: SQUARE( Wc2 - Uc2 - Vc2) where Wc is the wind module climatology, Uc & Vc are the wind speed climatology
E - Description of variables:
- LIM 3.0 variables output:
! ! ! title of the field ! name ! units ! save ! multipl. ! additive ! ! ! ! ! ! or not ! factor ! factor ! field_1 = 'Ice concentration ', 'iiceconc', '% ', 1 , 1.0 , 0.0 field_2 = 'Ice thickness ', 'iicethic', 'm ', 1 , 1.0 , 0.0 field_3 = 'Snow thickness ', 'isnowthi', 'm ', 1 , 1.0 , 0.0 field_4 = 'Daily bottom thermo ice production ', 'iicebopr', 'cm/day ', 1 , 100. , 0.0 field_5 = 'Daily dynamic ice production ', 'iicedypr', 'cm/day ', 1 , 100. , 0.0 field_6 = 'Oceanic flux at the ice base ', 'ioceflxb', 'w/m2 ', 1 , 1.0 , 0.0 field_7 = 'Ice velocity u ', 'iicevelu', 'm/s ', 1 , 1.0 , 0.0 field_8 = 'Ice velocity v ', 'iicevelv', 'm/s ', 1 , 1.0 , 0.0 field_9 = 'Sea surface temperature ', 'isstempe', 'C ', 1 , 1.0 , -273.15 field_10 = 'Sea surface salinity ', 'isssalin', 'PSU ', 1 , 1.0 , 0.0 field_11 = 'Total flux at ocean surface ', 'iocetflx', 'w/m2 ', 1 , 1.0 , 0.0 field_12 = 'Solar flux at ocean surface ', 'iocesflx', 'w/m2 ', 1 , 1.0 , 0.0 field_13 = 'Non-solar flux at ocean surface ', 'iocwnsfl', 'w/m2 ', 1 , 1.0 , 0.0 field_14 = 'Heat flux due to brine release ', 'iocehebr', 'w/m2 ', 1 , 1.0 , 0.0 field_15 = 'Wind stress u ', 'iocestru', 'Pa ', 1 , 1.0 , 0.0 field_16 = 'Wind stress v ', 'iocestrv', 'Pa ', 1 , 1.0 , 0.0 field_17 = 'Solar flux at ice/ocean surface ', 'iicesflx', 'w/m2 ', 1 , 1.0 , 0.0 field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2 ', 1 , 1.0 , 0.0 field_19 = 'Snow precipitation ', 'isnowpre', 'kg/day ', 1 , 1.0 , 0.0 field_20 = 'Mean ice salinity ', 'iicesali', 'psu ', 1 , 1.0 , 0.0 field_21 = 'Mean ice age ', 'iiceages', 'years ', 1 , 0.002739, 0.0 field_22 = 'Daily lateral thermo ice prod. ', 'iicelapr', 'cm/day ', 1 ,100. , 0.0 field_23 = 'Daily snowice ice production ', 'iicesipr', 'cm/day ', 1 ,100. , 0.0 field_24 = 'Mean ice temperature ', 'iicetemp', 'C ', 1 , 1.0 , -273.15 field_25 = 'Ice total heat content ', 'iiceheco', '10^9 J ', 1 , 1.0 , 0.0 field_26 = 'Ice surface temperature ', 'iicesurt', 'C ', 1 , 1.0 , -273.15 field_27 = 'Snow temperature ', 'isnotem2', 'C ', 1 , 1.0 , -273.15 field_28 = 'Fsbri - brine salt flux ', 'iicfsbri', 'kg/m2/s ', 1 , 1.0 , 0.0 field_29 = 'Fseqv - equivalent FW salt flux ', 'iicfseqv', 'kg/m2/s ', 1 , 1.0 , 0.0 field_30 = 'Brine volume ', 'ibrinvol', '% ', 1 , 100.0 , 0.0 field_31 = 'Frazil ice collection thickness ', 'iicecolf', 'm ', 1 , 1.0 , 0.0 field_32 = 'Ice strength ', 'iicestre', 'N/m ', 1 , 0.001 , 0.0 field_33 = 'Ice velocity ', 'iicevelo', 'm/s ', 1 , 1.0 , 0.0 field_34 = 'Surface melt ', 'iicesume', 'cm/day ', 1 ,100. , 0.0 field_35 = 'Bottom melt ', 'iicebome', 'cm/day ', 1 ,100. , 0.0 field_36 = 'Divergence ', 'iicedive', '10-8s-1 ', 1 , 1.0e8 , 0.0 field_37 = 'Shear ', 'iiceshea', '10-8s-1 ', 1 , 1.0e8 , 0.0