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 2902 – NEMO

Changeset 2902


Ignore:
Timestamp:
2011-10-11T15:13:47+02:00 (13 years ago)
Author:
cbricaud
Message:

improve documentation for diadct

File:
1 edited

Legend:

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

    r2871 r2902  
    684684the off-line use of this marvellous diagnostic tool. 
    685685 
     686 
     687% ------------------------------------------------------------------------------------------------------------- 
     688%       Sections transports 
     689% ------------------------------------------------------------------------------------------------------------- 
     690\section{Transports across sections (\key{diadct}) } 
     691\label{DIA_diag_dct} 
     692 
     693A module is available to compute volume, heat and salt tranports through sections. This diagnostic 
     694is actived with \key{diadct}. 
     695 
     696Each section is defined by the coordinates of its 2 extremities. The pathway between them is contructed 
     697with the tools \texttt{SECTIONS\_DIADCT} and is written in binary file \texttt{section\_ijglobal.diadct\_ORCA2\_LIM }  
     698which might be read by NEMO to compute on-line transports. 
     699 
     700On-line transports module will create three ascii files for outputs:  
     701 
     702- \texttt{volume\_transport} for volume transports (  unit: $10^{6} m^{3} s^{-1}$ ) 
     703 
     704- \texttt{heat\_transport}   for heat transports   (  unit: $10^{15} W $ ) 
     705 
     706- \texttt{salt\_transport}   for salt transports   (  unit: $10^{9}g s^{-1}$ ) 
     707 
     708 
     709%------------------------------------------namdct---------------------------------------------------- 
     710\namdisplay{namdct} 
     711%------------------------------------------------------------------------------------------------------------- 
     712 
     713Concerning the on-line transports computing, some parameters are available in namelist: 
     714 
     715- \texttt{nn\_dct} is the frequency of instantaneous transports computing 
     716 
     717- \texttt{nn\_dctwri} is the frequency of writing ( mean of instantaneous transports ) 
     718 
     719- \texttt{nn\_debug} to debug a section 
     720 
     721 
     722 
     723\subsubsection{ How to compute sections pathway and create the binary file? } 
     724 
     725In \texttt{NEMOGCM/TOOLS/SECTIONS\_DIADCT/run}, \texttt{ {list\_sections.ascii\_global}} 
     726contains all the sections we want to compute; this list of sections is inspired from MERSEA project metrics. 
     727 
     728Another file is available for GYRE configuration (\texttt{ {list\_sections.ascii\_GYRE}}).  
     729 
     730Each section is defined by: 
     731 
     732\noindent \texttt{ long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name }\\ 
     733with: 
     734 
     735- \texttt{long1 lat1} , coordinates of the first extremity of the section; 
     736 
     737- \texttt{long2 lat2} , coordinates of the second extremity of the section; 
     738 
     739- \texttt{nclass} the number of bounds of your classes ( 3 bounds for 2 classes, for example); 
     740 
     741- \texttt{okstrpond} to compute heat and salt transport, \texttt{nostrpond} if no; 
     742 
     743- \texttt{ice}  to compute surface and volume ice transports, \texttt{noice} if no. \\ 
     744 
     745 
     746\noindent Results of transports computing don't depend of the order of the 2 extremities in this file.\\  
     747 
     748 
     749\noindent If nclass =/ 0,the next lines will containes the class type and the nclass bounds: 
     750 
     751\texttt{long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name} 
     752 
     753\texttt{classtype} 
     754 
     755\texttt{zbound1} 
     756 
     757\texttt{zbound2} 
     758 
     759\texttt{.} 
     760 
     761\texttt{.} 
     762 
     763\texttt{nclass-1} 
     764 
     765\texttt{nclass} 
     766 
     767\noindent where \texttt{classtype} can be: 
     768 
     769- \texttt{zsal} for salinity classes 
     770 
     771- \texttt{ztem} for temperature classes 
     772 
     773- \texttt{zlay} for depth classes 
     774 
     775- \texttt{zsigi} for insitu density classes 
     776 
     777- \texttt{zsigp} for potential density classes \\ 
     778 
     779   
     780job.ksh compute the pathway for each section and create a binary file  
     781\texttt{section\_ijglobal.diadct\_ORCA2\_LIM } which will be read by NEMO. \\ 
     782 
     783 
     784Example of two sections: 
     785 
     786 
     787\noindent \texttt{ -68.    -54.5   -60.    -64.7  00 okstrpond noice ACC\_Drake\_Passage} 
     788 
     789\noindent \texttt{ -80.5    22.5   -80.5    25.5  05 nostrpond noice ATL\_Cuba\_Florida} 
     790 
     791\noindent \texttt{ztem} 
     792 
     793\noindent \texttt{-2.0} 
     794 
     795\noindent \texttt{ 4.5} 
     796 
     797\noindent \texttt{ 7.0} 
     798 
     799\noindent \texttt{12.0} 
     800 
     801\noindent \texttt{40.0} 
     802 
     803 
     804\subsubsection{ How to read the output files ? } 
     805 
     806The output format is : 
     807  
     808\small{\texttt{date, time-step number, section number, section name, section slope coefficient, class number,  
     809class name, class bound 1 , classe bound2, transport\_sens1 ,  transport\_sens2, transport\_total}}\\ 
     810 
     811 
     812If some sections have classes, the first lines correspond to the transport for each class  
     813and the last line corresponds to the total transport for all classes. If the other case, class number 
     814\texttt{ 1 } corresponds to \texttt{ total class } and this class is called  \texttt{N}, like \texttt{none}.\\ 
     815 
     816 
     817\noindent \texttt{ transport\_sens1 } is the positive part of the transport ( \texttt{ > = 0 } ). 
     818 
     819\noindent \texttt{ transport\_sens2 } is the negative part of the transport ( \texttt{ < = 0 } ).\\ 
     820 
     821 
     822\noindent \texttt{section slope coefficient} gives some informations about the significations of transports signs and sens:\\ 
     823 
     824 
     825 
     826\begin{tabular}{|c|c|c|c|p{4cm}|} 
     827\hline 
     828section slope coefficient & section type & sens 1 & sens 2 & total transport \\ \hline 
     8290.    &  horizontal & northward & southward & postive: northward  \\ \hline 
     8301000. &  vertical   & eastward  & westward  & postive: eastward  \\ \hline                 
     831\texttt{=/0, =/ 1000.}   &  diagonal   & eastward  & westward  & postive: eastward  \\ \hline                 
     832\end{tabular} 
     833 
     834 
     835 
    686836% ------------------------------------------------------------------------------------------------------------- 
    687837%       Other Diagnostics 
     
    735885\\ 
    736886\\ 
    737 A module is available to compute volume, heat and salt tranports through sections. This diagnostic 
    738 is actived with \key{diadct}. 
    739  
    740  
    741 %expliquer comment st ecrite les sorties, comment comprendre le sens 
    742  
    743 %------------------------------------------namdct---------------------------------------------------- 
    744 \namdisplay{namdct} 
    745 %------------------------------------------------------------------------------------------------------------- 
    746  
    747 In \texttt{NEMOGCM/TOOLS/SECTIONS\_DIADCT/run},\\ 
    748 the file \texttt{ {list\_sections.ascii\_keep}}\\ 
    749 contains all the sections we want to compute.\\ 
    750 Each section is defined by:\\ 
    751 \texttt{ long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name }\\ 
    752 with:\\ 
    753 - \texttt{long1 lat1} , coordinates of the first extremity of the section;\\ 
    754 - \texttt{long2 lat2} , coordinates of the second extremity of the section;\\ 
    755 - \texttt{nclass} the number of bounds of your classes ( 3 bounds for 2 classes, for example);\\ 
    756 - \texttt{okstrpond} to compute heat and salt transport, \texttt{nostrpond} if no;\\ 
    757 - \texttt{ice}  to compute surface and volume ice transports, \texttt{noice} if no.\\ 
    758 \\ 
    759 If nclass =/ 0,the next lines will containes the class type and the nclass bounds:\\ 
    760 \texttt{long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name}\\ 
    761 \texttt{classtype}\\ 
    762 \texttt{zbound1}\\ 
    763 \texttt{zbound2}\\ 
    764 \texttt{.}\\ 
    765 \texttt{.}\\ 
    766 \texttt{nclass-1}\\ 
    767 \texttt{nclass}\\ 
    768 where \texttt{classtype} can be:\\ 
    769 - \texttt{zsal} for salinity classes\\ 
    770 - \texttt{ztem} for temperature classes\\ 
    771 - \texttt{zlay} for depth classes\\ 
    772 - \texttt{zsigi} for insitu density classes\\ 
    773 - \texttt{zsigp} for potential density classes\\ 
    774 \\   
    775 job.ksh compute the pathway for each section and create a binary file  
    776 \texttt{section\_ijglobal.diadct\_ORCA2\_LIM } which will be read by NEMO. \\ 
    777 \\ 
    778 Example of two sections:\\ 
    779 \\ 
    780 \texttt{ -68.    -54.5   -60.    -64.7  00 okstrpond noice ACC\_Drake\_Passage}\\ 
    781 \texttt{ -80.5    22.5   -80.5    25.5  05 nostrpond noice ATL\_Cuba\_Florida}\\ 
    782 \texttt{ztem}\\ 
    783 \texttt{-2.0}\\ 
    784 \texttt{ 4.5}\\ 
    785 \texttt{ 7.0}\\ 
    786 \texttt{12.0}\\ 
    787 \texttt{40.0}\\ 
     887 
    788888 
    789889 
Note: See TracChangeset for help on using the changeset viewer.