Context Navigation

jthermcalc_e107.F @ 2489

Last change on this file since 2489 was 2429, checked in by emillour, 5 years ago
Mars GCM: Bug fix in jthermcalc_e107.F: correctly account for the actual Sun-Mars distance. FGG
File size: 37.3 KB

Rev	Line
[705]	1	c**********************************************************************
	2
	3	subroutine jthermcalc_e107
[1266]	4	$ (ig,nlayer,chemthermod,rm,nesptherm,tx,iz,zenit,zday)
[705]	5
	6
	7	c feb 2002 fgg first version
	8	c nov 2002 fgg second version
	9	c
	10	c modified from paramhr.F
	11	c MAC July 2003
	12	c**********************************************************************
	13
[1266]	14	use param_v4_h, only: jfotsout,crscabsi2,
	15	. c1_16,c17_24,c25_29,c30_31,c32,c33,c34,c35,c36,
	16	. co2crsc195,co2crsc295,t0,
	17	. jabsifotsintpar,ninter,nz2,
	18	. nabs,e107,date_e107,e107_tab,
	19	. coefit0,coefit1,coefit2,coefit3,coefit4
[2429]	20	use comsaison_h, only: dist_sol
[1266]	21
[705]	22	implicit none
	23
[1684]	24	include "callkeys.h"
	25
[705]	26	c input and output variables
	27
[2429]	28	integer,intent(in) :: ig,nlayer
	29	integer,intent(in) :: chemthermod
	30	integer,intent(in) :: nesptherm !Number of species considered
	31	real,intent(in) :: rm(nlayer,nesptherm) !Densities (cm-3)
	32	real,intent(in) :: tx(nlayer) !temperature
	33	real,intent(in) :: zenit !SZA
	34	real,intent(in) :: iz(nlayer) !Local altitude
	35	real,intent(in) :: zday !Martian day after Ls=0
[705]	36
[2429]	37	! NB: no output variable! (computed jfotsout() is in module param_v4_h)
[705]	38
	39	c local parameters and variables
	40
[1266]	41	real co2colx(nlayer) !column density of CO2 (cm^-2)
	42	real o2colx(nlayer) !column density of O2(cm^-2)
	43	real o3pcolx(nlayer) !column density of O(3P)(cm^-2)
	44	real h2colx(nlayer) !H2 column density (cm-2)
	45	real h2ocolx(nlayer) !H2O column density (cm-2)
	46	real h2o2colx(nlayer) !column density of H2O2(cm^-2)
	47	real o3colx(nlayer) !O3 column density (cm-2)
	48	real n2colx(nlayer) !N2 column density (cm-2)
	49	real ncolx(nlayer) !N column density (cm-2)
	50	real nocolx(nlayer) !NO column density (cm-2)
	51	real cocolx(nlayer) !CO column density (cm-2)
	52	real hcolx(nlayer) !H column density (cm-2)
	53	real no2colx(nlayer) !NO2 column density (cm-2)
	54	real t2(nlayer)
	55	real coltemp(nlayer)
	56	real sigma(ninter,nlayer)
	57	real alfa(ninter,nlayer)
[705]	58	real realday
	59
	60	integer i,j,k,indexint !indexes
	61	character dn
	62	integer tinf,tsup
	63
	64
	65
	66	c variables used in interpolation
	67
	68	real*8 auxcoltab(nz2)
	69	real*8 auxjco2(nz2)
	70	real*8 auxjo2(nz2)
	71	real*8 auxjo3p(nz2)
	72	real*8 auxjh2o(nz2)
	73	real*8 auxjh2(nz2)
	74	real*8 auxjh2o2(nz2)
	75	real*8 auxjo3(nz2)
	76	real*8 auxjn2(nz2)
	77	real*8 auxjn(nz2)
	78	real*8 auxjno(nz2)
	79	real*8 auxjco(nz2)
	80	real*8 auxjh(nz2)
	81	real*8 auxjno2(nz2)
[1266]	82	real*8 wp(nlayer),wm(nlayer)
	83	real*8 auxcolinp(nlayer)
	84	integer auxind(nlayer)
[705]	85	integer auxi
	86	integer ind
[1266]	87	real*8 cortemp(nlayer)
[705]	88
	89	real*8 limdown !limits for interpolation
	90	real*8 limup ! ""
	91
	92	!!!ATTENTION. Here i_co2 has to have the same value than in chemthermos.F90
	93	!!!If the value is changed there, if has to be changed also here !!!!
	94	integer,parameter :: i_co2=1
	95
	96
	97	c***********************PROGRAM STARTS*****************************
	98
	99	if(zenit.gt.140.) then
	100	dn='n'
	101	else
	102	dn='d'
	103	end if
	104	if(dn.eq.'n') then
	105	return
	106	endif
	107
	108	!Initializing the photoabsorption coefficients
	109	jfotsout(:,:,:)=0.
	110
	111	!Auxiliar temperature to take into account the temperature dependence
	112	!of CO2 cross section
[1266]	113	do i=1,nlayer
[705]	114	t2(i)=tx(i)
	115	if(t2(i).lt.195.0) t2(i)=195.0
	116	if(t2(i).gt.295.0) t2(i)=295.0
	117	end do
	118
	119	!Calculation of column amounts
[1266]	120	call column(ig,nlayer,chemthermod,rm,nesptherm,tx,iz,zenit,
[705]	121	$ co2colx,o2colx,o3pcolx,h2colx,h2ocolx,
	122	$ h2o2colx,o3colx,n2colx,ncolx,nocolx,cocolx,hcolx,no2colx)
	123
	124	!Auxiliar column to include the temperature dependence
	125	!of CO2 cross section
[1266]	126	coltemp(nlayer)=co2colx(nlayer)*abs(t2(nlayer)-t0(nlayer))
	127	do i=nlayer-1,1,-1
[705]	128	coltemp(i)=!coltemp(i+1)+ PQ SE ELIMINA? REVISAR
	129	$ ( rm(i,i_co2) + rm(i+1,i_co2) ) * 0.5
	130	$ * 1e5 * (iz(i+1)-iz(i)) * abs(t2(i)-t0(i))
	131	end do
	132
	133	!Calculation of CO2 cross section at temperature t0(i)
