1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> |
---|
2 | <html><head> |
---|
3 | <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM |
---|
4 | chapter 3.5.3</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;14001660"> <meta name="CHANGED" content="20041112;14483720"> <meta name="KEYWORDS" content="parallel LES model"> <style> |
---|
5 | <!-- |
---|
6 | @page { size: 21cm 29.7cm } |
---|
7 | --> |
---|
8 | </style></head> |
---|
9 | <body style="direction: ltr;" lang="en-US"><h3 style="line-height: 100%;">3.5.3 Definition of user-defined |
---|
10 | subdomains</h3> |
---|
11 | <p style="line-height: 100%;">By default, the values of |
---|
12 | the timeseries quantities and the horizontally averaged vertical |
---|
13 | profiles (saved in local files |
---|
14 | <a href="chapter_3.4.html#PLOT1D_DATA">PLOT1D_DATA</a> |
---|
15 | and <a href="chapter_3.4.html#LIST_PROFIL">LIST_PROFIL</a>) |
---|
16 | always refer to the total model domain. Independently, up to 9 time |
---|
17 | series or profiles for different user-defined subdomains can be |
---|
18 | computed and plotted additionally. Steering in principle is done |
---|
19 | via the initialization parameter <a href="chapter_4.1.html#statistic_regions">statistic_regions</a>. |
---|
20 | </p> |
---|
21 | <p style="line-height: 100%;">The exact definition of |
---|
22 | these subdomains |
---|
23 | has to be made by the user within the user-defined subroutine |
---|
24 | <tt><font style="font-size: 10pt;" size="2">init_user</font></tt>. |
---|
25 | The |
---|
26 | subdomains are defined with a mask array named <tt><font style="font-size: 10pt;" size="2">rmask</font></tt>, |
---|
27 | which has to be given the value 1.0 for all horizontal grid points |
---|
28 | belonging to |
---|
29 | the user-defined subdomain and the value 0.0, where grid points do not |
---|
30 | belong |
---|
31 | to the user-defined subdomain. In the model <tt><font size="2">rmask</font></tt> |
---|
32 | is |
---|
33 | declared as: </p> |
---|
34 | <ul> <p style="line-height: 100%;"><tt><font style="font-size: 10pt;" size="2">REAL :: rmask |
---|
35 | (nys-1:nyn+1,nxl-1,nxr+1,0:9) .</font></tt></p> |
---|
36 | </ul><p style="line-height: 100%;">The first two |
---|
37 | indices are the grid point |
---|
38 | indices in y and x-direction. With parallel model runs <tt><font style="font-size: 10pt;" size="2">nxl</font></tt><font style="font-size: 10pt;" size="2">, |
---|
39 | </font><tt><font style="font-size: 10pt;" size="2">nxr</font></tt><font style="font-size: 10pt;" size="2">, |
---|
40 | </font><tt><font style="font-size: 10pt;" size="2">nys</font></tt><font style="font-size: 10pt;" size="2"> |
---|
41 | <font size="3">and</font> </font><tt><font style="font-size: 10pt;" size="2">nyn</font></tt> |
---|
42 | are the array bounds of the respective subdomain (don't confuse this |
---|
43 | with the user-defined subdomain!) on the |
---|
44 | respective processor. With runs on one processor <tt><font style="font-size: 10pt;" size="2">nys</font></tt><font style="font-size: 10pt;" size="2"> |
---|
45 | = </font><tt><font style="font-size: 10pt;" size="2">nxl</font></tt><font style="font-size: 10pt;" size="2"> |
---|
46 | = 0</font> and <font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">nxr |
---|
47 | =</font> </font><a href="chapter_4.1.html#nx">nx</a> |
---|
48 | and <font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">nyn |
---|
49 | =</font></font> <a href="chapter_4.1.html#ny">ny</a>. |
---|
50 | The third index determines the user-defined subdomain. The total model |
---|
51 | domain |
---|
52 | carries the index 0, the user-defined subdomains have the values 1 to |
---|
53 | 9. </p> |
---|
54 | <p style="line-height: 100%;">The following example should |
---|
55 | illustrate |
---|
56 | this. Two subdomains are defined by the user. The first is determined |
---|
57 | by all grid |
---|
58 | points which lie within a circle whose center is equal to the |
---|
59 | (horizontal) center of the model domain and whose diameter is equal |
---|
60 | to half of the total horizontal domain size (square total domain |
---|
61 | assumed). |
---|
62 | The second subdomain should be defined by all points outside of this |
---|
63 | domain. This may be obtained by the following lines of code in <tt><font style="font-size: 10pt;" size="2">user_init</font></tt>: |
---|
64 | </p> |
---|
65 | <ul> <p style="line-height: 100%;"><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">USE |
---|
66 | grid_variables</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">USE indices</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">USE statistics</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">.</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">.</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">.</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">disc_center_x = dx * (nx + 1)/2</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">disc_center_y = dy * (ny + 1)/2</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">disc_radius = 0.5 * |
---|
67 | disc_center_x</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">DO i = nxl-1, nxr+1</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> x = i * dx</font></font></tt><br> |
---|
68 | <tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
69 | DO j = nys-1, nyn+1</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
70 | y = j * dy</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
71 | radial_distance = SQRT( ( x - disc_center_x )**2 + &</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
72 | ( y - disc_center_y )**2 )</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
73 | IF ( radial_distance |
---|
74 | > disc_radius ) THEN</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
75 | rmask(j,i,1) = 0.0</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
76 | rmask(j,i,2) = 1.0</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
77 | ELSE</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
78 | rmask(j,i,1) = 1.0</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
79 | rmask(j,i,2) = 0.0</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> |
---|
80 | ENDIF</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> <br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"> ENDDO</font></font></tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace"><br> </font></font><tt><font style="font-size: 10pt;" size="2"><font face="Cumberland, monospace">ENDDO</font></font></tt></p> |
---|
81 | </ul><p style="line-height: 100%;">The module <span style="font-family: monospace;">statistics</span> must |
---|
82 | be used, |
---|
83 | because it contains <span style="font-family: monospace;">rmask</span> |
---|
84 | and the modules <span style="font-family: monospace;">grid_variables</span> |
---|
85 | and <span style="font-family: monospace;">indices</span> |
---|
86 | are |
---|
87 | necessary in this example, because grid spacing and indices are used. |
---|
88 | All array elements of <span style="font-family: monospace;">rmask</span> |
---|
89 | (<span style="font-family: monospace;">rmask(:,:,:)</span>) |
---|
90 | are preset |
---|
91 | by the model with 1.0. In no case this assignment must be |
---|
92 | changed for the total domain (<span style="font-family: monospace;">rmask(:,:,0)</span>)! |
---|
93 | Computations and output for the user-defined subdomains only take place |
---|
94 | if |
---|
95 | the user sets <a href="chapter_4.1.html#statistic_regions">statistic_regions</a> |
---|
96 | ≥ <i>1</i>. Beyond that, names for the user-defined |
---|
97 | subdomains can be |
---|
98 | assigned |
---|
99 | via the initialization parameter <a href="chapter_4.3.html#region">region</a>. |
---|
100 | Output of the names of the selected user-defined subdomains happens in |
---|
101 | the local files <a href="chapter_3.4.html#HEADER">HEADER</a> |
---|
102 | and <a href="chapter_3.4.html#RUN_CONTROL">RUN_CONTROL</a> |
---|
103 | within the user-defined subroutine <tt><font style="font-size: 10pt;" size="2">user_header</font></tt>. |
---|
104 | <br> |
---|
105 | </p> |
---|
106 | <hr><p style="line-height: 100%;"><br> |
---|
107 | <font color="#000080"><font color="#000080"><a href="chapter_3.5.2.html"><font color="#000080"><img src="left.gif" name="Grafik1" align="bottom" border="2" height="32" width="32"></font></a><a href="index.html"><font color="#000080"><img src="up.gif" name="Grafik2" align="bottom" border="2" height="32" width="32"></font></a><a href="chapter_3.5.4.html"><font color="#000080"><img src="right.gif" name="Grafik3" align="bottom" border="2" height="32" width="32"></font></a></font></font></p><p style="line-height: 100%;"><i>Last change: </i> |
---|
108 | $Id: chapter_3.5.3.html 62 2007-03-13 02:52:40Z heinze $</p> |
---|
109 | </body></html> |
---|