New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3049 for branches/2011/dev_MERCATOR_2011_MERGE/DOC/TexFiles – NEMO

Ignore:
Timestamp:
2011-11-07T14:22:59+01:00 (13 years ago)
Author:
cbricaud
Message:

add changes from dev_r2802_MERCATOR10_diadct

Location:
branches/2011/dev_MERCATOR_2011_MERGE/DOC/TexFiles
Files:
2 edited
1 copied

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_MERCATOR_2011_MERGE/DOC/TexFiles/Chapters/Chap_DIA.tex

    r3046 r3049  
    684684the off-line use of this marvellous diagnostic tool. 
    685685 
     686 
    686687% ------------------------------------------------------------------------------------------------------------- 
    687688%       Harmonic analysis of tidal constituents 
     
    726727 
    727728We obtain in output $C_{j}$ and $S_{j}$ for each tidal wave. 
     729 
     730% ------------------------------------------------------------------------------------------------------------- 
     731%       Sections transports 
     732% ------------------------------------------------------------------------------------------------------------- 
     733\section{Transports across sections (\key{diadct}) } 
     734\label{DIA_diag_dct} 
     735 
     736A module is available to compute the transport of volume, heat and salt through sections. This diagnostic 
     737is actived with \key{diadct}. 
     738 
     739Each section is defined by the coordinates of its 2 extremities. The pathways between them are contructed 
     740using tools which can be found in  \texttt{NEMOGCM/TOOLS/SECTIONS\_DIADCT} and are written in a binary file 
     741 \texttt{section\_ijglobal.diadct\_ORCA2\_LIM} which is later read in by NEMO to compute on-line transports. 
     742 
     743The on-line transports module creates three output ascii files:  
     744 
     745- \texttt{volume\_transport} for volume transports (  unit: $10^{6} m^{3} s^{-1}$ ) 
     746 
     747- \texttt{heat\_transport}   for heat transports   (  unit: $10^{15} W $ ) 
     748 
     749- \texttt{salt\_transport}   for salt transports   (  unit: $10^{9}g s^{-1}$ )\\ 
     750 
     751 
     752Namelist parameters control how frequently the flows are summed and the time scales over which 
     753 they are averaged, as well as the level of output for debugging: 
     754 
     755%------------------------------------------namdct---------------------------------------------------- 
     756\namdisplay{namdct} 
     757%------------------------------------------------------------------------------------------------------------- 
     758 
     759\texttt{nn\_dct}: frequency of instantaneous transports computing 
     760 
     761\texttt{nn\_dctwri}: frequency of writing ( mean of instantaneous transports ) 
     762 
     763\texttt{nn\_debug}: debugging of the section 
     764 
     765\subsubsection{ To create a binary file containing the pathway of each section } 
     766 
     767In \texttt{NEMOGCM/TOOLS/SECTIONS\_DIADCT/run}, the file \texttt{ {list\_sections.ascii\_global}} 
     768contains a list of all the sections that are to be computed (this list of sections is based on MERSEA project metrics). 
     769 
     770Another file is available for the GYRE configuration (\texttt{ {list\_sections.ascii\_GYRE}}).  
     771 
     772Each section is defined by: 
     773 
     774\noindent \texttt{ long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name }\\ 
     775with: 
     776 
     777- \texttt{long1 lat1} , coordinates of the first extremity of the section; 
     778 
     779- \texttt{long2 lat2} , coordinates of the second extremity of the section; 
     780 
     781- \texttt{nclass} the number of bounds of your classes (e.g. 3 bounds for 2 classes); 
     782 
     783- \texttt{okstrpond} to compute heat and salt transport, \texttt{nostrpond} if no; 
     784 
     785- \texttt{ice}  to compute surface and volume ice transports, \texttt{noice} if no. \\ 
     786 
     787 
     788\noindent The results of the computing of transports, and the directions of positive 
     789 and negative flow do not depend on the order of the 2 extremities in this file.\\  
     790 
     791 
     792\noindent If nclass =/ 0,the next lines contain the class type and the nclass bounds: 
     793 
     794\texttt{long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name} 
     795 
     796\texttt{classtype} 
     797 
     798\texttt{zbound1} 
     799 
     800\texttt{zbound2} 
     801 
     802\texttt{.} 
     803 
     804\texttt{.} 
     805 
     806\texttt{nclass-1} 
     807 
     808\texttt{nclass} 
     809 
     810\noindent where \texttt{classtype} can be: 
     811 
     812- \texttt{zsal} for salinity classes 
     813 
     814- \texttt{ztem} for temperature classes 
     815 
     816- \texttt{zlay} for depth classes 
     817 
     818- \texttt{zsigi} for insitu density classes 
     819 
     820- \texttt{zsigp} for potential density classes \\ 
     821 
     822   
     823The script \texttt{job.ksh} computes the pathway for each section and creates a binary file  
     824\texttt{section\_ijglobal.diadct\_ORCA2\_LIM} which is read by NEMO. \\ 
     825 
     826It is possible to use this tools for new configuations: \texttt{job.ksh} has to be updated  
     827with the coordinates file name and path. \\ 
     828 
     829 
     830Examples of two sections, the ACC\_Drake\_Passage with no classes, and the 
     831 ATL\_Cuba\_Florida with 4 temperature clases (5 class bounds), are shown: 
     832 
     833\noindent \texttt{ -68.    -54.5   -60.    -64.7  00 okstrpond noice ACC\_Drake\_Passage} 
     834 
     835\noindent \texttt{ -80.5    22.5   -80.5    25.5  05 nostrpond noice ATL\_Cuba\_Florida} 
     836 
     837\noindent \texttt{ztem} 
     838 
     839\noindent \texttt{-2.0} 
     840 
     841\noindent \texttt{ 4.5} 
     842 
     843\noindent \texttt{ 7.0} 
     844 
     845\noindent \texttt{12.0} 
     846 
     847\noindent \texttt{40.0} 
     848 
     849 
     850\subsubsection{ To read the output files } 
     851 
     852The output format is : 
     853  
     854\small{\texttt{date, time-step number, section number, section name, section slope coefficient, class number,  
     855class name, class bound 1 , classe bound2, transport\_direction1 ,  transport\_direction2, transport\_total}}\\ 
     856 
     857 
     858For sections with classes, the first \texttt{nclass-1} lines correspond to the transport for each class  
     859and the last line corresponds to the total transport summed over all classes. For sections with no classes, class number 
     860\texttt{ 1 } corresponds to \texttt{ total class } and this class is called  \texttt{N}, meaning \texttt{none}.\\ 
     861 
     862 
     863\noindent \texttt{ transport\_direction1 } is the positive part of the transport ( \texttt{ > = 0 } ). 
     864 
     865\noindent \texttt{ transport\_direction2 } is the negative part of the transport ( \texttt{ < = 0 } ).\\ 
     866 
     867 
     868\noindent The \texttt{section slope coefficient} gives information about the significance of transports signs and direction:\\ 
     869 
     870 
     871 
     872\begin{tabular}{|c|c|c|c|p{4cm}|} 
     873\hline 
     874section slope coefficient & section type & direction 1 & direction 2 & total transport \\ \hline 
     8750.    &  horizontal & northward & southward & postive: northward  \\ \hline 
     8761000. &  vertical   & eastward  & westward  & postive: eastward  \\ \hline                 
     877\texttt{=/0, =/ 1000.}   &  diagonal   & eastward  & westward  & postive: eastward  \\ \hline                 
     878\end{tabular} 
     879 
     880 
    728881 
    729882% ------------------------------------------------------------------------------------------------------------- 
     
    769922are removed from the sub-basins. Note also that the Arctic Ocean has been split  
    770923into Atlantic and Pacific basins along the North fold line.  } 
    771 \end{center}   \end{figure} 
     924\end{center}   \end{figure}   
    772925%>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    773926 
     
    776929(see Section \ref{MISC_steric} below for one of them).  
    777930Activating those outputs requires to define the \key{diaar5} CPP key. 
     931\\ 
     932\\ 
     933 
    778934 
    779935 
Note: See TracChangeset for help on using the changeset viewer.