[1266]	134	do i=1,nlayer
[705]	135	do indexint=24,32
	136	sigma(indexint,i)=co2crsc195(indexint-23)
	137	alfa(indexint,i)=((co2crsc295(indexint-23)
	138	$ /sigma(indexint,i))-1.)/(295.-t0(i))
	139	end do
	140	end do
	141
[1684]	142	if (solvarmod==0) then
	143	e107=fixed_euv_value
	144	else
	145	!E10.7 for the day: linear interpolation to tabulated values
	146	realday=mod(zday,669.)
	147	if(realday.lt.date_e107(1)) then
[705]	148	e107=e107_tab(1)
[1684]	149	else if(realday.ge.date_e107(669)) then
[762]	150	e107=e107_tab(669)
[1684]	151	else if(realday.ge.date_e107(1).and.
[762]	152	$ realday.lt.date_e107(669)) then
	153	do i=1,668
[705]	154	if(realday.ge.date_e107(i).and.
	155	$ realday.lt.date_e107(i+1)) then
	156	tinf=i
	157	tsup=i+1
[762]	158	e107=e107_tab(tinf)+(realday-date_e107(tinf))*
[705]	159	$ (e107_tab(tsup)-e107_tab(tinf))
	160	endif
	161	enddo
[1684]	162	endif
	163	endif ! of if (solvarmod==0)
[762]	164
[705]	165	!Photoabsorption coefficients at TOA as a function of E10.7
	166	do j=1,nabs
	167	do indexint=1,ninter
[1266]	168	jfotsout(indexint,j,nlayer)=coefit0(indexint,j)+
[705]	169	$ coefit1(indexint,j)e107+coefit2(indexint,j)e107**2+
	170	$ coefit3(indexint,j)e1073+coefit4(indexint,j)e107**4
	171	enddo
	172	enddo
	173	! Interpolation to the tabulated photoabsorption coefficients for each species
	174	! in each spectral interval
	175
	176
	177	c auxcolinp-> Actual atmospheric column
	178	c auxj*-> Tabulated photoabsorption coefficients
	179	c auxcoltab-> Tabulated atmospheric columns
	180
	181	ccccccccccccccccccccccccccccccc
	182	c 0.1,5.0 (int 1)
	183	c
	184	c Absorption by:
	185	c CO2, O2, O, H2, N
	186	ccccccccccccccccccccccccccccccc
	187
	188	c Input atmospheric column
	189	indexint=1
[1266]	190	do i=1,nlayer
	191	auxcolinp(nlayer-i+1) = co2colx(i)*crscabsi2(1,indexint) +
[705]	192	$ o2colx(i)*crscabsi2(2,indexint) +
	193	$ o3pcolx(i)*crscabsi2(3,indexint) +
	194	$ h2colx(i)*crscabsi2(5,indexint) +
	195	$ ncolx(i)*crscabsi2(9,indexint)
	196	end do
	197	limdown=1.e-20
	198	limup=1.e26
	199
	200
	201	c Interpolations
	202
	203	do i=1,nz2
	204	auxi = nz2-i+1
	205	!CO2 tabulated coefficient
	206	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	207	!O2 tabulated coefficient
	208	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	209	!O3p tabulated coefficient
	210	auxjo3p(i) = jabsifotsintpar(auxi,3,indexint)
	211	!H2 tabulated coefficient
	212	auxjh2(i) = jabsifotsintpar(auxi,5,indexint)
	213	!Tabulated column
	214	auxcoltab(i) = c1_16(auxi,indexint)
	215	enddo
	216	!Only if chemthermod.ge.2
	217	!N tabulated coefficient
	218	if(chemthermod.ge.2) then
	219	do i=1,nz2
	220	auxjn(i) = jabsifotsintpar(nz2-i+1,9,indexint)
	221	enddo
	222	endif
	223
	224	call interfast
[1266]	225	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	226	do i=1,nlayer
[705]	227	ind=auxind(i)
[1266]	228	auxi=nlayer-i+1
[705]	229	!CO2 interpolated coefficient
[1266]	230	jfotsout(indexint,1,auxi) = jfotsout(indexint,1,nlayer) *
[705]	231	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind))
	232	!O2 interpolated coefficient
[1266]	233	jfotsout(indexint,2,auxi) = jfotsout(indexint,2,nlayer) *
[705]	234	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind))
	235	!O3p interpolated coefficient
[1266]	236	jfotsout(indexint,3,auxi) = jfotsout(indexint,3,nlayer) *
[705]	237	$ (wm(i)auxjo3p(ind+1) + wp(i)auxjo3p(ind))
	238	!H2 interpolated coefficient
	239	jfotsout(indexint,5,auxi) = jfotsout(indexint,5,auxi) *
	240	$ (wm(i)auxjh2(ind+1) + wp(i)auxjh2(ind))
	241	enddo
	242	!Only if chemthermod.ge.2
	243	!N interpolated coefficient
	244	if(chemthermod.ge.2) then
[1266]	245	do i=1,nlayer
[705]	246	ind=auxind(i)
[1266]	247	jfotsout(indexint,9,nlayer-i+1) =
	248	$ jfotsout(indexint,9,nlayer) *
[705]	249	$ (wm(i)auxjn(ind+1) + wp(i)auxjn(ind))
	250	enddo
	251	endif
	252
	253
	254	c End interval 1
	255
	256
	257	ccccccccccccccccccccccccccccccc
	258	c 5-80.5nm (int 2-15)
	259	c
	260	c Absorption by:
	261	c CO2, O2, O, H2, N2, N,
	262	c NO, CO, H, NO2
	263	ccccccccccccccccccccccccccccccc
	264
	265	c Input atmospheric column
	266	do indexint=2,15
[1266]	267	do i=1,nlayer
	268	auxcolinp(nlayer-i+1) = co2colx(i)*crscabsi2(1,indexint)+
