Home

Context Navigation

source: palm/trunk/LIB/rrtmg/rrtmg_lw_read_nc.f90 @ 1618

Last change on this file since 1618 was 1585, checked in by maronga, 10 years ago
Added support for RRTMG radiation code
File size: 43.3 KB

Rev	Line
[1585]	1	! path: $Source: /storm/rc1/cvsroot/rc/rrtmg_lw/src/rrtmg_lw_read_nc.f90,v $
	2	! author: $Author: mike $
	3	! revision: $Revision: 1.1 $
	4	! created: $Date: 2009/05/22 21:02:13 $
	5	!
	6
	7	!===============================================================================
	8	! rrtmg_lw_read_nc.f90
	9	!
	10	! Description: This program reads all of the RRTM longwave data from a NetCDF
	11	! file in band by band subroutines, as a replacement for the
	12	! rrtmg_lw_k_g.f90 data statements.
	13	!
	14	! Written By: Patrick Hofmann
	15	! Last Update: 1/23/2009
	16	!===============================================================================
	17
	18	!*******************************************************************************
	19	subroutine lw_kgb01
	20	use rrlw_kg01, only : fracrefao, fracrefbo, kao, kbo, kao_mn2, kbo_mn2, &
	21	selfrefo, forrefo, no1
	22	use rrlw_ncpar
	23	use netcdf
	24
	25	implicit none
	26	save
	27
	28	integer(kind=im) :: ab
	29	integer(kind=im), parameter :: bandNumber = 1, numGPoints = no1
	30	integer(kind=im), parameter :: gPointSetNumber = 1
	31	integer(kind=im) :: ncid, varID
	32
	33	status(:) = nf90_NoErr
	34	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	35
	36	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	37	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	38	start = (/1,1,bandNumber,gPointSetNumber/), &
	39	count = (/numGPoints,1,1,1/))
	40
	41	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	42	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	43	start = (/1,1,bandNumber,gPointSetNumber/), &
	44	count = (/numGPoints,1,1,1/))
	45
	46	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	47	status(7) = nf90_get_var(ncid, varID, kao, &
	48	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	49	count = (/1,Tdiff,plower,numGPoints,1,1/))
	50
	51	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	52	status(9) = nf90_get_var(ncid, varID, kbo, &
	53	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	54	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	55
	56	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	57	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	58	start = (/1,1,bandNumber,gPointSetNumber/), &
	59	count = (/Tself,numGPoints,1,1/))
	60
	61	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	62	status(13) = nf90_get_var(ncid, varID, forrefo, &
	63	start = (/1,1,bandNumber,gPointSetNumber/), &
	64	count = (/Tforeign,numGPoints,1,1/))
	65
	66	!Get absorber index for N2
	67	call getAbsorberIndex('N2',ab)
	68	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	69	status(15) = nf90_get_var(ncid, varID, kao_mn2, &
	70	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	71	count = (/1,T,numGPoints,1,1,1/))
	72
	73	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	74	status(17) = nf90_get_var(ncid, varID, kbo_mn2, &
	75	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	76	count = (/1,T,numGPoints,1,1,1/))
	77
	78	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	79
	80	status(1) = nf90_close(ncid)
	81
	82	end subroutine lw_kgb01
	83	!*******************************************************************************
	84
	85	!*******************************************************************************
	86	subroutine lw_kgb02
	87	use rrlw_kg02, only : fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, no2
	88	use rrlw_ncpar
	89	use netcdf
	90
	91	implicit none
	92	save
	93
	94	integer(kind=im), parameter :: bandNumber = 2
	95	integer(kind=im), parameter :: numGPoints = no2
	96	integer(kind=im), parameter :: gPointSetNumber = 1
	97	integer(kind=im) :: ncid, varID
	98
	99	status(:) = nf90_NoErr
	100	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	101
	102	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	103	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	104	start = (/1,1,bandNumber,gPointSetNumber/), &
	105	count = (/numGPoints,1,1,1/))
	106
	107	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	108	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	109	start = (/1,1,bandNumber,gPointSetNumber/), &
	110	count = (/numGPoints,1,1,1/))
	111
	112	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	113	status(7) = nf90_get_var(ncid, varID, kao, &
	114	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	115	count = (/1,Tdiff,plower,numGPoints,1,1/))
	116
	117	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	118	status(9) = nf90_get_var(ncid, varID, kbo, &
	119	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	120	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	121
	122	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	123	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	124	start = (/1,1,bandNumber,gPointSetNumber/), &
	125	count = (/Tself,numGPoints,1,1/))
	126
	127	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	128	status(13) = nf90_get_var(ncid, varID, forrefo, &
	129	start = (/1,1,bandNumber,gPointSetNumber/), &
	130	count = (/Tforeign,numGPoints,1,1/))
	131
	132	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	133
	134	status(1) = nf90_close(ncid)
	135
	136	end subroutine lw_kgb02
	137	!*******************************************************************************
	138
	139	!*******************************************************************************
	140	subroutine lw_kgb03
	141	use rrlw_kg03, only : fracrefao, fracrefbo, kao, kbo, kao_mn2o, &
	142	kbo_mn2o, selfrefo, forrefo, no3
	143	use rrlw_ncpar
	144	use netcdf
	145
	146	implicit none
	147	save
	148
	149	integer(kind=im) :: ab
	150	integer(kind=im), parameter :: bandNumber = 3
	151	integer(kind=im), parameter :: numGPoints = no3
	152	integer(kind=im), parameter :: gPointSetNumber = 1
	153	integer(kind=im) :: ncid, varID
	154
	155	status(:) = nf90_NoErr
	156	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	157
	158	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	159	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	160	start = (/1,1,bandNumber,gPointSetNumber/), &
	161	count = (/numGPoints,keylower,1,1/))
	162
	163	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	164	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	165	start = (/1,1,bandNumber,gPointSetNumber/), &
	166	count = (/numGPoints,keyupper,1,1/))
	167
	168	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	169	status(7) = nf90_get_var(ncid, varID, kao, &
	170	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	171	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	172
	173	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	174	status(9) = nf90_get_var(ncid, varID, kbo, &
	175	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	176	count = (/keyupper,Tdiff,pupper,numGPoints,1,1/))
	177
	178	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	179	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	180	start = (/1,1,bandNumber,gPointSetNumber/), &
	181	count = (/Tself,numGPoints,1,1/))
	182
	183	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	184	status(13) = nf90_get_var(ncid, varID, forrefo, &
	185	start = (/1,1,bandNumber,gPointSetNumber/), &
	186	count = (/Tforeign,numGPoints,1,1/))
	187
	188	!Get absorber index for N2
	189	call getAbsorberIndex('N2O',ab)
	190	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	191	status(15) = nf90_get_var(ncid, varID, kao_mn2o, &
	192	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	193	count = (/keylower,T,numGPoints,1,1,1/))
	194
	195	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	196	status(17) = nf90_get_var(ncid, varID, kbo_mn2o, &
	197	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	198	count = (/keyupper,T,numGPoints,1,1,1/))
	199
	200	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	201
	202	status(1) = nf90_close(ncid)
	203
	204	end subroutine lw_kgb03
	205	!*******************************************************************************
	206
	207	!*******************************************************************************
	208	subroutine lw_kgb04
	209	use rrlw_kg04, only : fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, no4
	210	use rrlw_ncpar
	211	use netcdf
	212
	213	implicit none
	214	save
	215
	216	integer(kind=im), parameter :: bandNumber = 4
	217	integer(kind=im), parameter :: numGPoints = no4
	218	integer(kind=im), parameter :: gPointSetNumber = 1
	219	integer(kind=im) :: ncid, varID
	220
	221	status(:) = nf90_NoErr
	222	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	223
	224	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	225	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	226	start = (/1,1,bandNumber,gPointSetNumber/), &
	227	count = (/numGPoints,keylower,1,1/))
	228
	229	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	230	status(5) = nf90_get_var(ncid, varID, fracrefbo(:,1:5), &
	231	start = (/1,1,bandNumber,gPointSetNumber/), &
	232	count = (/numGPoints,keyupper,1,1/))
	233
	234	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	235	status(7) = nf90_get_var(ncid, varID, kao, &
	236	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	237	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	238
	239	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	240	status(9) = nf90_get_var(ncid, varID, kbo, &
	241	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	242	count = (/keyupper,Tdiff,pupper,numGPoints,1,1/))
	243
	244	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	245	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	246	start = (/1,1,bandNumber,gPointSetNumber/), &
	247	count = (/Tself,numGPoints,1,1/))
	248
	249	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	250	status(13) = nf90_get_var(ncid, varID, forrefo, &
	251	start = (/1,1,bandNumber,gPointSetNumber/), &
	252	count = (/Tforeign,numGPoints,1,1/))
	253
	254	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	255
	256	status(1) = nf90_close(ncid)
	257
	258	end subroutine lw_kgb04
	259	!*******************************************************************************
	260
	261	!*******************************************************************************
	262	subroutine lw_kgb05
	263	use rrlw_kg05, only : fracrefao, fracrefbo, kao, kbo, kao_mo3, &
	264	selfrefo, forrefo, ccl4o, no5
	265	use rrlw_ncpar
	266	use netcdf
	267
	268	implicit none
	269	save
	270
	271	integer(kind=im) :: ab
	272	integer(kind=im), parameter :: bandNumber = 5
	273	integer(kind=im), parameter :: numGPoints = no5
	274	integer(kind=im), parameter :: gPointSetNumber = 1
	275	integer(kind=im) :: ncid, varID
	276
	277	status(:) = nf90_NoErr
	278	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	279
	280	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	281	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	282	start = (/1,1,bandNumber,gPointSetNumber/), &
	283	count = (/numGPoints,keylower,1,1/))
	284
	285	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	286	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	287	start = (/1,1,bandNumber,gPointSetNumber/), &
	288	count = (/numGPoints,keyupper,1,1/))
	289
	290	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	291	status(7) = nf90_get_var(ncid, varID, kao, &
	292	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	293	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	294
	295	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	296	status(9) = nf90_get_var(ncid, varID, kbo, &
	297	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	298	count = (/keyupper,Tdiff,pupper,numGPoints,1,1/))
	299
	300	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	301	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	302	start = (/1,1,bandNumber,gPointSetNumber/), &
	303	count = (/Tself,numGPoints,1,1/))
	304
	305	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	306	status(13) = nf90_get_var(ncid, varID, forrefo, &
	307	start = (/1,1,bandNumber,gPointSetNumber/), &
	308	count = (/Tforeign,numGPoints,1,1/))
	309
	310	!Get absorber index for O3
	311	call getAbsorberIndex('O3',ab)
	312	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	313	status(15) = nf90_get_var(ncid, varID, kao_mo3, &
	314	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	315	count = (/keylower,T,numGPoints,1,1,1/))
	316
	317	!Get absorber index for CCL4
	318	call getAbsorberIndex('CCL4',ab)
	319	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	320	status(17) = nf90_get_var(ncid, varID, ccl4o, &
	321	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	322	count = (/1,1,numGPoints,1,1,1/))
	323
	324	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	325
	326	status(1) = nf90_close(ncid)
	327
	328	end subroutine lw_kgb05
	329	!*******************************************************************************
	330
	331	!*******************************************************************************
	332	subroutine lw_kgb06
	333	use rrlw_kg06, only : fracrefao, kao, kao_mco2, selfrefo, forrefo, &
	334	cfc11adjo, cfc12o, no6
	335	use rrlw_ncpar
	336	use netcdf
	337
	338	implicit none
	339	save
	340
	341	integer(kind=im) :: ab
	342	integer(kind=im), parameter :: bandNumber = 6
	343	integer(kind=im), parameter :: numGPoints = no6
	344	integer(kind=im), parameter :: gPointSetNumber = 1
	345	integer(kind=im) :: ncid, varID
	346
	347	status(:) = nf90_NoErr
	348	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	349
	350	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	351	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	352	start = (/1,1,bandNumber,gPointSetNumber/), &
	353	count = (/numGPoints,1,1,1/))
	354
	355	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	356	status(7) = nf90_get_var(ncid, varID, kao, &
	357	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	358	count = (/1,Tdiff,plower,numGPoints,1,1/))
	359
	360	status(8) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	361	status(9) = nf90_get_var(ncid, varID, selfrefo, &
	362	start = (/1,1,bandNumber,gPointSetNumber/), &
	363	count = (/Tself,numGPoints,1,1/))
	364
	365	status(10) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	366	status(11) = nf90_get_var(ncid, varID, forrefo, &
	367	start = (/1,1,bandNumber,gPointSetNumber/), &
	368	count = (/Tforeign,numGPoints,1,1/))
	369
	370	!Get absorber index for CO2
	371	call getAbsorberIndex('CO2',ab)
	372	status(12) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	373	status(13) = nf90_get_var(ncid, varID, kao_mco2, &
	374	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	375	count = (/1,T,numGPoints,1,1,1/))
	376
	377	!Get absorber index for CFC11
	378	call getAbsorberIndex('CFC11',ab)
	379	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	380	status(15) = nf90_get_var(ncid, varID, cfc11adjo, &
	381	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	382	count = (/1,1,numGPoints,1,1,1/))
	383
	384	!Get absorber index for CFC12
	385	call getAbsorberIndex('CFC12',ab)
	386	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	387	status(17) = nf90_get_var(ncid, varID, cfc12o, &
	388	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	389	count = (/1,1,numGPoints,1,1,1/))
	390
	391	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	392
	393	status(1) = nf90_close(ncid)
	394
	395	end subroutine lw_kgb06
	396	!*******************************************************************************
	397
	398	!*******************************************************************************
	399	subroutine lw_kgb07
	400	use rrlw_kg07, only : fracrefao, fracrefbo, kao, kbo, kao_mco2, &
	401	kbo_mco2, selfrefo, forrefo, no7
	402	use rrlw_ncpar
	403	use netcdf
	404
	405	implicit none
	406	save
	407
	408	integer(kind=im) :: ab
	409	integer(kind=im), parameter :: bandNumber = 7
	410	integer(kind=im), parameter :: numGPoints = no7
	411	integer(kind=im), parameter :: gPointSetNumber = 1
	412	integer(kind=im) :: ncid, varID
	413
	414	status(:) = nf90_NoErr
	415	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	416
	417	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	418	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	419	start = (/1,1,bandNumber,gPointSetNumber/), &
	420	count = (/numGPoints,keylower,1,1/))
	421
	422	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	423	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	424	start = (/1,1,bandNumber,gPointSetNumber/), &
	425	count = (/numGPoints,1,1,1/))
	426
	427	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	428	status(7) = nf90_get_var(ncid, varID, kao, &
	429	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	430	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	431
	432	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	433	status(9) = nf90_get_var(ncid, varID, kbo, &
	434	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	435	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	436
	437	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	438	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	439	start = (/1,1,bandNumber,gPointSetNumber/), &
	440	count = (/Tself,numGPoints,1,1/))
	441
	442	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	443	status(13) = nf90_get_var(ncid, varID, forrefo, &
	444	start = (/1,1,bandNumber,gPointSetNumber/), &
	445	count = (/Tforeign,numGPoints,1,1/))
	446
	447	!Get absorber index for CO2
	448	call getAbsorberIndex('CO2',ab)
	449	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	450	status(15) = nf90_get_var(ncid, varID, kao_mco2, &
	451	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	452	count = (/keylower,T,numGPoints,1,1,1/))
	453
	454	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	455	status(17) = nf90_get_var(ncid, varID, kbo_mco2, &
	456	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	457	count = (/1,T,numGPoints,1,1,1/))
	458
	459	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	460
	461	status(1) = nf90_close(ncid)
	462
	463	end subroutine lw_kgb07
	464	!*******************************************************************************
	465
	466	!*******************************************************************************
	467	subroutine lw_kgb08
	468	use rrlw_kg08, only : fracrefao, fracrefbo, kao, kao_mco2, kao_mn2o, &
	469	kao_mo3, kbo, kbo_mco2, kbo_mn2o, selfrefo, forrefo, &
	470	cfc12o, cfc22adjo, no8
	471	use rrlw_ncpar
	472	use netcdf
	473
	474	implicit none
	475	save
	476
	477	integer(kind=im) :: ab
	478	integer(kind=im), parameter :: bandNumber = 8
	479	integer(kind=im), parameter :: numGPoints = no8
	480	integer(kind=im), parameter :: gPointSetNumber = 1
	481	integer(kind=im) :: ncid, varID
	482
	483	status(:) = nf90_NoErr
	484	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	485
	486	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	487	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	488	start = (/1,1,bandNumber,gPointSetNumber/), &
	489	count = (/numGPoints,1,1,1/))
	490
	491	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	492	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	493	start = (/1,1,bandNumber,gPointSetNumber/), &
	494	count = (/numGPoints,1,1,1/))
	495
	496	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	497	status(7) = nf90_get_var(ncid, varID, kao, &
	498	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	499	count = (/1,Tdiff,plower,numGPoints,1,1/))
	500
	501	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	502	status(9) = nf90_get_var(ncid, varID, kbo, &
	503	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	504	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	505
	506	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	507	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	508	start = (/1,1,bandNumber,gPointSetNumber/), &
	509	count = (/Tself,numGPoints,1,1/))
	510
	511	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	512	status(13) = nf90_get_var(ncid, varID, forrefo, &
	513	start = (/1,1,bandNumber,gPointSetNumber/), &
	514	count = (/Tforeign,numGPoints,1,1/))
	515
	516	!Get absorber index for O3
	517	call getAbsorberIndex('O3',ab)
	518	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	519	status(15) = nf90_get_var(ncid, varID, kao_mo3, &
	520	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	521	count = (/1,T,numGPoints,1,1,1/))
	522
	523	!Get absorber index for CO2
	524	call getAbsorberIndex('CO2',ab)
	525	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	526	status(17) = nf90_get_var(ncid, varID, kao_mco2, &
	527	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	528	count = (/1,T,numGPoints,1,1,1/))
	529
	530	status(18) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	531	status(19) = nf90_get_var(ncid, varID, kbo_mco2, &
	532	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	533	count = (/1,T,numGPoints,1,1,1/))
	534
	535	!Get absorber index for N2O
	536	call getAbsorberIndex('N2O',ab)
	537	status(20) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	538	status(21) = nf90_get_var(ncid, varID, kao_mn2o, &
	539	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	540	count = (/1,T,numGPoints,1,1,1/))
	541
	542	status(22) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	543	status(23) = nf90_get_var(ncid, varID, kbo_mn2o, &
	544	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	545	count = (/1,T,numGPoints,1,1,1/))
	546
	547	!Get absorber index for CFC12
	548	call getAbsorberIndex('CFC12',ab)
	549	status(24) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	550	status(25) = nf90_get_var(ncid, varID, cfc12o, &
	551	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	552	count = (/1,1,numGPoints,1,1,1/))
	553
	554	!Get absorber index for CFC22
	555	call getAbsorberIndex('CFC22',ab)
	556	status(26) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	557	status(27) = nf90_get_var(ncid, varID, cfc22adjo, &
	558	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	559	count = (/1,1,numGPoints,1,1,1/))
	560
	561	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	562
	563	status(1) = nf90_close(ncid)
	564
	565	end subroutine lw_kgb08
	566	!*******************************************************************************
	567
	568	!*******************************************************************************
	569	subroutine lw_kgb09
	570	use rrlw_kg09, only : fracrefao, fracrefbo, kao, kbo, kao_mn2o, &
	571	kbo_mn2o, selfrefo, forrefo, no9
	572	use rrlw_ncpar
	573	use netcdf
	574
	575	implicit none
	576	save
	577
	578	integer(kind=im) :: ab
	579	integer(kind=im), parameter :: bandNumber = 9
	580	integer(kind=im), parameter :: numGPoints = no9
	581	integer(kind=im), parameter :: gPointSetNumber = 1
	582	integer(kind=im) :: ncid, varID
	583
	584	status(:) = nf90_NoErr
	585	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	586
	587	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	588	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	589	start = (/1,1,bandNumber,gPointSetNumber/), &
	590	count = (/numGPoints,keylower,1,1/))
	591
	592	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	593	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	594	start = (/1,1,bandNumber,gPointSetNumber/), &
	595	count = (/numGPoints,1,1,1/))
	596
	597	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	598	status(7) = nf90_get_var(ncid, varID, kao, &
	599	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	600	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	601
	602	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	603	status(9) = nf90_get_var(ncid, varID, kbo, &
	604	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	605	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	606
	607	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	608	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	609	start = (/1,1,bandNumber,gPointSetNumber/), &
	610	count = (/Tself,numGPoints,1,1/))
	611
	612	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	613	status(13) = nf90_get_var(ncid, varID, forrefo, &
	614	start = (/1,1,bandNumber,gPointSetNumber/), &
	615	count = (/Tforeign,numGPoints,1,1/))
	616
	617	!Get absorber index for N2O
	618	call getAbsorberIndex('N2O',ab)
	619	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	620	status(15) = nf90_get_var(ncid, varID, kao_mn2o, &
	621	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	622	count = (/keylower,T,numGPoints,1,1,1/))
	623
	624	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	625	status(17) = nf90_get_var(ncid, varID, kbo_mn2o, &
	626	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	627	count = (/1,T,numGPoints,1,1,1/))
	628
	629	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	630
	631	status(1) = nf90_close(ncid)
	632
	633	end subroutine lw_kgb09
	634	!*******************************************************************************
	635
	636	!*******************************************************************************
	637	subroutine lw_kgb10
	638	use rrlw_kg10, only : fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, no10
	639	use rrlw_ncpar
	640	use netcdf
	641
	642	implicit none
	643	save
	644
	645	integer(kind=im), parameter :: bandNumber = 10
	646	integer(kind=im), parameter :: numGPoints = no10
	647	integer(kind=im), parameter :: gPointSetNumber = 1
	648	integer(kind=im) :: ncid, varID
	649
	650	status(:) = nf90_NoErr
	651	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	652
	653	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	654	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	655	start = (/1,1,bandNumber,gPointSetNumber/), &
	656	count = (/numGPoints,1,1,1/))
	657
	658	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	659	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	660	start = (/1,1,bandNumber,gPointSetNumber/), &
	661	count = (/numGPoints,1,1,1/))
	662
	663	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	664	status(7) = nf90_get_var(ncid, varID, kao, &
	665	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	666	count = (/1,Tdiff,plower,numGPoints,1,1/))
	667
	668	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	669	status(9) = nf90_get_var(ncid, varID, kbo, &
	670	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	671	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	672
	673	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	674	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	675	start = (/1,1,bandNumber,gPointSetNumber/), &
	676	count = (/Tself,numGPoints,1,1/))
	677
	678	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	679	status(13) = nf90_get_var(ncid, varID, forrefo, &
	680	start = (/1,1,bandNumber,gPointSetNumber/), &
	681	count = (/Tforeign,numGPoints,1,1/))
	682
	683	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	684
	685	status(1) = nf90_close(ncid)
	686
	687	end subroutine lw_kgb10
	688	!*******************************************************************************
	689
	690	!*******************************************************************************
	691	subroutine lw_kgb11
	692	use rrlw_kg11, only : fracrefao, fracrefbo, kao, kbo, kao_mo2, &
	693	kbo_mo2, selfrefo, forrefo, no11
	694	use rrlw_ncpar
	695	use netcdf
	696
	697	implicit none
	698	save
	699
	700	integer(kind=im) :: ab
	701	integer(kind=im), parameter :: bandNumber = 11
	702	integer(kind=im), parameter :: numGPoints = no11
	703	integer(kind=im), parameter :: gPointSetNumber = 1
	704	integer(kind=im) :: ncid, varID
	705
	706	status(:) = nf90_NoErr
	707	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	708
	709	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	710	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	711	start = (/1,1,bandNumber,gPointSetNumber/), &
	712	count = (/numGPoints,1,1,1/))
	713
	714	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	715	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	716	start = (/1,1,bandNumber,gPointSetNumber/), &
	717	count = (/numGPoints,1,1,1/))
	718
	719	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	720	status(7) = nf90_get_var(ncid, varID, kao, &
	721	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	722	count = (/1,Tdiff,plower,numGPoints,1,1/))
	723
	724	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	725	status(9) = nf90_get_var(ncid, varID, kbo, &
	726	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	727	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	728
	729	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	730	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	731	start = (/1,1,bandNumber,gPointSetNumber/), &
	732	count = (/Tself,numGPoints,1,1/))
	733
	734	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	735	status(13) = nf90_get_var(ncid, varID, forrefo, &
	736	start = (/1,1,bandNumber,gPointSetNumber/), &
	737	count = (/Tforeign,numGPoints,1,1/))
	738
	739	!Get absorber index for O2
	740	call getAbsorberIndex('O2',ab)
	741	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	742	status(15) = nf90_get_var(ncid, varID, kao_mo2, &
	743	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	744	count = (/1,T,numGPoints,1,1,1/))
	745
	746	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	747	status(17) = nf90_get_var(ncid, varID, kbo_mo2, &
	748	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	749	count = (/1,T,numGPoints,1,1,1/))
	750
	751	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	752
	753	status(1) = nf90_close(ncid)
	754
	755	end subroutine lw_kgb11
	756	!*******************************************************************************
	757
	758	!*******************************************************************************
	759	subroutine lw_kgb12
	760	use rrlw_kg12, only : fracrefao, kao, selfrefo, forrefo, no12
	761	use rrlw_ncpar
	762	use netcdf
	763
	764	implicit none
	765	save
	766
	767	integer(kind=im), parameter :: bandNumber = 12
	768	integer(kind=im), parameter :: numGPoints = no12
	769	integer(kind=im), parameter :: gPointSetNumber = 1
	770	integer(kind=im) :: ncid, varID
	771
	772	status(:) = nf90_NoErr
	773	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	774
	775	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	776	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	777	start = (/1,1,bandNumber,gPointSetNumber/), &
	778	count = (/numGPoints,keylower,1,1/))
	779
	780	status(4) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	781	status(5) = nf90_get_var(ncid, varID, kao, &
	782	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	783	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	784
	785	status(6) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	786	status(7) = nf90_get_var(ncid, varID, selfrefo, &
	787	start = (/1,1,bandNumber,gPointSetNumber/), &
	788	count = (/Tself,numGPoints,1,1/))
	789
	790	status(8) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	791	status(9) = nf90_get_var(ncid, varID, forrefo, &
	792	start = (/1,1,bandNumber,gPointSetNumber/), &
	793	count = (/Tforeign,numGPoints,1,1/))
	794
	795	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	796
	797	status(1) = nf90_close(ncid)
	798
	799	end subroutine lw_kgb12
	800	!*******************************************************************************
	801
	802	!*******************************************************************************
	803	subroutine lw_kgb13
	804	use rrlw_kg13, only : fracrefao, fracrefbo, kao, kao_mco2, kao_mco, &
	805	kbo_mo3, selfrefo, forrefo, no13
	806	use rrlw_ncpar
	807	use netcdf
	808
	809	implicit none
	810	save
	811
	812	integer(kind=im) :: ab
	813	integer(kind=im), parameter :: bandNumber = 13
	814	integer(kind=im), parameter :: numGPoints = no13
	815	integer(kind=im), parameter :: gPointSetNumber = 1
	816	integer(kind=im) :: ncid, varID
	817
	818	status(:) = nf90_NoErr
	819	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	820
	821	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	822	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	823	start = (/1,1,bandNumber,gPointSetNumber/), &
	824	count = (/numGPoints,keylower,1,1/))
	825
	826	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	827	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	828	start = (/1,1,bandNumber,gPointSetNumber/), &
	829	count = (/numGPoints,1,1,1/))
	830
	831	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	832	status(7) = nf90_get_var(ncid, varID, kao, &
	833	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	834	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	835
	836	status(8) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	837	status(9) = nf90_get_var(ncid, varID, selfrefo, &
	838	start = (/1,1,bandNumber,gPointSetNumber/), &
	839	count = (/Tself,numGPoints,1,1/))
	840
	841	status(10) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	842	status(11) = nf90_get_var(ncid, varID, forrefo, &
	843	start = (/1,1,bandNumber,gPointSetNumber/), &
	844	count = (/Tforeign,numGPoints,1,1/))
	845
	846	!Get absorber index for O3
	847	call getAbsorberIndex('O3',ab)
	848	status(12) = nf90_inq_varid(ncid,"AbsorptionCoefficientsUpperAtmos",varID)
	849	status(13) = nf90_get_var(ncid, varID, kbo_mo3, &
	850	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	851	count = (/1,T,numGPoints,1,1,1/))
	852
	853	!Get absorber index for CO2
	854	call getAbsorberIndex('CO2',ab)
	855	status(14) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	856	status(15) = nf90_get_var(ncid, varID, kao_mco2, &
	857	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	858	count = (/keylower,T,numGPoints,1,1,1/))
	859
	860	!Get absorber index for CO
	861	call getAbsorberIndex('CO',ab)
	862	status(16) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	863	status(17) = nf90_get_var(ncid, varID, kao_mco, &
	864	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	865	count = (/keylower,T,numGPoints,1,1,1/))
	866
	867	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	868
	869	status(1) = nf90_close(ncid)
	870
	871	end subroutine lw_kgb13
	872	!*******************************************************************************
	873
	874	!*******************************************************************************
	875	subroutine lw_kgb14
	876	use rrlw_kg14, only : fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, no14
	877	use rrlw_ncpar
	878	use netcdf
	879
	880	implicit none
	881	save
	882
	883	integer(kind=im), parameter :: bandNumber = 14
	884	integer(kind=im), parameter :: numGPoints = no14
	885	integer(kind=im), parameter :: gPointSetNumber = 1
	886	integer(kind=im) :: ncid, varID
	887
	888	status(:) = nf90_NoErr
	889	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	890
	891	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	892	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	893	start = (/1,1,bandNumber,gPointSetNumber/), &
	894	count = (/numGPoints,1,1,1/))
	895
	896	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	897	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	898	start = (/1,1,bandNumber,gPointSetNumber/), &
	899	count = (/numGPoints,1,1,1/))
	900
	901	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	902	status(7) = nf90_get_var(ncid, varID, kao, &
	903	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	904	count = (/1,Tdiff,plower,numGPoints,1,1/))
	905
	906	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	907	status(9) = nf90_get_var(ncid, varID, kbo, &
	908	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	909	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	910
	911	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	912	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	913	start = (/1,1,bandNumber,gPointSetNumber/), &
	914	count = (/Tself,numGPoints,1,1/))
	915
	916	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	917	status(13) = nf90_get_var(ncid, varID, forrefo, &
	918	start = (/1,1,bandNumber,gPointSetNumber/), &
	919	count = (/Tforeign,numGPoints,1,1/))
	920
	921	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	922
	923	status(1) = nf90_close(ncid)
	924
	925	end subroutine lw_kgb14
	926	!*******************************************************************************
	927
	928	!*******************************************************************************
	929	subroutine lw_kgb15
	930	use rrlw_kg15, only : fracrefao, kao, kao_mn2, selfrefo, forrefo, no15
	931	use rrlw_ncpar
	932	use netcdf
	933
	934	implicit none
	935	save
	936
	937	integer(kind=im) :: ab
	938	integer(kind=im), parameter :: bandNumber = 15
	939	integer(kind=im), parameter :: numGPoints = no15
	940	integer(kind=im), parameter :: gPointSetNumber = 1
	941	integer(kind=im) :: ncid, varID
	942
	943	status(:) = nf90_NoErr
	944	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	945
	946	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	947	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	948	start = (/1,1,bandNumber,gPointSetNumber/), &
	949	count = (/numGPoints,keylower,1,1/))
	950
	951	status(4) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	952	status(5) = nf90_get_var(ncid, varID, kao, &
	953	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	954	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	955
	956	status(6) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	957	status(7) = nf90_get_var(ncid, varID, selfrefo, &
	958	start = (/1,1,bandNumber,gPointSetNumber/), &
	959	count = (/Tself,numGPoints,1,1/))
	960
	961	status(8) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	962	status(9) = nf90_get_var(ncid, varID, forrefo, &
	963	start = (/1,1,bandNumber,gPointSetNumber/), &
	964	count = (/Tforeign,numGPoints,1,1/))
	965
	966	!Get absorber index for N2
	967	call getAbsorberIndex('N2',ab)
	968	status(10) = nf90_inq_varid(ncid,"AbsorptionCoefficientsLowerAtmos",varID)
	969	status(11) = nf90_get_var(ncid, varID, kao_mn2, &
	970	start = (/1,1,1,ab,bandNumber,gPointSetNumber/), &
	971	count = (/keylower,T,numGPoints,1,1,1/))
	972
	973	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	974
	975	status(1) = nf90_close(ncid)
	976
	977	end subroutine lw_kgb15
	978	!*******************************************************************************
	979
	980	!*******************************************************************************
	981	subroutine lw_kgb16
	982	use rrlw_kg16, only : fracrefao, fracrefbo, kao, kbo, selfrefo, forrefo, no16
	983	use rrlw_ncpar
	984	use netcdf
	985
	986	implicit none
	987	save
	988
	989	integer(kind=im), parameter :: bandNumber = 16
	990	integer(kind=im), parameter :: numGPoints = no16
	991	integer(kind=im), parameter :: gPointSetNumber = 1
	992	integer(kind=im) :: ncid, varID
	993
	994	status(:) = nf90_NoErr
	995	status(1) = nf90_open('rrtmg_lw.nc',nf90_nowrite,ncid)
	996
	997	status(2) = nf90_inq_varid(ncid,"PlanckFractionLowerAtmos",varID)
	998	status(3) = nf90_get_var(ncid, varID, fracrefao, &
	999	start = (/1,1,bandNumber,gPointSetNumber/), &
	1000	count = (/numGPoints,keylower,1,1/))
	1001
	1002	status(4) = nf90_inq_varid(ncid,"PlanckFractionUpperAtmos",varID)
	1003	status(5) = nf90_get_var(ncid, varID, fracrefbo, &
	1004	start = (/1,1,bandNumber,gPointSetNumber/), &
	1005	count = (/numGPoints,1,1,1/))
	1006
	1007	status(6) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsLowerAtmos",varID)
	1008	status(7) = nf90_get_var(ncid, varID, kao, &
	1009	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	1010	count = (/keylower,Tdiff,plower,numGPoints,1,1/))
	1011
	1012	status(8) = nf90_inq_varid(ncid,"KeySpeciesAbsorptionCoefficientsUpperAtmos",varID)
	1013	status(9) = nf90_get_var(ncid, varID, kbo, &
	1014	start = (/1,1,1,1,bandNumber,gPointSetNumber/), &
	1015	count = (/1,Tdiff,pupper,numGPoints,1,1/))
	1016
	1017	status(10) = nf90_inq_varid(ncid,"H20SelfAbsorptionCoefficients",varID)
	1018	status(11) = nf90_get_var(ncid, varID, selfrefo, &
	1019	start = (/1,1,bandNumber,gPointSetNumber/), &
	1020	count = (/Tself,numGPoints,1,1/))
	1021
	1022	status(12) = nf90_inq_varid(ncid,"H20ForeignAbsorptionCoefficients",varID)
	1023	status(13) = nf90_get_var(ncid, varID, forrefo, &
	1024	start = (/1,1,bandNumber,gPointSetNumber/), &
	1025	count = (/Tforeign,numGPoints,1,1/))
	1026
	1027	if(any(status(:) /= nf90_NoErr)) stop "Error reading variables from file"
	1028
	1029	status(1) = nf90_close(ncid)
	1030
	1031	end subroutine lw_kgb16
	1032	!*******************************************************************************

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |