Changeset 6396


Ignore:
Timestamp:
2016-03-18T14:56:37+01:00 (6 years ago)
Author:
mathiot
Message:

Modification of the ISF documentation based on the comments from N. Jourdain

Location:
branches/2015/nemo_v3_6_STABLE/DOC/TexFiles
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/DOC/TexFiles/Chapters/Chap_DOM.tex

    r6317 r6396  
    887887By default, misfdep(:,:)=1 and no cells are masked. 
    888888 
    889 In case of ice shelf cavities, modifications of the model bathymetry and ice shelf draft into  
    890 the cavities are performed in the \textit{zgr\_isf} routine. The compatibility between ice shelf draft and bathymetry is checked.  
    891 If only one cell on the water column is opened at $t$-, $u$- or $v$-points, the bathymetry or the ice shelf draft is dug to fit this constrain. 
    892 If the incompatibility is too strong (need to dig more than 1 cell), the cell is masked.\\  
     889In case of ice shelf cavities (\np{ln\_isfcav}~=~true), modifications of the model bathymetry and ice shelf draft in  
     890the cavities are performed through the \textit{zgr\_isf} routine. The compatibility between ice shelf draft and bathymetry is checked:  
     891if only one cell on the water column is opened at $t$-, $u$- or $v$-points, the bathymetry or the ice shelf draft is dug to have a 2-level water column  
     892(i.e. two unmasked levels). If the incompatibility is too strong (i.e. need to dig more than one cell), the entire water column is masked.\\  
    893893 
    894894From the \textit{mbathy} array, the mask fields are defined as follows: 
  • branches/2015/nemo_v3_6_STABLE/DOC/TexFiles/Chapters/Chap_DYN.tex

    r6317 r6396  
    655655\label{DYN_hpg_isf} 
    656656Beneath an ice shelf, the total pressure gradient is the sum of the pressure gradient due to the ice shelf load and 
    657  the pressure gradient due to the ocean load. If cavity opened (\np{ln\_isfcav}~=~true) these 2 terms can be 
    658  calculated by setting \np{ln\_dynhpg\_isf}~=~true. No other scheme are working with the ice shelf.\\ 
    659  
    660 $\bullet$ The main hypothesis to compute the ice shelf load is that the ice shelf is in an isostatic equilibrium. 
    661  The top pressure is computed integrating from surface to the base of the ice shelf a reference density profile  
    662 (prescribed as density of a water at 34.4 PSU and -1.9$\degres C$) and corresponds to the water replaced by the ice shelf.  
    663 This top pressure is constant over time. A detailed description of this method is described in \citet{Losch2008}.\\ 
     657 the pressure gradient due to the ocean load. If cavities are present (\np{ln\_isfcav}~=~true) these two terms can be 
     658 calculated by setting \np{ln\_dynhpg\_isf}~=~true. No other scheme is working with ice shelves.\\ 
     659 
     660$\bullet$ The main hypothesis to compute the ice shelf load is that the ice shelf is in isostatic equilibrium. 
     661 The top pressure is computed integrating a reference density profile (prescribed as density of a water at 34.4  
     662PSU and -1.9$\degres C$) from the sea surface to the ice shelf base, which corresponds to the load of the water 
     663column in which the ice shelf is floatting. This top pressure is constant over time. A detailed description of  
     664this method is described in \citet{Losch2008}.\\ 
    664665 
    665666$\bullet$ The ocean load is computed using the expression \eqref{Eq_dynhpg_sco} described in \ref{DYN_hpg_sco}.  
    666 A treatment of the partial cell for top and bottom similar to the one described in \ref{DYN_hpg_zps} is done  
     667A treatment of the top and bottom partial cells similar to the one described in \ref{DYN_hpg_zps} is done  
    667668to reduce the residual circulation generated by the top partial cell.  
    668669 
  • branches/2015/nemo_v3_6_STABLE/DOC/TexFiles/Chapters/Chap_SBC.tex

    r6350 r6396  
    948948\namdisplay{namsbc_isf} 
    949949%-------------------------------------------------------------------------------------------------------- 
    950 Namelist variable in \ngn{namsbc}, \np{nn\_isf}, controls the ice shelf representation used.  
     950Namelist variable in \ngn{namsbc}, \np{nn\_isf}, controls the ice shelf representation used (Fig. \ref{Fig_SBC_isf}):  
     951 
     952%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     953\begin{figure}[!h]    \begin{center} 
     954\includegraphics[width=0.8\textwidth]{./TexFiles/Figures/Fig_SBC_isf.pdf} 
     955\caption{ \label{Fig_SBC_isf} 
     956Schematic for all the options available trough \np{nn\_isf}.} 
     957\end{center}   \end{figure} 
     958%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     959 
    951960\begin{description} 
     961\item[\np{nn\_isf}~=~0] 
     962The ice shelf routines are not used. The ice shelf melting is not computed or prescribed, the cavity have to be closed.  
     963If needed, the ice shelf melting should be added to the runoff or the precipitation file. 
     964 
    952965\item[\np{nn\_isf}~=~1] 
    953966The ice shelf cavity is represented. The fwf and heat flux are computed. Two different bulk formula are available: 
     
    10021015$\bullet$ \np{nn\_isf}~=~1 and \np{nn\_isf}~=~2 compute a melt rate based on the water mass properties, ocean velocities and depth. 
    10031016 This flux is thus highly dependent of the model resolution (horizontal and vertical), realism of the water masses onto the shelf ...\\ 
    1004  
    10051017 
    10061018$\bullet$ \np{nn\_isf}~=~3 and \np{nn\_isf}~=~4 read the melt rate from a file. You have total control of the fwf forcing. 
Note: See TracChangeset for help on using the changeset viewer.