[705]	269	$ o2colx(i)*crscabsi2(2,indexint)+
	270	$ o3pcolx(i)*crscabsi2(3,indexint)+
	271	$ h2colx(i)*crscabsi2(5,indexint)+
	272	$ n2colx(i)*crscabsi2(8,indexint)+
	273	$ ncolx(i)*crscabsi2(9,indexint)+
	274	$ nocolx(i)*crscabsi2(10,indexint)+
	275	$ cocolx(i)*crscabsi2(11,indexint)+
	276	$ hcolx(i)*crscabsi2(12,indexint)+
	277	$ no2colx(i)*crscabsi2(13,indexint)
	278	end do
	279
	280	c Interpolations
	281
	282	do i=1,nz2
	283	auxi = nz2-i+1
	284	!O2 tabulated coefficient
	285	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	286	!O3p tabulated coefficient
	287	auxjo3p(i) = jabsifotsintpar(auxi,3,indexint)
	288	!CO2 tabulated coefficient
	289	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	290	!H2 tabulated coefficient
	291	auxjh2(i) = jabsifotsintpar(auxi,5,indexint)
	292	!N2 tabulated coefficient
	293	auxjn2(i) = jabsifotsintpar(auxi,8,indexint)
	294	!CO tabulated coefficient
	295	auxjco(i) = jabsifotsintpar(auxi,11,indexint)
	296	!H tabulated coefficient
	297	auxjh(i) = jabsifotsintpar(auxi,12,indexint)
	298	!tabulated column
	299	auxcoltab(i) = c1_16(auxi,indexint)
	300	enddo
	301	!Only if chemthermod.ge.2
	302	if(chemthermod.ge.2) then
	303	do i=1,nz2
	304	auxi = nz2-i+1
	305	!N tabulated coefficient
	306	auxjn(i) = jabsifotsintpar(auxi,9,indexint)
	307	!NO tabulated coefficient
	308	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	309	!NO2 tabulated coefficient
	310	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	311	enddo
	312	endif
	313
[1266]	314	call interfast(wm,wp,auxind,auxcolinp,nlayer,
[705]	315	$ auxcoltab,nz2,limdown,limup)
[1266]	316	do i=1,nlayer
[705]	317	ind=auxind(i)
[1266]	318	auxi = nlayer-i+1
[705]	319	!O2 interpolated coefficient
	320	jfotsout(indexint,2,auxi) =
[1266]	321	$ jfotsout(indexint,2,nlayer) *
[705]	322	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind))
	323	!O3p interpolated coefficient
	324	jfotsout(indexint,3,auxi) =
[1266]	325	$ jfotsout(indexint,3,nlayer) *
[705]	326	$ (wm(i)auxjo3p(ind+1) + wp(i)auxjo3p(ind))
	327	!CO2 interpolated coefficient
	328	jfotsout(indexint,1,auxi) =
[1266]	329	$ jfotsout(indexint,1,nlayer) *
[705]	330	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind))
	331	!H2 interpolated coefficient
	332	jfotsout(indexint,5,auxi) =
[1266]	333	$ jfotsout(indexint,5,nlayer) *
[705]	334	$ (wm(i)auxjh2(ind+1) + wp(i)auxjh2(ind))
	335	!N2 interpolated coefficient
	336	jfotsout(indexint,8,auxi) =
[1266]	337	$ jfotsout(indexint,8,nlayer) *
[705]	338	$ (wm(i)auxjn2(ind+1) + wp(i)auxjn2(ind))
	339	!CO interpolated coefficient
	340	jfotsout(indexint,11,auxi) =
[1266]	341	$ jfotsout(indexint,11,nlayer) *
[705]	342	$ (wm(i)auxjco(ind+1) + wp(i)auxjco(ind))
	343	!H interpolated coefficient
	344	jfotsout(indexint,12,auxi) =
[1266]	345	$ jfotsout(indexint,12,nlayer) *
[705]	346	$ (wm(i)auxjh(ind+1) + wp(i)auxjh(ind))
	347	enddo
	348	!Only if chemthermod.ge.2
	349	if(chemthermod.ge.2) then
[1266]	350	do i=1,nlayer
[705]	351	ind=auxind(i)
[1266]	352	auxi = nlayer-i+1
[705]	353	!N interpolated coefficient
	354	jfotsout(indexint,9,auxi) =
[1266]	355	$ jfotsout(indexint,9,nlayer) *
[705]	356	$ (wm(i)auxjn(ind+1) + wp(i)auxjn(ind))
	357	!NO interpolated coefficient
	358	jfotsout(indexint,10,auxi)=
[1266]	359	$ jfotsout(indexint,10,nlayer) *
[705]	360	$ (wm(i)auxjno(ind+1) + wp(i)auxjno(ind))
	361	!NO2 interpolated coefficient
	362	jfotsout(indexint,13,auxi)=
[1266]	363	$ jfotsout(indexint,13,nlayer) *
[705]	364	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	365	enddo
	366	endif
	367	end do
	368
	369	c End intervals 2-15
	370
	371
	372	ccccccccccccccccccccccccccccccc
	373	c 80.6-90.8nm (int16)
	374	c
	375	c Absorption by:
	376	c CO2, O2, O, N2, N, NO,
	377	c CO, H, NO2
	378	ccccccccccccccccccccccccccccccc
	379
	380	c Input atmospheric column
	381	indexint=16
[1266]	382	do i=1,nlayer
	383	auxcolinp(nlayer-i+1) = co2colx(i)*crscabsi2(1,indexint)+
