Changeset 148 for palm/trunk/SOURCE/user_interface.f90

Timestamp:

Feb 13, 2008 4:06:08 PM (17 years ago)

Author:

letzel

Message:

Small bugfixes in the user_interface.

File:

• palm/trunk/SOURCE/user_interface.f90 (modified) (5 diffs)

palm/trunk/SOURCE/user_interface.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! small bugfixes for sample code (comments):
+! - initialize ustvst with 0.0 as it is now computed only until nxr and nyn
+! - two ALLOCATE statements moved from user_read_restart_data back to user_init
+! - remove 'READ (13) u2_av' statement in user_read_restart_data
 ! +routines user_read_restart_data, user_spectra
 !
@@ -214 +218 @@
 !-- Here the reading of user-defined restart data follows:
 !-- Sample for user-defined output
-!    ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
-!    ALLOCATE( ustvst(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
 !
 !    IF ( initializing_actions == 'read_restart_data' )  THEN
@@ -228 +230 @@
 !                ENDIF
 !                READ ( 13 )  tmp_3d
-!                READ ( 13 )  u2_av
 !                u2_av(:,nysc-1:nync+1,nxlc-1:nxrc+1) = &
 !                                          tmp_3d(:,nysf-1:nynf+1,nxlf-1:nxrf+1)
@@ -269 +270 @@
 !
 !-- Here the user-defined initializing actions follow:
+!-- Sample for user-defined output
+!    ALLOCATE( u2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
+!    ALLOCATE( ustvst(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) );  ustvst = 0.0
+
 !-- Sample for user-defined time series
 !-- For each time series quantity you have to give a label and a unit,
@@ -502 +507 @@
 !                ENDDO
 !             ENDDO
-!             DO  i = nxl-1, nxr+1
-!                DO  j = nys-1, nyn+1
+!             DO  i = nxl-1, nxr
+!                DO  j = nys-1, nyn
 !                   DO  k = nzb, nzt+1
 !                      ustvst(k,j,i) =  &