Last edited [[Timestamp]] [[BR]] [[PageOutline]] = '''ORCA2_LIM reference runs''' = [[BR]] == A - Purposes: == [[BR]] Document the changes from ORCA2-LIM based on OPA8.2 and ORCA2-LIM base on last NEMO version [[BR]] Key changes to be documented:[[BR]] # Existing one : * partial step representation of bottom topography * updated TKE physics (surface/internal wave breaking) * new sea-ice component (LIM-3) * test (again) free-slip option # 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) * New TKE physics #2 (MERCATOR) * New advection scheme on the vertical COMPACT (Laurent) * Time splitting over the vertical for the advection (Gurvan) * Freshwater input under Ice Shelves (JMM DRAKKAR) # Forth coming one : * Fox-Kemper vertical velocities around the MLD parameterization (Gurvan & Xavier) * Advective BBL * bound Aeiv as a function of mean velocity * add Aeiv = advection-diffusion equation Eden et Greatbatch 2008-2009 * add 3D source for ice-shelf melting, runoff and cross-land adv * Code and test Gent coef. vertical dependence to N2 * Background Kz as a function of latitutde * mesoscale eddy parametrisation : Griffies scheme + Visbeck coefficient (MetOffice, NOCS) * sensitivity of the ediss parameter (MERCATOR) * npc reactivation (Claude/Gurvan) * conservative leap-frog (Matthieu) * vvl + conservative leap-frog (Matthieu) == B - Ongoing simulations: == # [wiki:private/ORCA2Runs/LIM2 ORCA2 + LIM 2.0] # [wiki:private/ORCA2Runs/LIM3 ORCA2 + LIM 3.0] == C - 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:[[BR]] * Damping term in the Red sea is not clear, strong differences with Levitus in this area, since it is reduced in dtatem.F90 ''' 2 check''' * 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. '' == D - 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 * == 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 }}}