Changes between Version 58 and Version 59 of doc/app/userint/output


Ignore:
Timestamp:
Oct 11, 2010 11:22:33 AM (14 years ago)
Author:
heinze
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • doc/app/userint/output

    v58 v59  
    7777   3. The data array has to be allocated in subroutine [../int#user_init user_init]: \\
    7878
    79    {{{ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) }}} \\
    80 
    81    In case that output of time averaged data is requested, the array containing the sum has possibly to be read from the restart file (local filename [../../iofiles#BININ BININ]) by executing the following code in [../int#user_init user_init]:
    82 
    83    {{{IF ( initializing_actions == 'read_restart_data' )  THEN }}} \\
    84       {{{    READ ( 13 )  field_char }}} \\
    85       {{{    DO  WHILE ( TRIM( field_char ) /= '*** end user ***' ) }}} \\
    86              
    87          {{{    SELECT CASE ( TRIM( field_char ) ) }}} \\
    88 
    89            {{{     CASE ( 'u2_av' ) }}} \\
    90               {{{     ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) }}} \\
    91               {{{     READ ( 13 )  u2_av }}} \\
    92 
    93            {{{     CASE DEFAULT }}} \\
    94               {{{     PRINT*, '+++ user_init: unknown variable named "', & }}}
    95                  [[span(   {{{TRIM( field_char ), '" found in' }}},style=color: black)]] \\
    96               {{{     PRINT*, '               data from prior run on PE ', myid }}} \\
    97               {{{     CALL local_stop }}} \\
    98           {{{    END SELECT  }}} \\
    99       {{{    ENDDO }}} \\
    100    {{{    ENDIF }}} \\
     79   {{{ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) }}}
    10180
    10281   4. The quantity has to be given a unit (subroutine [../int#user_check_data_output user_check_data_output]): \\
     
    240219  Otherwise, the calculated time-average may be wrong. In the restart run, this quantity has to be read from the restart file by including the following code in subroutine [../#user_read_restart_data user_read_restart_data]:
    241220
    242   {{{   IF ( initializing_actions == 'read_restart_data' )  THEN }}}\\
    243       {{{       READ ( 13 )  field_char }}}\\
    244       {{{       DO  WHILE ( TRIM( field_char ) /= '*** end user ***' ) }}}\\
    245          {{{       SELECT CASE ( TRIM( field_char ) ) }}}\\
    246             {{{       CASE ( 'u2_av' ) }}}\\
    247                {{{       IF ( .NOT. ALLOCATED( u2_av ) ) THEN }}}\\
    248                   {{{       ALLOCATE( u2_av(nzb:nzt+1,nysd:nynd,nxld:nxrd) ) }}}\\
    249                {{{       ENDIF }}}\\
    250                {{{       READ ( 13 )  tmp_3d }}}\\
    251                   {{{    u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) }}}\\
    252 
    253             {{{       CASE DEFAULT }}}\\
    254                {{{       PRINT*, '+++ user_init: unknown variable named "', & }}}\\
    255                {{{               TRIM( field_char ), '" found in' }}}\\
    256                {{{       PRINT*, '               data from prior run on PE ', myid }}}\\
    257                {{{       CALL local_stop }}}\\
    258          {{{       END SELECT }}}\\
    259          {{{       READ ( 13 )  field_char }}}\\
    260       {{{    ENDDO }}}\\
     221  {{{    IF ( initializing_actions == 'read_restart_data' )  THEN }}}\\
     222      {{{      READ ( 13 )  field_char }}}\\
     223      {{{      DO  WHILE ( TRIM( field_char ) /= '*** end user ***' ) }}}\\
     224         {{{      DO  k = 1, overlap_count }}}\\
     225             {{{       nxlf = nxlfa(i,k) }}}\\
     226             {{{       nxlc = nxlfa(i,k) + offset_xa(i,k)}}}\\
     227             {{{       nxrf = nxrfa(i,k)}}}\\
     228             {{{       nxrc = nxrfa(i,k) + offset_xa(i,k)}}}\\
     229             {{{       nysf = nysfa(i,k)}}}\\
     230             {{{       nysc = nysfa(i,k) + offset_ya(i,k)}}}\\
     231             {{{       nynf = nynfa(i,k)}}}\\
     232             {{{       nync = nynfa(i,k) + offset_ya(i,k)}}}\\\\
     233             {{{       SELECT CASE ( TRIM( field_char ) ) }}}\\
     234                 {{{        CASE ( 'u2_av' ) }}}\\
     235                    {{{        IF ( .NOT. ALLOCATED( u2_av ) ) THEN }}}\\
     236                       {{{        ALLOCATE( u2_av(nzb:nzt+1,nysd:nynd,nxld:nxrd) ) }}}\\
     237                    {{{        ENDIF }}}\\
     238                    {{{        IF ( k == 1 ) READ ( 13 )  tmp_3d }}}\\
     239                    {{{        u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = & }}}\\
     240                    {{{            tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1) }}}\\
     241 
     242                 {{{        CASE DEFAULT }}}\\
     243                    {{{        WRITE( message_string, * ) 'unknown variable named "',  &}}}\\
     244                    {{{             TRIM( field_char ), '" found in', &}}}\\
     245                    {{{             '&data from prior run on PE ', myid}}}\\
     246                    {{{         CALL message('user_read_restart_data', 'UI0012', 1, 2, 0, 6, 0) }}}\\
     247             {{{       END SELECT }}}\\
     248         {{{      ENDDO          }}}\\
     249         {{{      READ ( 13 )  field_char }}}\\
     250      {{{      ENDDO }}}\\
    261251  {{{    ENDIF }}}
    262252\\\\