[[ParentLinkage(=)]] [[BR]] [[PageOutline]] = '''ORCA2_LIM reference runs''' = [[BR]] [[BR]] [[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) # 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 - Ongoing simulations: == [[BR]] Abstract for the simulations performed over 10 years (latest ones are the first presented):[[BR]] * for [wiki:private/ORCA2Runs/PHC16 PHC_16] based on PHC_15 but :[[BR]] {{{ - apply -9% on precipitation & runoff instead of -8% in PHC_15 experiment }}} * for [wiki:private/ORCA2Runs/PHC15 PHC_15]: based on PHC_11 but :[[BR]] {{{ just uses "good" input fields for K1 and M2 waves }}} * for PHC_11: based on PHC_10 but includes both:[[BR]] {{{ it seems that input fields were not the good ones - global Kz tides - increased Kz tides in the ITF area }}} * for [wiki:private/ORCA2Runs/PHC10 PHC_10]: based on PHC_9 but includes both:[[BR]] {{{ - 3 waves lenghts penetration (RGB) through the traqsr.F90 module - ocean color: climatology of a 2D chlorophyll field }}} * for [wiki:private/ORCA2Runs/PHC9 PHC_9]: based on the PHC_8 experiment[[BR]] {{{ - change the way to compute wind stress & wind module: use the climatology of the field to be considered instead of the climatology of fields based on climatology components. says climatology of ||u||.u & ||u||.v fields instead using the wind module of climatological wind speed and climatological wind speed components u/v - in doing this we don't take at all the relative velocity of the ocean/sea-ice in the computation of the wind stress }}} * from PHC_1 to PHC_8 : {{{ - freshwater budget adjustement in playing with new namelist parameters alpha_precip, alpha_runoff - Gibraltar, Bab El Mandeb & Sound sraits transports adjustement ( based on scale factors or/and bondary condition fmask ) - computation of the wind module is modified: we add the difference between the variance of the climatological module (based on 6hourly components u/v) and wind module (based on climatological components u&v) }}} [[BR]] Results from the 10th year:[[BR]] ||'''Runs''' || '''PHC_1''' || '''PHC_2''' || '''PHC_3''' || '''PHC_5''' || '''PHC_6''' || '''PHC_7''' || '''PHC_8''' || '''PHC_9''' || '''PHC_10''' || '''PHC_11''' || '''PHC_15''' || '''PHC_16''' || || '''Precip''' || / || -10% || -12% || / || -8% || -10% || -8% || -8% || -8% || -8% || -8% || -9% || || '''Evap''' || / || / || / || +10% || +5% || +10% || +10% || +10% || +10% || +10% || +10% || +10% || || '''Empold''' (Sv) ||+1.88 || +0.26 || -0.06 || +1.37 || +0.31 || -0.25 || +0.07 || +0.13 || +0.1 || +0.12 || +0.13 || -0.02 || || '''SST/rt PHC3''' || / || / || / || ++ || + || ++ || ++ || ?? || ?? || ?? || ?? || ?? || || '''SSH drift''' (m) || / || / || -0.006 || +0.12 || +0.028 || -0.022 || +0.006 || +0.01 || +0.008 || +0.01 || +0.01 || -0.002 || [[BR]] == 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]] * 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 - Last results: == Plot of the differences between: wind module climatology and wind module computed with climatology wind speed: [[Image(Diff_wndm_strong_smooth_ymean.png, 60%)]] Plot: SQUAREROOT[ Wc^2 - (Uc^2 + Vc^2) ] where Wc is the wind module climatology, Uc & Vc are the wind speed climatology [[Image(Diff_wndm_variance_ymean.png, 60%)]] == F - 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 }}}