sbuv2%code_ibsl19$ diff anclv6_x.f ../code_ibsl1/anclv6_x.f 34,45d33 < c-h < c-h Add a 2nd set of tables for N19 to switch Channel-12 wavelength from 340 nm to 343 nm. < c-h This will allow the code to process the two different wavelength settings in the same < c-h day cotinuously. To run this switching, change the nml file with the following < c-h IOPTS(18)=T < c-h xancilfile1='ancln19n.seq', < c-h ancilfilewv0='auxn19n.dat', < c-h xancilfile2='' < c-h The switching subroutine is in chnglam.f, which is called by lamchk. The Day-1 calcons < c-h subroutines, which are also modified to load parameters for proper wavelength, are called < c-h in lamchk. < c-h ----- April 2013 54d41 < include "multx.inc" 71c58 < INTEGER*4 DATA(179,40),ianclx1(3),ianclx2(3),i2xancf --- > INTEGER*4 DATA(179,40),ianclx1(3),ianclx2(3) 117,125d103 < c-h If given only ONE extra ancilfile, then, it is not for the N16 TOMS wavelength scan mode. < c-h This is used to get the extra ancilfile for the N19 SBUV/2 wavelength change (340nm -> 343 nm) < i2xancf=index(xancilfile2,' ') < if(index(xancilfile1,' ').ne.0 .and. index(ancilfilewv0,' ').ne.0 < 1 .and. index(xancilfile2,' ').eq. 0)then < lamswtch=.TRUE. < else < lamswtch=.FALSE. < endif 149,161d126 < c-h < if(lamswtch) then < open(18,file=ancilfilewv0) < do i=1,12 < read(18,1400) x1,x2,x3,x4,x5,x6 < alfac1(i)=x2 < calfa11(i)=x3 < calfa21(i)=x4 < betac1(i)=x5 < oobrcr1(i)=x6 < enddo < close(18) < endif 194,226d158 < c-h < if(lamswtch) then < OPEN (UNIT=LUXANCL2,FILE=xancilfile1,FORM='UNFORMATTED' < + ,STATUS='OLD') < c-h To get all infor for N19 changing wavelength < 205 continue < call rdswap(LUXANCL2,B_REC1B,id_rec,*900) < IF (REC1B(1).NE.465) GO TO 205 < CALL FNMOVE (WLEN0x(1),13*4,REC1B(2)) < CALL FNMOVE (ALFA0x(1),13*4,REC1B(15)) < CALL FNMOVE (BETAx(1),13*4,REC1B(28)) < K = 1 < 210 continue < call rdswap(LUXANCL2,B_REC1B,id_rec,*900) < IF (REC1B(1).NE.475) GO TO 900 < CALL FNMOVE(DATA(1,K),179*4,REC1B(2)) < K = K + 1 < IF (K.LT.41) GO TO 210 < CALL FNMOVE (QLOGx(1,1,1,1),7130*4,DATA(1,1)) < K = 1 < 231 continue < call rdswap(LUXANCL2,B_REC1B,id_rec,*900) < IF (REC1B(1).NE.485 .and. K .eq. 1) GO TO 231 < CALL FNMOVE (DATA(1,K),179*4,REC1B(2)) < K = K + 1 < IF (K.LT.34) GO TO 231 < close(LUXANCL2) < CALL FNMOVE (XLOGI0X(1,1,1,1),5796*4,DATA(1,1)) < write(6,6665) xancilfile1 < write(6,6666) wlen0x,alfa0x,betax < 6665 format(1x,'Loaded ',a120) < endif < sbuv2%code_ibsl19$ diff lamchk.f ../code_ibsl1/lamchk.f 114,122d113 < c-h N19 Channel-12 wavelength change < if(SCNO.eq.N19 .and. i .eq. 12)then < if(((abs(GRTWAV(I)-3430.22) .lt. 5) .and. < 1 (abs(WLEN0(I)-3398.77) .lt. 5)) .or. < 1 ((abs(GRTWAV(I)-3398.77) .lt. 5) .and. < 1 (abs(WLEN0(I)-3430.22) .lt. 5)))then < call chnglam < endif < endif sbuv2%code_ibsl19$ diff prelch_62x.f ../code_ibsl1/prelch_62x.f 18c18 < real*4 rirr12(12),rirr23(12) !,acf(13) --- > real*4 rirr12(12),rirr23(12),acf(13) sbuv2%code_ibsl19$ diff params.inc ../code_ibsl1/params.inc 18,21c18,21 < 3 alfac(13),calfa1(13),calfa2(13),betac(13),oobrcr(12), < 4 alfac0(13),calfa10(13),calfa20(13),betac0(13), < 5 alfac1(13),calfa11(13),calfa21(13),betac1(13),oobrcr1(12), < 6 wlen0x(13),alfa0x(13),betax(13),RSQR --- > 3 alfac(13),calfa1(13),calfa2(13), > 4 betac(13),oobrcr(12), > 5 alfac0(13),calfa10(13),calfa20(13), > 6 betac0(13) sbuv2%code_ibsl19$ diff day1cals.inc ../code_ibsl1/day1cals.inc 6c6 < + ebert,ccr_os,acf --- > + ebert,ccr_os 8c8 < + rirr12_value,rirr23_day1(12),ebert(3),ccr_os(12),acf(13) --- > + rirr12_value,rirr23_day1(12),ebert(3),ccr_os(12)