! module rdstdin ! ! Uars solstice flux data: integer :: i,j,ii,l,ll,nbins integer,parameter :: nflxbins=170 real,dimension(nflxbins) :: | wv_low, wv_high, sminflx, smedflx, smaxflx ! ! From sftrf.h: integer,parameter :: nsftrf=158 real,dimension(nsftrf) :: | vl1 ,vl2 ,sfx ,sgry ,sgo2hz ,sgo31 ,sgo32 integer,parameter :: nsfdim1=17, nsfdima=9, nsfdimb=5 real,dimension(nsfdim1,nsfdima) :: aij real,dimension(nsfdim1,nsfdimb) :: bij real,dimension(nsfdima) :: ad00, ad10 real,dimension(nsfdimb) :: bd00, bd10 ! ! From atmosp.h: integer, parameter :: natm=71,nfatm=8 real :: atm(natm,nfatm) ! contains !----------------------------------------------------------------------- subroutine readstdin ! ! Read uars solar fluxes (see readsf.F and uarsflx.h): read(5,"(85x,i4)") nbins ! (smin,smed,smax uars days are in the 85x) if (nbins.ne.nflxbins) then write(6,"('>>> mkinputdat: read nbins=',i4,' but nflxbins=', | i4)") nbins,nflxbins stop 'nbins' endif do i=1,nflxbins read(5,"(5x,2f12.2,3e15.4)") wv_low(i),wv_high(i), | sminflx(i),smedflx(i),smaxflx(i) ! write(6,"('bin ',i3,': wv_low,high=',f8.2,',',f8.2, ! | ' sminflx=',1pe12.4,' smedflx=',1pe12.4,' smaxflx=', ! | 1pe12.4)") i,wv_low(i),wv_high(i),sminflx(i),smedflx(i), ! | smaxflx(i) enddo ! ! Read vl1,vl2, etc (see sftrf.h) do l=1,32 read 100,ll,vl1(l),vl2(l),sfx(l),sgry(l),sgo2hz(l),sgo31(l), | sgo32(l) enddo do l=33,60 read 101,ll,vl1(l),vl2(l),sfx(l),sgry(l),sgo31(l),sgo32(l) enddo do l=61,158 read 102,ll,vl1(l),vl2(l),sfx(l),sgry(l),sgo31(l) enddo DO I=1,17 READ 103,II,(AIJ(I,J),J=1,4) enddo DO I=1,17 READ 103,II,(AIJ(I,J),J=5,8) enddo DO I=1,17 READ 103,II,AIJ(I,9) enddo DO I=1,17 READ 103,II,(BIJ(I,J),J=1,4) enddo DO I=1,17 READ 103,II,BIJ(I,5) enddo READ 103,II,(AD00(J),J=1,4) READ 103,II,(AD00(J),J=5,8) READ 103,II,AD00(9) READ 103,II,(AD10(J),J=1,4) READ 103,II,(AD10(J),J=5,8) READ 103,II,AD10(9) READ 103,II,(BD00(J),J=1,4) READ 103,II,BD00(5) READ 103,II,(BD10(J),J=1,4) READ 103,II,BD10(5) 100 FORMAT(I3,F9.3,F9.3,E10.2,E10.2,E10.2,E10.2,E10.2) 101 FORMAT(I3,F9.3,F9.3,E10.2,E10.2,E10.2,E10.2) 102 FORMAT(I3,F7.1,F7.1,E10.2,E10.2,E10.2) 103 FORMAT(I3,E15.6,E15.6,E15.6,E15.6) 104 FORMAT(1X,I3,F9.3,F9.3,E10.2,E10.2,E10.2,E10.2,E10.2) 105 FORMAT(1X,I3,F9.3,F9.3,E10.2,E10.2,E10.2,E10.2) 106 FORMAT(1X,I3,F7.1,F7.1,E10.2,E10.2,E10.2) 107 FORMAT(1X,I3,E15.6,E15.6,E15.6,E15.6) ! real,dimension(158) :: ! | vl1 ,vl2 ,sfx ,sgry ,sgo2hz ,sgo31 ,sgo32 ! ! write(6,"('vl1=',/,(6e12.4))") vl1 ! write(6,"('vl2=',/,(6e12.4))") vl2 ! write(6,"('sfx=',/,(6e12.4))") sfx ! write(6,"('sgry=',/,(6e12.4))") sgry ! write(6,"('sgo2hz=',/,(6e12.4))") sgo2hz ! write(6,"('sgo31=',/,(6e12.4))") sgo31 ! write(6,"('sgo32=',/,(6e12.4))") sgo32 ! real,dimension(17,9) :: aij ! real,dimension(17,5) :: bij ! real,dimension(9) :: ad00, ad10 ! real,dimension(5) :: bd00, bd10 ! do i=1,9 ! write(6,"('i=',i3,' aij(:,i)=',/,(6e12.4))") ! | i,aij(:,i) ! enddo ! do i=1,5 ! write(6,"('i=',i3,' bij(:,i)=',/,(6e12.4))") ! | i,bij(:,i) ! enddo ! write(6,"('ad00=',/,(6e12.4))") ad00 ! write(6,"('ad10=',/,(6e12.4))") ad10 ! write(6,"('bd00=',/,(6e12.4))") bd00 ! write(6,"('bd10=',/,(6e12.4))") bd10 ! ! Read Fomichev data (as in glbmean.F): ! integer, parameter :: natm=71,nfatm=8 ! real :: atm(natm,nfatm) ! do i=1,natm read(5,*) (atm(i,ii),ii=1,nfatm) enddo do i=1,nfatm write(6,"('readstdin: i=',i3,' atm(:,i)=',/,(6e12.4))") | i,atm(:,i) enddo end subroutine readstdin end module rdstdin