Changes between Version 7 and Version 8 of doc/app/userint/subd

Timestamp:

Sep 15, 2010 9:30:49 AM (14 years ago)

Author:

suehring

Comment:

--

doc/app/userint/subd

@@ -6 +6 @@
 REAL :: rmask(nysd:nynd,nxld:nxrd,0:9) .
 }}}
-The first two indices denote the array bounds (including the ghost points) in y and x-direction for each subdomain (see the corresponding description of [../../parallel#parallelization Parallelization]; don't confuse this with the user-defined subdomain!). The third index determines the user-defined subdomain., where 0 indicates the total model domain and ''1'' to ''9'' the user-defined subdomains.\\\\
+The first two indices denote the array bounds (including the ghost points) in y and x-direction for each subdomain (see the corresponding description of [../../parallel#parallelization Parallelization]; don't confuse this with the user-defined subdomain!). The third index determines the user-defined subdomain, where 0 indicates the total model domain and ''1'' to ''9'' the user-defined subdomains.\\\\
 The following example should illustrate this. Two subdomains are defined by the user. The first is determined by all grid points which lie within a circle whose center is equal to the (horizontal) center of the model domain and whose diameter is equal to half of the total horizontal domain size (square total domain assumed). The second subdomain should be defined by all points outside of this domain. This may be obtained by the following lines of code in [..\int#user_init user_init]:
 {{{
@@ -18 +18 @@
 disc_center_y = dy * (ny + 1)/2
 disc_radius = 0.5 * disc_center_x
-DO  i = nxl-1, nxr+1
+DO  i = nxld, nxrd
    x = i * dx
-   DO  j = nys-1, nyn+1
+   DO  j = nysd, nynd
       y = j * dy
       radial_distance = SQRT( ( x - disc_center_x )**2 + &