Changes between Version 14 and Version 15 of doc/app/userint/subd
- Timestamp:
- Jan 3, 2011 12:38:45 PM (14 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
doc/app/userint/subd
v14 v15 5 5 The exact definition of these domains has to be made by the user within the user-defined routine [../int#user_init user_init]. The domains are defined with a mask array named {{{rmask}}}, which has to be given the value ''1.0'' for all horizontal grid points belonging to the user-defined domain and the value ''0.0'', where grid points do not belong to the user-defined domain. In the model {{{rmask}}} is declared as: 6 6 {{{ 7 REAL :: rmask(nys d:nynd,nxld:nxrd,0:9) .7 REAL :: rmask(nysg:nyng,nxlg:nxrg,0:9) . 8 8 }}} 9 9 The first two indices denote the array bounds (including the ghost points) in y and x-direction for each domain (see the corresponding description of [../../parallel#parallelization Parallelization]; don't confuse this with the user-defined domain!). The third index determines the user-defined domain, where 0 indicates the total model domain and ''1'' to ''9'' the user-defined domains.\\\\ … … 19 19 disc_center_y = dy * (ny + 1)/2 20 20 disc_radius = 0.5 * disc_center_x 21 DO i = nxl d, nxrd21 DO i = nxlg, nxrg 22 22 x = i * dx 23 DO j = nys d, nynd23 DO j = nysg, nyng 24 24 y = j * dy 25 25 radial_distance = SQRT( ( x - disc_center_x )**2 + &