[705]	384	$ o2colx(i)*crscabsi2(2,indexint)+
	385	$ o3pcolx(i)*crscabsi2(3,indexint)+
	386	$ n2colx(i)*crscabsi2(8,indexint)+
	387	$ ncolx(i)*crscabsi2(9,indexint)+
	388	$ nocolx(i)*crscabsi2(10,indexint)+
	389	$ cocolx(i)*crscabsi2(11,indexint)+
	390	$ hcolx(i)*crscabsi2(12,indexint)+
	391	$ no2colx(i)*crscabsi2(13,indexint)
	392	end do
	393
	394	c Interpolations
	395
	396	do i=1,nz2
	397	auxi = nz2-i+1
	398	!O2 tabulated coefficient
	399	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	400	!CO2 tabulated coefficient
	401	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	402	!O3p tabulated coefficient
	403	auxjo3p(i) = jabsifotsintpar(auxi,3,indexint)
	404	!N2 tabulated coefficient
	405	auxjn2(i) = jabsifotsintpar(auxi,8,indexint)
	406	!CO tabulated coefficient
	407	auxjco(i) = jabsifotsintpar(auxi,11,indexint)
	408	!H tabulated coefficient
	409	auxjh(i) = jabsifotsintpar(auxi,12,indexint)
	410	!NO2 tabulated coefficient
	411	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	412	!Tabulated column
	413	auxcoltab(i) = c1_16(auxi,indexint)
	414	enddo
	415	!Only if chemthermod.ge.2
	416	if(chemthermod.ge.2) then
	417	do i=1,nz2
	418	auxi = nz2-i+1
	419	!N tabulated coefficient
	420	auxjn(i) = jabsifotsintpar(auxi,9,indexint)
	421	!NO tabulated coefficient
	422	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	423	!NO2 tabulated coefficient
	424	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	425	enddo
	426	endif
	427
	428	call interfast
[1266]	429	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	430	do i=1,nlayer
[705]	431	ind=auxind(i)
[1266]	432	auxi = nlayer-i+1
[705]	433	!O2 interpolated coefficient
[1266]	434	jfotsout(indexint,2,auxi) = jfotsout(indexint,2,nlayer) *
[705]	435	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind))
	436	!CO2 interpolated coefficient
[1266]	437	jfotsout(indexint,1,auxi) = jfotsout(indexint,1,nlayer) *
[705]	438	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind))
	439	!O3p interpolated coefficient
[1266]	440	jfotsout(indexint,3,auxi) = jfotsout(indexint,3,nlayer) *
[705]	441	$ (wm(i)auxjo3p(ind+1) + wp(i)auxjo3p(ind))
	442	!N2 interpolated coefficient
[1266]	443	jfotsout(indexint,8,auxi) = jfotsout(indexint,8,nlayer) *
[705]	444	$ (wm(i)auxjn2(ind+1) + wp(i)auxjn2(ind))
	445	!CO interpolated coefficient
	446	jfotsout(indexint,11,auxi) =
[1266]	447	$ jfotsout(indexint,11,nlayer) *
[705]	448	$ (wm(i)auxjco(ind+1) + wp(i)auxjco(ind))
	449	!H interpolated coefficient
	450	jfotsout(indexint,12,auxi) =
[1266]	451	$ jfotsout(indexint,12,nlayer) *
[705]	452	$ (wm(i)auxjh(ind+1) + wp(i)auxjh(ind))
	453	enddo
	454	!Only if chemthermod.ge.2
	455	if(chemthermod.ge.2) then
[1266]	456	do i=1,nlayer
[705]	457	ind=auxind(i)
[1266]	458	auxi = nlayer-i+1
[705]	459	!N interpolated coefficient
	460	jfotsout(indexint,9,auxi) =
[1266]	461	$ jfotsout(indexint,9,nlayer) *
[705]	462	$ (wm(i)auxjn(ind+1) + wp(i)auxjn(ind))
	463	!NO interpolated coefficient
	464	jfotsout(indexint,10,auxi) =
[1266]	465	$ jfotsout(indexint,10,nlayer) *
[705]	466	$ (wm(i)auxjno(ind+1) + wp(i)auxjno(ind))
	467	!NO2 interpolated coefficient
	468	jfotsout(indexint,13,auxi) =
[1266]	469	$ jfotsout(indexint,13,nlayer) *
[705]	470	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	471	enddo
	472	endif
	473	c End interval 16
	474
	475
	476	ccccccccccccccccccccccccccccccc
	477	c 90.9-119.5nm (int 17-24)
	478	c
	479	c Absorption by:
	480	c CO2, O2, N2, NO, CO, NO2
	481	ccccccccccccccccccccccccccccccc
	482
	483	c Input column
	484
[1266]	485	do i=1,nlayer
	486	auxcolinp(nlayer-i+1) = co2colx(i) + o2colx(i) + n2colx(i) +
[705]	487	$ nocolx(i) + cocolx(i) + no2colx(i)
	488	end do
	489
	490	do indexint=17,24
	491
	492	c Interpolations
	493
	494	do i=1,nz2
	495	auxi = nz2-i+1
	496	!CO2 tabulated coefficient
	497	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	498	!O2 tabulated coefficient
	499	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	500	!N2 tabulated coefficient
	501	auxjn2(i) = jabsifotsintpar(auxi,8,indexint)
	502	!CO tabulated coefficient
	503	auxjco(i) = jabsifotsintpar(auxi,11,indexint)
	504	!Tabulated column
	505	auxcoltab(i) = c17_24(auxi)
	506	enddo
	507	!Only if chemthermod.ge.2
	508	if(chemthermod.ge.2) then
	509	do i=1,nz2
	510	auxi = nz2-i+1
	511	!NO tabulated coefficient
	512	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	513	!NO2 tabulated coefficient
	514	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	515	enddo
	516	endif
	517
	518	call interfast
[1266]	519	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
[705]	520	!Correction to include T variation of CO2 cross section
	521	if(indexint.eq.24) then
[1266]	522	do i=1,nlayer
	523	auxi = nlayer-i+1
[705]	524	if(sigma(indexint,auxi)*
	525	$ alfa(indexint,auxi)*coltemp(auxi)
	526	$ .lt.60.) then
	527	cortemp(i)=exp(-sigma(indexint,auxi)*
	528	$ alfa(indexint,auxi)*coltemp(auxi))
	529	else
	530	cortemp(i)=0.
	531	end if
	532	enddo
	533	else
[1266]	534	do i=1,nlayer
[705]	535	cortemp(i)=1.
	536	enddo
	537	end if
[1266]	538	do i=1,nlayer
[705]	539	ind=auxind(i)
[1266]	540	auxi = nlayer-i+1
[705]	541	!O2 interpolated coefficient
	542	jfotsout(indexint,2,auxi) =
