c----------------------------------------------------------------------- c c ----------- c ssmitdrtb 2.0 c ----------- c c Copy a Ta file produced by ssmitdrta, converting the antenna c temperatures to brightness temperatures. c c----------------------------------------------------------------------- c c Written 11/94 by Don Moss (University of Alabama in Huntsville) c c Rev 1 2/16/95 Don Moss DOR95-065 Release Level 1.1 c Use single quotes instead of double quotes in output strings. c c Rev 2 5/28/96 Don Moss DOR96-101 Release Level 2.0 c Use new Ta-Tb conversion subroutine. c c----------------------------------------------------------------------- c c Standard Input (keyboard): c Input file name c Output file name c c File Input: c SSM/I TDR antenna temperature file produced by ssmitdrta. c For file formats, see the "File Output" description in the c companion program, ssmitdrta.c. c c File Output: c Same as input, except with Ta's converted to Tb's. c c----------------------------------------------------------------------- c c Note: Variables used by the Ta-Tb conversion subroutines are c spelled with all capital letters throughout. c c----------------------------------------------------------------------- c c *-*-*-*-*-*-*-*-*-*-*-* c c I M P O R T A N T c c *-*-*-*-*-*-*-*-*-*-*-* c c NOTE: When compiling this program, use the "static" compiler c option if available. It must be compiled in such a way that c all variables are treated as "static" rather than "automatic". c This may be the default mode with your compiler - consult your c manual to be sure. Static means that subprogram variables c retain their values between subprogram calls. Automatic means c that subprogram variables are created and destroyed each time c the subprogram is called. c c----------------------------------------------------------------------- character*48 infile, outfile ! file names integer ids/11/ ! input device number integer ods/21/ ! output device number integer irev ! revolution # integer nscan/0/ ! scan counter c---variables for Wentz Tb conversion subroutines--- REAL*8 REV,XTIME INTEGER*4 ITIME,ITIMSC,IVOLT,IAGC,ICOLDA,IHOTA,ICOLDB,IHOTB,IASCTM INTEGER*4 IATOIL,IBTOIL,ISAT REAL*4 XLATSC,XLONSC,ALTSC,THT,HLTEMP,RFTEMP,FRTEMP REAL*4 PERIOD,ASCLOC,ANGINC,AXIS,ECC,ANGPER,SPACER REAL*4 ALAT,ALON,BLAT,BLON,TALO,ATAHI,BTAHI COMMON /OUTDAT/ REV,XTIME,ITIME,ITIMSC,XLATSC,XLONSC,ALTSC,THT, 1 HLTEMP(3),IVOLT(2),RFTEMP,FRTEMP,IAGC(6), 2 IASCTM,PERIOD,ASCLOC,ANGINC,AXIS,ECC,ANGPER,ISAT,SPACER(6), 2 ICOLDA(5,7),IHOTA(5,7),ICOLDB(5,2),IHOTB(5,2), 3 ALAT(128),ALON(128),BLAT(128),BLON(128), 4 TALO(5,64),ATAHI(2,128),BTAHI(2,128),IATOIL(128),IBTOIL(128) COMMON /TBDATA/ TBLO(5,64),ATBHI(2,128),BTBHI(2,128) c---show program name and version--- print *, '===== ssmitdrtb 2.0 =====' c---get input file name--- print *, 'Enter input file name' read (*,555) infile 555 format (a48) c---get output file name--- print *, 'Enter output file name' read (*,555) outfile c---show intention--- print *, 'Ta''s will be converted to Tb''s' c---open input file--- print *, 'Opening input file: ', infile open (ids,file=infile) c---open output file--- print *, 'Opening output file: ', outfile open (ods,file=outfile) c====================main loop==================== c---read a scan (16 records)--- 2 read (ids,101,err=7,end=8) ISAT, XTIME, irev read (ids,102) ALAT read (ids,102) BLAT read (ids,102) ALON read (ids,102) BLON read (ids,103) IATOIL read (ids,103) IBTOIL read (ids,104) ((TALO(i,j),j=1,64),i=1,5) read (ids,105) (ATAHI(1,j),j=1,128) read (ids,105) (BTAHI(1,j),j=1,128) read (ids,105) (ATAHI(2,j),j=1,128) read (ids,105) (BTAHI(2,j),j=1,128) 101 format (i3, f16.4, i9) 102 format (128f7.2) 103 format (128i2) 104 format (64f7.2) 105 format (128f7.2) c---convert Ta's to Tb's--- call FDTATB(1) c---write the scan (16 records)--- write (ods,101) ISAT, XTIME, irev write (ods,102) ALAT write (ods,102) BLAT write (ods,102) ALON write (ods,102) BLON write (ods,103) IATOIL write (ods,103) IBTOIL write (ods,104) ((TBLO(i,j),j=1,64),i=1,5) write (ods,105) (ATBHI(1,j),j=1,128) write (ods,105) (BTBHI(1,j),j=1,128) write (ods,105) (ATBHI(2,j),j=1,128) write (ods,105) (BTBHI(2,j),j=1,128) nscan = nscan + 1 go to 2 c================================================= 7 print *, '* read error *' stop 8 print *, '- eof -' print *, nscan, ' scan(s) converted' stop end c======================================================================= c c The following subroutine, FDTATB, is part of the NASA "Third c Precipitation Intercomparison Project" (PIP-3) package. FDTATB c was supplied by Frank J. Wentz, Remote Sensing Systems, Santa c Rosa, CA, and is included here by permission. c c======================================================================= SUBROUTINE FDTATB(I85GHZ) C C THIS SUBROUTINE ADJUST THE ANTENNA TEMPS AND FIND BRIGHTNESS TEMPS C REAL*4 DELTA(4),CHI(2,4),AVGTA(5),XTA(5) REAL*4 GCOEF10(5),HCOEF10(5),GCOEF11(5),HCOEF11(5) REAL*4 SBIAS08(64,5),SBIAS10(64,5),SBIAS11(64,5) REAL*4 AVV(4),AHV(4),AOV(4),AHH(4),AVH(4),AOH(4) C C SPECIFY COMMON /OUTDAT/ C REAL*8 REV,XTIME INTEGER*4 ITIME,ITIMSC,IVOLT,IAGC,ICOLDA,IHOTA,ICOLDB,IHOTB,IASCTM INTEGER*4 IATOIL,IBTOIL,ISAT REAL*4 XLATSC,XLONSC,ALTSC,THT,HLTEMP,RFTEMP,FRTEMP REAL*4 PERIOD,ASCLOC,ANGINC,AXIS,ECC,ANGPER,SPACER REAL*4 ALAT,ALON,BLAT,BLON,TALO,ATAHI,BTAHI COMMON/OUTDAT/ REV,XTIME,ITIME,ITIMSC,XLATSC,XLONSC,ALTSC,THT, 1 HLTEMP(3),IVOLT(2),RFTEMP,FRTEMP,IAGC(6), 2 IASCTM,PERIOD,ASCLOC,ANGINC,AXIS,ECC,ANGPER,ISAT,SPACER(6), 2 ICOLDA(5,7),IHOTA(5,7),ICOLDB(5,2),IHOTB(5,2), 3 ALAT(128),ALON(128),BLAT(128),BLON(128), 4 TALO(5,64),ATAHI(2,128),BTAHI(2,128),IATOIL(128),IBTOIL(128) C COMMON /TBDATA/ TBLO(5,64),ATBHI(2,128),BTBHI(2,128) C C DATA INITIALIZATION C DATA ISTART/1/ DATA DELTA/0.03199,0.02685,0.01434,0.01186/ DATA CHI/.00379,.00525,.00983,0.0,.02136,.02664,.01387,.01967/ DATA AVGTA/190.93,130.14,215.42,211.39,158.16/ C C GCOEF10,HCOEF10 ARE F10-F08 INTERCALIBRATION COEFS. C GCOEF11,HCOEF11 ARE F11-F08 INTERCALIBRATION COEFS. C DATA GCOEF10/0.08, 0.35,-0.33,-0.01, 0.44/ DATA HCOEF10/.221132E-02,.786968E-03,.161037E-02,.335131E-02, 1 .165331E-02/ C DATA GCOEF11/0.44, -0.16, 0.30, -0.01, -0.03/ DATA HCOEF11/.221132E-02,.786968E-03,.161037E-02,.335131E-02, 1 .165331E-02/ C C F08 X-SCAN BIASES C DATA(SBIAS08(I,1),I=1,64)/ -.03, .01, .04, .06, + .09, .10, .12, .12, .13, .13, .13, .13, .13, .14, + .14, .14, .14, .14, .15, .15, .14, .14, .13, .12, + .12, .11, .11, .11, .11, .11, .12, .11, .10, .10, + .09, .07, .06, .05, .05, .04, .03, .03, .02, .01, + -.01, -.02, -.04, -.05, -.06, -.08, -.09, -.09, -.11, -.12, + -.14, -.18, -.21, -.26, -.30, -.36, -.43, -.52, -.63, -.75/ C DATA(SBIAS08(I,2),I=1,64)/ -.03, -.01, .00, .01, + .01, .01, .01, .01, .01, .01, .01, .02, .03, .04, + .04, .05, .05, .05, .05, .05, .04, .04, .04, .04, + .04, .04, .04, .04, .04, .05, .06, .07, .07, .08, + .07, .07, .06, .07, .07, .06, .06, .06, .06, .06, + .05, .04, .04, .04, .04, .02, .00, -.01, -.02, -.03, + -.04, -.06, -.09, -.12, -.16, -.20, -.26, -.32, -.39, -.48/ C DATA(SBIAS08(I,3),I=1,64)/ -.05, .00, .04, .07, + .08, .09, .10, .11, .11, .12, .12, .13, .14, .14, + .15, .15, .17, .17, .16, .16, .15, .14, .14, .14, + .13, .13, .13, .13, .13, .13, .14, .13, .14, .13, + .12, .11, .10, .09, .09, .08, .07, .05, .04, .03, + .02, .00, -.01, -.02, -.04, -.05, -.07, -.10, -.12, -.14, + -.17, -.20, -.25, -.31, -.36, -.42, -.50, -.63, -.79, -.97/ C DATA(SBIAS08(I,4),I=1,64)/ -.06, -.03, .00, .01, + .03, .04, .05, .06, .06, .07, .09, .10, .11, .13, + .14, .15, .15, .16, .16, .16, .16, .16, .16, .15, + .14, .14, .14, .13, .13, .13, .14, .13, .13, .13, + .12, .11, .10, .10, .09, .08, .08, .06, .06, .05, + .04, .03, .02, .01, .00, -.02, -.05, -.07, -.09, -.12, + -.16, -.18, -.22, -.27, -.33, -.42, -.50, -.60, -.74, -.90/ C DATA(SBIAS08(I,5),I=1,64)/ -.23, -.20, -.17, -.15, + -.14, -.13, -.12, -.10, -.09, -.10, -.09, -.07, -.05, -.04, + -.03, .00, .01, .03, .04, .04, .04, .05, .05, .05, + .04, .04, .04, .05, .05, .07, .08, .09, .10, .11, + .10, .10, .09, .10, .10, .09, .09, .09, .09, .08, + .07, .07, .07, .09, .09, .08, .08, .08, .08, .08, + .06, .04, .00, -.03, -.07, -.11, -.16, -.22, -.29, -.41/ C C F10 XSCAN BIASES C DATA(SBIAS10(I,1),I=1,64)/ -.01, .03, .05, .07, + .09, .10, .11, .12, .12, .14, .14, .14, .14, .15, + .16, .15, .15, .14, .14, .14, .14, .14, .13, .13, + .14, .14, .14, .14, .14, .14, .13, .13, .12, .12, + .11, .10, .08, .07, .07, .07, .06, .05, .04, .02, + .00, -.01, -.03, -.04, -.05, -.06, -.08, -.09, -.11, -.14, + -.15, -.18, -.21, -.26, -.31, -.37, -.46, -.57, -.69, -.82/ C DATA(SBIAS10(I,2),I=1,64)/ .00, .02, .02, .02, + .01, .01, .01, .01, .01, .00, .00, .00, .00, .00, + .00, .00, .00, -.01, -.02, -.02, -.02, -.02, -.02, -.01, + .00, .01, .02, .03, .04, .05, .06, .06, .06, .06, + .05, .05, .05, .06, .07, .08, .08, .08, .07, .08, + .08, .08, .08, .08, .08, .07, .06, .04, .03, .02, + .01, -.01, -.04, -.07, -.11, -.15, -.21, -.27, -.36, -.45/ C DATA(SBIAS10(I,3),I=1,64)/ -.04, .01, .04, .06, + .08, .09, .11, .11, .12, .13, .14, .15, .15, .15, + .15, .15, .15, .14, .13, .14, .14, .14, .14, .14, + .14, .15, .15, .16, .16, .18, .18, .16, .15, .14, + .14, .14, .14, .13, .11, .10, .09, .08, .07, .05, + .03, .01, -.01, -.03, -.05, -.07, -.09, -.12, -.14, -.15, + -.17, -.20, -.24, -.28, -.34, -.42, -.52, -.64, -.77, -.95/ C DATA(SBIAS10(I,4),I=1,64)/ -.09, -.06, -.03, -.02, + -.01, .01, .02, .03, .05, .06, .07, .08, .09, .10, + .11, .12, .12, .12, .13, .13, .13, .14, .14, .14, + .15, .16, .16, .16, .16, .17, .16, .16, .15, .15, + .14, .13, .12, .12, .12, .11, .11, .09, .08, .07, + .06, .05, .04, .03, .01, -.01, -.04, -.05, -.07, -.09, + -.12, -.16, -.21, -.25, -.31, -.38, -.46, -.56, -.69, -.84/ C DATA(SBIAS10(I,5),I=1,64)/ -.20, -.18, -.16, -.15, + -.15, -.15, -.14, -.13, -.11, -.11, -.10, -.09, -.08, -.07, + -.07, -.05, -.04, -.05, -.05, -.04, -.04, -.03, -.02, -.01, + .01, .02, .04, .05, .06, .08, .09, .09, .10, .10, + .09, .09, .10, .11, .12, .12, .13, .13, .13, .13, + .12, .13, .12, .13, .13, .13, .12, .12, .12, .11, + .10, .08, .05, .01, -.02, -.05, -.11, -.18, -.27, -.39/ C C F11 XSCAN BIASES C DATA(SBIAS11(I,1),I=1,64)/ -.11, -.08, -.07, -.05, + -.02, .00, -.01, -.02, -.02, .00, .01, .00, .00, .00, + .02, .03, .02, .02, .02, .04, .04, .04, .03, .03, + .05, .06, .06, .05, .06, .08, .07, .07, .07, .07, + .09, .09, .08, .08, .08, .09, .10, .09, .08, .08, + .10, .10, .09, .08, .08, .10, .10, .09, .06, .04, + .03, .00, -.03, -.08, -.14, -.20, -.32, -.43, -.58, -.75/ C DATA(SBIAS11(I,2),I=1,64)/ -.02, .00, .00, .01, + .01, .01, .02, .01, .02, .02, .02, .03, .04, .05, + .05, .06, .07, .07, .07, .06, .06, .07, .07, .07, + .07, .07, .07, .07, .08, .07, .06, .06, .07, .07, + .06, .06, .06, .07, .08, .08, .08, .08, .08, .09, + .09, .10, .09, .08, .06, .04, .03, .01, -.01, -.03, + -.05, -.09, -.12, -.16, -.21, -.29, -.37, -.46, -.56, -.67/ C DATA(SBIAS11(I,3),I=1,64)/ -.13, -.10, -.06, -.04, + -.02, -.01, .00, .01, .02, .03, .03, .03, .04, .05, + .06, .06, .07, .07, .07, .08, .07, .08, .08, .09, + .09, .10, .10, .11, .12, .13, .13, .13, .13, .13, + .13, .13, .13, .13, .13, .12, .11, .11, .10, .10, + .09, .08, .08, .08, .07, .07, .06, .05, .03, .01, + -.02, -.07, -.13, -.19, -.25, -.34, -.47, -.64, -.81, -.99/ C DATA(SBIAS11(I,4),I=1,64)/ -.16, -.09, -.06, -.03, + -.05, -.04, -.02, -.01, .01, -.01, .00, .02, .04, .06, + .04, .05, .07, .09, .10, .07, .07, .09, .11, .11, + .09, .09, .11, .13, .14, .11, .09, .11, .12, .13, + .10, .10, .11, .12, .13, .10, .10, .10, .11, .12, + .10, .10, .11, .11, .10, .07, .07, .06, .04, .02, + -.03, -.07, -.10, -.15, -.23, -.36, -.50, -.60, -.73, -.89/ C DATA(SBIAS11(I,5),I=1,64)/ -.17, -.15, -.14, -.12, + -.10, -.10, -.09, -.09, -.09, -.07, -.06, -.05, -.04, -.03, + -.01, .01, .01, .02, .03, .03, .03, .04, .04, .04, + .05, .06, .05, .06, .06, .06, .05, .05, .06, .06, + .06, .07, .08, .08, .09, .09, .10, .09, .08, .09, + .11, .12, .12, .12, .11, .11, .12, .12, .10, .08, + .06, .04, .01, -.03, -.08, -.14, -.22, -.31, -.43, -.58/ C C BEGIN EXECUTION C IF(ISTART.EQ.1) THEN ISTART=0 DO 10 IFREQ=1,4 IF(IFREQ.EQ.2) GO TO 10 XFAC=(1.-CHI(1,IFREQ)*CHI(2,IFREQ))*(1.-DELTA(IFREQ)) AVV(IFREQ)=(1.+CHI(1,IFREQ))/XFAC AHV(IFREQ)=-CHI(1,IFREQ)*(1.+CHI(2,IFREQ))/XFAC AOV(IFREQ)=(1.-AVV(IFREQ)-AHV(IFREQ))*2.7 AHH(IFREQ)=(1.+CHI(2,IFREQ))/XFAC AVH(IFREQ)=-CHI(2,IFREQ)*(1.+CHI(1,IFREQ))/XFAC AOH(IFREQ)=(1.-AHH(IFREQ)-AVH(IFREQ))*2.7 10 CONTINUE C DO 20 ICH=1,5 DO 20 ICEL=1,64 SBIAS08(ICEL,ICH)= 1-SBIAS08(ICEL,ICH)/AVGTA(ICH) SBIAS10(ICEL,ICH)=(1-SBIAS10(ICEL,ICH)/AVGTA(ICH))* 1 (1-HCOEF10(ICH)) SBIAS11(ICEL,ICH)=(1-SBIAS11(ICEL,ICH)/AVGTA(ICH))* 1 (1-HCOEF11(ICH)) 20 CONTINUE ENDIF C DO 100 ICEL=1,64 C IBAD1=0 IBAD2=0 IBAD3=0 IF(TALO(1,ICEL).LT.55 .OR. TALO(1,ICEL).GT.320) IBAD1=1 IF(TALO(2,ICEL).LT.55 .OR. TALO(2,ICEL).GT.320) IBAD1=1 IF(TALO(3,ICEL).LT.55 .OR. TALO(3,ICEL).GT.320) IBAD2=1 IF(TALO(4,ICEL).LT.55 .OR. TALO(4,ICEL).GT.320) IBAD3=1 IF(TALO(5,ICEL).LT.55 .OR. TALO(5,ICEL).GT.320) IBAD3=1 C DO 30 ICH=1,5 XTA(ICH)= TALO(ICH,ICEL)*SBIAS08(ICEL,ICH) IF(ISAT.EQ.10) 1 XTA(ICH)=TALO(ICH,ICEL)*SBIAS10(ICEL,ICH)-GCOEF10(ICH) IF(ISAT.EQ.11) 1 XTA(ICH)=TALO(ICH,ICEL)*SBIAS11(ICEL,ICH)-GCOEF11(ICH) 30 CONTINUE C IF(IBAD1.EQ.0) THEN TALO(1,ICEL)=XTA(1) TALO(2,ICEL)=XTA(2) TBLO(1,ICEL)=AVV(1)*XTA(1) + AHV(1)*XTA(2) + AOV(1) TBLO(2,ICEL)=AHH(1)*XTA(2) + AVH(1)*XTA(1) + AOH(1) ELSE TBLO(1,ICEL)=TALO(1,ICEL) TBLO(2,ICEL)=TALO(2,ICEL) ENDIF C IF(IBAD2.EQ.0) THEN TALO(3,ICEL)=XTA(3) TBLO(3,ICEL)=1.01993*XTA(3) + 1.994 ELSE TBLO(3,ICEL)=TALO(3,ICEL) ENDIF C IF(IBAD3.EQ.0) THEN TALO(4,ICEL)=XTA(4) TALO(5,ICEL)=XTA(5) TBLO(4,ICEL)=AVV(3)*XTA(4) + AHV(3)*XTA(5) + AOV(3) TBLO(5,ICEL)=AHH(3)*XTA(5) + AVH(3)*XTA(4) + AOH(3) ELSE TBLO(4,ICEL)=TALO(4,ICEL) TBLO(5,ICEL)=TALO(5,ICEL) ENDIF C 100 CONTINUE C IF(I85GHZ.EQ.0) RETURN C DO 200 ICEL=1,128 IF(ATAHI(1,ICEL).GE.55 .AND. ATAHI(1,ICEL).LE.320 .AND. 1 ATAHI(2,ICEL).GE.55 .AND. ATAHI(2,ICEL).LE.320) THEN ATBHI(1,ICEL)=AVV(4)*ATAHI(1,ICEL)+AHV(4)*ATAHI(2,ICEL)+AOV(4) ATBHI(2,ICEL)=AHH(4)*ATAHI(2,ICEL)+AVH(4)*ATAHI(1,ICEL)+AOH(4) ELSE ATBHI(1,ICEL)=ATAHI(1,ICEL) ATBHI(2,ICEL)=ATAHI(2,ICEL) ENDIF C IF(BTAHI(1,ICEL).GE.55 .AND. BTAHI(1,ICEL).LE.320 .AND. 1 BTAHI(2,ICEL).GE.55 .AND. BTAHI(2,ICEL).LE.320) THEN BTBHI(1,ICEL)=AVV(4)*BTAHI(1,ICEL)+AHV(4)*BTAHI(2,ICEL)+AOV(4) BTBHI(2,ICEL)=AHH(4)*BTAHI(2,ICEL)+AVH(4)*BTAHI(1,ICEL)+AOH(4) ELSE BTBHI(1,ICEL)=BTAHI(1,ICEL) BTBHI(2,ICEL)=BTAHI(2,ICEL) ENDIF 200 CONTINUE RETURN END