[211] | 1 | SUBROUTINE user_init |
---|
| 2 | |
---|
| 3 | !------------------------------------------------------------------------------! |
---|
[484] | 4 | ! Current revisions: |
---|
[211] | 5 | ! ----------------- |
---|
[226] | 6 | ! |
---|
[211] | 7 | ! |
---|
| 8 | ! Former revisions: |
---|
| 9 | ! ----------------- |
---|
| 10 | ! $Id: user_init.f90 484 2010-02-05 07:36:54Z raasch $ |
---|
| 11 | ! |
---|
[226] | 12 | ! 211 2008-11-11 04:46:24Z raasch |
---|
| 13 | ! Former file user_interface.f90 split into one file per subroutine |
---|
| 14 | ! |
---|
[211] | 15 | ! Description: |
---|
| 16 | ! ------------ |
---|
| 17 | ! Execution of user-defined initializing actions |
---|
| 18 | !------------------------------------------------------------------------------! |
---|
| 19 | |
---|
| 20 | USE control_parameters |
---|
| 21 | USE indices |
---|
| 22 | USE netcdf_control |
---|
| 23 | USE pegrid |
---|
| 24 | USE user |
---|
| 25 | |
---|
| 26 | IMPLICIT NONE |
---|
| 27 | |
---|
| 28 | CHARACTER (LEN=20) :: field_char |
---|
| 29 | ! |
---|
| 30 | !-- Here the user-defined initializing actions follow: |
---|
| 31 | !-- Sample for user-defined output |
---|
| 32 | ! ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) |
---|
| 33 | ! ALLOCATE( ustvst(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ); ustvst = 0.0 |
---|
| 34 | |
---|
| 35 | !-- Sample for user-defined time series |
---|
| 36 | !-- For each time series quantity you have to give a label and a unit, |
---|
| 37 | !-- which will be used for the NetCDF file. They must not contain more than |
---|
| 38 | !-- seven characters. The value of dots_num has to be increased by the |
---|
| 39 | !-- number of new time series quantities. Its old value has to be store in |
---|
| 40 | !-- dots_num_palm. See routine user_statistics on how to output calculate |
---|
| 41 | !-- and output these quantities. |
---|
| 42 | ! dots_label(dots_num+1) = 'abs_umx' |
---|
| 43 | ! dots_unit(dots_num+1) = 'm/s' |
---|
| 44 | ! dots_label(dots_num+2) = 'abs_vmx' |
---|
| 45 | ! dots_unit(dots_num+2) = 'm/s' |
---|
| 46 | ! |
---|
| 47 | ! dots_num_palm = dots_num |
---|
| 48 | ! dots_num = dots_num + 2 |
---|
| 49 | |
---|
| 50 | END SUBROUTINE user_init |
---|
| 51 | |
---|