[1266]	543	$ jfotsout(indexint,2,nlayer) *
[705]	544	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind)) *
	545	$ cortemp(i)
	546	!CO2 interpolated coefficient
	547	jfotsout(indexint,1,auxi) =
[1266]	548	$ jfotsout(indexint,1,nlayer) *
[705]	549	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind))
	550	$ * cortemp(i)
	551	if(indexint.eq.24) jfotsout(indexint,1,auxi)=
	552	$ jfotsout(indexint,1,auxi)*
	553	$ (1+alfa(indexint,auxi)*
	554	$ (t2(auxi)-t0(auxi)))
	555	!N2 interpolated coefficient
	556	jfotsout(indexint,8,auxi) =
[1266]	557	$ jfotsout(indexint,8,nlayer) *
[705]	558	$ (wm(i)auxjn2(ind+1) + wp(i)auxjn2(ind)) *
	559	$ cortemp(i)
	560	!CO interpolated coefficient
	561	jfotsout(indexint,11,auxi) =
[1266]	562	$ jfotsout(indexint,11,nlayer) *
[705]	563	$ (wm(i)auxjco(ind+1) + wp(i)auxjco(ind)) *
	564	$ cortemp(i)
	565	enddo
	566	!Only if chemthermod.ge.2
	567	if(chemthermod.ge.2) then
[1266]	568	do i=1,nlayer
[705]	569	ind=auxind(i)
[1266]	570	auxi = nlayer-i+1
[705]	571	!NO interpolated coefficient
	572	jfotsout(indexint,10,auxi)=
[1266]	573	$ jfotsout(indexint,10,nlayer) *
[705]	574	$ (wm(i)auxjno(ind+1) + wp(i)auxjno(ind)) *
	575	$ cortemp(i)
	576	!NO2 interpolated coefficient
	577	jfotsout(indexint,13,auxi)=
[1266]	578	$ jfotsout(indexint,13,nlayer) *
[705]	579	$ (wm(i)auxjno2(ind+1)+ wp(i)auxjno2(ind)) *
	580	$ cortemp(i)
	581	enddo
	582	endif
	583	end do
	584	c End intervals 17-24
	585
	586
	587	ccccccccccccccccccccccccccccccc
	588	c 119.6-167.0nm (int 25-29)
	589	c
	590	c Absorption by:
	591	c CO2, O2, H2O, H2O2, NO,
	592	c CO, NO2
	593	ccccccccccccccccccccccccccccccc
	594
	595	c Input atmospheric column
	596
[1266]	597	do i=1,nlayer
	598	auxcolinp(nlayer-i+1) = co2colx(i) + o2colx(i) + h2ocolx(i) +
[705]	599	$ h2o2colx(i) + nocolx(i) + cocolx(i) + no2colx(i)
	600	end do
	601
	602	do indexint=25,29
	603
	604	c Interpolations
	605
	606	do i=1,nz2
	607	auxi = nz2-i+1
	608	!CO2 tabulated coefficient
	609	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	610	!O2 tabulated coefficient
	611	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	612	!H2O tabulated coefficient
	613	auxjh2o(i) = jabsifotsintpar(auxi,4,indexint)
	614	!H2O2 tabulated coefficient
	615	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	616	!CO tabulated coefficient
	617	auxjco(i) = jabsifotsintpar(auxi,11,indexint)
	618	!Tabulated column
	619	auxcoltab(i) = c25_29(auxi)
	620	enddo
	621	!Only if chemthermod.ge.2
	622	if(chemthermod.ge.2) then
	623	do i=1,nz2
	624	auxi = nz2-i+1
	625	!NO tabulated coefficient
	626	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	627	!NO2 tabulated coefficient
	628	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	629	enddo
	630	endif
	631	call interfast
[1266]	632	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	633	do i=1,nlayer
[705]	634	ind=auxind(i)
[1266]	635	auxi = nlayer-i+1
[705]	636	!Correction to include T variation of CO2 cross section
	637	if(sigma(indexint,auxi)alfa(indexint,auxi)
	638	$ coltemp(auxi).lt.60.) then
	639	cortemp(i)=exp(-sigma(indexint,auxi)*
	640	$ alfa(indexint,auxi)*coltemp(auxi))
	641	else
	642	cortemp(i)=0.
	643	end if
	644	!CO2 interpolated coefficient
	645	jfotsout(indexint,1,auxi) =
[1266]	646	$ jfotsout(indexint,1,nlayer) *
[705]	647	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind)) *
	648	$ cortemp(i) *
	649	$ (1+alfa(indexint,auxi)*
	650	$ (t2(auxi)-t0(auxi)))
	651	!O2 interpolated coefficient
	652	jfotsout(indexint,2,auxi) =
[1266]	653	$ jfotsout(indexint,2,nlayer) *
[705]	654	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind)) *
	655	$ cortemp(i)
	656	!H2O interpolated coefficient
	657	jfotsout(indexint,4,auxi) =
[1266]	658	$ jfotsout(indexint,4,nlayer) *
[705]	659	$ (wm(i)auxjh2o(ind+1) + wp(i)auxjh2o(ind)) *
	660	$ cortemp(i)
	661	!H2O2 interpolated coefficient
	662	jfotsout(indexint,6,auxi) =
[1266]	663	$ jfotsout(indexint,6,nlayer) *
[705]	664	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind)) *
	665	$ cortemp(i)
	666	!CO interpolated coefficient
	667	jfotsout(indexint,11,auxi) =
[1266]	668	$ jfotsout(indexint,11,nlayer) *
[705]	669	$ (wm(i)auxjco(ind+1) + wp(i)auxjco(ind)) *
	670	$ cortemp(i)
	671	enddo
	672	!Only if chemthermod.ge.2
	673	if(chemthermod.ge.2) then
[1266]	674	do i=1,nlayer
[705]	675	ind=auxind(i)
[1266]	676	auxi = nlayer-i+1
[705]	677	!NO interpolated coefficient
	678	jfotsout(indexint,10,auxi)=
