- Timestamp:
- 2015-10-13T16:20:38+02:00 (9 years ago)
- Location:
- branches/2015/dev_r5151_UKMO_ISF/DOC/TexFiles
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5151_UKMO_ISF/DOC/TexFiles/Biblio/Biblio.bib
r5120 r5792 276 276 author = {A. Beckmann and H. Goosse}, 277 277 title = {A parameterization of ice shelf-ocean interaction for climate models}, 278 journal = OM 279 year = {2003} 280 volume = {5} 278 journal = OM, 279 year = {2003}, 280 volume = {5}, 281 281 pages = {157--170} 282 282 } … … 1347 1347 } 1348 1348 1349 @ARTICLE{Holland1999, 1350 author = {D. Holland and A. Jenkins}, 1351 title = {Modeling Thermodynamic Ice-Ocean Interactions at the Base of an Ice Shelf}, 1352 journal = JPO, 1353 year = {1999}, 1354 volume = {29}, 1355 pages = {1787--1800}, 1356 } 1357 1349 1358 @ARTICLE{HollowayOM86, 1350 1359 author = {Greg Holloway}, … … 1493 1502 } 1494 1503 1504 @ARTICLE{Jenkins1991, 1505 author = {A. Jenkins}, 1506 title = {A one-dimensional model of ice shelf-ocean interaction}, 1507 journal = JGR, 1508 year = {1991}, 1509 volume = {96}, number = {C11}, 1510 pages = {2298--2312} 1511 } 1512 1513 @ARTICLE{Jenkins2010, 1514 author = {A. Jenkins}, 1515 title = {observation and parameterization of ablation at the base of Ronne Ice Shelf, Antarctica}, 1516 journal = JPO, 1517 year = {2010}, 1518 volume = {40}, number = {10}, 1519 pages = {2298--2312} 1520 } 1521 1495 1522 @BOOK{Jerlov_Bk68, 1496 1523 title = {Optical Oceanography}, … … 1898 1925 volume = {51}, 1899 1926 pages = {737--769} 1927 } 1928 1929 @ARTICLE{Losch2008, 1930 author = {M. Losch}, 1931 title = {Modeling ice shelf cavities in a z coordinate ocean general circulation model}, 1932 journal = JGR, 1933 year = {2008}, 1934 volume = {113}, number = {C13}, 1900 1935 } 1901 1936 -
branches/2015/dev_r5151_UKMO_ISF/DOC/TexFiles/Chapters/Chap_SBC.tex
r5120 r5792 49 49 (\np{nn\_ice}~=~0,1, 2 or 3); the addition of river runoffs as surface freshwater 50 50 fluxes or lateral inflow (\np{ln\_rnf}~=~true); the addition of isf melting as lateral inflow (parameterisation) 51 (\np{nn\_isf}~=~2 or 3 and \np{ln\_isfcav}~=~false) or as surface flux at the land-ice ocean interface 52 (\np{nn\_isf}~=~1 or 4 and \np{ln\_isfcav}~=~true); 51 or as surface flux at the land-ice ocean interface (\np{ln\_isf}=~true); 53 52 the addition of a freshwater flux adjustment in order to avoid a mean sea-level drift (\np{nn\_fwb}~=~0,~1~or~2); the 54 53 transformation of the solar radiation (if provided as daily mean) into a diurnal … … 958 957 \namdisplay{namsbc_isf} 959 958 %-------------------------------------------------------------------------------------------------------- 960 Namelist variable in \ngn{namsbc}, \np{nn\_isf}, 959 Namelist variable in \ngn{namsbc}, \np{nn\_isf}, control the kind of ice shelf representation used. 961 960 \begin{description} 962 961 \item[\np{nn\_isf}~=~1] 963 The ice shelf cavity is represented. The fwf and heat flux are computed. 962 The ice shelf cavity is represented. The fwf and heat flux are computed. 2 bulk formulations are available: the ISOMIP one (\np{nn\_isfblk = 1}) described in (\np{nn\_isfblk = 2}), the 3 equation formulation described in \citet{Jenkins1991}. In addition to this, 963 3 different way to compute the exchange coefficient are available. $\gamma\_{T/S}$ is constant (\np{nn\_gammablk = 0}), $\gamma\_{T/S}$ is velocity dependant \citep{Jenkins2010} (\np{nn\_gammablk = 1}) and $\gamma\_{T/S}$ is velocity dependant and stratification dependent \citep{Holland1999} (\np{nn\_gammablk = 2}). For each of them, the thermal/salt exchange coefficient (\np{rn\_gammat0} and \np{rn\_gammas0}) have to be specified (the default values are for the ISOMIP case). 964 964 Full description, sensitivity and validation in preparation. 965 965 … … 969 969 (\np{sn\_depmax\_isf}) and the base of the ice shelf along the calving front (\np{sn\_depmin\_isf}) as in (\np{nn\_isf}~=~3). 970 970 Furthermore the fwf is computed using the \citet{Beckmann2003} parameterisation of isf melting. 971 The effective melting length (\np{sn\_Leff\_isf}) is read from a file .971 The effective melting length (\np{sn\_Leff\_isf}) is read from a file and the exchange coefficients are set as (\np{rn\_gammat0}) and (\np{rn\_gammas0}). 972 972 973 973 \item[\np{nn\_isf}~=~3] … … 987 987 \np{nn\_isf}~=~3 and \np{nn\_isf}~=~4 read the melt rate and heat flux from a file. You have total control of the fwf scenario. 988 988 989 989 This can be usefull if the water masses on the shelf are not realistic or the resolution (horizontal/vertical) are too 990 990 coarse to have realistic melting or for sensitivity studies where you want to control your input. 991 991 Full description, sensitivity and validation in preparation. 992 992 993 There is 2 ways to apply the fwf to NEMO. The first possibility (\np{ln\_divisf}~=~false) applied the fwf 994 and heat flux directly on the salinity and temperature tendancy. The second possibility (\np{ln\_divisf}~=~true) 995 apply the fwf as for the runoff fwf (see \S\ref{SBC_rnf}). The mass/volume addition due to the ice shelf melting is, 996 at each relevant depth level, added to the horizontal divergence (\textit{hdivn}) in the subroutine \rou{sbc\_isf\_div} 997 (called from \mdl{divcur}). 993 \np{rn\_hisf\_tbl} is the top boundary layer (tbl) thickness used by the Losch parametrisation \citep{Losch2008} to compute the melt. if 0, temperature/salt/velocity in the top cell is used to compute the melt. 994 Otherwise, NEMO used the mean value into the tbl. 998 995 % 999 996 % ================================================================ -
branches/2015/dev_r5151_UKMO_ISF/DOC/TexFiles/Namelist/namsbc
r5120 r5792 19 19 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave 20 20 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) 21 nn_isf = 0 ! ice shelf melting/freezing (/=0 => fill namsbc_isf) 22 ! 0 =no isf 1 = presence of ISF 23 ! 2 = bg03 parametrisation 3 = rnf file for isf 24 ! 4 = ISF fwf specified 25 ! option 1 and 4 need ln_isfcav = .true. (domzgr) 21 ln_isf = .false. ! ice shelf melting/freezing (T => fill namsbc_isf) 26 22 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) 27 23 nn_fwb = 3 ! FreshWater Budget: =0 unchecked -
branches/2015/dev_r5151_UKMO_ISF/DOC/TexFiles/Namelist/namsbc_isf
r5120 r5792 2 2 &namsbc_isf ! Top boundary layer (ISF) 3 3 !----------------------------------------------------------------------- 4 nn_isf = -99 ! ice shelf melting/freezing (/=0 => fill namsbc_isf) 5 ! 1 = ISF explicit (cavity open) 6 ! 2 = bg03 parametrisation (cavity closed) 7 ! 3 = ISF specified in depth along the calving front (cavity closed) 8 ! 4 = ISF melting specified (cavity open) 9 ! option 1 and 4 need ln_isfcav = .true. (domzgr) 4 10 ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! 5 11 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! 6 12 ! nn_isf == 4 7 sn_qisf = 'rnfisf' , -12 ,'sohflisf', .false. , .true. , 'yearly' , '' , ''8 13 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' 9 14 ! nn_isf == 3 10 15 sn_rnfisf = 'runoffs' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' 11 ! nn_isf == 2 and316 ! nn_isf == 2 or 3 12 17 sn_depmax_isf = 'runoffs' , -12 ,'sozisfmax' , .false. , .true. , 'yearly' , '' , '' 13 18 sn_depmin_isf = 'runoffs' , -12 ,'sozisfmin' , .false. , .true. , 'yearly' , '' , '' 14 19 ! nn_isf == 2 15 20 sn_Leff_isf = 'rnfisf' , 0 ,'Leff' , .false. , .true. , 'yearly' , '' , '' 16 ! for all case17 ln_divisf = .true. ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?)18 21 ! only for nn_isf = 1 or 2 19 22 rn_gammat0 = 1.0e-4 ! gammat coefficient used in blk formula 20 23 rn_gammas0 = 1.0e-4 ! gammas coefficient used in blk formula 24 ! only for nn_isf = 1 or 4 25 nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 21 26 ! only for nn_isf = 1 22 nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 23 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 27 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 24 28 ! 0 => thickness of the tbl = thickness of the first wet cell 25 ln_conserve = .true. ! conservative case (take into account meltwater advection)26 29 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 27 30 ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010)
Note: See TracChangeset
for help on using the changeset viewer.