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

Timestamp:

Feb 1, 2008 12:41:46 PM (17 years ago)

Author:

raasch

Message:

further updates for turbulent inflow: reading input data of a precursor run using a smaller total domain is working

File:

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

palm/trunk/SOURCE/user_interface.f90

@@ -177 +177 @@
 
 
- SUBROUTINE user_read_restart_data
+ SUBROUTINE user_read_restart_data( nxlc, nxlf, nxl_on_file, nxrc, nxrf,  &
+                                    nxr_on_file, nync, nynf, nyn_on_file, &
+                                    nysc, nysf, nys_on_file, tmp_2d, tmp_3d )
 
 !------------------------------------------------------------------------------!
@@ -184 +186 @@
 ! ------------
 ! Reading restart data from file(s)
+! Subdomain index limits on file are given by nxl_on_file, etc.
+! Indices nxlc, etc. indicate the range of gridpoints to be mapped from the
+! subdomain on file (f) to the subdomain of the current PE (c). They have been
+! calculated in routine read_3d_binary.
 !------------------------------------------------------------------------------!
 
@@ -194 +200 @@
 
     CHARACTER (LEN=20) :: field_char
+
+    INTEGER ::  nxlc, nxlf, nxl_on_file, nxrc, nxrf, nxr_on_file, nync, nynf, &
+                nyn_on_file, nysc, nysf, nys_on_file
+
+    REAL, DIMENSION(nys_on_file-1:nyn_on_file+1,nxl_on_file-1:nxr_on_file+1) ::&
+          tmp_2d
+
+    REAL, DIMENSION(nzb:nzt+1,nys_on_file-1:nyn_on_file+1, &
+                    nxl_on_file-1:nxr_on_file+1) ::        &
+          tmp_3d
+
 !
 !-- Here the reading of user-defined restart data follows:
@@ -207 +224 @@
 !
 !             CASE ( 'u2_av' )
-!                ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
+!                IF ( .NOT. ALLOCATED( u2_av )  THEN
+!                   ALLOCATE( u2_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
+!                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)
 !
 !             CASE DEFAULT