[1266]	679	$ jfotsout(indexint,10,nlayer) *
[705]	680	$ (wm(i)auxjno(ind+1) + wp(i)auxjno(ind)) *
	681	$ cortemp(i)
	682	!NO2 interpolated coefficient
	683	jfotsout(indexint,13,auxi)=
[1266]	684	$ jfotsout(indexint,13,nlayer) *
[705]	685	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind)) *
	686	$ cortemp(i)
	687	enddo
	688	endif
	689
	690	end do
	691
	692	c End intervals 25-29
	693
	694
	695	cccccccccccccccccccccccccccccccc
	696	c 167.1-202.5nm (int 30-31)
	697	c
	698	c Absorption by:
	699	c CO2, O2, H2O, H2O2, NO,
	700	c NO2
	701	cccccccccccccccccccccccccccccccc
	702
	703	c Input atmospheric column
	704
[1266]	705	do i=1,nlayer
	706	auxcolinp(nlayer-i+1) = co2colx(i) + o2colx(i) + h2ocolx(i) +
[705]	707	$ h2o2colx(i) + nocolx(i) + no2colx(i)
	708	end do
	709
	710	c Interpolation
	711
	712	do indexint=30,31
	713
	714	do i=1,nz2
	715	auxi = nz2-i+1
	716	!CO2 tabulated coefficient
	717	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	718	!O2 tabulated coefficient
	719	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	720	!H2O tabulated coefficient
	721	auxjh2o(i) = jabsifotsintpar(auxi,4,indexint)
	722	!H2O2 tabulated coefficient
	723	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	724	!Tabulated column
	725	auxcoltab(i) = c30_31(auxi)
	726	enddo
	727	!Only if chemthermod.ge.2
	728	if(chemthermod.ge.2) then
	729	do i=1,nz2
	730	auxi = nz2-i+1
	731	!NO tabulated coefficient
	732	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	733	!NO2 tabulated coefficient
	734	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	735	enddo
	736	endif
	737
	738	call interfast
[1266]	739	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	740	do i=1,nlayer
[705]	741	ind=auxind(i)
[1266]	742	auxi = nlayer-i+1
[705]	743	!Correction to include T variation of CO2 cross section
	744	if(sigma(indexint,auxi)alfa(indexint,auxi)
	745	$ coltemp(auxi).lt.60.) then
	746	cortemp(i)=exp(-sigma(indexint,auxi)*
	747	$ alfa(indexint,auxi)*coltemp(auxi))
	748	else
	749	cortemp(i)=0.
	750	end if
	751	!CO2 interpolated coefficient
	752	jfotsout(indexint,1,auxi) =
[1266]	753	$ jfotsout(indexint,1,nlayer) *
[705]	754	$ (wm(i)auxjco2(ind+1) + wp(i)auxjco2(ind)) *
	755	$ cortemp(i) *
	756	$ (1+alfa(indexint,auxi)*
	757	$ (t2(auxi)-t0(auxi)))
	758	!O2 interpolated coefficient
	759	jfotsout(indexint,2,auxi) =
[1266]	760	$ jfotsout(indexint,2,nlayer) *
[705]	761	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind)) *
	762	$ cortemp(i)
	763	!H2O interpolated coefficient
	764	jfotsout(indexint,4,auxi) =
[1266]	765	$ jfotsout(indexint,4,nlayer) *
[705]	766	$ (wm(i)auxjh2o(ind+1) + wp(i)auxjh2o(ind)) *
	767	$ cortemp(i)
	768	!H2O2 interpolated coefficient
	769	jfotsout(indexint,6,auxi) =
[1266]	770	$ jfotsout(indexint,6,nlayer) *
[705]	771	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind)) *
	772	$ cortemp(i)
	773	enddo
	774	!Only if chemthermod.ge.2
	775	if(chemthermod.ge.2) then
[1266]	776	do i=1,nlayer
[705]	777	ind=auxind(i)
[1266]	778	auxi = nlayer-i+1
[705]	779	!NO interpolated coefficient
	780	jfotsout(indexint,10,auxi)=
[1266]	781	$ jfotsout(indexint,10,nlayer) *
[705]	782	$ (wm(i)auxjno(ind+1) +wp(i)auxjno(ind)) *
	783	$ cortemp(i)
	784	!NO2 interpolated coefficient
	785	jfotsout(indexint,13,auxi)=
	786	$ jfotsout(indexint,13,auxi) *
	787	$ (wm(i)auxjno2(ind+1)+wp(i)auxjno2(ind)) *
	788	$ cortemp(i)
	789	enddo
	790	endif
	791
	792	end do
	793
	794	c End intervals 30-31
	795
	796
	797	ccccccccccccccccccccccccccccccc
	798	c 202.6-210.0nm (int 32)
	799	c
	800	c Absorption by:
	801	c CO2, O2, H2O2, NO, NO2
	802	ccccccccccccccccccccccccccccccc
	803
	804	c Input atmospheric column
	805
	806	indexint=32
[1266]	807	do i=1,nlayer
	808	auxcolinp(nlayer-i+1) =co2colx(i) + o2colx(i) + h2o2colx(i) +
[705]	809	$ nocolx(i) + no2colx(i)
	810	end do
	811
	812	c Interpolation
	813
	814	do i=1,nz2
	815	auxi = nz2-i+1
	816	!CO2 tabulated coefficient
	817	auxjco2(i) = jabsifotsintpar(auxi,1,indexint)
	818	!O2 tabulated coefficient
	819	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	820	!H2O2 tabulated coefficient
	821	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	822	!Tabulated column
	823	auxcoltab(i) = c32(auxi)
	824	enddo
	825	!Only if chemthermod.ge.2
	826	if(chemthermod.ge.2) then
	827	do i=1,nz2
	828	auxi = nz2-i+1
	829	!NO tabulated coefficient
	830	auxjno(i) = jabsifotsintpar(auxi,10,indexint)
	831	!NO2 tabulated coefficient
	832	auxjno2(i) = jabsifotsintpar(auxi,13,indexint)
	833	enddo
	834	endif
	835	call interfast
