SUBROUTINE COPYSB(LUNIN,LUNOT,IRET) C$$$ SUBPROGRAM DOCUMENTATION BLOCK C C SUBPROGRAM: COPYSB C PRGMMR: WOOLLEN ORG: NP20 DATE: 1994-01-06 C C ABSTRACT: THIS SUBROUTINE COPIES A PACKED DATA SUBSET, INTACT, FROM C LOGICAL UNIT LUNIN, OPENED FOR INPUT VIA A PREVIOUS CALL TO BUFR C ARCHIVE LIBRARY SUBROUTINE OPENBF, TO LOGICAL UNIT LUNOT, OPENED C FOR OUTPUT VIA A PREVIOUS CALL TO OPENBF. THE BUFR MESSAGE MUST C HAVE BEEN PREVIOUSLY READ FROM UNIT LUNIT USING BUFR ARCHIVE C LIBRARY SUBROUTINE READMG OR READERME AND MAY BE EITHER COMPRESSED C OR UNCOMPRESSED. ALSO, BUFR ARCHIVE LIBRARY SUBROUTINE OPENMG OR C OPENMB MUST HAVE BEEN PREVIOUSLY CALLED TO OPEN AND INITIALIZE A C BUFR MESSAGE WITHIN MEMORY FOR UNIT LUNOT. EACH CALL TO COPYSB C ADVANCES THE POINTER TO THE BEGINNING OF THE NEXT SUBSET IN BOTH C THE INPUT AND OUTPUT FILES, UNLESS INPUT PARAMETER LUNOT IS .LE. C ZERO, IN WHICH CASE THE OUTPUT POINTER IS NOT ADVANCED. THE C COMPRESSION STATUS OF THE OUTPUT SUBSET/BUFR MESSAGE WILL ALWAYS C MATCH THAT OF THE INPUT SUBSET/BUFR MESSAGE {I.E., IF INPUT MESSAGE C IS UNCOMPRESSED(COMPRESSED) OUTPUT MESSAGE WILL BE UNCOMPRESSED C (COMPRESSED)}. C C PROGRAM HISTORY LOG: C 1994-01-06 J. WOOLLEN -- ORIGINAL AUTHOR C 1998-07-08 J. WOOLLEN -- REPLACED CALL TO CRAY LIBRARY ROUTINE C "ABORT" WITH CALL TO NEW INTERNAL BUFRLIB C ROUTINE "BORT" C 1999-11-18 J. WOOLLEN -- THE NUMBER OF BUFR FILES WHICH CAN BE C OPENED AT ONE TIME INCREASED FROM 10 TO 32 C (NECESSARY IN ORDER TO PROCESS MULTIPLE C BUFR FILES UNDER THE MPI) C 2000-09-19 J. WOOLLEN -- MAXIMUM MESSAGE LENGTH INCREASED FROM C 10,000 TO 20,000 BYTES C 2002-05-14 J. WOOLLEN -- REMOVED OLD CRAY COMPILER DIRECTIVES C 2003-11-04 S. BENDER -- ADDED REMARKS/BUFRLIB ROUTINE C INTERDEPENDENCIES C 2003-11-04 D. KEYSER -- UNIFIED/PORTABLE FOR WRF; ADDED C DOCUMENTATION (INCLUDING HISTORY); OUTPUTS C MORE COMPLETE DIAGNOSTIC INFO WHEN ROUTINE C TERMINATES ABNORMALLY C 2004-08-09 J. ATOR -- MAXIMUM MESSAGE LENGTH INCREASED FROM C 20,000 TO 50,000 BYTES C 2005-09-16 J. WOOLLEN -- NOW WRITES OUT COMPRESSED SUBSET/MESSAGE IF C INPUT SUBSET/MESSAGE IS COMPRESSED (BEFORE C COULD ONLY WRITE OUT UNCOMPRESSED SUBSET/ C MESSAGE REGARDLESS OF COMPRESSION STATUS OF C INPUT SUBSET/MESSAGE) C 2009-06-26 J. ATOR -- USE IOK2CPY C C USAGE: CALL COPYSB ( LUNIN, LUNOT, IRET ) C INPUT ARGUMENT LIST: C LUNIN - INTEGER: FORTRAN LOGICAL UNIT NUMBER FOR INPUT BUFR C FILE C LUNOT - INTEGER: FORTRAN LOGICAL UNIT NUMBER FOR OUTPUT BUFR C FILE C C OUTPUT ARGUMENT LIST: C IRET - INTEGER: RETURN CODE: C 0 = normal return C -1 = there are no more subsets in the input C BUFR message C C REMARKS: C THIS ROUTINE CALLS: BORT CMPMSG CPYUPD IOK2CPY C MESGBC READSB STATUS UFBCPY C UPB WRITSB C THIS ROUTINE IS CALLED BY: ICOPYSB C Also called by application programs. C C ATTRIBUTES: C LANGUAGE: FORTRAN 77 C MACHINE: PORTABLE TO ALL PLATFORMS C C$$$ INCLUDE 'bufrlib.prm' COMMON /MSGCWD/ NMSG(NFILES),NSUB(NFILES),MSUB(NFILES), . INODE(NFILES),IDATE(NFILES) COMMON /BITBUF/ MAXBYT,IBIT,IBAY(MXMSGLD4),MBYT(NFILES), . MBAY(MXMSGLD4,NFILES) COMMON /BTABLES/ MAXTAB,NTAB,TAG(MAXJL),TYP(MAXJL),KNT(MAXJL), . JUMP(MAXJL),LINK(MAXJL),JMPB(MAXJL), . IBT(MAXJL),IRF(MAXJL),ISC(MAXJL), . ITP(MAXJL),VALI(MAXJL),KNTI(MAXJL), . ISEQ(MAXJL,2),JSEQ(MAXJL) CHARACTER*10 TAG CHARACTER*3 TYP CHARACTER*128 BORT_STR C----------------------------------------------------------------------- C----------------------------------------------------------------------- IRET = 0 C CHECK THE FILE STATUSES C ----------------------- CALL STATUS(LUNIN,LIN,IL,IM) IF(IL.EQ.0) GOTO 900 IF(IL.GT.0) GOTO 901 IF(IM.EQ.0) GOTO 902 IF(LUNOT.GT.0) THEN CALL STATUS(LUNOT,LOT,IL,IM) IF(IL.EQ.0) GOTO 903 IF(IL.LT.0) GOTO 904 IF(IM.EQ.0) GOTO 905 IF(INODE(LIN).NE.INODE(LOT)) THEN IF( (TAG(INODE(LIN)).NE.TAG(INODE(LOT))) .OR. . (IOK2CPY(LIN,LOT).NE.1) ) GOTO 906 ENDIF ENDIF C SEE IF THERE IS ANOTHER SUBSET IN THE MESSAGE C --------------------------------------------- IF(NSUB(LIN).EQ.MSUB(LIN)) THEN IRET = -1 GOTO 100 ENDIF C CHECK COMPRESSION STATUS OF INPUT MESSAGE, OUTPUT MESSAGE WILL MATCH C -------------------------------------------------------------------- CALL MESGBC(-LUNIN,MEST,ICMP) IF(ICMP.EQ.1) THEN C ------------------------------------------------------- C THIS BRANCH IS FOR COMPRESSED INPUT/OUTPUT MESSAGES C ------------------------------------------------------- C READ IN AND UNCOMPRESS SUBSET, THEN COPY IT TO COMPRESSED OUTPUT MSG C -------------------------------------------------------------------- CALL READSB(LUNIN,IRET) IF(LUNOT.GT.0) THEN CALL UFBCPY(LUNIN,LUNOT) CALL CMPMSG('Y') CALL WRITSB(LUNOT) CALL CMPMSG('N') ENDIF GOTO 100 ELSE IF(ICMP.EQ.0) THEN C ------------------------------------------------------- C THIS BRANCH IS FOR UNCOMPRESSED INPUT/OUTPUT MESSAGES C ------------------------------------------------------- C COPY THE SUBSET TO THE OUTPUT MESSAGE AND/OR RESET THE POINTERS C --------------------------------------------------------------- IBIT = (MBYT(LIN))*8 CALL UPB(NBYT,16,MBAY(1,LIN),IBIT) IF(LUNOT.GT.0) CALL CPYUPD(LUNOT,LIN,LOT,NBYT) MBYT(LIN) = MBYT(LIN) + NBYT NSUB(LIN) = NSUB(LIN) + 1 ELSE GOTO 907 ENDIF C EXITS C ----- 100 RETURN 900 CALL BORT('BUFRLIB: COPYSB - INPUT BUFR FILE IS CLOSED, IT '// . 'MUST BE OPEN FOR INPUT') 901 CALL BORT('BUFRLIB: COPYSB - INPUT BUFR FILE IS OPEN FOR '// . 'OUTPUT, IT MUST BE OPEN FOR INPUT') 902 CALL BORT('BUFRLIB: COPYSB - A MESSAGE MUST BE OPEN IN INPUT '// . 'BUFR FILE, NONE ARE') 903 CALL BORT('BUFRLIB: COPYSB - OUTPUT BUFR FILE IS CLOSED, IT '// . 'MUST BE OPEN FOR OUTPUT') 904 CALL BORT('BUFRLIB: COPYSB - OUTPUT BUFR FILE IS OPEN FOR '// . 'INPUT, IT MUST BE OPEN FOR OUTPUT') 905 CALL BORT('BUFRLIB: COPYSB - A MESSAGE MUST BE OPEN IN OUTPUT '// . 'BUFR FILE, NONE ARE') 906 CALL BORT('BUFRLIB: COPYSB - INPUT AND OUTPUT BUFR FILES MUST '// . 'HAVE THE SAME INTERNAL TABLES, THEY ARE DIFFERENT HERE') 907 WRITE(BORT_STR,'("BUFRLIB: COPYSB - INVALID COMPRESSION '// . 'INDICATOR (ICMP=",I3," RETURNED FROM BUFR ARCHIVE LIBRARY '// . 'ROUTINE MESGBC")') ICMP CALL BORT(BORT_STR) END