Changeset 2902
- Timestamp:
- 2011-10-11T15:13:47+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR10_diadct/DOC/TexFiles/Chapters/Chap_DIA.tex
r2871 r2902 684 684 the off-line use of this marvellous diagnostic tool. 685 685 686 687 % ------------------------------------------------------------------------------------------------------------- 688 % Sections transports 689 % ------------------------------------------------------------------------------------------------------------- 690 \section{Transports across sections (\key{diadct}) } 691 \label{DIA_diag_dct} 692 693 A module is available to compute volume, heat and salt tranports through sections. This diagnostic 694 is actived with \key{diadct}. 695 696 Each section is defined by the coordinates of its 2 extremities. The pathway between them is contructed 697 with the tools \texttt{SECTIONS\_DIADCT} and is written in binary file \texttt{section\_ijglobal.diadct\_ORCA2\_LIM } 698 which might be read by NEMO to compute on-line transports. 699 700 On-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 713 Concerning 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 725 In \texttt{NEMOGCM/TOOLS/SECTIONS\_DIADCT/run}, \texttt{ {list\_sections.ascii\_global}} 726 contains all the sections we want to compute; this list of sections is inspired from MERSEA project metrics. 727 728 Another file is available for GYRE configuration (\texttt{ {list\_sections.ascii\_GYRE}}). 729 730 Each section is defined by: 731 732 \noindent \texttt{ long1 lat1 long2 lat2 nclass (ok/no)strpond (no)ice section\_name }\\ 733 with: 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 780 job.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 784 Example 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 806 The output format is : 807 808 \small{\texttt{date, time-step number, section number, section name, section slope coefficient, class number, 809 class name, class bound 1 , classe bound2, transport\_sens1 , transport\_sens2, transport\_total}}\\ 810 811 812 If some sections have classes, the first lines correspond to the transport for each class 813 and 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 828 section slope coefficient & section type & sens 1 & sens 2 & total transport \\ \hline 829 0. & horizontal & northward & southward & postive: northward \\ \hline 830 1000. & vertical & eastward & westward & postive: eastward \\ \hline 831 \texttt{=/0, =/ 1000.} & diagonal & eastward & westward & postive: eastward \\ \hline 832 \end{tabular} 833 834 835 686 836 % ------------------------------------------------------------------------------------------------------------- 687 837 % Other Diagnostics … … 735 885 \\ 736 886 \\ 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 788 888 789 889
Note: See TracChangeset
for help on using the changeset viewer.