[1266]	836	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	837	do i=1,nlayer
[705]	838	ind=auxind(i)
[1266]	839	auxi = nlayer-i+1
[705]	840	!Correction to include T variation of CO2 cross section
[1266]	841	if(sigma(indexint,nlayer-i+1)alfa(indexint,auxi)
[705]	842	$ coltemp(auxi).lt.60.) then
	843	cortemp(i)=exp(-sigma(indexint,auxi)*
	844	$ alfa(indexint,auxi)*coltemp(auxi))
	845	else
	846	cortemp(i)=0.
	847	end if
	848	!CO2 interpolated coefficient
	849	jfotsout(indexint,1,auxi) =
[1266]	850	$ jfotsout(indexint,1,nlayer) *
[705]	851	$ (wm(i)auxjco2(ind+1)+wp(i)auxjco2(ind)) *
	852	$ cortemp(i) *
	853	$ (1+alfa(indexint,auxi)*
	854	$ (t2(auxi)-t0(auxi)))
	855	!O2 interpolated coefficient
	856	jfotsout(indexint,2,auxi) =
[1266]	857	$ jfotsout(indexint,2,nlayer) *
[705]	858	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind)) *
	859	$ cortemp(i)
	860	!H2O2 interpolated coefficient
	861	jfotsout(indexint,6,auxi) =
[1266]	862	$ jfotsout(indexint,6,nlayer) *
[705]	863	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind)) *
	864	$ cortemp(i)
	865	enddo
	866	!Only if chemthermod.ge.2
	867	if(chemthermod.ge.2) then
[1266]	868	do i=1,nlayer
	869	auxi = nlayer-i+1
[705]	870	ind=auxind(i)
	871	!NO interpolated coefficient
	872	jfotsout(indexint,10,auxi) =
[1266]	873	$ jfotsout(indexint,10,nlayer) *
[705]	874	$ (wm(i)auxjno(ind+1) + wp(i)auxjno(ind)) *
	875	$ cortemp(i)
	876	!NO2 interpolated coefficient
	877	jfotsout(indexint,13,auxi) =
[1266]	878	$ jfotsout(indexint,13,nlayer) *
[705]	879	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind)) *
	880	$ cortemp(i)
	881	enddo
	882	endif
	883
	884	c End of interval 32
	885
	886
	887	ccccccccccccccccccccccccccccccc
	888	c 210.1-231.0nm (int 33)
	889	c
	890	c Absorption by:
	891	c O2, H2O2, NO2
	892	ccccccccccccccccccccccccccccccc
	893
	894	c Input atmospheric column
	895
	896	indexint=33
[1266]	897	do i=1,nlayer
	898	auxcolinp(nlayer-i+1) = o2colx(i) + h2o2colx(i) + no2colx(i)
[705]	899	end do
	900
	901	c Interpolation
	902
	903	do i=1,nz2
	904	auxi = nz2-i+1
	905	!O2 tabulated coefficient
	906	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	907	!H2O2 tabulated coefficient
	908	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	909	!Tabulated column
	910	auxcoltab(i) = c33(auxi)
	911	enddo
	912	!Only if chemthermod.ge.2
	913	if(chemthermod.ge.2) then
	914	do i=1,nz2
	915	!NO2 tabulated coefficient
	916	auxjno2(i) = jabsifotsintpar(nz2-i+1,13,indexint)
	917	enddo
	918	endif
	919	call interfast
[1266]	920	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	921	do i=1,nlayer
[705]	922	ind=auxind(i)
[1266]	923	auxi = nlayer-i+1
[705]	924	!O2 interpolated coefficient
[1266]	925	jfotsout(indexint,2,auxi) = jfotsout(indexint,2,nlayer) *
[705]	926	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind))
	927	!H2O2 interpolated coefficient
[1266]	928	jfotsout(indexint,6,auxi) = jfotsout(indexint,6,nlayer) *
[705]	929	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind))
	930	enddo
	931	!Only if chemthermod.ge.2
	932	if(chemthermod.ge.2) then
[1266]	933	do i=1,nlayer
[705]	934	ind=auxind(i)
	935	!NO2 interpolated coefficient
[1266]	936	jfotsout(indexint,13,nlayer-i+1) =
	937	$ jfotsout(indexint,13,nlayer) *
[705]	938	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	939	enddo
	940	endif
	941
	942	c End of interval 33
	943
	944
	945	ccccccccccccccccccccccccccccccc
	946	c 231.1-240.0nm (int 34)
	947	c
	948	c Absorption by:
	949	c O2, H2O2, O3, NO2
	950	ccccccccccccccccccccccccccccccc
	951
	952	c Input atmospheric column
	953
	954	indexint=34
[1266]	955	do i=1,nlayer
	956	auxcolinp(nlayer-i+1) = h2o2colx(i) + o2colx(i) + o3colx(i) +
[705]	957	$ no2colx(i)
	958	end do
	959
	960	c Interpolation
	961
	962	do i=1,nz2
	963	auxi = nz2-i+1
	964	!O2 tabulated coefficient
	965	auxjo2(i) = jabsifotsintpar(auxi,2,indexint)
	966	!H2O2 tabulated coefficient
	967	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	968	!O3 tabulated coefficient
	969	auxjo3(i) = jabsifotsintpar(auxi,7,indexint)
	970	!Tabulated column
	971	auxcoltab(i) = c34(nz2-i+1)
	972	enddo
	973	!Only if chemthermod.ge.2
	974	if(chemthermod.ge.2) then
	975	do i=1,nz2
	976	!NO2 tabulated coefficient
	977	auxjno2(i) = jabsifotsintpar(nz2-i+1,13,indexint)
	978	enddo
	979	endif
	980	call interfast
[1266]	981	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	982	do i=1,nlayer
[705]	983	ind=auxind(i)
[1266]	984	auxi = nlayer-i+1
[705]	985	!O2 interpolated coefficient
[1266]	986	jfotsout(indexint,2,auxi) = jfotsout(indexint,2,nlayer) *
[705]	987	$ (wm(i)auxjo2(ind+1) + wp(i)auxjo2(ind))
	988	!H2O2 interpolated coefficient
