|
Last change
on this file since 4 was
4,
checked in by xlvlod, 18 years ago
|
|
import initial des vraies codes.
|
|
File size:
1.0 KB
|
| Line | |
|---|
| 1 | subroutine user_BCs |
|---|
| 2 | #define DEBUG_LEVEL 1 |
|---|
| 3 | |
|---|
| 4 | !********** Nearly everything available via modules ***** |
|---|
| 5 | use grid_info |
|---|
| 6 | use dependent_variables |
|---|
| 7 | use dimensional_scales |
|---|
| 8 | use mpi_parameters |
|---|
| 9 | use user_parameters |
|---|
| 10 | use counters_flags_etc |
|---|
| 11 | use user_routines |
|---|
| 12 | |
|---|
| 13 | implicit none |
|---|
| 14 | |
|---|
| 15 | ! local variables |
|---|
| 16 | include '../input/problem_size.h' |
|---|
| 17 | real, dimension(nx,ny,locnz) :: vals |
|---|
| 18 | integer :: field_id |
|---|
| 19 | !********************************************************* |
|---|
| 20 | |
|---|
| 21 | if( user_BCs_flag == 'no' ) return |
|---|
| 22 | |
|---|
| 23 | #if DEBUG_LEVEL >=1 |
|---|
| 24 | if(myid == 0) write(0,*) 'hello world from user_BCs' |
|---|
| 25 | #endif |
|---|
| 26 | |
|---|
| 27 | if ( istep >= 1 ) then |
|---|
| 28 | !============== |
|---|
| 29 | field_id = 1 !!! imposed surface current only for "U" |
|---|
| 30 | |
|---|
| 31 | ! do i = 1,nx |
|---|
| 32 | call user_forcing( Grid(0)%x(:,1)*L_scale , Grid(0)%y(:)*L_scale , Grid(0)%z(1,:)*L_scale , & |
|---|
| 33 | t*t_scale , vals , field_id , & |
|---|
| 34 | nx , ny , locnz ) |
|---|
| 35 | ! u(i,:,:) = u(i,:,:) + vals/U_scale |
|---|
| 36 | ! enddo |
|---|
| 37 | u = u + vals/U_scale |
|---|
| 38 | |
|---|
| 39 | !============== |
|---|
| 40 | |
|---|
| 41 | endif |
|---|
| 42 | |
|---|
| 43 | |
|---|
| 44 | end subroutine user_BCs |
|---|
Note: See
TracBrowser
for help on using the repository browser.
