OB.DAAC Logo
NASA Logo
Ocean Color Science Software

ocssw V2022
outcrfty.f
Go to the documentation of this file.
1  subroutine outcrfty(tmcfd,tmcfu,ocrftzu,ocrftzd)
2 
3 c***********************************************************************
4 c subroutine outcrfty creates output datasets for upwelling and
5 c downwelling diffused radiation at the aircraft height when the
6 c lower surface is rough ocean
7 c***********************************************************************
8 c.....include the common blocks
9  implicit real*8 (a-h,o-z)
10  include 'afrt_rt2.cmn'
11  real*8 tma(nsz),tmb(nsz),tmc(nsz),tmfd(nsz),tmfu(nsz),tms(nsz)
12  real*8 tmcfd(nsz),tmcfu(nsz),ocrftzu(nph,nth,nsz),ocrftzd(nph,nth,nsz)
13 c***********************************************************************
14 c
15  lsza=bfr1(35)+0.001
16  do i=1,lsza,1
17  m=msza(i)
18  tma(i)=the0in(m)
19  tmb(i)=fdirc(m)
20  tmc(i)=sbarz(m)
21  tmfd(i)=fdown(m)
22  tmfu(i)=fup(m)
23  tms(i)=oalb(m)
24  tmcfd(i)=cfdown(m)
25  tmcfu(i)=cfup(m)
26  enddo
27 c
28  iprn=14
29  call headr(iprn)
30 c
31  write(iprn,6659)
32  do 5567 i=1,lsza
33  it=msza(i)
34  write(iprn,6660)tma(i),pi
35  write(iprn,6661)tmb(i),tmfd(i)
36  write(iprn,6662)tmfu(i),tms(i)
37  write(iprn,6674)tmcfd(i),tmcfu(i)
38  write(iprn,6668)
39  write(iprn,6664)(the(ir),ir=1,(nx-1))
40  do is=1,jpart
41  write(iprn,6667)jphi(is),(crftzu(it,ir,is)*pi,ir=1,(nx-1))
42  do ir=1,(nx-1)
43  ocrftzu(is,ir,it) = crftzu(it,ir,is)*pi
44  enddo
45  enddo
46 5567 continue
47 c
48  iprn=13
49  call headr(iprn)
50 c
51  write(iprn,6659)
52  do 5568 i=1,lsza
53  it=msza(i)
54  write(iprn,6660)tma(i),pi
55  write(iprn,6661)tmb(i),tmfd(i)
56  write(iprn,6662)tmfu(i),tms(i)
57  write(iprn,6674)tmcfd(i),tmcfu(i)
58  write(iprn,6668)
59  write(iprn,6664)(the(ir),ir=1,(nx-1))
60  do is=1,jpart
61  write(iprn,6667)jphi(is),(crftzd(it,ir,is)*pi,ir=1,(nx-1))
62  do ir=1,(nx-1)
63  ocrftzd(it,ir,is) = crftzd(it,ir,is)*pi
64  enddo
65  enddo
66 5568 continue
67 c*****format statements*************************************************
68 c
69 6659 format(t1,'fluxes')
70 6660 format(t8,'sza',t32,f6.1,t40,'f0_top',t64,1pe12.4)
71 6661 format(t8,'fdir_btm',t26,1pe12.4,t40,'fdif_btm',t64,1pe12.4)
72 6662 format(t8,'fdif_up',t26,1pe12.4,t40,'hem. ref(ocean)',t64,1pe12.4)
73 6664 format('phi theta',7(2x,f5.1,3x)/3(9x,7(2x,f5.1,3x)/))
74 6667 format(i3,6x,1p7e10.3/3(9x,1p7e10.3/))
75 6668 format(t1,'radiances (f0=pi)')
76 6672 format(t8,'fdif_up',t26,1pe12.4,t40,'sbar',t64,1pe12.4)
77 6674 format(t8,'acrft_fdif_dn',t26,1pe12.4,t40,'acrft_fdif_up',
78  1 t64,1pe12.4)
79 c***********************************************************************
80  return
81  end
82 c************************************************************************
#define real
Definition: DbAlgOcean.cpp:26
#define pi
Definition: vincenty.c:23
subroutine headr(iprn)
Definition: headr.f:2
subroutine outcrfty(tmcfd, tmcfu, ocrftzu, ocrftzd)
Definition: outcrfty.f:2