Context Navigation

photo_O1D.f90 @ 3481cc1

10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bug

Last change on this file since 3481cc1 was 3481cc1, checked in by Ignacio Pisso <ip@…>, 4 years ago
move license from headers to a different file
Property mode set to `100644`
File size: 2.4 KB

Line
1	real function photo_O1D(sza)
2
3	!*****************************************************************************
4	! *
5	! *
6	! Author: A. Stohl *
7	! *
8	! Nov 2014 *
9	! *
10	! *
11	!*****************************************************************************
12	! *
13	! INPUT: *
14	! sza solar zenith angle (degrees) *
15	! *
16	! OUTPUT: *
17	! photo_O1D J(O1D) photoylsis rate *
18	! *
19	!*****************************************************************************
20
21	implicit none
22
23	integer :: iz,ik
24	real :: sza
25	real :: z1,z2,zg,f1,f2,dummy
26	real :: photo_NO2
27	integer, parameter :: nzenith=11
28	real, parameter :: pi=3.1415927
29	real, dimension(nzenith) :: zangle,fact_photo
30
31	! zangle: zenith angles for which fact_photo is tabulated
32	! fact_photo: conversion of photolysis rate of NO2 to photolysis
33	! rate of O3 into O1D as a function of solar zenith angle
34
35	zangle=(/0.,10.,20.,30.,40.,50.,60.,70.,78.,86.,90.0001/)
36	fact_photo=(/0.4616E-02,0.4478E-02,0.4131E-02,0.3583E-02,0.2867E-02,&
37	&0.2081E-02,0.1235E-02,0.5392E-03,0.2200E-03,0.1302E-03,0.0902E-03/)
38
39	if (sza.lt.90.) then
40	do iz=1,nzenith-1
41	if(sza.ge.zangle(iz)) ik=iz
42	end do
43	z1=1./cos(zangle(ik)*pi/180.)
44	z2=1./cos(zangle(ik+1)*pi/180.)
45	zg=1./cos(sza*pi/180.)
46	dummy=(zg-z1)/(z2-z1)
47	f1=alog(fact_photo(ik))
48	f2=alog(fact_photo(ik+1))
49	photo_NO2=1.45e-2exp(-0.4/cos(szapi/180.))
50	photo_O1D=photo_NO2exp(f1+(f2-f1)dummy)
51	else
52	photo_O1D=0.
53	endif
54
55	return
56
57	end function photo_O1D
58
59

Note: See TracBrowser for help on using the repository browser.