!WRF:MEDIATION:IO
SUBROUTINE wrf_ext_read_field_arr(DataHandle,DateStr,Var & 4,6
,Field &
,idx4, idx5, idx6, idx7 &
,nx4 , nx5 , nx6 &
,TypeSizeInBytes &
,FieldType,Comm,IOComm &
,DomainDesc &
,bdy_mask &
,MemoryOrder &
,Stagger &
,debug_message &
,ds1, de1, ds2, de2, ds3, de3 &
,ms1, me1, ms2, me2, ms3, me3 &
,ps1, pe1, ps2, pe2, ps3, pe3, Status )
USE module_io
USE module_wrf_error
USE module_state_description
USE module_timing
IMPLICIT NONE
INTEGER, INTENT(IN) :: idx4, idx5, idx6, idx7
INTEGER, INTENT(IN) :: nx4 , nx5 , nx6
INTEGER, INTENT(IN) :: TypeSizeInBytes
INTEGER ,INTENT(IN ) :: DataHandle
CHARACTER*(*) ,INTENT(IN ) :: DateStr
CHARACTER*(*) ,INTENT(IN ) :: Var
INTEGER ,INTENT(INOUT) :: Field(*)
INTEGER ,INTENT(IN ) :: FieldType
INTEGER ,INTENT(IN ) :: Comm
INTEGER ,INTENT(IN ) :: IOComm
INTEGER ,INTENT(IN ) :: DomainDesc
CHARACTER*(*) ,INTENT(IN ) :: MemoryOrder
LOGICAL, DIMENSION(4) ,INTENT(IN ) :: bdy_mask
CHARACTER*(*) ,INTENT(IN ) :: Stagger
CHARACTER*(*) ,INTENT(IN ) :: debug_message
INTEGER , INTENT(IN ) :: ds1, de1, ds2, de2, ds3, de3, &
ms1, me1, ms2, me2, ms3, me3, &
ps1, pe1, ps2, pe2, ps3, pe3
INTEGER , INTENT(INOUT) :: Status
! Local
INTEGER tsfac ! Type size factor
CHARACTER*256 mess
tsfac = TypeSizeInBytes / IWORDSIZE
IF ( tsfac .LE. 0 ) THEN
CALL wrf_message('wrf_ext_read_field_arr')
WRITE(mess,*)'Internal error: email this message to wrfhelp@ucar.edu ',TypeSizeInBytes,IWORDSIZE
CALL wrf_error_fatal(mess)
ENDIF
CALL wrf_ext_read_field( DataHandle,DateStr,Var &
,Field(1 &
+tsfac*(0 &
+(idx4-1)*(me3-ms3+1)*(me2-ms2+1)*(me1-ms1+1) &
+(idx5-1)*nx4*(me3-ms3+1)*(me2-ms2+1)*(me1-ms1+1) &
+(idx6-1)*nx5*nx4*(me3-ms3+1)*(me2-ms2+1)*(me1-ms1+1) &
+(idx7-1)*nx6*nx5*nx4*(me3-ms3+1)*(me2-ms2+1)*(me1-ms1+1))) &
,FieldType,Comm,IOComm &
,DomainDesc &
,bdy_mask &
,MemoryOrder &
,Stagger &
,debug_message &
,ds1, de1, ds2, de2, ds3, de3 &
,ms1, me1, ms2, me2, ms3, me3 &
,ps1, pe1, ps2, pe2, ps3, pe3, Status )
END SUBROUTINE wrf_ext_read_field_arr
SUBROUTINE wrf_ext_read_field( DataHandle,DateStr,Var,Field,FieldType,Comm,IOComm, & 19,5
DomainDesc, bdy_mask, MemoryOrder,Stagger, &
debug_message , &
ds1, de1, ds2, de2, ds3, de3, &
ms1, me1, ms2, me2, ms3, me3, &
ps1, pe1, ps2, pe2, ps3, pe3, Status )
USE module_io
USE module_wrf_error
IMPLICIT NONE
integer :: DataHandle
character*(*) :: DateStr
character*(*) :: Var
integer :: Field(*)
integer :: FieldType
integer :: Comm
integer :: IOComm
integer :: DomainDesc
logical, dimension(4) :: bdy_mask
character*(*) :: MemoryOrder
character*(*) :: Stagger
character*(*) :: debug_message
INTEGER , INTENT(IN ) :: ds1, de1, ds2, de2, ds3, de3, &
ms1, me1, ms2, me2, ms3, me3, &
ps1, pe1, ps2, pe2, ps3, pe3
INTEGER itrace
INTEGER , DIMENSION(3) :: domain_start , domain_end
INTEGER , DIMENSION(3) :: memory_start , memory_end
INTEGER , DIMENSION(3) :: patch_start , patch_end
CHARACTER*80 , DIMENSION(3) :: dimnames
integer ,intent(inout) :: Status
domain_start(1) = ds1 ; domain_end(1) = de1 ;
patch_start(1) = ps1 ; patch_end(1) = pe1 ;
memory_start(1) = ms1 ; memory_end(1) = me1 ;
domain_start(2) = ds2 ; domain_end(2) = de2 ;
patch_start(2) = ps2 ; patch_end(2) = pe2 ;
memory_start(2) = ms2 ; memory_end(2) = me2 ;
domain_start(3) = ds3 ; domain_end(3) = de3 ;
patch_start(3) = ps3 ; patch_end(3) = pe3 ;
memory_start(3) = ms3 ; memory_end(3) = me3 ;
CALL debug_io_wrf ( debug_message,DateStr, &
domain_start,domain_end,patch_start,patch_end, &
memory_start,memory_end )
#if 0
Status = 1
if ( de1 - ds1 < 0 ) return
if ( de2 - ds2 < 0 ) return
if ( de3 - ds3 < 0 ) return
if ( pe1 - ps1 < 0 ) return
if ( pe2 - ps2 < 0 ) return
if ( pe3 - ps3 < 0 ) return
if ( me1 - ms1 < 0 ) return
if ( me2 - ms2 < 0 ) return
if ( me3 - ms3 < 0 ) return
#endif
Status = 0
CALL wrf_read_field ( &
DataHandle & ! DataHandle
,DateStr & ! DateStr
,Var & ! Data Name
,Field & ! Field
,FieldType & ! FieldType
,Comm & ! Comm
,IOComm & ! IOComm
,DomainDesc & ! DomainDesc
,bdy_mask & ! bdy_mask
,MemoryOrder & ! MemoryOrder
,Stagger & ! Stagger
,dimnames & ! JMMOD 1109
,domain_start & ! DomainStart
,domain_end & ! DomainEnd
,memory_start & ! MemoryStart
,memory_end & ! MemoryEnd
,patch_start & ! PatchStart
,patch_end & ! PatchEnd
,Status )
IF ( wrf_at_debug_level(300) ) THEN
WRITE(wrf_err_message,*) debug_message,' Status = ',Status
CALL wrf_message ( TRIM(wrf_err_message) )
ENDIF
END SUBROUTINE wrf_ext_read_field