Ignore:
Timestamp:
Feb 13, 2008 4:06:08 PM (14 years ago)
Author:
letzel
Message:

Small bugfixes in the user_interface.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/user_interface.f90

    r147 r148  
    44! Actual revisions:
    55! -----------------
     6! small bugfixes for sample code (comments):
     7! - initialize ustvst with 0.0 as it is now computed only until nxr and nyn
     8! - two ALLOCATE statements moved from user_read_restart_data back to user_init
     9! - remove 'READ (13) u2_av' statement in user_read_restart_data
    610! +routines user_read_restart_data, user_spectra
    711!
     
    214218!-- Here the reading of user-defined restart data follows:
    215219!-- Sample for user-defined output
    216 !    ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    217 !    ALLOCATE( ustvst(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    218220!
    219221!    IF ( initializing_actions == 'read_restart_data' )  THEN
     
    228230!                ENDIF
    229231!                READ ( 13 )  tmp_3d
    230 !                READ ( 13 )  u2_av
    231232!                u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
    232233!                                          tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
     
    269270!
    270271!-- Here the user-defined initializing actions follow:
     272!-- Sample for user-defined output
     273!    ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
     274!    ALLOCATE( ustvst(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) );  ustvst = 0.0
     275
    271276!-- Sample for user-defined time series
    272277!-- For each time series quantity you have to give a label and a unit,
     
    502507!                ENDDO
    503508!             ENDDO
    504 !             DO  i = nxl-1, nxr+1
    505 !                DO  j = nys-1, nyn+1
     509!             DO  i = nxl-1, nxr
     510!                DO  j = nys-1, nyn
    506511!                   DO  k = nzb, nzt+1
    507512!                      ustvst(k,j,i) =  &
Note: See TracChangeset for help on using the changeset viewer.