[1266]	989	jfotsout(indexint,6,auxi) = jfotsout(indexint,6,nlayer) *
[705]	990	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind))
	991	!O3 interpolated coefficient
[1266]	992	jfotsout(indexint,7,auxi) = jfotsout(indexint,7,nlayer) *
[705]	993	$ (wm(i)auxjo3(ind+1) + wp(i)auxjo3(ind))
	994	enddo
	995	!Only if chemthermod.ge.2
	996	if(chemthermod.ge.2) then
[1266]	997	do i=1,nlayer
[705]	998	ind=auxind(i)
	999	!NO2 interpolated coefficient
[1266]	1000	jfotsout(indexint,13,nlayer-i+1) =
	1001	$ jfotsout(indexint,13,nlayer) *
[705]	1002	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	1003	enddo
	1004	endif
	1005
	1006	c End of interval 34
	1007
	1008
	1009	ccccccccccccccccccccccccccccccc
	1010	c 240.1-337.7nm (int 35)
	1011	c
	1012	c Absorption by:
	1013	c H2O2, O3, NO2
	1014	ccccccccccccccccccccccccccccccc
	1015
	1016	c Input atmospheric column
	1017
	1018	indexint=35
[1266]	1019	do i=1,nlayer
	1020	auxcolinp(nlayer-i+1) = h2o2colx(i) + o3colx(i) + no2colx(i)
[705]	1021	end do
	1022
	1023	c Interpolation
	1024
	1025	do i=1,nz2
	1026	auxi = nz2-i+1
	1027	!H2O2 tabulated coefficient
	1028	auxjh2o2(i) = jabsifotsintpar(auxi,6,indexint)
	1029	!O3 tabulated coefficient
	1030	auxjo3(i) = jabsifotsintpar(auxi,7,indexint)
	1031	!Tabulated column
	1032	auxcoltab(i) = c35(auxi)
	1033	enddo
	1034	!Only if chemthermod.ge.2
	1035	if(chemthermod.ge.2) then
	1036	do i=1,nz2
	1037	!NO2 tabulated coefficient
	1038	auxjno2(i) = jabsifotsintpar(nz2-i+1,13,indexint)
	1039	enddo
	1040	endif
	1041	call interfast
[1266]	1042	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	1043	do i=1,nlayer
[705]	1044	ind=auxind(i)
[1266]	1045	auxi = nlayer-i+1
[705]	1046	!H2O2 interpolated coefficient
[1266]	1047	jfotsout(indexint,6,auxi) = jfotsout(indexint,6,nlayer) *
[705]	1048	$ (wm(i)auxjh2o2(ind+1) + wp(i)auxjh2o2(ind))
	1049	!O3 interpolated coefficient
[1266]	1050	jfotsout(indexint,7,auxi) = jfotsout(indexint,7,nlayer) *
[705]	1051	$ (wm(i)auxjo3(ind+1) + wp(i)auxjo3(ind))
	1052	enddo
	1053	if(chemthermod.ge.2) then
[1266]	1054	do i=1,nlayer
[705]	1055	ind=auxind(i)
	1056	!NO2 interpolated coefficient
[1266]	1057	jfotsout(indexint,13,nlayer-i+1) =
	1058	$ jfotsout(indexint,13,nlayer) *
[705]	1059	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	1060	enddo
	1061	endif
	1062
	1063	c End of interval 35
	1064
	1065	ccccccccccccccccccccccccccccccc
	1066	c 337.8-800.0 nm (int 36)
	1067	c
	1068	c Absorption by:
	1069	c O3, NO2
	1070	ccccccccccccccccccccccccccccccc
	1071
	1072	c Input atmospheric column
	1073
	1074	indexint=36
[1266]	1075	do i=1,nlayer
	1076	auxcolinp(nlayer-i+1) = o3colx(i) + no2colx(i)
[705]	1077	end do
	1078
	1079	c Interpolation
	1080
	1081	do i=1,nz2
	1082	auxi = nz2-i+1
	1083	!O3 tabulated coefficient
	1084	auxjo3(i) = jabsifotsintpar(auxi,7,indexint)
	1085	!Tabulated column
	1086	auxcoltab(i) = c36(auxi)
	1087	enddo
	1088	!Only if chemthermod.ge.2
	1089	if(chemthermod.ge.2) then
	1090	do i=1,nz2
	1091	!NO2 tabulated coefficient
	1092	auxjno2(i) = jabsifotsintpar(nz2-i+1,13,indexint)
	1093	enddo
	1094	endif
	1095	call interfast
[1266]	1096	$ (wm,wp,auxind,auxcolinp,nlayer,auxcoltab,nz2,limdown,limup)
	1097	do i=1,nlayer
[705]	1098	ind=auxind(i)
	1099	!O3 interpolated coefficient
[1266]	1100	jfotsout(indexint,7,nlayer-i+1) =
	1101	$ jfotsout(indexint,7,nlayer) *
[705]	1102	$ (wm(i)auxjo3(ind+1) + wp(i)auxjo3(ind))
	1103	enddo
	1104	!Only if chemthermod.ge.2
	1105	if(chemthermod.ge.2) then
[1266]	1106	do i=1,nlayer
[705]	1107	ind=auxind(i)
	1108	!NO2 interpolated coefficient
[1266]	1109	jfotsout(indexint,13,nlayer-i+1) =
	1110	$ jfotsout(indexint,13,nlayer) *
[705]	1111	$ (wm(i)auxjno2(ind+1) + wp(i)auxjno2(ind))
	1112	enddo
	1113	endif
	1114
	1115	c End of interval 36
	1116
	1117	c End of interpolation to obtain photoabsorption rates
	1118
[2429]	1119	c Correction to the actual Sun-Mars distance
[705]	1120
[2429]	1121	jfotsout(:,:,:)=jfotsout(:,:,:)(1.52/dist_sol)*2
[705]	1122
	1123	end
	1124
	1125
	1126
	1127
	1128
	1129

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

Download in other formats:

Original Format