Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

cosp.f90 @ 5440

Last change on this file since 5440 was 5268, checked in by abarral, 2 months ago
.f90 <-> .F90 depending on cpp key use
File size: 215.5 KB

Rev	Line
[3491]	1	! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	2	! Copyright (c) 2015, Regents of the University of Colorado
	3	! All rights reserved.
	4	!
	5	! Redistribution and use in source and binary forms, with or without modification, are
	6	! permitted provided that the following conditions are met:
	7	!
	8	! 1. Redistributions of source code must retain the above copyright notice, this list of
	9	! conditions and the following disclaimer.
	10	!
	11	! 2. Redistributions in binary form must reproduce the above copyright notice, this list
	12	! of conditions and the following disclaimer in the documentation and/or other
	13	! materials provided with the distribution.
	14	!
	15	! 3. Neither the name of the copyright holder nor the names of its contributors may be
	16	! used to endorse or promote products derived from this software without specific prior
	17	! written permission.
	18	!
	19	! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
	20	! EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
	21	! MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL
	22	! THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
	23	! SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
	24	! OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
	25	! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
	26	! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
	27	! OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
	28	!
	29	! History:
	30	! May 2015- D. Swales - Original version
	31	! Mar 2018- R. Guzman - Added OPAQ diagnostics and GLID simulator
	32	! Apr 2018- R. Guzman - Added ATLID simulator
	33	!
	34	! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	35
	36
	37	MODULE MOD_COSP
	38	USE COSP_KINDS, ONLY: wp
	39	USE MOD_COSP_CONFIG, ONLY: R_UNDEF,PARASOL_NREFL,LIDAR_NCAT,LIDAR_NTYPE, SR_BINS,&
	40	N_HYDRO,RTTOV_MAX_CHANNELS,numMISRHgtBins, &
	41	cloudsat_DBZE_BINS,LIDAR_NTEMP,calipso_histBsct,&
	42	use_vgrid,Nlvgrid,vgrid_zu,vgrid_zl,vgrid_z, &
	43	numMODISTauBins,numMODISPresBins, &
	44	numMODISReffIceBins,numMODISReffLiqBins, &
	45	numISCCPTauBins,numISCCPPresBins,numMISRTauBins,&
	46	ntau,modis_histTau,tau_binBounds, &
	47	modis_histTauEdges,tau_binEdges,nCloudsatPrecipClass,&
	48	modis_histTauCenters,tau_binCenters, &
	49	! cloudsat_preclvl,grLidar532_histBsct,atlid_histBsct !PREC_BUG
	50	grLidar532_histBsct,atlid_histBsct
	51	USE MOD_COSP_MODIS_INTERFACE, ONLY: cosp_modis_init, modis_IN
	52	! USE MOD_COSP_RTTOV_INTERFACE, ONLY: cosp_rttov_init, rttov_IN
	53	USE MOD_COSP_MISR_INTERFACE, ONLY: cosp_misr_init, misr_IN
	54	USE MOD_COSP_ISCCP_INTERFACE, ONLY: cosp_isccp_init, isccp_IN
	55	USE MOD_COSP_CALIPSO_INTERFACE, ONLY: cosp_calipso_init, calipso_IN
	56	USE MOD_COSP_ATLID_INTERFACE, ONLY: cosp_atlid_init, atlid_IN
	57	USE MOD_COSP_GRLIDAR532_INTERFACE, ONLY: cosp_grLidar532_init, grLidar532_IN
	58	USE MOD_COSP_PARASOL_INTERFACE, ONLY: cosp_parasol_init, parasol_in
	59	USE MOD_COSP_CLOUDSAT_INTERFACE, ONLY: cosp_cloudsat_init, cloudsat_IN
	60	USE quickbeam, ONLY: quickbeam_subcolumn, quickbeam_column, radar_cfg
	61	USE MOD_ICARUS, ONLY: icarus_subcolumn, icarus_column
	62	USE MOD_MISR_SIMULATOR, ONLY: misr_subcolumn, misr_column
	63	USE MOD_LIDAR_SIMULATOR, ONLY: lidar_subcolumn, lidar_column
	64	USE MOD_MODIS_SIM, ONLY: modis_subcolumn, modis_column
	65	USE MOD_PARASOL, ONLY: parasol_subcolumn, parasol_column
	66	! use mod_cosp_rttov, ONLY: rttov_column
	67	USE MOD_COSP_STATS, ONLY: COSP_LIDAR_ONLY_CLOUD,COSP_CHANGE_VERTICAL_GRID
	68
	69	IMPLICIT NONE
	70
	71	logical :: linitialization ! Initialization flag
	72
	73	! ######################################################################################
	74	! TYPE cosp_column_inputs
	75	! ######################################################################################
	76	type cosp_column_inputs
	77	integer :: &
	78	Npoints, & ! Number of gridpoints.
	79	Ncolumns, & ! Number of columns.
	80	Nlevels ! Number of levels.
	81
	82	integer,allocatable,dimension(:) :: &
	83	sunlit ! Sunlit flag (0-1)
	84
	85	real(wp),allocatable,dimension(:,:) :: &
	86	at, & ! Temperature (K)
	87	pfull, & ! Pressure (Pa)
	88	phalf, & ! Pressure at half-levels (Pa)
	89	qv, & ! Specific humidity (kg/kg)
	90	hgt_matrix, & ! Height of hydrometeors (km)
	91	hgt_matrix_half ! Height of hydrometeors at half levels (km)
	92
	93	real(wp),allocatable,dimension(:) :: &
	94	land, & ! Land/Sea mask (0-1)
	95	skt, & ! Surface temperature (K)
	96	surfelev ! Surface Elevation (m)
	97	! Fields used ONLY by RTTOV
	98	integer :: &
	99	month ! Month for surface emissivty atlas (1-12)
	100	real(wp) :: &
	101	zenang, & ! Satellite zenith angle for RTTOV (deg)
	102	co2, & ! CO2 (kg/kg)
	103	ch4, & ! Methane (kg/kg)
	104	n2o, & ! N2O (kg/kg)
	105	co ! CO (kg/kg)
	106	real(wp),allocatable,dimension(:) :: &
	107	emis_sfc, & ! Surface emissivity (1)
	108	u_sfc, & ! Surface u-wind (m/s)
	109	v_sfc, & ! Surface v-wind (m/s)
	110	seaice, & ! Sea-ice fraction (0-1)
	111	lat, & ! Latitude (deg)
	112	lon ! Longitude (deg)
	113	real(wp),allocatable,dimension(:,:) :: &
	114	o3, & ! Ozone (kg/kg)
	115	tca, & ! Total column cloud fraction (0-1)
	116	cloudIce, & ! Cloud ice water mixing ratio (kg/kg)
	117	cloudLiq, & ! Cloud liquid water mixing ratio (kg/kg)
	118	fl_rain, & ! Precipitation (rain) flux (kg/m2/s)
	119	fl_snow ! Precipitation (snow) flux (kg/m2/s)
	120	end type cosp_column_inputs
	121
	122	! ######################################################################################
	123	! TYPE cosp_optical_inputs
	124	! ######################################################################################
	125	type cosp_optical_inputs
	126	integer :: &
	127	Npoints, & ! Number of gridpoints.
	128	Ncolumns, & ! Number of columns.
	129	Nlevels, & ! Number of levels.
	130	Npart, & ! Number of cloud meteors for LIDAR simulators.
	131	Nrefl ! Number of reflectances for PARASOL simulator
	132	real(wp) :: &
	133	emsfc_lw ! Surface emissivity @ 11micron
	134	real(wp),allocatable,dimension(:,:,:) :: &
	135	frac_out, & ! Cloud fraction
	136	tau_067, & ! Optical depth @ 0.67micron
	137	emiss_11, & ! Emissivity @ 11 micron
	138	fracLiq, & ! Fraction of optical-depth due to liquid (MODIS)
	139	asym, & ! Assymetry parameter @ 3.7micron (MODIS)
	140	ss_alb, & ! Single-scattering albedo @ 3.7micron (MODIS)
	141	betatot_calipso, & ! Lidar backscatter coefficient (calipso @ 532nm)
	142	betatot_grLidar532, & ! Lidar backscatter coefficient (ground-lidar @ 532nm)
	143	betatot_atlid, & ! Lidar backscatter coefficient (atlid @ 355nm)
	144	betatot_ice_calipso, & ! Lidar backscatter coefficient ICE (calipso @ 532nm)
	145	betatot_liq_calipso, & ! Lidar backscatter coefficient LIQUID (calipso @ 532nm)
	146	tautot_calipso, & ! Lidar Optical thickness (calipso @ 532nm)
	147	tautot_grLidar532, & ! Lidar Optical thickness (ground-lidar @ 532nm)
	148	tautot_atlid, & ! Lidar Optical thickness (atlid @ 355nm)
	149	tautot_ice_calipso, & ! Lidar Ice Optical thickness (calipso @ 532nm)
	150	tautot_liq_calipso, & ! Lidar Liquid Optical thickness (calipso @ 532nm)
	151	z_vol_cloudsat, & ! Effective reflectivity factor (mm^6/m^3)
	152	kr_vol_cloudsat, & ! Attenuation coefficient hydro (dB/km)
	153	g_vol_cloudsat ! Attenuation coefficient gases (dB/km)
	154	real(wp),allocatable,dimension(:,:) :: &
	155	beta_mol_calipso, & ! Lidar molecular backscatter coefficient (calipso @ 532nm)
	156	beta_mol_grLidar532, & ! Lidar molecular backscatter coefficient (ground-lidar @ 532nm)
	157	beta_mol_atlid, & ! Lidar molecular backscatter coefficient (atlid @ 355nm)
	158	tau_mol_calipso, & ! Lidar molecular optical depth (calipso @ 532nm)
	159	tau_mol_grLidar532, & ! Lidar molecular optical depth (ground-lidar @ 532nm)
	160	tau_mol_atlid, & ! Lidar molecular optical depth (atlid @ 355nm)
	161	tautot_S_liq, & ! Parasol Liquid water optical thickness, from TOA to SFC
	162	tautot_S_ice, & ! Parasol Ice water optical thickness, from TOA to SFC
	163	fracPrecipIce ! Fraction of precipitation which is frozen (1).
	164	type(radar_cfg) :: &
	165	rcfg_cloudsat ! Radar configuration information (CLOUDSAT)
	166	end type cosp_optical_inputs
	167
	168	! ######################################################################################
	169	! TYPE cosp_outputs
	170	! ######################################################################################
	171	type cosp_outputs
	172
	173	! CALIPSO outputs
	174	real(wp),dimension(:,:,:),pointer :: &
	175	calipso_betaperp_tot => null(), & ! Total backscattered signal
	176	calipso_beta_tot => null(), & ! Total backscattered signal
	177	calipso_tau_tot => null(), & ! Optical thickness integrated from top to level z
	178	calipso_lidarcldphase => null(), & ! 3D "lidar" phase cloud fraction
	179	calipso_lidarcldtype => null(), & ! 3D "lidar" OPAQ type fraction
	180	calipso_cldlayerphase => null(), & ! low, mid, high-level lidar phase cloud cover
	181	calipso_lidarcldtmp => null(), & ! 3D "lidar" phase cloud temperature
	182	calipso_cfad_sr => null() ! CFAD of scattering ratio
	183	real(wp), dimension(:,:),pointer :: &
	184	calipso_lidarcld => null(), & ! 3D "lidar" cloud fraction
	185	calipso_cldlayer => null(), & ! low, mid, high-level, total lidar cloud cover
	186	calipso_cldtype => null(), & ! opaque and thin lidar cloud cover + z_opaque altitude
	187	calipso_cldtypetemp => null(), & ! opaque and thin cloud temperature
	188	calipso_cldtypemeanz => null(), & ! opaque and thin cloud altitude
	189	calipso_cldtypemeanzse => null(),& ! same as just above with respect to SE
	190	calipso_beta_mol => null(), & ! Molecular backscatter
	191	calipso_temp_tot => null()
	192	real(wp), dimension(:),pointer :: &
	193	calipso_cldthinemis => null(), & ! thin cloud emissivity
	194	calipso_srbval => null() ! SR bins in cfad_sr
	195
	196	! GROUND LIDAR outputs
	197	real(wp),dimension(:,:,:),pointer :: &
	198	grLidar532_beta_tot => null(), & ! Total GROUND LIDAR backscattered signal
	199	grLidar532_cfad_sr => null() ! CFAD of GROUND LIDAR scattering ratio
	200	real(wp), dimension(:,:),pointer :: &
	201	grLidar532_lidarcld => null(), & ! 3D GROUND "lidar" cloud fraction
	202	grLidar532_cldlayer => null(), & ! low, mid, high-level, total GROUND lidar cloud cover
	203	grLidar532_beta_mol => null() ! GROUND LIDAR Molecular backscatter
	204	real(wp), dimension(:),pointer :: &
	205	grLidar532_srbval => null() ! SR bins in cfad_sr
	206
	207	! ATLID outputs
	208	real(wp),dimension(:,:,:),pointer :: &
	209	atlid_beta_tot => null(), & ! Total ATLID backscattered signal
	210	atlid_cfad_sr => null() ! CFAD of ATLID scattering ratio
	211	real(wp), dimension(:,:),pointer :: &
	212	atlid_lidarcld => null(), & ! 3D ATLID cloud fraction
	213	atlid_cldlayer => null(), & ! low, mid, high-level, total ATLID cloud cover
	214	atlid_beta_mol => null() ! ATLID Molecular backscatter
	215	real(wp), dimension(:),pointer :: &
	216	atlid_srbval => null() ! SR bins in cfad_sr
	217
	218
	219	! PARASOL outputs
	220	real(wp),dimension(:,:,:),pointer :: &
	221	parasolPix_refl => null() ! PARASOL reflectances (subcolumn)
	222	real(wp),dimension(:,:),pointer :: &
	223	parasolGrid_refl => null() ! PARASOOL reflectances (column)
	224
	225	! CLOUDSAT outputs
	226	real(wp),dimension(:,:,:),pointer :: &
	227	cloudsat_Ze_tot => null(), & ! Effective reflectivity factor (Npoints,Ncolumns,Nlevels)
	228	cloudsat_cfad_ze => null() ! Ze CFAD(Npoints,dBZe_bins,Nlevels)
	229	real(wp), dimension(:,:),pointer :: &
	230	lidar_only_freq_cloud => null(), & ! (Npoints,Nlevels)
	231	cloudsat_precip_cover => null() ! Radar total cloud amount by CloudSat precip flag (Npoints,dBZe_bins)
	232	real(wp),dimension(:),pointer :: &
	233	cloudsat_tcc => null(), &
	234	cloudsat_tcc2 => null(), &
	235	radar_lidar_tcc => null(), & ! Radar&lidar total cloud amount, grid-box scale (Npoints)
	236	cloudsat_pia => null() ! Radar path integrated attenuation (Npoints)
	237
	238	! ISCCP outputs
	239	real(wp),dimension(:),pointer :: &
	240	isccp_totalcldarea => null(), & ! The fraction of model grid box columns with cloud
	241	! somewhere in them. (%)
	242	isccp_meantb => null(), & ! Mean all-sky 10.5 micron brightness temperature. (K)
	243	isccp_meantbclr => null(), & ! Mean clear-sky 10.5 micron brightness temperature. (K)
	244	isccp_meanptop => null(), & ! Mean cloud top pressure (mb).
	245	isccp_meantaucld => null(), & ! Mean optical thickness. (1)
	246	isccp_meanalbedocld => null() ! Mean cloud albedo. (1)
	247	real(wp),dimension(:,:),pointer ::&
	248	isccp_boxtau => null(), & ! Optical thickness in each column. (1)
	249	isccp_boxptop => null() ! Cloud top pressure in each column. (mb)
	250	real(wp),dimension(:,:,:),pointer :: &
	251	isccp_fq => null() ! The fraction of the model grid box covered by each of
	252	! the 49 ISCCP D level cloud types. (%)
	253
	254	! MISR outptus
	255	real(wp),dimension(:,:,:),pointer :: & !
	256	misr_fq => null() ! Fraction of the model grid box covered by each of the MISR
	257	! cloud types
	258	real(wp),dimension(:,:),pointer :: & !
	259	misr_dist_model_layertops => null() !
	260	real(wp),dimension(:),pointer :: & !
	261	misr_meanztop => null(), & ! Mean MISR cloud top height
	262	misr_cldarea => null() ! Mean MISR cloud cover area
	263
	264	! MODIS outptus
	265	real(wp),pointer,dimension(:) :: & !
	266	modis_Cloud_Fraction_Total_Mean => null(), & ! L3 MODIS retrieved cloud fraction (total)
	267	modis_Cloud_Fraction_Water_Mean => null(), & ! L3 MODIS retrieved cloud fraction (liq)
	268	modis_Cloud_Fraction_Ice_Mean => null(), & ! L3 MODIS retrieved cloud fraction (ice)
	269	modis_Cloud_Fraction_High_Mean => null(), & ! L3 MODIS retrieved cloud fraction (high)
	270	modis_Cloud_Fraction_Mid_Mean => null(), & ! L3 MODIS retrieved cloud fraction (middle)
	271	modis_Cloud_Fraction_Low_Mean => null(), & ! L3 MODIS retrieved cloud fraction (low )
	272	modis_Optical_Thickness_Total_Mean => null(), & ! L3 MODIS retrieved optical thickness (tot)
	273	modis_Optical_Thickness_Water_Mean => null(), & ! L3 MODIS retrieved optical thickness (liq)
	274	modis_Optical_Thickness_Ice_Mean => null(), & ! L3 MODIS retrieved optical thickness (ice)
	275	modis_Optical_Thickness_Total_LogMean => null(), & ! L3 MODIS retrieved log10 optical thickness
	276	modis_Optical_Thickness_Water_LogMean => null(), & ! L3 MODIS retrieved log10 optical thickness
	277	modis_Optical_Thickness_Ice_LogMean => null(), & ! L3 MODIS retrieved log10 optical thickness
	278	modis_Cloud_Particle_Size_Water_Mean => null(), & ! L3 MODIS retrieved particle size (liquid)
	279	modis_Cloud_Particle_Size_Ice_Mean => null(), & ! L3 MODIS retrieved particle size (ice)
	280	modis_Cloud_Top_Pressure_Total_Mean => null(), & ! L3 MODIS retrieved cloud top pressure
	281	modis_Liquid_Water_Path_Mean => null(), & ! L3 MODIS retrieved liquid water path
	282	modis_Ice_Water_Path_Mean => null() ! L3 MODIS retrieved ice water path
	283	real(wp),pointer,dimension(:,:,:) :: &
	284	modis_Optical_Thickness_vs_Cloud_Top_Pressure => null(), & ! Tau/Pressure joint histogram
	285	modis_Optical_Thickness_vs_ReffICE => null(), & ! Tau/ReffICE joint histogram
	286	modis_Optical_Thickness_vs_ReffLIQ => null() ! Tau/ReffLIQ joint histogram
	287
	288	! RTTOV outputs
	289	real(wp),pointer :: &
	290	rttov_tbs(:,:) => null() ! Brightness Temperature
	291
	292	end type cosp_outputs
	293
	294	CONTAINS
	295	! ######################################################################################
	296	! FUNCTION cosp_simulator
	297	! ######################################################################################
	298	function COSP_SIMULATOR(cospIN,cospgridIN,cospOUT,start_idx,stop_idx,debug)
	299	type(cosp_optical_inputs),intent(in),target :: cospIN ! Optical inputs to COSP simulator
	300	type(cosp_column_inputs), intent(in),target :: cospgridIN ! Host model inputs to COSP
	301
	302	! Inputs into the simulators
	303	type(isccp_IN) :: isccpIN ! Input to the ISCCP simulator
	304	type(misr_IN) :: misrIN ! Input to the LIDAR simulator
	305	type(calipso_IN) :: calipsoIN ! Input to the LIDAR simulator
	306	type(grLidar532_IN) :: grLidar532IN ! Input to the GROUND LIDAR simulator
	307	type(atlid_IN) :: atlidIN ! Input to the ATLID simulator
	308	type(parasol_IN) :: parasolIN ! Input to the PARASOL simulator
	309	type(cloudsat_IN) :: cloudsatIN ! Input to the CLOUDSAT radar simulator
	310	type(modis_IN) :: modisIN ! Input to the MODIS simulator
	311	! type(rttov_IN) :: rttovIN ! Input to the RTTOV simulator
	312	integer,optional :: start_idx,stop_idx
	313	logical,optional :: debug
	314
	315	! Outputs from the simulators (nested simulator output structure)
	316	type(cosp_outputs), intent(inout) :: cospOUT
	317	character(len=256),dimension(100) :: cosp_simulator
	318
	319	! Local variables
	320	integer :: &
	321	i,icol,ij,ik,nError
	322	integer,target :: &
	323	Npoints
	324	logical :: &
	325	Lisccp_subcolumn, & ! On/Off switch for subcolumn ISCCP simulator
	326	Lmisr_subcolumn, & ! On/Off switch for subcolumn MISR simulator
	327	Lcalipso_subcolumn, & ! On/Off switch for subcolumn CALIPSO simulator
	328	LgrLidar532_subcolumn,& ! On/Off switch for subcolumn GROUND LIDAR simulator
	329	Latlid_subcolumn, & ! On/Off switch for subcolumn ATLID simulator
	330	Lparasol_subcolumn, & ! On/Off switch for subcolumn PARASOL simulator
	331	Lcloudsat_subcolumn, & ! On/Off switch for subcolumn CLOUDSAT simulator
	332	Lmodis_subcolumn, & ! On/Off switch for subcolumn MODIS simulator
	333	Lrttov_subcolumn, & ! On/Off switch for subcolumn RTTOV simulator
	334	Lisccp_column, & ! On/Off switch for column ISCCP simulator
	335	Lmisr_column, & ! On/Off switch for column MISR simulator
	336	Lcalipso_column, & ! On/Off switch for column CALIPSO simulator
	337	LgrLidar532_column, & ! On/Off switch for column GROUND LIDAR simulator
	338	Latlid_column, & ! On/Off switch for column ATLID simulator
	339	Lparasol_column, & ! On/Off switch for column PARASOL simulator
	340	Lcloudsat_column, & ! On/Off switch for column CLOUDSAT simulator
	341	Lmodis_column, & ! On/Off switch for column MODIS simulator
	342	Lrttov_column, & ! On/Off switch for column RTTOV simulator (not used)
	343	Lradar_lidar_tcc, & ! On/Off switch from joint Calipso/Cloudsat product
	344	Lcloudsat_tcc, & !
	345	Lcloudsat_tcc2, & !
	346	Llidar_only_freq_cloud ! On/Off switch from joint Calipso/Cloudsat product
	347	logical :: &
	348	ok_lidar_cfad = .false., &
	349	ok_lidar_cfad_grLidar532 = .false., &
	350	ok_lidar_cfad_atlid = .false., &
	351	lrttov_cleanUp = .false.
	352
	353	integer, dimension(:,:),allocatable :: &
	354	modisRetrievedPhase,isccpLEVMATCH
	355	real(wp), dimension(:), allocatable :: &
	356	modisCfTotal,modisCfLiquid,modisMeanIceWaterPath, isccp_meantbclr, &
	357	modisCfIce, modisCfHigh, modisCfMid, modisCfLow,modisMeanTauTotal, &
	358	modisMeanTauLiquid, modisMeanTauIce, modisMeanLogTauTotal, &
	359	modisMeanLogTauLiquid, modisMeanLogTauIce, modisMeanSizeLiquid, &
	360	modisMeanSizeIce, modisMeanCloudTopPressure, modisMeanLiquidWaterPath, &
	361	radar_lidar_tcc, cloudsat_tcc, cloudsat_tcc2
	362	REAL(WP), dimension(:,:),allocatable :: &
	363	modisRetrievedCloudTopPressure,modisRetrievedTau,modisRetrievedSize, &
	364	misr_boxtau,misr_boxztop,misr_dist_model_layertops,isccp_boxtau, &
	365	isccp_boxttop,isccp_boxptop,calipso_beta_mol,lidar_only_freq_cloud, &
	366	grLidar532_beta_mol,atlid_beta_mol
	367	REAL(WP), dimension(:,:,:),allocatable :: &
	368	modisJointHistogram,modisJointHistogramIce,modisJointHistogramLiq, &
	369	calipso_beta_tot,calipso_betaperp_tot, cloudsatDBZe,parasolPix_refl, &
	370	grLidar532_beta_tot,atlid_beta_tot,cloudsatZe_non
	371	real(wp),dimension(:),allocatable,target :: &
	372	out1D_1,out1D_2,out1D_3,out1D_4,out1D_5,out1D_6,out1D_7,out1D_8, &
	373	out1D_9,out1D_10,out1D_11,out1D_12
	374	real(wp),dimension(:,:,:),allocatable :: &
	375	betamol_in,betamoli,pnormi,ze_toti,ze_noni
	376
	377	! Initialize error reporting for output
	378	cosp_simulator(:)=''
	379
	380	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	381	! 1) Determine if using full inputs or subset
	382	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	383	if (present(start_idx) .and. present(stop_idx)) then
	384	ij=start_idx
	385	ik=stop_idx
	386	else
	387	ij=1
	388	ik=cospIN%Npoints
	389	endif
	390	Npoints = ik-ij+1
	391
	392	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	393	! 2a) Determine which simulators to run and which statistics to compute
	394	! - If any of the subcolumn fields are allocated, then run the subcolumn simulators.
	395	! - If any of the column fields are allocated, then compute the statistics for that
	396	! simulator, but only save the requested fields.
	397	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	398	! Start with all simulators and joint-diagnostics off
	399	Lisccp_subcolumn = .false.
	400	Lmisr_subcolumn = .false.
	401	Lcalipso_subcolumn = .false.
	402	LgrLidar532_subcolumn = .false.
	403	Latlid_subcolumn = .false.
	404	Lparasol_subcolumn = .false.
	405	Lcloudsat_subcolumn = .false.
	406	Lmodis_subcolumn = .false.
	407	Lrttov_subcolumn = .false.
	408	Lisccp_column = .false.
	409	Lmisr_column = .false.
	410	Lcalipso_column = .false.
	411	LgrLidar532_column = .false.
	412	Latlid_column = .false.
	413	Lparasol_column = .false.
	414	Lcloudsat_column = .false.
	415	Lmodis_column = .false.
	416	Lrttov_column = .false.
	417	Lradar_lidar_tcc = .false.
	418	Llidar_only_freq_cloud = .false.
	419	Lcloudsat_tcc = .false.
	420	Lcloudsat_tcc2 = .false.
	421
	422	! CLOUDSAT subcolumn
	423	if (associated(cospOUT%cloudsat_Ze_tot)) Lcloudsat_subcolumn = .true.
	424
	425	! MODIS subcolumn
	426	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean) .or. &
	427	associated(cospOUT%modis_Cloud_Fraction_Total_Mean) .or. &
	428	associated(cospOUT%modis_Cloud_Fraction_Ice_Mean) .or. &
	429	associated(cospOUT%modis_Cloud_Fraction_High_Mean) .or. &
	430	associated(cospOUT%modis_Cloud_Fraction_Mid_Mean) .or. &
	431	associated(cospOUT%modis_Cloud_Fraction_Low_Mean) .or. &
	432	associated(cospOUT%modis_Optical_Thickness_Total_Mean) .or. &
	433	associated(cospOUT%modis_Optical_Thickness_Water_Mean) .or. &
	434	associated(cospOUT%modis_Optical_Thickness_Ice_Mean) .or. &
	435	associated(cospOUT%modis_Optical_Thickness_Total_LogMean) .or. &
	436	associated(cospOUT%modis_Optical_Thickness_Water_LogMean) .or. &
	437	associated(cospOUT%modis_Optical_Thickness_Ice_LogMean) .or. &
	438	associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean) .or. &
	439	associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean) .or. &
	440	associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean) .or. &
	441	associated(cospOUT%modis_Liquid_Water_Path_Mean) .or. &
	442	associated(cospOUT%modis_Ice_Water_Path_Mean) .or. &
	443	associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	444	Lmodis_subcolumn = .true.
	445
	446	! ISCCP subcolumn
	447	if (associated(cospOUT%isccp_boxtau) .or. &
	448	associated(cospOUT%isccp_boxptop)) &
	449	Lisccp_subcolumn = .true.
	450
	451	! MISR subcolumn
	452	if (associated(cospOUT%misr_dist_model_layertops)) &
	453	Lmisr_subcolumn = .true.
	454
	455	! CALIPOSO subcolumn
	456	if (associated(cospOUT%calipso_tau_tot) .or. &
	457	associated(cospOUT%calipso_beta_mol) .or. &
	458	associated(cospOUT%calipso_temp_tot) .or. &
	459	associated(cospOUT%calipso_betaperp_tot) .or. &
	460	associated(cospOUT%calipso_beta_tot)) &
	461	Lcalipso_subcolumn = .true.
	462
	463	! GROUND LIDAR subcolumn
	464	if (associated(cospOUT%grLidar532_beta_mol) .or. &
	465	associated(cospOUT%grLidar532_beta_tot)) &
	466	LgrLidar532_subcolumn = .true.
	467
	468	! ATLID subcolumn
	469	if (associated(cospOUT%atlid_beta_mol) .or. &
	470	associated(cospOUT%atlid_beta_tot)) &
	471	Latlid_subcolumn = .true.
	472
	473	! PARASOL subcolumn
	474	if (associated(cospOUT%parasolPix_refl)) &
	475	Lparasol_subcolumn = .true.
	476
	477	! RTTOV column
	478	if (associated(cospOUT%rttov_tbs)) &
	479	Lrttov_column = .true.
	480
	481	! Set flag to deallocate rttov types (only done on final call to simulator)
	482	if (size(cospOUT%isccp_meantb) .eq. stop_idx) lrttov_cleanUp = .true.
	483
	484	! ISCCP column
	485	if (associated(cospOUT%isccp_fq) .or. &
	486	associated(cospOUT%isccp_meanalbedocld) .or. &
	487	associated(cospOUT%isccp_meanptop) .or. &
	488	associated(cospOUT%isccp_meantaucld) .or. &
	489	associated(cospOUT%isccp_totalcldarea) .or. &
	490	associated(cospOUT%isccp_meantb)) then
	491	Lisccp_column = .true.
	492	Lisccp_subcolumn = .true.
	493	endif
	494
	495	! MISR column
	496	if (associated(cospOUT%misr_cldarea) .or. &
	497	associated(cospOUT%misr_meanztop) .or. &
	498	associated(cospOUT%misr_fq)) then
	499	Lmisr_column = .true.
	500	Lmisr_subcolumn = .true.
	501	endif
	502
	503	! CALIPSO column
	504	if (associated(cospOUT%calipso_cfad_sr) .or. &
	505	associated(cospOUT%calipso_lidarcld) .or. &
	506	associated(cospOUT%calipso_lidarcldphase) .or. &
	507	associated(cospOUT%calipso_lidarcldtype) .or. &
	508	associated(cospOUT%calipso_cldlayer) .or. &
	509	associated(cospOUT%calipso_cldtype) .or. &
	510	associated(cospOUT%calipso_cldtypetemp) .or. &
	511	associated(cospOUT%calipso_cldtypemeanz) .or. &
	512	associated(cospOUT%calipso_cldtypemeanzse) .or. &
	513	associated(cospOUT%calipso_cldthinemis) .or. &
	514	associated(cospOUT%calipso_cldlayerphase) .or. &
	515	associated(cospOUT%calipso_lidarcldtmp)) then
	516	Lcalipso_column = .true.
	517	Lcalipso_subcolumn = .true.
	518	endif
	519
	520	! GROUND LIDAR column
	521	if (associated(cospOUT%grLidar532_cfad_sr) .or. &
	522	associated(cospOUT%grLidar532_lidarcld) .or. &
	523	associated(cospOUT%grLidar532_cldlayer)) then
	524	LgrLidar532_column = .true.
	525	LgrLidar532_subcolumn = .true.
	526	endif
	527
	528	! ATLID column
	529	if (associated(cospOUT%atlid_cfad_sr) .or. &
	530	associated(cospOUT%atlid_lidarcld) .or. &
	531	associated(cospOUT%atlid_cldlayer)) then
	532	Latlid_column = .true.
	533	Latlid_subcolumn = .true.
	534	endif
	535
	536	! PARASOL column
	537	if (associated(cospOUT%parasolGrid_refl)) then
	538	Lparasol_column = .true.
	539	Lparasol_subcolumn = .true.
	540	endif
	541
	542	! CLOUDSAT column
	543	if (associated(cospOUT%cloudsat_cfad_ze)) then
	544	Lcloudsat_column = .true.
	545	Lcloudsat_subcolumn = .true.
	546	endif
	547
	548	! MODIS column
	549	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean) .or. &
	550	associated(cospOUT%modis_Cloud_Fraction_Water_Mean) .or. &
	551	associated(cospOUT%modis_Cloud_Fraction_Ice_Mean) .or. &
	552	associated(cospOUT%modis_Cloud_Fraction_High_Mean) .or. &
	553	associated(cospOUT%modis_Cloud_Fraction_Mid_Mean) .or. &
	554	associated(cospOUT%modis_Cloud_Fraction_Low_Mean) .or. &
	555	associated(cospOUT%modis_Optical_Thickness_Total_Mean) .or. &
	556	associated(cospOUT%modis_Optical_Thickness_Water_Mean) .or. &
	557	associated(cospOUT%modis_Optical_Thickness_Ice_Mean) .or. &
	558	associated(cospOUT%modis_Optical_Thickness_Total_LogMean) .or. &
	559	associated(cospOUT%modis_Optical_Thickness_Water_LogMean) .or. &
	560	associated(cospOUT%modis_Optical_Thickness_Ice_LogMean) .or. &
	561	associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean) .or. &
	562	associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean) .or. &
	563	associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean) .or. &
	564	associated(cospOUT%modis_Liquid_Water_Path_Mean) .or. &
	565	associated(cospOUT%modis_Ice_Water_Path_Mean) .or. &
	566	associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) then
	567	Lmodis_column = .true.
	568	Lmodis_subcolumn = .true.
	569	endif
	570
	571	! Joint simulator products
	572	if (associated(cospOUT%lidar_only_freq_cloud) .or. associated(cospOUT%radar_lidar_tcc) .or. &
	573	associated(cospOUT%cloudsat_tcc) .or. associated(cospOUT%cloudsat_tcc2)) then
	574	Lcalipso_column = .true.
	575	Lcalipso_subcolumn = .true.
	576	Lcloudsat_column = .true.
	577	Lcloudsat_subcolumn = .true.
	578	Lradar_lidar_tcc = .true.
	579	Llidar_only_freq_cloud = .true.
	580	Lcloudsat_tcc = .true.
	581	Lcloudsat_tcc2 = .true.
	582	endif
	583
	584	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	585	! 2b) Error Checking
	586	! Enforce bounds on input fields. If input field is out-of-bounds, report error
	587	! and turn off simulator
	588	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	589	call cosp_errorCheck(cospgridIN, cospIN, Lisccp_subcolumn, Lisccp_column, &
	590	Lmisr_subcolumn, Lmisr_column, Lmodis_subcolumn, Lmodis_column, &
	591	Lcloudsat_subcolumn, Lcloudsat_column, Lcalipso_subcolumn, Lcalipso_column, &
	592	Latlid_subcolumn, Latlid_column, LgrLidar532_subcolumn, LgrLidar532_column, &
	593	Lrttov_subcolumn, Lrttov_column, Lparasol_subcolumn, Lparasol_column, &
	594	Lradar_lidar_tcc, Llidar_only_freq_cloud, Lcloudsat_tcc,Lcloudsat_tcc2, cospOUT,&
	595	cosp_simulator, nError)
	596
	597	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	598	! 3) Populate instrument simulator inputs
	599	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	600	if (Lisccp_subcolumn .or. Lmodis_subcolumn) then
	601	isccpIN%Npoints => Npoints
	602	isccpIN%Ncolumns => cospIN%Ncolumns
	603	isccpIN%Nlevels => cospIN%Nlevels
	604	isccpIN%emsfc_lw => cospIN%emsfc_lw
	605	isccpIN%skt => cospgridIN%skt
	606	isccpIN%qv => cospgridIN%qv
	607	isccpIN%at => cospgridIN%at
	608	isccpIN%frac_out => cospIN%frac_out
	609	isccpIN%dtau => cospIN%tau_067
	610	isccpIN%dem => cospIN%emiss_11
	611	isccpIN%phalf => cospgridIN%phalf
	612	isccpIN%sunlit => cospgridIN%sunlit
	613	isccpIN%pfull => cospgridIN%pfull
	614	endif
	615
	616	if (Lmisr_subcolumn) then
	617	misrIN%Npoints => Npoints
	618	misrIN%Ncolumns => cospIN%Ncolumns
	619	misrIN%Nlevels => cospIN%Nlevels
	620	misrIN%dtau => cospIN%tau_067
	621	misrIN%sunlit => cospgridIN%sunlit
	622	misrIN%zfull => cospgridIN%hgt_matrix
	623	misrIN%at => cospgridIN%at
	624	endif
	625
	626	if (Lcalipso_subcolumn) then
	627	calipsoIN%Npoints => Npoints
	628	calipsoIN%Ncolumns => cospIN%Ncolumns
	629	calipsoIN%Nlevels => cospIN%Nlevels
	630	calipsoIN%beta_mol => cospIN%beta_mol_calipso
	631	calipsoIN%betatot => cospIN%betatot_calipso
	632	calipsoIN%betatot_liq => cospIN%betatot_liq_calipso
	633	calipsoIN%betatot_ice => cospIN%betatot_ice_calipso
	634	calipsoIN%tau_mol => cospIN%tau_mol_calipso
	635	calipsoIN%tautot => cospIN%tautot_calipso
	636	calipsoIN%tautot_liq => cospIN%tautot_liq_calipso
	637	calipsoIN%tautot_ice => cospIN%tautot_ice_calipso
	638	endif
	639
	640	if (LgrLidar532_subcolumn) then
	641	grLidar532IN%Npoints => Npoints
	642	grLidar532IN%Ncolumns => cospIN%Ncolumns
	643	grLidar532IN%Nlevels => cospIN%Nlevels
	644	grLidar532IN%beta_mol => cospIN%beta_mol_grLidar532
	645	grLidar532IN%betatot => cospIN%betatot_grLidar532
	646	grLidar532IN%tau_mol => cospIN%tau_mol_grLidar532
	647	grLidar532IN%tautot => cospIN%tautot_grLidar532
	648	endif
	649
	650	if (Latlid_subcolumn) then
	651	atlidIN%Npoints => Npoints
	652	atlidIN%Ncolumns => cospIN%Ncolumns
	653	atlidIN%Nlevels => cospIN%Nlevels
	654	atlidIN%beta_mol_atlid => cospIN%beta_mol_atlid
	655	atlidIN%betatot_atlid => cospIN%betatot_atlid
	656	atlidIN%tau_mol_atlid => cospIN%tau_mol_atlid
	657	atlidIN%tautot_atlid => cospIN%tautot_atlid
	658	endif
	659
	660	if (Lparasol_subcolumn) then
	661	parasolIN%Npoints => Npoints
	662	parasolIN%Nlevels => cospIN%Nlevels
	663	parasolIN%Ncolumns => cospIN%Ncolumns
	664	parasolIN%Nrefl => cospIN%Nrefl
	665	parasolIN%tautot_S_liq => cospIN%tautot_S_liq
	666	parasolIN%tautot_S_ice => cospIN%tautot_S_ice
	667	endif
	668
	669	if (Lcloudsat_subcolumn) then
	670	cloudsatIN%Npoints => Npoints
	671	cloudsatIN%Nlevels => cospIN%Nlevels
	672	cloudsatIN%Ncolumns => cospIN%Ncolumns
	673	cloudsatIN%z_vol => cospIN%z_vol_cloudsat
	674	cloudsatIN%kr_vol => cospIN%kr_vol_cloudsat
	675	cloudsatIN%g_vol => cospIN%g_vol_cloudsat
	676	cloudsatIN%rcfg => cospIN%rcfg_cloudsat
	677	cloudsatIN%hgt_matrix => cospgridIN%hgt_matrix
	678	endif
	679
	680	if (Lmodis_subcolumn) then
	681	modisIN%Ncolumns => cospIN%Ncolumns
	682	modisIN%Nlevels => cospIN%Nlevels
	683	modisIN%Npoints => Npoints
	684	modisIN%liqFrac => cospIN%fracLiq
	685	modisIN%tau => cospIN%tau_067
	686	modisIN%g => cospIN%asym
	687	modisIN%w0 => cospIN%ss_alb
	688	modisIN%Nsunlit = count(cospgridIN%sunlit > 0)
	689	if (modisIN%Nsunlit .gt. 0) then
	690	allocate(modisIN%sunlit(modisIN%Nsunlit),modisIN%pres(modisIN%Nsunlit,cospIN%Nlevels+1))
	691	modisIN%sunlit = pack((/ (i, i = 1, Npoints ) /),mask = cospgridIN%sunlit > 0)
	692	modisIN%pres = cospgridIN%phalf(int(modisIN%sunlit(:)),:)
	693	endif
	694	if (count(cospgridIN%sunlit <= 0) .gt. 0) then
	695	allocate(modisIN%notSunlit(count(cospgridIN%sunlit <= 0)))
	696	modisIN%notSunlit = pack((/ (i, i = 1, Npoints ) /),mask = .not. cospgridIN%sunlit > 0)
	697	endif
	698	endif
	699
	700	! if (Lrttov_column) then
	701	! rttovIN%nPoints => Npoints
	702	! rttovIN%nLevels => cospIN%nLevels
	703	! rttovIN%nSubCols => cospIN%nColumns
	704	! rttovIN%zenang => cospgridIN%zenang
	705	! rttovIN%co2 => cospgridIN%co2
	706	! rttovIN%ch4 => cospgridIN%ch4
	707	! rttovIN%n2o => cospgridIN%n2o
	708	! rttovIN%co => cospgridIN%co
	709	! rttovIN%surfem => cospgridIN%emis_sfc
	710	! rttovIN%h_surf => cospgridIN%hgt_matrix_half(:,cospIN%Nlevels+1)
	711	! rttovIN%u_surf => cospgridIN%u_sfc
	712	! rttovIN%v_surf => cospgridIN%v_sfc
	713	! rttovIN%t_skin => cospgridIN%skt
	714	! rttovIN%p_surf => cospgridIN%phalf(:,cospIN%Nlevels+1)
	715	! rttovIN%q2m => cospgridIN%qv(:,cospIN%Nlevels)
	716	! rttovIN%t2m => cospgridIN%at(:,cospIN%Nlevels)
	717	! rttovIN%lsmask => cospgridIN%land
	718	! rttovIN%latitude => cospgridIN%lat
	719	! rttovIN%longitude => cospgridIN%lon
	720	! rttovIN%seaice => cospgridIN%seaice
	721	! rttovIN%p => cospgridIN%pfull
	722	! rttovIN%ph => cospgridIN%phalf
	723	! rttovIN%t => cospgridIN%at
	724	! rttovIN%q => cospgridIN%qv
	725	! rttovIN%o3 => cospgridIN%o3
	726	! Below only needed for all-sky RTTOV calculation
	727	! rttovIN%month => cospgridIN%month
	728	! rttovIN%tca => cospgridIN%tca
	729	! rttovIN%cldIce => cospgridIN%cloudIce
	730	! rttovIN%cldLiq => cospgridIN%cloudLiq
	731	! rttovIN%fl_rain => cospgridIN%fl_rain
	732	! rttovIN%fl_snow => cospgridIN%fl_snow
	733	! endif
	734
	735	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	736	! 4) Call subcolumn simulators
	737	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	738
	739	! ISCCP (icarus) subcolumn simulator
	740	if (Lisccp_subcolumn .or. Lmodis_subcolumn) then
	741	! Allocate space for local variables
	742	allocate(isccpLEVMATCH(Npoints,isccpIN%Ncolumns), &
	743	isccp_boxttop(Npoints,isccpIN%Ncolumns), &
	744	isccp_boxptop(Npoints,isccpIN%Ncolumns), &
	745	isccp_boxtau(Npoints,isccpIN%Ncolumns), isccp_meantbclr(Npoints))
	746	! Call simulator
	747	call icarus_subcolumn(isccpIN%npoints,isccpIN%ncolumns,isccpIN%nlevels, &
	748	isccpIN%sunlit,isccpIN%dtau,isccpIN%dem,isccpIN%skt, &
	749	isccpIN%emsfc_lw,isccpIN%qv,isccpIN%at,isccpIN%pfull, &
	750	isccpIN%phalf,isccpIN%frac_out,isccpLEVMATCH, &
	751	isccp_boxtau(:,:),isccp_boxptop(:,:), &
	752	isccp_boxttop(:,:),isccp_meantbclr(:))
	753	! Store output (if requested)
	754	if (associated(cospOUT%isccp_boxtau)) then
	755	cospOUT%isccp_boxtau(ij:ik,:) = isccp_boxtau
	756	endif
	757	if (associated(cospOUT%isccp_boxptop)) then
	758	cospOUT%isccp_boxptop(ij:ik,:) = isccp_boxptop
	759	endif
	760	if (associated(cospOUT%isccp_meantbclr)) then
	761	cospOUT%isccp_meantbclr(ij:ik) = isccp_meantbclr
	762	endif
	763	endif
	764
	765	! MISR subcolumn simulator
	766	if (Lmisr_subcolumn) then
	767	! Allocate space for local variables
	768	allocate(misr_boxztop(Npoints,misrIN%Ncolumns), &
	769	misr_boxtau(Npoints,misrIN%Ncolumns), &
	770	misr_dist_model_layertops(Npoints,numMISRHgtBins))
	771	! Call simulator
	772	call misr_subcolumn(misrIN%Npoints,misrIN%Ncolumns,misrIN%Nlevels,misrIN%dtau, &
	773	misrIN%zfull,misrIN%at,misrIN%sunlit,misr_boxtau, &
	774	misr_dist_model_layertops,misr_boxztop)
	775	! Store output (if requested)
	776	if (associated(cospOUT%misr_dist_model_layertops)) then
	777	cospOUT%misr_dist_model_layertops(ij:ik,:) = misr_dist_model_layertops
	778	endif
	779	endif
	780
	781	! Calipso subcolumn simulator
	782	if (Lcalipso_subcolumn) then
	783	! Allocate space for local variables
	784	allocate(calipso_beta_mol(calipsoIN%Npoints,calipsoIN%Nlevels), &
	785	calipso_beta_tot(calipsoIN%Npoints,calipsoIN%Ncolumns,calipsoIN%Nlevels),&
	786	calipso_betaperp_tot(calipsoIN%Npoints,calipsoIN%Ncolumns,calipsoIN%Nlevels))
	787	! Call simulator
	788	call lidar_subcolumn(calipsoIN%npoints, calipsoIN%ncolumns, calipsoIN%nlevels, .false., &
	789	calipsoIN%beta_mol, calipsoIN%tau_mol, calipsoIN%betatot, calipsoIN%tautot, &
	790	calipso_beta_mol(:,:), calipso_beta_tot(:,:,:), calipsoIN%betatot_ice, &
	791	calipsoIN%tautot_ice, calipsoIN%betatot_liq, calipsoIN%tautot_liq, &
	792	calipso_betaperp_tot(:,:,:))
	793	! Store output (if requested)
	794	if (associated(cospOUT%calipso_beta_mol)) &
	795	cospOUT%calipso_beta_mol(ij:ik,calipsoIN%Nlevels:1:-1) = calipso_beta_mol
	796	if (associated(cospOUT%calipso_beta_tot)) &
	797	cospOUT%calipso_beta_tot(ij:ik,:,calipsoIN%Nlevels:1:-1) = calipso_beta_tot
	798	if (associated(cospOUT%calipso_betaperp_tot)) &
	799	cospOUT%calipso_betaperp_tot(ij:ik,:,:) = calipso_betaperp_tot
	800	endif
	801
	802	! GROUND LIDAR subcolumn simulator
	803	if (LgrLidar532_subcolumn) then
	804	! Allocate space for local variables
	805	allocate(grLidar532_beta_mol(grLidar532IN%Npoints,grLidar532IN%Nlevels), &
	806	grLidar532_beta_tot(grLidar532IN%Npoints,grLidar532IN%Ncolumns,grLidar532IN%Nlevels))
	807	! Call simulator
	808	call lidar_subcolumn(grLidar532IN%npoints, grLidar532IN%ncolumns, grLidar532IN%nlevels,&
	809	.true., grLidar532IN%beta_mol, grLidar532IN%tau_mol, grLidar532IN%betatot,&
	810	grLidar532IN%tautot, grLidar532_beta_mol(:,:), grLidar532_beta_tot(:,:,:))
	811	! Store output (if requested)
	812	if (associated(cospOUT%grLidar532_beta_mol)) &
	813	cospOUT%grLidar532_beta_mol(ij:ik,grLidar532IN%Nlevels:1:-1) = grLidar532_beta_mol
	814	if (associated(cospOUT%grLidar532_beta_tot)) &
	815	cospOUT%grLidar532_beta_tot(ij:ik,:,grLidar532IN%Nlevels:1:-1) = grLidar532_beta_tot
	816	endif
	817
	818	! ATLID subcolumn simulator
	819	if (Latlid_subcolumn) then
	820	! Allocate space for local variables
	821	allocate(atlid_beta_mol(atlidIN%Npoints,atlidIN%Nlevels), &
	822	atlid_beta_tot(atlidIN%Npoints,atlidIN%Ncolumns,atlidIN%Nlevels))
	823	! Call simulator
	824	call lidar_subcolumn(atlidIN%npoints, atlidIN%ncolumns, atlidIN%nlevels,&
	825	.false., atlidIN%beta_mol_atlid, atlidIN%tau_mol_atlid, atlidIN%betatot_atlid,&
	826	atlidIN%tautot_atlid, atlid_beta_mol(:,:), atlid_beta_tot(:,:,:))
	827	! Store output (if requested)
	828	if (associated(cospOUT%atlid_beta_mol)) &
	829	cospOUT%atlid_beta_mol(ij:ik,atlidIN%Nlevels:1:-1) = atlid_beta_mol
	830	if (associated(cospOUT%atlid_beta_tot)) &
	831	cospOUT%atlid_beta_tot(ij:ik,:,atlidIN%Nlevels:1:-1) = atlid_beta_tot
	832	endif
	833
	834	! PARASOL subcolumn simulator
	835	if (Lparasol_subcolumn) then
	836	! Allocate space for local variables
	837	allocate(parasolPix_refl(parasolIN%Npoints,parasolIN%Ncolumns,PARASOL_NREFL))
	838	! Call simulator
	839	do icol=1,parasolIN%Ncolumns
	840	call parasol_subcolumn(parasolIN%npoints, PARASOL_NREFL, &
	841	parasolIN%tautot_S_liq(1:parasolIN%Npoints,icol), &
	842	parasolIN%tautot_S_ice(1:parasolIN%Npoints,icol), &
	843	parasolPix_refl(:,icol,1:PARASOL_NREFL))
	844	! Store output (if requested)
	845	if (associated(cospOUT%parasolPix_refl)) then
	846	cospOUT%parasolPix_refl(ij:ik,icol,1:PARASOL_NREFL) = &
	847	parasolPix_refl(:,icol,1:PARASOL_NREFL)
	848	endif
	849	enddo
	850	endif
	851
	852	! Cloudsat (quickbeam) subcolumn simulator
	853	if (Lcloudsat_subcolumn) then
	854	! Allocate space for local variables
	855	allocate(cloudsatDBZe(cloudsatIN%Npoints,cloudsatIN%Ncolumns,cloudsatIN%Nlevels), &
	856	cloudsatZe_non(cloudsatIN%Npoints,cloudsatIN%Ncolumns,cloudsatIN%Nlevels))
	857	do icol=1,cloudsatIN%ncolumns
	858	call quickbeam_subcolumn(cloudsatIN%rcfg,cloudsatIN%Npoints,cloudsatIN%Nlevels,&
	859	cloudsatIN%hgt_matrix/1000._wp, &
	860	cloudsatIN%z_vol(:,icol,:), &
	861	cloudsatIN%kr_vol(:,icol,:), &
	862	cloudsatIN%g_vol(:,1,:),cloudsatDBze(:,icol,:),cloudsatZe_non(:,icol,:))
	863	enddo
	864	! Store output (if requested)
	865	if (associated(cospOUT%cloudsat_Ze_tot)) then
	866	cospOUT%cloudsat_Ze_tot(ij:ik,:,:) = cloudsatDBZe(:,:,1:cloudsatIN%Nlevels)
	867	endif
	868	endif
	869
	870	if (Lmodis_subcolumn) then
	871	if(modisiN%nSunlit > 0) then
	872	! Allocate space for local variables
	873	allocate(modisRetrievedTau(modisIN%nSunlit,modisIN%nColumns), &
	874	modisRetrievedSize(modisIN%nSunlit,modisIN%nColumns), &
	875	modisRetrievedPhase(modisIN%nSunlit,modisIN%nColumns), &
	876	modisRetrievedCloudTopPressure(modisIN%nSunlit,modisIN%nColumns))
	877	! Call simulator
	878	do i = 1, modisIN%nSunlit
	879	call modis_subcolumn(modisIN%Ncolumns,modisIN%Nlevels,modisIN%pres(i,:), &
	880	modisIN%tau(int(modisIN%sunlit(i)),:,:), &
	881	modisIN%liqFrac(int(modisIN%sunlit(i)),:,:), &
	882	modisIN%g(int(modisIN%sunlit(i)),:,:), &
	883	modisIN%w0(int(modisIN%sunlit(i)),:,:), &
	884	isccp_boxptop(int(modisIN%sunlit(i)),:), &
	885	modisRetrievedPhase(i,:), &
	886	modisRetrievedCloudTopPressure(i,:), &
	887	modisRetrievedTau(i,:),modisRetrievedSize(i,:))
	888	end do
	889	endif
	890	endif
	891
	892	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	893	! 5) Call column simulators
	894	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	895
	896	! ISCCP
	897	if (Lisccp_column) then
	898	! Check to see which outputs are requested. If not requested, use a local dummy array
	899	if(.not. associated(cospOUT%isccp_meanalbedocld)) then
	900	allocate(out1D_1(Npoints))
	901	cospOUT%isccp_meanalbedocld(ij:ik) => out1D_1
	902	endif
	903	if(.not. associated(cospOUT%isccp_meanptop)) then
	904	allocate(out1D_2(Npoints))
	905	cospOUT%isccp_meanptop(ij:ik) => out1D_2
	906	endif
	907	if(.not. associated(cospOUT%isccp_meantaucld)) then
	908	allocate(out1D_3(Npoints))
	909	cospOUT%isccp_meantaucld(ij:ik) => out1D_3
	910	endif
	911	if(.not. associated(cospOUT%isccp_totalcldarea)) then
	912	allocate(out1D_4(Npoints))
	913	cospOUT%isccp_totalcldarea(ij:ik) => out1D_4
	914	endif
	915	if(.not. associated(cospOUT%isccp_meantb)) then
	916	allocate(out1D_5(Npoints))
	917	cospOUT%isccp_meantb(ij:ik) => out1D_5
	918	endif
	919	if(.not. associated(cospOUT%isccp_fq)) then
	920	allocate(out1D_6(NpointsnumISCCPTauBinsnumISCCPPresBins))
	921	cospOUT%isccp_fq(ij:ik,1:numISCCPTauBins,1:numISCCPPresBins) => out1D_6
	922	endif
	923
	924	! Call simulator
	925	call icarus_column(isccpIN%npoints, isccpIN%ncolumns,isccp_boxtau(:,:), &
	926	isccp_boxptop(:,:)/100._wp, isccpIN%sunlit,isccp_boxttop, &
	927	cospOUT%isccp_fq(ij:ik,:,:), &
	928	cospOUT%isccp_meanalbedocld(ij:ik), &
	929	cospOUT%isccp_meanptop(ij:ik),cospOUT%isccp_meantaucld(ij:ik), &
	930	cospOUT%isccp_totalcldarea(ij:ik),cospOUT%isccp_meantb(ij:ik))
	931	cospOUT%isccp_fq(ij:ik,:,:) = cospOUT%isccp_fq(ij:ik,:,7:1:-1)
	932
	933	! Check if there is any value slightly greater than 1
	934	where ((cospOUT%isccp_totalcldarea > 1.0-1.e-5) .and. &
	935	(cospOUT%isccp_totalcldarea < 1.0+1.e-5))
	936	cospOUT%isccp_totalcldarea = 1.0
	937	endwhere
	938
	939	! Clear up memory (if necessary)
	940	if (allocated(isccp_boxttop)) deallocate(isccp_boxttop)
	941	if (allocated(isccp_boxptop)) deallocate(isccp_boxptop)
	942	if (allocated(isccp_boxtau)) deallocate(isccp_boxtau)
	943	if (allocated(isccp_meantbclr)) deallocate(isccp_meantbclr)
	944	if (allocated(isccpLEVMATCH)) deallocate(isccpLEVMATCH)
	945	if (allocated(out1D_1)) then
	946	deallocate(out1D_1)
	947	nullify(cospOUT%isccp_meanalbedocld)
	948	endif
	949	if (allocated(out1D_2)) then
	950	deallocate(out1D_2)
	951	nullify(cospOUT%isccp_meanptop)
	952	endif
	953	if (allocated(out1D_3)) then
	954	deallocate(out1D_3)
	955	nullify(cospOUT%isccp_meantaucld)
	956	endif
	957	if (allocated(out1D_4)) then
	958	deallocate(out1D_4)
	959	nullify(cospOUT%isccp_totalcldarea)
	960	endif
	961	if (allocated(out1D_5)) then
	962	deallocate(out1D_5)
	963	nullify(cospOUT%isccp_meantb)
	964	endif
	965	if (allocated(out1D_6)) then
	966	deallocate(out1D_6)
	967	nullify(cospOUT%isccp_fq)
	968	endif
	969	endif
	970
	971	! MISR
	972	if (Lmisr_column) then
	973	! Check to see which outputs are requested. If not requested, use a local dummy array
	974	if (.not. associated(cospOUT%misr_cldarea)) then
	975	allocate(out1D_1(Npoints))
	976	cospOUT%misr_cldarea(ij:ik) => out1D_1
	977	endif
	978	if (.not. associated(cospOUT%misr_meanztop)) then
	979	allocate(out1D_2(Npoints))
	980	cospOUT%misr_meanztop(ij:ik) => out1D_2
	981	endif
	982	if (.not. associated(cospOUT%misr_fq)) then
	983	allocate(out1D_3(NpointsnumMISRTauBinsnumMISRHgtBins))
	984	cospOUT%misr_fq(ij:ik,1:numMISRTauBins,1:numMISRHgtBins) => out1D_3
	985	endif
	986
	987	! Call simulator
	988	call misr_column(misrIN%Npoints,misrIN%Ncolumns,misr_boxztop,misrIN%sunlit,&
	989	misr_boxtau,cospOUT%misr_cldarea(ij:ik), &
	990	cospOUT%misr_meanztop(ij:ik),cospOUT%misr_fq(ij:ik,:,:))
	991
	992	! Clear up memory
	993	if (allocated(misr_boxtau)) deallocate(misr_boxtau)
	994	if (allocated(misr_boxztop)) deallocate(misr_boxztop)
	995	if (allocated(misr_dist_model_layertops)) deallocate(misr_dist_model_layertops)
	996	if (allocated(out1D_1)) then
	997	deallocate(out1D_1)
	998	nullify(cospOUT%misr_cldarea)
	999	endif
	1000	if (allocated(out1D_2)) then
	1001	deallocate(out1D_2)
	1002	nullify(cospOUT%misr_meanztop)
	1003	endif
	1004	if (allocated(out1D_3)) then
	1005	deallocate(out1D_3)
	1006	nullify(cospOUT%misr_fq)
	1007	endif
	1008	endif
	1009
	1010	! CALIPSO LIDAR Simulator
	1011	if (Lcalipso_column) then
	1012	! Check to see which outputs are requested. If not requested, use a local dummy array
	1013	if (.not. associated(cospOUT%calipso_cfad_sr)) then
	1014	allocate(out1D_1(NpointsSR_BINSNlvgrid))
	1015	cospOUT%calipso_cfad_sr(ij:ik,1:SR_BINS,1:Nlvgrid) => out1D_1
	1016	endif
	1017	if (.not. associated(cospOUT%calipso_lidarcld)) then
	1018	allocate(out1D_2(Npoints*Nlvgrid))
	1019	cospOUT%calipso_lidarcld(ij:ik,1:Nlvgrid) => out1D_2
	1020	endif
	1021	if (.not. associated(cospOUT%calipso_lidarcldphase)) then
	1022	allocate(out1D_3(NpointsNlvgrid6))
	1023	cospOUT%calipso_lidarcldphase(ij:ik,1:Nlvgrid,1:6) => out1D_3
	1024	endif
	1025	if (.not. associated(cospOUT%calipso_cldlayer)) then
	1026	allocate(out1D_4(Npoints*LIDAR_NCAT))
	1027	cospOUT%calipso_cldlayer(ij:ik,1:LIDAR_NCAT) => out1D_4
	1028	endif
	1029	if (.not. associated(cospOUT%calipso_cldlayerphase)) then
	1030	allocate(out1D_5(NpointsLIDAR_NCAT6))
	1031	cospOUT%calipso_cldlayerphase(ij:ik,1:LIDAR_NCAT,1:6) => out1D_5
	1032	endif
	1033	if (.not. associated(cospOUT%calipso_lidarcldtmp)) then
	1034	allocate(out1D_6(Npoints405))
	1035	cospOUT%calipso_lidarcldtmp(ij:ik,1:40,1:5) => out1D_6
	1036	endif
	1037	if (.not. associated(cospOUT%calipso_lidarcldtype)) then
	1038	allocate(out1D_7(NpointsNlvgrid4))
	1039	cospOUT%calipso_lidarcldtype(ij:ik,1:Nlvgrid,1:4) => out1D_7
	1040	endif
	1041	if (.not. associated(cospOUT%calipso_cldtype)) then
	1042	allocate(out1D_8(Npoints*LIDAR_NTYPE))
	1043	cospOUT%calipso_cldtype(ij:ik,1:LIDAR_NTYPE) => out1D_8
	1044	endif
	1045	if (.not. associated(cospOUT%calipso_cldtypetemp)) then
	1046	allocate(out1D_9(Npoints*LIDAR_NTYPE))
	1047	cospOUT%calipso_cldtypetemp(ij:ik,1:LIDAR_NTYPE) => out1D_9
	1048	endif
	1049	if (.not. associated(cospOUT%calipso_cldtypemeanz)) then
	1050	allocate(out1D_10(Npoints*2))
	1051	cospOUT%calipso_cldtypemeanz(ij:ik,1:2) => out1D_10
	1052	endif
	1053	if (.not. associated(cospOUT%calipso_cldtypemeanzse)) then
	1054	allocate(out1D_12(Npoints*3))
	1055	cospOUT%calipso_cldtypemeanzse(ij:ik,1:3) => out1D_12
	1056	endif
	1057	if (.not. associated(cospOUT%calipso_cldthinemis)) then
	1058	allocate(out1D_11(Npoints))
	1059	cospOUT%calipso_cldthinemis(ij:ik) => out1D_11
	1060	endif
	1061
	1062	! Call simulator
	1063	ok_lidar_cfad=.true.
	1064	call lidar_column(calipsoIN%Npoints, calipsoIN%Ncolumns, calipsoIN%Nlevels, &
	1065	Nlvgrid, SR_BINS, LIDAR_NTYPE, 'calipso',calipso_beta_tot(:,:,:), calipso_beta_mol(:,:),&
	1066	cospgridIN%phalf(:,2:calipsoIN%Nlevels),cospgridIN%hgt_matrix, &
	1067	cospgridIN%hgt_matrix_half, vgrid_z(:), ok_lidar_cfad, LIDAR_NCAT, &
	1068	cospOUT%calipso_cfad_sr(ij:ik,:,:), cospOUT%calipso_lidarcld(ij:ik,:), &
	1069	cospOUT%calipso_cldlayer(ij:ik,:), &
	1070	cospgridIN%at(:,:), calipso_betaperp_tot(:,:,:), cospgridIN%surfelev, &
	1071	cospOUT%calipso_lidarcldphase(ij:ik,:,:), &
	1072	cospOUT%calipso_lidarcldtype(ij:ik,:,:), cospOUT%calipso_cldtype(ij:ik,:), &
	1073	cospOUT%calipso_cldtypetemp(ij:ik,:), cospOUT%calipso_cldtypemeanz(ij:ik,:), &
	1074	cospOUT%calipso_cldtypemeanzse(ij:ik,:), cospOUT%calipso_cldthinemis(ij:ik), &
	1075	cospOUT%calipso_cldlayerphase(ij:ik,:,:), cospOUT%calipso_lidarcldtmp(ij:ik,:,:))
	1076
	1077	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval = calipso_histBsct
	1078
	1079	! Free up memory (if necessary)
	1080	if (allocated(out1D_1)) then
	1081	deallocate(out1D_1)
	1082	nullify(cospOUT%calipso_cfad_sr)
	1083	endif
	1084	if (allocated(out1D_2)) then
	1085	deallocate(out1D_2)
	1086	nullify(cospOUT%calipso_lidarcld)
	1087	endif
	1088	if (allocated(out1D_3)) then
	1089	deallocate(out1D_3)
	1090	nullify(cospOUT%calipso_lidarcldphase)
	1091	endif
	1092	if (allocated(out1D_4)) then
	1093	deallocate(out1D_4)
	1094	nullify(cospOUT%calipso_cldlayer)
	1095	endif
	1096	if (allocated(out1D_5)) then
	1097	deallocate(out1D_5)
	1098	nullify(cospOUT%calipso_cldlayerphase)
	1099	endif
	1100	if (allocated(out1D_6)) then
	1101	deallocate(out1D_6)
	1102	nullify(cospOUT%calipso_lidarcldtmp)
	1103	endif
	1104	if (allocated(out1D_7)) then
	1105	deallocate(out1D_7)
	1106	nullify(cospOUT%calipso_lidarcldtype)
	1107	endif
	1108	if (allocated(out1D_8)) then
	1109	deallocate(out1D_8)
	1110	nullify(cospOUT%calipso_cldtype)
	1111	endif
	1112	if (allocated(out1D_9)) then
	1113	deallocate(out1D_9)
	1114	nullify(cospOUT%calipso_cldtypetemp)
	1115	endif
	1116	if (allocated(out1D_10)) then
	1117	deallocate(out1D_10)
	1118	nullify(cospOUT%calipso_cldtypemeanz)
	1119	endif
	1120	if (allocated(out1D_12)) then
	1121	deallocate(out1D_12)
	1122	nullify(cospOUT%calipso_cldtypemeanzse)
	1123	endif
	1124	if (allocated(out1D_11)) then
	1125	deallocate(out1D_11)
	1126	nullify(cospOUT%calipso_cldthinemis)
	1127	endif
	1128
	1129	endif
	1130
	1131	! GROUND LIDAR Simulator
	1132	if (LgrLidar532_column) then
	1133	! Check to see which outputs are requested. If not requested, use a local dummy array
	1134	if (.not. associated(cospOUT%grLidar532_cfad_sr)) then
	1135	allocate(out1D_1(NpointsSR_BINSNlvgrid))
	1136	cospOUT%grLidar532_cfad_sr(ij:ik,1:SR_BINS,1:Nlvgrid) => out1D_1
	1137	endif
	1138	if (.not. associated(cospOUT%grLidar532_lidarcld)) then
	1139	allocate(out1D_2(Npoints*Nlvgrid))
	1140	cospOUT%grLidar532_lidarcld(ij:ik,1:Nlvgrid) => out1D_2
	1141	endif
	1142	if (.not. associated(cospOUT%grLidar532_cldlayer)) then
	1143	allocate(out1D_3(Npoints*LIDAR_NCAT))
	1144	cospOUT%grLidar532_cldlayer(ij:ik,1:LIDAR_NCAT) => out1D_3
	1145	endif
	1146
	1147	! Call simulator
	1148	ok_lidar_cfad_grLidar532=.true.
	1149	call lidar_column(grLidar532IN%Npoints, grLidar532IN%Ncolumns, grLidar532IN%Nlevels, &
	1150	Nlvgrid, SR_BINS, LIDAR_NTYPE, 'grlidar532',grLidar532_beta_tot(:,:,:), grLidar532_beta_mol(:,:),&
	1151	cospgridIN%phalf(:,2:grLidar532IN%Nlevels),cospgridIN%hgt_matrix, &
	1152	cospgridIN%hgt_matrix_half, vgrid_z(:), ok_lidar_cfad_grLidar532, LIDAR_NCAT, &
	1153	cospOUT%grLidar532_cfad_sr(ij:ik,:,:), cospOUT%grLidar532_lidarcld(ij:ik,:), &
	1154	cospOUT%grLidar532_cldlayer(ij:ik,:))
	1155
	1156	if (associated(cospOUT%grLidar532_srbval)) cospOUT%grLidar532_srbval = grLidar532_histBsct
	1157
	1158	! Free up memory (if necessary)
	1159	if (allocated(out1D_1)) then
	1160	deallocate(out1D_1)
	1161	nullify(cospOUT%grLidar532_cfad_sr)
	1162	endif
	1163	if (allocated(out1D_2)) then
	1164	deallocate(out1D_2)
	1165	nullify(cospOUT%grLidar532_lidarcld)
	1166	endif
	1167	if (allocated(out1D_3)) then
	1168	deallocate(out1D_3)
	1169	nullify(cospOUT%grLidar532_cldlayer)
	1170	endif
	1171
	1172	endif
	1173
	1174	! ATLID Simulator
	1175	if (Latlid_column) then
	1176	! Check to see which outputs are requested. If not requested, use a local dummy array
	1177	if (.not. associated(cospOUT%atlid_cfad_sr)) then
	1178	allocate(out1D_1(NpointsSR_BINSNlvgrid))
	1179	cospOUT%atlid_cfad_sr(ij:ik,1:SR_BINS,1:Nlvgrid) => out1D_1
	1180	endif
	1181	if (.not. associated(cospOUT%atlid_lidarcld)) then
	1182	allocate(out1D_2(Npoints*Nlvgrid))
	1183	cospOUT%atlid_lidarcld(ij:ik,1:Nlvgrid) => out1D_2
	1184	endif
	1185	if (.not. associated(cospOUT%atlid_cldlayer)) then
	1186	allocate(out1D_3(Npoints*LIDAR_NCAT))
	1187	cospOUT%atlid_cldlayer(ij:ik,1:LIDAR_NCAT) => out1D_3
	1188	endif
	1189
	1190	! Call simulator
	1191	ok_lidar_cfad_atlid=.true.
	1192	call lidar_column(atlidIN%Npoints, atlidIN%Ncolumns, atlidIN%Nlevels, &
	1193	Nlvgrid, SR_BINS, LIDAR_NTYPE, 'atlid',atlid_beta_tot(:,:,:), &
	1194	atlid_beta_mol(:,:), cospgridIN%phalf(:,2:atlidIN%Nlevels), &
	1195	cospgridIN%hgt_matrix, cospgridIN%hgt_matrix_half, vgrid_z(:), &
	1196	ok_lidar_cfad_atlid, LIDAR_NCAT, cospOUT%atlid_cfad_sr(ij:ik,:,:), &
	1197	cospOUT%atlid_lidarcld(ij:ik,:), cospOUT%atlid_cldlayer(ij:ik,:))
	1198
	1199	if (associated(cospOUT%atlid_srbval)) cospOUT%atlid_srbval = atlid_histBsct
	1200
	1201	! Free up memory (if necessary)
	1202	if (allocated(out1D_1)) then
	1203	deallocate(out1D_1)
	1204	nullify(cospOUT%atlid_cfad_sr)
	1205	endif
	1206	if (allocated(out1D_2)) then
	1207	deallocate(out1D_2)
	1208	nullify(cospOUT%atlid_lidarcld)
	1209	endif
	1210	if (allocated(out1D_3)) then
	1211	deallocate(out1D_3)
	1212	nullify(cospOUT%atlid_cldlayer)
	1213	endif
	1214
	1215	endif
	1216
	1217	! PARASOL
	1218	if (Lparasol_column) then
	1219	call parasol_column(parasolIN%Npoints,PARASOL_NREFL,parasolIN%Ncolumns, &
	1220	cospgridIN%land(:),parasolPix_refl(:,:,:), &
	1221	cospOUT%parasolGrid_refl(ij:ik,:))
	1222	if (allocated(parasolPix_refl)) deallocate(parasolPix_refl)
	1223	endif
	1224
	1225	! CLOUDSAT
	1226	if (Lcloudsat_column) then
	1227	! Check to see which outputs are requested. If not requested, use a local dummy array
	1228	if (.not. associated(cospOUT%cloudsat_cfad_ze)) then
	1229	allocate(out1D_1(Npointscloudsat_DBZE_BINSNlvgrid))
	1230	cospOUT%cloudsat_cfad_ze(ij:ik,1:cloudsat_DBZE_BINS,1:Nlvgrid) => out1D_1
	1231	endif
	1232
	1233	if (.not. associated(cospOUT%cloudsat_pia)) then
	1234	allocate(out1D_2(Npoints))
	1235	cospOUT%cloudsat_pia(ij:ik) => out1D_2
	1236	endif
	1237	if (.not. associated(cospOUT%cloudsat_precip_cover)) then
	1238	allocate(out1D_3(Npoints*nCloudsatPrecipClass))
	1239	cospOUT%cloudsat_precip_cover(ij:ik,1:nCloudsatPrecipClass) => out1D_3
	1240	endif
	1241
	1242	! Call simulator
	1243	call quickbeam_column(cloudsatIN%Npoints, cloudsatIN%Ncolumns, cloudsatIN%Nlevels,&
	1244	Nlvgrid, cloudsat_DBZE_BINS, 'cloudsat', cloudsatDBZe, cloudsatZe_non, &
	1245	! cospgridIN%land(:), cospgridIN%at(:,cospIN%Nlevels), cospIN%fracPrecipIce, & !PREC_BUG
	1246	! cospgridIN%hgt_matrix, cospgridIN%hgt_matrix_half, & !PREC_BUG
	1247	cospgridIN%land(:), cospgridIN%surfelev(:), cospgridIN%at(:,cospIN%Nlevels), & !PREC_BUG
	1248	cospIN%fracPrecipIce, cospgridIN%hgt_matrix, cospgridIN%hgt_matrix_half, & !PREC_BUG
	1249	cospOUT%cloudsat_cfad_ze(ij:ik,:,:), cospOUT%cloudsat_precip_cover, &
	1250	cospOUT%cloudsat_pia)
	1251	! Free up memory (if necessary)
	1252	if (allocated(out1D_1)) then
	1253	deallocate(out1D_1)
	1254	nullify(cospOUT%cloudsat_cfad_ze)
	1255	endif
	1256	if (allocated(out1D_2)) then
	1257	deallocate(out1D_2)
	1258	nullify(cospOUT%cloudsat_pia)
	1259	endif
	1260	if (allocated(out1D_3)) then
	1261	deallocate(out1D_3)
	1262	nullify(cospOUT%cloudsat_precip_cover)
	1263	endif
	1264	endif
	1265
	1266	! MODIS
	1267	if (Lmodis_column) then
	1268	if(modisiN%nSunlit > 0) then
	1269	! Allocate space for local variables
	1270	allocate(modisCftotal(modisIN%nSunlit), modisCfLiquid(modisIN%nSunlit), &
	1271	modisCfIce(modisIN%nSunlit),modisCfHigh(modisIN%nSunlit), &
	1272	modisCfMid(modisIN%nSunlit),modisCfLow(modisIN%nSunlit), &
	1273	modisMeanTauTotal(modisIN%nSunlit), &
	1274	modisMeanTauLiquid(modisIN%nSunlit),modisMeanTauIce(modisIN%nSunlit), &
	1275	modisMeanLogTauTotal(modisIN%nSunlit), &
	1276	modisMeanLogTauLiquid(modisIN%nSunlit), &
	1277	modisMeanLogTauIce(modisIN%nSunlit), &
	1278	modisMeanSizeLiquid(modisIN%nSunlit), &
	1279	modisMeanSizeIce(modisIN%nSunlit), &
	1280	modisMeanCloudTopPressure(modisIN%nSunlit), &
	1281	modisMeanLiquidWaterPath(modisIN%nSunlit), &
	1282	modisMeanIceWaterPath(modisIN%nSunlit), &
	1283	modisJointHistogram(modisIN%nSunlit,numMODISTauBins,numMODISPresBins),&
	1284	modisJointHistogramIce(modisIN%nSunlit,numModisTauBins,numMODISReffIceBins),&
	1285	modisJointHistogramLiq(modisIN%nSunlit,numModisTauBins,numMODISReffLiqBins))
	1286	! Call simulator
	1287	call modis_column(modisIN%nSunlit, modisIN%Ncolumns,modisRetrievedPhase, &
	1288	modisRetrievedCloudTopPressure,modisRetrievedTau, &
	1289	modisRetrievedSize, modisCfTotal, modisCfLiquid, modisCfIce,&
	1290	modisCfHigh, modisCfMid, modisCfLow, modisMeanTauTotal, &
	1291	modisMeanTauLiquid, modisMeanTauIce, modisMeanLogTauTotal, &
	1292	modisMeanLogTauLiquid, modisMeanLogTauIce, &
	1293	modisMeanSizeLiquid, modisMeanSizeIce, &
	1294	modisMeanCloudTopPressure, modisMeanLiquidWaterPath, &
	1295	modisMeanIceWaterPath, modisJointHistogram, &
	1296	modisJointHistogramIce,modisJointHistogramLiq)
	1297	! Store data (if requested)
	1298	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) then
	1299	cospOUT%modis_Cloud_Fraction_Total_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1300	modisCfTotal
	1301	endif
	1302	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) then
	1303	cospOUT%modis_Cloud_Fraction_Water_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1304	modisCfLiquid
	1305	endif
	1306	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) then
	1307	cospOUT%modis_Cloud_Fraction_Ice_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1308	modisCfIce
	1309	endif
	1310	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) then
	1311	cospOUT%modis_Cloud_Fraction_High_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1312	modisCfHigh
	1313	endif
	1314	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) then
	1315	cospOUT%modis_Cloud_Fraction_Mid_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1316	modisCfMid
	1317	endif
	1318	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) then
	1319	cospOUT%modis_Cloud_Fraction_Low_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1320	modisCfLow
	1321	endif
	1322	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) then
	1323	cospOUT%modis_Optical_Thickness_Total_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1324	modisMeanTauTotal
	1325	endif
	1326	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) then
	1327	cospOUT%modis_Optical_Thickness_Water_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1328	modisMeanTauLiquid
	1329	endif
	1330	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) then
	1331	cospOUT%modis_Optical_Thickness_Ice_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1332	modisMeanTauIce
	1333	endif
	1334	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) then
	1335	cospOUT%modis_Optical_Thickness_Total_LogMean(ij+int(modisIN%sunlit(:))-1)= &
	1336	modisMeanLogTauTotal
	1337	endif
	1338	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) then
	1339	cospOUT%modis_Optical_Thickness_Water_LogMean(ij+int(modisIN%sunlit(:))-1) = &
	1340	modisMeanLogTauLiquid
	1341	endif
	1342	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) then
	1343	cospOUT%modis_Optical_Thickness_Ice_LogMean(ij+int(modisIN%sunlit(:))-1) = &
	1344	modisMeanLogTauIce
	1345	endif
	1346	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) then
	1347	cospOUT%modis_Cloud_Particle_Size_Water_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1348	modisMeanSizeLiquid
	1349	endif
	1350	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) then
	1351	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1352	modisMeanSizeIce
	1353	endif
	1354	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) then
	1355	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1356	modisMeanCloudTopPressure
	1357	endif
	1358	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) then
	1359	cospOUT%modis_Liquid_Water_Path_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1360	modisMeanLiquidWaterPath
	1361	endif
	1362	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) then
	1363	cospOUT%modis_Ice_Water_Path_Mean(ij+int(modisIN%sunlit(:))-1) = &
	1364	modisMeanIceWaterPath
	1365	endif
	1366	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) then
	1367	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(ij+ &
	1368	int(modisIN%sunlit(:))-1, 1:numModisTauBins, :) = modisJointHistogram(:, :, :)
	1369	! Reorder pressure bins in joint histogram to go from surface to TOA
	1370	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(ij:ik,:,:) = &
	1371	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(ij:ik,:,numMODISPresBins:1:-1)
	1372	endif
	1373	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffIce)) then
	1374	cospOUT%modis_Optical_Thickness_vs_ReffIce(ij+int(modisIN%sunlit(:))-1, 1:numMODISTauBins,:) = &
	1375	modisJointHistogramIce(:,:,:)
	1376	endif
	1377	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLiq)) then
	1378	cospOUT%modis_Optical_Thickness_vs_ReffLiq(ij+int(modisIN%sunlit(:))-1, 1:numMODISTauBins,:) = &
	1379	modisJointHistogramLiq(:,:,:)
	1380	endif
	1381
	1382	if(modisIN%nSunlit < modisIN%Npoints) then
	1383	! Where it's night and we haven't done the retrievals the values are undefined
	1384	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	1385	cospOUT%modis_Cloud_Fraction_Total_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1386	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	1387	cospOUT%modis_Cloud_Fraction_Water_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1388	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	1389	cospOUT%modis_Cloud_Fraction_Ice_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1390	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	1391	cospOUT%modis_Cloud_Fraction_High_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1392	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	1393	cospOUT%modis_Cloud_Fraction_Mid_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1394	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	1395	cospOUT%modis_Cloud_Fraction_Low_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1396	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	1397	cospOUT%modis_Optical_Thickness_Total_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1398	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	1399	cospOUT%modis_Optical_Thickness_Water_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1400	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	1401	cospOUT%modis_Optical_Thickness_Ice_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1402	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	1403	cospOUT%modis_Optical_Thickness_Total_LogMean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1404	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	1405	cospOUT%modis_Optical_Thickness_Water_LogMean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1406	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	1407	cospOUT%modis_Optical_Thickness_Ice_LogMean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1408	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	1409	cospOUT%modis_Cloud_Particle_Size_Water_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1410	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	1411	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1412	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	1413	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1414	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	1415	cospOUT%modis_Liquid_Water_Path_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1416	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	1417	cospOUT%modis_Ice_Water_Path_Mean(ij+int(modisIN%notSunlit(:))-1) = R_UNDEF
	1418	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	1419	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(ij+int(modisIN%notSunlit(:))-1, :, :) = R_UNDEF
	1420	end if
	1421	else
	1422	! It's nightime everywhere - everything is undefined
	1423	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	1424	cospOUT%modis_Cloud_Fraction_Total_Mean(ij:ik) = R_UNDEF
	1425	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	1426	cospOUT%modis_Cloud_Fraction_Water_Mean(ij:ik) = R_UNDEF
	1427	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	1428	cospOUT%modis_Cloud_Fraction_Ice_Mean(ij:ik) = R_UNDEF
	1429	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	1430	cospOUT%modis_Cloud_Fraction_High_Mean(ij:ik) = R_UNDEF
	1431	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	1432	cospOUT%modis_Cloud_Fraction_Mid_Mean(ij:ik) = R_UNDEF
	1433	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	1434	cospOUT%modis_Cloud_Fraction_Low_Mean(ij:ik) = R_UNDEF
	1435	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	1436	cospOUT%modis_Optical_Thickness_Total_Mean(ij:ik) = R_UNDEF
	1437	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	1438	cospOUT%modis_Optical_Thickness_Water_Mean(ij:ik) = R_UNDEF
	1439	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	1440	cospOUT%modis_Optical_Thickness_Ice_Mean(ij:ik) = R_UNDEF
	1441	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	1442	cospOUT%modis_Optical_Thickness_Total_LogMean(ij:ik) = R_UNDEF
	1443	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	1444	cospOUT%modis_Optical_Thickness_Water_LogMean(ij:ik) = R_UNDEF
	1445	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	1446	cospOUT%modis_Optical_Thickness_Ice_LogMean(ij:ik) = R_UNDEF
	1447	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	1448	cospOUT%modis_Cloud_Particle_Size_Water_Mean(ij:ik) = R_UNDEF
	1449	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	1450	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(ij:ik) = R_UNDEF
	1451	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	1452	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(ij:ik) = R_UNDEF
	1453	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	1454	cospOUT%modis_Liquid_Water_Path_Mean(ij:ik) = R_UNDEF
	1455	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	1456	cospOUT%modis_Ice_Water_Path_Mean(ij:ik) = R_UNDEF
	1457	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	1458	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(ij:ik, :, :) = R_UNDEF
	1459	endif
	1460	! Free up memory (if necessary)
	1461	if (allocated(modisRetrievedTau)) deallocate(modisRetrievedTau)
	1462	if (allocated(modisRetrievedSize)) deallocate(modisRetrievedSize)
	1463	if (allocated(modisRetrievedPhase)) deallocate(modisRetrievedPhase)
	1464	if (allocated(modisRetrievedCloudTopPressure)) deallocate(modisRetrievedCloudTopPressure)
	1465	if (allocated(modisCftotal)) deallocate(modisCftotal)
	1466	if (allocated(modisCfLiquid)) deallocate(modisCfLiquid)
	1467	if (allocated(modisCfIce)) deallocate(modisCfIce)
	1468	if (allocated(modisCfHigh)) deallocate(modisCfHigh)
	1469	if (allocated(modisCfMid)) deallocate(modisCfMid)
	1470	if (allocated(modisCfLow)) deallocate(modisCfLow)
	1471	if (allocated(modisMeanTauTotal)) deallocate(modisMeanTauTotal)
	1472	if (allocated(modisMeanTauLiquid)) deallocate(modisMeanTauLiquid)
	1473	if (allocated(modisMeanTauIce)) deallocate(modisMeanTauIce)
	1474	if (allocated(modisMeanLogTauTotal)) deallocate(modisMeanLogTauTotal)
	1475	if (allocated(modisMeanLogTauLiquid)) deallocate(modisMeanLogTauLiquid)
	1476	if (allocated(modisMeanLogTauIce)) deallocate(modisMeanLogTauIce)
	1477	if (allocated(modisMeanSizeLiquid)) deallocate(modisMeanSizeLiquid)
	1478	if (allocated(modisMeanSizeIce)) deallocate(modisMeanSizeIce)
	1479	if (allocated(modisMeanCloudTopPressure)) deallocate(modisMeanCloudTopPressure)
	1480	if (allocated(modisMeanLiquidWaterPath)) deallocate(modisMeanLiquidWaterPath)
	1481	if (allocated(modisMeanIceWaterPath)) deallocate(modisMeanIceWaterPath)
	1482	if (allocated(modisJointHistogram)) deallocate(modisJointHistogram)
	1483	if (allocated(modisJointHistogramIce)) deallocate(modisJointHistogramIce)
	1484	if (allocated(modisJointHistogramLiq)) deallocate(modisJointHistogramLiq)
	1485	if (allocated(isccp_boxttop)) deallocate(isccp_boxttop)
	1486	if (allocated(isccp_boxptop)) deallocate(isccp_boxptop)
	1487	if (allocated(isccp_boxtau)) deallocate(isccp_boxtau)
	1488	if (allocated(isccp_meantbclr)) deallocate(isccp_meantbclr)
	1489	if (allocated(isccpLEVMATCH)) deallocate(isccpLEVMATCH)
	1490	endif
	1491
	1492	! RTTOV
	1493	! if (lrttov_column) then
	1494	! call rttov_column(rttovIN%nPoints,rttovIN%nLevels,rttovIN%nSubCols,rttovIN%q, &
	1495	! rttovIN%p,rttovIN%t,rttovIN%o3,rttovIN%ph,rttovIN%h_surf, &
	1496	! rttovIN%u_surf,rttovIN%v_surf,rttovIN%p_surf,rttovIN%t_skin, &
	1497	! rttovIN%t2m,rttovIN%q2m,rttovIN%lsmask,rttovIN%longitude, &
	1498	! rttovIN%latitude,rttovIN%seaice,rttovIN%co2,rttovIN%ch4, &
	1499	! rttovIN%n2o,rttovIN%co,rttovIN%zenang,lrttov_cleanUp, &
	1500	! cospOUT%rttov_tbs(ij:ik,:),cosp_simulator(nError+1), &
	1501	! ! Optional arguments for surface emissivity calculation
	1502	! month=rttovIN%month)
	1503	! ! Optional arguments to rttov for all-sky calculation
	1504	! ! rttovIN%month, rttovIN%tca,rttovIN%cldIce,rttovIN%cldLiq, &
	1505	! ! rttovIN%fl_rain,rttovIN%fl_snow)
	1506	! endif
	1507
	1508	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1509	! 6) Compute multi-instrument products
	1510	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1511
	1512	! CLOUDSAT/CALIPSO products
	1513	if (Lradar_lidar_tcc .or. Llidar_only_freq_cloud .or. Lcloudsat_tcc .or. Lcloudsat_tcc2) then
	1514
	1515	if (use_vgrid) then
	1516	allocate(lidar_only_freq_cloud(cloudsatIN%Npoints,Nlvgrid), &
	1517	radar_lidar_tcc(cloudsatIN%Npoints), cloudsat_tcc(cloudsatIN%Npoints), &
	1518	cloudsat_tcc2(cloudsatIN%Npoints))
	1519	allocate(betamol_in(cloudsatIN%Npoints,1,cloudsatIN%Nlevels), &
	1520	betamoli(cloudsatIN%Npoints,1,Nlvgrid), &
	1521	pnormI(cloudsatIN%Npoints,cloudsatIN%Ncolumns,Nlvgrid), &
	1522	Ze_totI(cloudsatIN%Npoints,cloudsatIN%Ncolumns,Nlvgrid))
	1523
	1524	! Regrid in the vertical (NOTE This routine requires SFC-2-TOA ordering, so flip
	1525	! inputs and outputs to maintain TOA-2-SFC ordering convention in COSP2.)
	1526	betamol_in(:,1,:) = calipso_beta_mol(:,cloudsatIN%Nlevels:1:-1)
	1527	call cosp_change_vertical_grid(cloudsatIN%Npoints,1,cloudsatIN%Nlevels, &
	1528	cospgridIN%hgt_matrix(:,cloudsatIN%Nlevels:1:-1), &
	1529	cospgridIN%hgt_matrix_half(:,cloudsatIN%Nlevels:1:-1),betamol_in, &
	1530	Nlvgrid,vgrid_zl(Nlvgrid:1:-1),vgrid_zu(Nlvgrid:1:-1), &
	1531	betamolI(:,1,Nlvgrid:1:-1))
	1532
	1533	call cosp_change_vertical_grid(cloudsatIN%Npoints,cloudsatIN%Ncolumns, &
	1534	cloudsatIN%Nlevels,cospgridIN%hgt_matrix(:,cloudsatIN%Nlevels:1:-1), &
	1535	cospgridIN%hgt_matrix_half(:,cloudsatIN%Nlevels:1:-1), &
	1536	calipso_beta_tot(:,:,cloudsatIN%Nlevels:1:-1),Nlvgrid, &
	1537	vgrid_zl(Nlvgrid:1:-1),vgrid_zu(Nlvgrid:1:-1),pnormI(:,:,Nlvgrid:1:-1))
	1538
	1539	call cosp_change_vertical_grid(cloudsatIN%Npoints,cloudsatIN%Ncolumns, &
	1540	cloudsatIN%Nlevels,cospgridIN%hgt_matrix(:,cloudsatIN%Nlevels:1:-1), &
	1541	cospgridIN%hgt_matrix_half(:,cloudsatIN%Nlevels:1:-1), &
	1542	cloudsatDBZe(:,:,cloudsatIN%Nlevels:1:-1),Nlvgrid,vgrid_zl(Nlvgrid:1:-1), &
	1543	vgrid_zu(Nlvgrid:1:-1),Ze_totI(:,:,Nlvgrid:1:-1),log_units=.true.)
	1544
	1545	call cosp_lidar_only_cloud(cloudsatIN%Npoints, cloudsatIN%Ncolumns, Nlvgrid, &
	1546	pnormI, betamolI, Ze_totI, lidar_only_freq_cloud, radar_lidar_tcc, &
	1547	cloudsat_tcc, cloudsat_tcc2)
	1548
	1549	deallocate(betamol_in,betamolI,pnormI,ze_totI)
	1550	else
	1551	allocate(lidar_only_freq_cloud(cloudsatIN%Npoints,cloudsatIN%Nlevels), &
	1552	radar_lidar_tcc(cloudsatIN%Npoints), cloudsat_tcc(cloudsatIN%Npoints), &
	1553	cloudsat_tcc2(cloudsatIN%Npoints))
	1554	call cosp_lidar_only_cloud(cloudsatIN%Npoints,cloudsatIN%Ncolumns, &
	1555	cospIN%Nlevels,calipso_beta_tot(:,:,cloudsatIN%Nlevels:1:-1), &
	1556	calipso_beta_mol(:,cloudsatIN%Nlevels:1:-1), &
	1557	cloudsatDBZe(:,:,cloudsatIN%Nlevels:1:-1),lidar_only_freq_cloud, &
	1558	radar_lidar_tcc, cloudsat_tcc, cloudsat_tcc2)
	1559	endif
	1560
	1561	! Store, when necessary
	1562	if (associated(cospOUT%lidar_only_freq_cloud)) then
	1563	cospOUT%lidar_only_freq_cloud(ij:ik,:) = lidar_only_freq_cloud
	1564	endif
	1565	if (associated(cospOUT%radar_lidar_tcc)) then
	1566	cospOUT%radar_lidar_tcc(ij:ik) = radar_lidar_tcc
	1567	endif
	1568	if (associated(cospOUT%cloudsat_tcc)) then
	1569	cospOUT%cloudsat_tcc(ij:ik) = cloudsat_tcc
	1570	endif
	1571	if (associated(cospOUT%cloudsat_tcc2)) then
	1572	cospOUT%cloudsat_tcc2(ij:ik) = cloudsat_tcc2
	1573	endif
	1574	endif
	1575
	1576	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1577	! 7) Cleanup
	1578	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1579	if (Lisccp_subcolumn .or. Lmodis_subcolumn) then
	1580	nullify(isccpIN%Npoints,isccpIN%Ncolumns,isccpIN%Nlevels,isccpIN%emsfc_lw, &
	1581	isccpIN%skt,isccpIN%qv,isccpIN%at,isccpIN%frac_out,isccpIN%dtau, &
	1582	isccpIN%dem,isccpIN%phalf,isccpIN%sunlit,isccpIN%pfull)
	1583	endif
	1584
	1585	if (Lmisr_subcolumn) then
	1586	nullify(misrIN%Npoints,misrIN%Ncolumns,misrIN%Nlevels,misrIN%dtau,misrIN%sunlit, &
	1587	misrIN%zfull,misrIN%at)
	1588	endif
	1589
	1590	if (Lcalipso_subcolumn) then
	1591	nullify(calipsoIN%Npoints,calipsoIN%Ncolumns,calipsoIN%Nlevels,calipsoIN%beta_mol,&
	1592	calipsoIN%betatot,calipsoIN%betatot_liq,calipsoIN%betatot_ice, &
	1593	calipsoIN%tau_mol,calipsoIN%tautot,calipsoIN%tautot_liq,calipsoIN%tautot_ice)
	1594	endif
	1595
	1596	if (LgrLidar532_subcolumn) then
	1597	nullify(grLidar532IN%Npoints,grLidar532IN%Ncolumns,grLidar532IN%Nlevels,grLidar532IN%beta_mol, &
	1598	grLidar532IN%betatot,grLidar532IN%tau_mol,grLidar532IN%tautot)
	1599	endif
	1600
	1601	if (Latlid_subcolumn) then
	1602	nullify(atlidIN%Npoints,atlidIN%Ncolumns,atlidIN%Nlevels,atlidIN%beta_mol_atlid, &
	1603	atlidIN%betatot_atlid,atlidIN%tau_mol_atlid,atlidIN%tautot_atlid)
	1604	endif
	1605
	1606	if (Lparasol_subcolumn) then
	1607	nullify(parasolIN%Npoints,parasolIN%Nlevels,parasolIN%Ncolumns,parasolIN%Nrefl, &
	1608	parasolIN%tautot_S_liq,parasolIN%tautot_S_ice)
	1609	endif
	1610
	1611
	1612	if (Lcloudsat_subcolumn) then
	1613	nullify(cloudsatIN%Npoints,cloudsatIN%Nlevels,cloudsatIN%Ncolumns,cloudsatIN%rcfg,&
	1614	cloudsatIN%kr_vol,cloudsatIN%g_vol,cloudsatIN%z_vol,cloudsatIN%hgt_matrix)
	1615	endif
	1616
	1617	if (Lmodis_subcolumn) then
	1618	nullify(modisIN%Npoints,modisIN%Ncolumns,modisIN%Nlevels,modisIN%tau,modisIN%g, &
	1619	modisIN%liqFrac,modisIN%w0)
	1620	if (allocated(modisIN%sunlit)) deallocate(modisIN%sunlit)
	1621	if (allocated(modisIN%notSunlit)) deallocate(modisIN%notSunlit)
	1622	if (allocated(modisIN%pres)) deallocate(modisIN%pres)
	1623	endif
	1624
	1625	if (allocated(calipso_beta_tot)) deallocate(calipso_beta_tot)
	1626	if (allocated(grLidar532_beta_tot)) deallocate(grLidar532_beta_tot)
	1627	if (allocated(atlid_beta_tot)) deallocate(atlid_beta_tot)
	1628	if (allocated(calipso_beta_mol)) deallocate(calipso_beta_mol)
	1629	if (allocated(grLidar532_beta_mol)) deallocate(grLidar532_beta_mol)
	1630	if (allocated(atlid_beta_mol)) deallocate(atlid_beta_mol)
	1631	if (allocated(calipso_betaperp_tot)) deallocate(calipso_betaperp_tot)
	1632	if (allocated(cloudsatDBZe)) deallocate(cloudsatDBZe)
	1633	if (allocated(lidar_only_freq_cloud)) deallocate(lidar_only_freq_cloud)
	1634	if (allocated(radar_lidar_tcc)) deallocate(radar_lidar_tcc)
	1635	if (allocated(cloudsat_tcc)) deallocate(cloudsat_tcc)
	1636	if (allocated(cloudsat_tcc2)) deallocate(cloudsat_tcc2)
	1637
	1638	end function COSP_SIMULATOR
	1639	! ######################################################################################
	1640	! SUBROUTINE cosp_init
	1641	! ######################################################################################
	1642	SUBROUTINE COSP_INIT(Lisccp, Lmodis, Lmisr, Lcloudsat, Lcalipso, LgrLidar532, Latlid, Lparasol, Lrttov, &
	1643	cloudsat_radar_freq, cloudsat_k2, cloudsat_use_gas_abs, cloudsat_do_ray, &
	1644	isccp_top_height, isccp_top_height_direction, surface_radar, rcfg, lusevgrid, &
	1645	luseCSATvgrid, Nvgrid, Nlevels, cloudsat_micro_scheme)
	1646
	1647	! INPUTS
	1648	logical,intent(in) :: Lisccp,Lmodis,Lmisr,Lcloudsat,Lcalipso,LgrLidar532,Latlid,Lparasol,Lrttov
	1649	integer,intent(in) :: &
	1650	cloudsat_use_gas_abs, & !
	1651	cloudsat_do_ray, & !
	1652	isccp_top_height, & !
	1653	isccp_top_height_direction, & !
	1654	Nlevels, & !
	1655	Nvgrid, & ! Number of levels for new L3 grid
	1656	surface_radar !
	1657	real(wp),intent(in) :: &
	1658	cloudsat_radar_freq, & !
	1659	cloudsat_k2 !
	1660	logical,intent(in) :: &
	1661	lusevgrid, & ! Switch to use different vertical grid
	1662	luseCSATvgrid ! Switch to use CLOUDSAT grid spacing for new
	1663	! vertical grid
	1664	character(len=64),intent(in) :: &
	1665	cloudsat_micro_scheme ! Microphysical scheme used by CLOUDSAT
	1666
	1667	! OUTPUTS
	1668	type(radar_cfg) :: rcfg
	1669
	1670	! Local variables
	1671	integer :: i
	1672	real(wp) :: zstep
	1673
	1674	! Initialize MODIS optical-depth bin boundaries for joint-histogram. (defined in cosp_config.F90)
	1675	if (.not. allocated(modis_histTau)) then
	1676	allocate(modis_histTau(ntau+1),modis_histTauEdges(2,ntau),modis_histTauCenters(ntau))
	1677	numMODISTauBins = ntau
	1678	modis_histTau = tau_binBounds
	1679	modis_histTauEdges = tau_binEdges
	1680	modis_histTauCenters = tau_binCenters
	1681	endif
	1682
	1683	! Set up vertical grid used by CALIPSO and CLOUDSAT L3
	1684	use_vgrid = lusevgrid
	1685
	1686	if (use_vgrid) then
	1687	Nlvgrid = Nvgrid
	1688	allocate(vgrid_zl(Nlvgrid),vgrid_zu(Nlvgrid),vgrid_z(Nlvgrid))
	1689	! CloudSat grid requested
	1690	if (luseCSATvgrid) zstep = 480._wp
	1691	! Other grid requested. Constant vertical spacing with top at 20 km
	1692	if (.not. luseCSATvgrid) zstep = 20000._wp/Nvgrid
	1693	do i=1,Nvgrid
	1694	vgrid_zl(Nlvgrid-i+1) = (i-1)*zstep
	1695	vgrid_zu(Nlvgrid-i+1) = i*zstep
	1696	enddo
	1697	vgrid_z = (vgrid_zl+vgrid_zu)/2._wp
	1698	else
	1699	Nlvgrid = Nlevels
	1700	allocate(vgrid_zl(Nlvgrid),vgrid_zu(Nlvgrid),vgrid_z(Nlvgrid))
	1701	endif
	1702
	1703	! Initialize simulators
	1704	if (Lisccp) call cosp_isccp_init(isccp_top_height,isccp_top_height_direction)
	1705	if (Lmodis) call cosp_modis_init()
	1706	if (Lmisr) call cosp_misr_init()
	1707	! if (Lrttov) call cosp_rttov_init()
	1708	if (Lcloudsat) call cosp_cloudsat_init(cloudsat_radar_freq,cloudsat_k2, &
	1709	cloudsat_use_gas_abs,cloudsat_do_ray,R_UNDEF,N_HYDRO, surface_radar, &
	1710	rcfg,cloudsat_micro_scheme)
	1711	if (Lcalipso) call cosp_calipso_init()
	1712	if (LgrLidar532) call cosp_grLidar532_init()
	1713	if (Latlid) call cosp_atlid_init()
	1714	if (Lparasol) call cosp_parasol_init()
	1715
	1716	linitialization = .FALSE.
	1717	END SUBROUTINE COSP_INIT
	1718
	1719	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1720	! SUBROUTINE cosp_cleanUp
	1721	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1722	subroutine cosp_cleanUp()
	1723	deallocate(vgrid_zl,vgrid_zu,vgrid_z)
	1724	end subroutine cosp_cleanUp
	1725
	1726	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1727	! SUBROUTINE cosp_errorCheck
	1728	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1729	subroutine cosp_errorCheck(cospgridIN, cospIN, Lisccp_subcolumn, Lisccp_column, &
	1730	Lmisr_subcolumn, Lmisr_column, Lmodis_subcolumn, Lmodis_column, Lcloudsat_subcolumn, &
	1731	Lcloudsat_column, Lcalipso_subcolumn, Lcalipso_column, Latlid_subcolumn, &
	1732	Latlid_column, LgrLidar532_subcolumn, LgrLidar532_column, Lrttov_subcolumn, &
	1733	Lrttov_column, Lparasol_subcolumn, Lparasol_column, Lradar_lidar_tcc, &
	1734	Llidar_only_freq_cloud, Lcloudsat_tcc, Lcloudsat_tcc2, cospOUT, errorMessage, nError)
	1735
	1736	! Inputs
	1737	type(cosp_column_inputs),intent(in) :: &
	1738	cospgridIN ! Model grid inputs to COSP
	1739	type(cosp_optical_inputs),intent(in) :: &
	1740	cospIN ! Derived (optical) input to COSP
	1741
	1742	! Outputs
	1743	logical,intent(inout) :: &
	1744	Lisccp_subcolumn, & ! ISCCP subcolumn simulator on/off switch
	1745	Lisccp_column, & ! ISCCP column simulator on/off switch
	1746	Lmisr_subcolumn, & ! MISR subcolumn simulator on/off switch
	1747	Lmisr_column, & ! MISR column simulator on/off switch
	1748	Lmodis_subcolumn, & ! MODIS subcolumn simulator on/off switch
	1749	Lmodis_column, & ! MODIS column simulator on/off switch
	1750	Lcloudsat_subcolumn, & ! CLOUDSAT subcolumn simulator on/off switch
	1751	Lcloudsat_column, & ! CLOUDSAT column simulator on/off switch
	1752	Lcalipso_subcolumn, & ! CALIPSO subcolumn simulator on/off switch
	1753	Lcalipso_column, & ! CALIPSO column simulator on/off switch
	1754	Latlid_subcolumn, & ! EarthCare subcolumn simulator on/off switch
	1755	Latlid_column, & ! EarthCare column simulator on/off switch
	1756	LgrLidar532_subcolumn, & ! Ground Lidar subcolumn simulator on/off switch
	1757	LgrLidar532_column, & ! Ground Lidar column simulator on/off switch
	1758	Lparasol_subcolumn, & ! PARASOL subcolumn simulator on/off switch
	1759	Lparasol_column, & ! PARASOL column simulator on/off switch
	1760	Lrttov_subcolumn, & ! RTTOV subcolumn simulator on/off switch
	1761	Lrttov_column, & ! RTTOV column simulator on/off switch
	1762	Lcloudsat_tcc, & !
	1763	Lcloudsat_tcc2, & !
	1764	Lradar_lidar_tcc, & ! On/Off switch for joint Calipso/Cloudsat product
	1765	Llidar_only_freq_cloud ! On/Off switch for joint Calipso/Cloudsat product
	1766	type(cosp_outputs),intent(inout) :: &
	1767	cospOUT ! COSP Outputs
	1768	character(len=256),dimension(100) :: errorMessage
	1769	integer,intent(out) :: nError
	1770
	1771	! Local variables
	1772	character(len=100) :: parasolErrorMessage
	1773	logical :: alloc_status
	1774
	1775	nError = 0
	1776	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1777	! PART 0: Ensure that the inputs needed by the requested simulators are allocated.
	1778	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	1779	! ISCCP simulator
	1780	if (Lisccp_subcolumn .or. Lisccp_column) then
	1781	alloc_status = .true.
	1782	if (.not. allocated(cospgridIN%skt)) then
	1783	nError=nError+1
	1784	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%skt has not been allocated'
	1785	alloc_status = .false.
	1786	endif
	1787	if (.not. allocated(cospgridIN%qv)) then
	1788	nError=nError+1
	1789	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%qv has not been allocated'
	1790	alloc_status = .false.
	1791	endif
	1792	if (.not. allocated(cospgridIN%at)) then
	1793	nError=nError+1
	1794	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%at has not been allocated'
	1795	alloc_status = .false.
	1796	endif
	1797	if (.not. allocated(cospIN%frac_out)) then
	1798	nError=nError+1
	1799	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospIN%frac_out has not been allocated'
	1800	alloc_status = .false.
	1801	endif
	1802	if (.not. allocated(cospIN%tau_067)) then
	1803	nError=nError+1
	1804	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospIN%tau_067 has not been allocated'
	1805	alloc_status = .false.
	1806	endif
	1807	if (.not. allocated(cospIN%emiss_11)) then
	1808	nError=nError+1
	1809	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospIN%emiss_11 has not been allocated'
	1810	alloc_status = .false.
	1811	endif
	1812	if (.not. allocated(cospgridIN%phalf)) then
	1813	nError=nError+1
	1814	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%phalf has not been allocated'
	1815	alloc_status = .false.
	1816	endif
	1817	if (.not. allocated(cospgridIN%sunlit)) then
	1818	nError=nError+1
	1819	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%sunlit has not been allocated'
	1820	alloc_status = .false.
	1821	endif
	1822	if (.not. allocated(cospgridIN%pfull)) then
	1823	nError=nError+1
	1824	errorMessage(nError) = 'ERROR: COSP input variable (ISSCP simulator): cospgridIN%pfull has not been allocated'
	1825	alloc_status = .false.
	1826	endif
	1827	if (.not. alloc_status) then
	1828	Lisccp_subcolumn = .false.
	1829	Lisccp_column = .false.
	1830	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	1831	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	1832	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	1833	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	1834	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	1835	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	1836	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	1837	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	1838	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	1839	endif
	1840	endif
	1841
	1842	! MISR simulator
	1843	if (Lmisr_subcolumn .or. Lmisr_column) then
	1844	alloc_status = .true.
	1845	if (.not. allocated(cospIN%tau_067)) then
	1846	nError=nError+1
	1847	errorMessage(nError) = 'ERROR: COSP input variable (MISR simulator): cospIN%tau_067 has not been allocated'
	1848	alloc_status = .false.
	1849	endif
	1850	if (.not. allocated(cospgridIN%sunlit)) then
	1851	nError=nError+1
	1852	errorMessage(nError) = 'ERROR: COSP input variable (MISR simulator): cospgridIN%sunlit has not been allocated'
	1853	alloc_status = .false.
	1854	endif
	1855	if (.not. allocated(cospgridIN%hgt_matrix)) then
	1856	nError=nError+1
	1857	errorMessage(nError) = 'ERROR: COSP input variable (MISR simulator): cospgridIN%hgt_matrix has not been allocated'
	1858	alloc_status = .false.
	1859	endif
	1860	if (.not. allocated(cospgridIN%at)) then
	1861	nError=nError+1
	1862	errorMessage(nError) = 'ERROR: COSP input variable (MISR simulator): cospgridIN%at has not been allocated'
	1863	alloc_status = .false.
	1864	endif
	1865	if (.not. alloc_status) then
	1866	Lmisr_subcolumn = .false.
	1867	Lmisr_column = .false.
	1868	if (associated(cospOUT%misr_fq)) cospOUT%misr_fq(:,:,:) = R_UNDEF
	1869	if (associated(cospOUT%misr_dist_model_layertops)) cospOUT%misr_dist_model_layertops(:,:) = R_UNDEF
	1870	if (associated(cospOUT%misr_meanztop)) cospOUT%misr_meanztop(:) = R_UNDEF
	1871	if (associated(cospOUT%misr_cldarea)) cospOUT%misr_cldarea(:) = R_UNDEF
	1872	endif
	1873	endif
	1874
	1875	! EarthCare Lidar simulator.
	1876	if (Latlid_subcolumn .or. Latlid_column) then
	1877	alloc_status = .true.
	1878	if (.not. allocated(cospIN%beta_mol_atlid)) then
	1879	nError=nError+1
	1880	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospIN%beta_mol_atlid has not been allocated'
	1881	alloc_status = .false.
	1882	endif
	1883	if (.not. allocated(cospIN%betatot_atlid)) then
	1884	nError=nError+1
	1885	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospIN%betatot_atlid has not been allocated'
	1886	alloc_status = .false.
	1887	endif
	1888	if (.not. allocated(cospIN%tau_mol_atlid)) then
	1889	nError=nError+1
	1890	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospIN%tau_mol_atlid has not been allocated'
	1891	alloc_status = .false.
	1892	endif
	1893	if (.not. allocated(cospIN%tautot_atlid)) then
	1894	nError=nError+1
	1895	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospIN%tautot_atlid has not been allocated'
	1896	alloc_status = .false.
	1897	endif
	1898	if (.not. alloc_status) then
	1899	Latlid_subcolumn = .false.
	1900	Latlid_column = .false.
	1901	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	1902	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	1903	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	1904	if (associated(cospOUT%atlid_beta_mol)) cospOUT%atlid_beta_mol(:,:) = R_UNDEF
	1905	if (associated(cospOUT%atlid_beta_tot)) cospOUT%atlid_beta_tot(:,:,:) = R_UNDEF
	1906	endif
	1907
	1908	! EarthCare column simulator requires additional inputs not required by the subcolumn simulator.
	1909	alloc_status = .true.
	1910	if (.not. allocated(cospgridIN%hgt_matrix)) then
	1911	nError=nError+1
	1912	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospgridIN%hgt_matrix has not been allocated'
	1913	alloc_status = .false.
	1914	endif
	1915	if (.not. allocated(cospgridIN%hgt_matrix_half)) then
	1916	nError=nError+1
	1917	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospgridIN%hgt_matrix_half has not been allocated'
	1918	alloc_status = .false.
	1919	endif
	1920	if (.not. allocated(cospgridIN%phalf)) then
	1921	nError=nError+1
	1922	errorMessage(nError) = 'ERROR: COSP input variable (EarthCare Lidar simulator): cospgridIN%phalf has not been allocated'
	1923	alloc_status = .false.
	1924	endif
	1925	if (.not. alloc_status) then
	1926	Latlid_column = .false.
	1927	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	1928	endif
	1929	endif
	1930
	1931	! Ground Lidar simulator.
	1932	if (LgrLidar532_subcolumn .or. LgrLidar532_column) then
	1933	alloc_status = .true.
	1934	if (.not. allocated(cospIN%beta_mol_grLidar532)) then
	1935	nError=nError+1
	1936	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospIN%beta_mol_grLidar532 has not been allocated'
	1937	alloc_status = .false.
	1938	endif
	1939	if (.not. allocated(cospIN%betatot_grLidar532)) then
	1940	nError=nError+1
	1941	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospIN%betatot_grLidar532 has not been allocated'
	1942	alloc_status = .false.
	1943	endif
	1944	if (.not. allocated(cospIN%tau_mol_grLidar532)) then
	1945	nError=nError+1
	1946	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospIN%tau_mol_grLidar532 has not been allocated'
	1947	alloc_status = .false.
	1948	endif
	1949	if (.not. allocated(cospIN%tautot_grLidar532)) then
	1950	nError=nError+1
	1951	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospIN%tautot_grLidar532 has not been allocated'
	1952	alloc_status = .false.
	1953	endif
	1954	if (.not. alloc_status) then
	1955	LgrLidar532_subcolumn = .false.
	1956	LgrLidar532_column = .false.
	1957	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	1958	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	1959	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	1960	if (associated(cospOUT%grLidar532_beta_mol)) cospOUT%grLidar532_beta_mol(:,:) = R_UNDEF
	1961	if (associated(cospOUT%grLidar532_beta_tot)) cospOUT%grLidar532_beta_tot(:,:,:) = R_UNDEF
	1962	endif
	1963
	1964	! Ground Lidar column simulator requires additional inputs not required by the subcolumn simulator.
	1965	alloc_status = .true.
	1966	if (.not. allocated(cospgridIN%hgt_matrix)) then
	1967	nError=nError+1
	1968	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospgridIN%hgt_matrix has not been allocated'
	1969	alloc_status = .false.
	1970	endif
	1971	if (.not. allocated(cospgridIN%hgt_matrix_half)) then
	1972	nError=nError+1
	1973	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospgridIN%hgt_matrix_half has not been allocated'
	1974	alloc_status = .false.
	1975	endif
	1976	if (.not. allocated(cospgridIN%phalf)) then
	1977	nError=nError+1
	1978	errorMessage(nError) = 'ERROR: COSP input variable (Ground Lidar simulator): cospgridIN%phalf has not been allocated'
	1979	alloc_status = .false.
	1980	endif
	1981	if (.not. alloc_status) then
	1982	LgrLidar532_column = .false.
	1983	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	1984	endif
	1985	endif
	1986
	1987	! Calipso Lidar simulator
	1988	if (Lcalipso_subcolumn .or. Lcalipso_column) then
	1989	alloc_status = .true.
	1990	if (.not. allocated(cospIN%beta_mol_calipso)) then
	1991	nError=nError+1
	1992	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospIN%beta_mol_calipso has not been allocated'
	1993	alloc_status = .false.
	1994	endif
	1995	if (.not. allocated(cospIN%betatot_calipso)) then
	1996	nError=nError+1
	1997	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospIN%betatot_calipso has not been allocated'
	1998	alloc_status = .false.
	1999	endif
	2000	if (.not. allocated(cospIN%betatot_liq_calipso)) then
	2001	nError=nError+1
	2002	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2003	' cospIN%betatot_liq_calipso has not been allocated'
	2004	alloc_status = .false.
	2005	endif
	2006	if (.not. allocated(cospIN%betatot_ice_calipso)) then
	2007	nError=nError+1
	2008	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2009	' cospIN%betatot_ice_calipso has not been allocated'
	2010	alloc_status = .false.
	2011	endif
	2012	if (.not. allocated(cospIN%tau_mol_calipso)) then
	2013	nError=nError+1
	2014	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospIN%tau_mol_calipso has not been allocated'
	2015	alloc_status = .false.
	2016	endif
	2017	if (.not. allocated(cospIN%tautot_calipso)) then
	2018	nError=nError+1
	2019	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospIN%tautot_calipso has not been allocated'
	2020	alloc_status = .false.
	2021	endif
	2022	if (.not. allocated(cospIN%tautot_liq_calipso)) then
	2023	nError=nError+1
	2024	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2025	' cospIN%tautot_liq has not been allocated'
	2026	alloc_status = .false.
	2027	endif
	2028	if (.not. allocated(cospIN%tautot_ice_calipso)) then
	2029	nError=nError+1
	2030	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2031	' cospIN%tautot_ice has not been allocated'
	2032	alloc_status = .false.
	2033	endif
	2034	if (.not. alloc_status) then
	2035	Lcalipso_subcolumn = .false.
	2036	Lcalipso_column = .false.
	2037	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2038	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2039	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2040	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2041	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2042	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2043	if (associated(cospOUT%calipso_beta_mol)) cospOUT%calipso_beta_mol(:,:) = R_UNDEF
	2044	if (associated(cospOUT%calipso_beta_tot)) cospOUT%calipso_beta_tot(:,:,:) = R_UNDEF
	2045	if (associated(cospOUT%calipso_betaperp_tot)) cospOUT%calipso_betaperp_tot(:,:,:) = R_UNDEF
	2046	! Also, turn-off joint-products
	2047	if (Lradar_lidar_tcc) then
	2048	Lradar_lidar_tcc = .false.
	2049	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2050	endif
	2051	if (Lcloudsat_tcc) then
	2052	Lcloudsat_tcc = .false.
	2053	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2054	endif
	2055	if (Lcloudsat_tcc2) then
	2056	Lcloudsat_tcc2 = .false.
	2057	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2058	endif
	2059	if (Llidar_only_freq_cloud) then
	2060	Llidar_only_freq_cloud = .false.
	2061	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2062	endif
	2063	endif
	2064
	2065	! Calipso column simulator requires additional inputs not required by the subcolumn simulator.
	2066	alloc_status = .true.
	2067	if (.not. allocated(cospgridIN%hgt_matrix)) then
	2068	nError=nError+1
	2069	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2070	' cospgridIN%hgt_matrix has not been allocated'
	2071	alloc_status = .false.
	2072	endif
	2073	if (.not. allocated(cospgridIN%hgt_matrix_half)) then
	2074	nError=nError+1
	2075	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator):'//&
	2076	' cospgridIN%hgt_matrix_half has not been allocated'
	2077	alloc_status = .false.
	2078	endif
	2079	if (.not. allocated(cospgridIN%at)) then
	2080	nError=nError+1
	2081	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospgridIN%at has not been allocated'
	2082	alloc_status = .false.
	2083	endif
	2084	if (.not. allocated(cospgridIN%phalf)) then
	2085	nError=nError+1
	2086	errorMessage(nError) = 'ERROR: COSP input variable (Calipso Lidar simulator): cospgridIN%phalf has not been allocated'
	2087	alloc_status = .false.
	2088	endif
	2089	if (.not. alloc_status) then
	2090	Lcalipso_column = .false.
	2091	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2092	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2093	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2094	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2095	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2096	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2097	if (Lcloudsat_tcc) then
	2098	Lcloudsat_tcc = .false.
	2099	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2100	endif
	2101	if (Lcloudsat_tcc2) then
	2102	Lcloudsat_tcc2 = .false.
	2103	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2104	endif
	2105	! Also, turn-off joint-products
	2106	if (Lradar_lidar_tcc) then
	2107	Lradar_lidar_tcc = .false.
	2108	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2109	endif
	2110	if (Llidar_only_freq_cloud) then
	2111	Llidar_only_freq_cloud = .false.
	2112	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2113	endif
	2114	endif
	2115	endif
	2116
	2117	! PARASOL simulator
	2118	if (Lparasol_subcolumn .or. Lparasol_column) then
	2119	alloc_status = .true.
	2120	if (.not. allocated(cospIN%tautot_S_liq)) then
	2121	nError=nError+1
	2122	errorMessage(nError) = 'ERROR: COSP input variable (PARASOL simulator): cospIN%tautot_S_liq has not been allocated'
	2123	alloc_status = .false.
	2124	endif
	2125	if (.not. allocated(cospIN%tautot_S_ice)) then
	2126	nError=nError+1
	2127	errorMessage(nError) = 'ERROR: COSP input variable (PARASOL simulator): cospIN%tautot_S_ice has not been allocated'
	2128	alloc_status = .false.
	2129	endif
	2130	if (.not. alloc_status) then
	2131	Lparasol_subcolumn = .false.
	2132	Lparasol_column = .false.
	2133	if (associated(cospOUT%parasolPix_refl)) cospOUT%parasolPix_refl(:,:,:) = R_UNDEF
	2134	if (associated(cospOUT%parasolGrid_refl)) cospOUT%parasolGrid_refl(:,:) = R_UNDEF
	2135	endif
	2136
	2137	! PARASOL column simulator requires additional inputs not required by the subcolumn simulator.
	2138	alloc_status = .true.
	2139	if (.not. allocated(cospgridIN%land)) then
	2140	nError=nError+1
	2141	errorMessage(nError) = 'ERROR: COSP input variable (PARASOL simulator): cospgridIN%land has not been allocated'
	2142	alloc_status = .false.
	2143	endif
	2144	if (.not. alloc_status) then
	2145	Lparasol_column = .false.
	2146	if (associated(cospOUT%parasolGrid_refl)) cospOUT%parasolGrid_refl(:,:) = R_UNDEF
	2147	endif
	2148	endif
	2149
	2150	! Cloudsat radar simulator
	2151	if (Lcloudsat_subcolumn .or. Lcloudsat_column) then
	2152	alloc_status = .true.
	2153	if (.not. allocated(cospIN%z_vol_cloudsat)) then
	2154	nError=nError+1
	2155	errorMessage(nError) = 'ERROR: COSP input variable (Cloudsat radar simulator):'//&
	2156	' cospIN%z_vol_cloudsat has not been allocated'
	2157	alloc_status = .false.
	2158	endif
	2159	if (.not. allocated(cospIN%kr_vol_cloudsat)) then
	2160	nError=nError+1
	2161	errorMessage(nError) = 'ERROR: COSP input variable (Cloudsat radar simulator):'//&
	2162	' cospIN%kr_vol_cloudsat has not been allocated'
	2163	alloc_status = .false.
	2164	endif
	2165	if (.not. allocated(cospIN%g_vol_cloudsat)) then
	2166	nError=nError+1
	2167	errorMessage(nError) = 'ERROR: COSP input variable (Cloudsat radar simulator):'//&
	2168	' cospIN%g_vol_cloudsat has not been allocated'
	2169	alloc_status = .false.
	2170	endif
	2171	if (.not. allocated(cospgridIN%hgt_matrix)) then
	2172	nError=nError+1
	2173	errorMessage(nError) = 'ERROR: COSP input variable (Cloudsat radar simulator):'//&
	2174	' cospgridIN%hgt_matrix has not been allocated'
	2175	alloc_status = .false.
	2176	endif
	2177	if (.not. alloc_status) then
	2178	Lcloudsat_subcolumn = .false.
	2179	Lcloudsat_column = .false.
	2180	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	2181	if (associated(cospOUT%cloudsat_Ze_tot)) cospOUT%cloudsat_Ze_tot(:,:,:) = R_UNDEF
	2182	if (Lcloudsat_tcc) then
	2183	Lcloudsat_tcc = .false.
	2184	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2185	endif
	2186	if (Lcloudsat_tcc2) then
	2187	Lcloudsat_tcc2 = .false.
	2188	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2189	endif
	2190	! Also, turn-off joint-products
	2191	if (Lradar_lidar_tcc) then
	2192	Lradar_lidar_tcc = .false.
	2193	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2194	endif
	2195	if (Llidar_only_freq_cloud) then
	2196	Llidar_only_freq_cloud = .false.
	2197	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2198	endif
	2199	endif
	2200
	2201	! Cloudsat column simulator requires additional inputs not required by the subcolumn simulator.
	2202	alloc_status = .true.
	2203	if (.not. allocated(cospgridIN%hgt_matrix_half)) then
	2204	nError=nError+1
	2205	errorMessage(nError) = 'ERROR: COSP input variable (Cloudsat radar simulator):'//&
	2206	' cospgridIN%hgt_matrix_half has not been allocated'
	2207	alloc_status = .false.
	2208	endif
	2209	if (.not. alloc_status) then
	2210	Lcloudsat_column = .false.
	2211	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	2212	if (Lcloudsat_tcc) then
	2213	Lcloudsat_tcc = .false.
	2214	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2215	endif
	2216	if (Lcloudsat_tcc2) then
	2217	Lcloudsat_tcc2 = .false.
	2218	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2219	endif
	2220	! Also, turn-off joint-products
	2221	if (Lradar_lidar_tcc) then
	2222	Lradar_lidar_tcc = .false.
	2223	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2224	endif
	2225	if (Llidar_only_freq_cloud) then
	2226	Llidar_only_freq_cloud = .false.
	2227	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2228	endif
	2229	endif
	2230	endif
	2231
	2232	! MODIS simulator
	2233	if (Lmodis_subcolumn .or. Lmodis_column) then
	2234	alloc_status = .true.
	2235	if (.not. allocated(cospIN%fracLiq)) then
	2236	nError=nError+1
	2237	errorMessage(nError) = 'ERROR: COSP input variable (MODIS simulator): cospIN%fracLiq has not been allocated'
	2238	alloc_status = .false.
	2239	endif
	2240	if (.not. allocated(cospIN%tau_067)) then
	2241	nError=nError+1
	2242	errorMessage(nError) = 'ERROR: COSP input variable (MODIS simulator): cospIN%tau_067 has not been allocated'
	2243	alloc_status = .false.
	2244	endif
	2245	if (.not. allocated(cospIN%asym)) then
	2246	nError=nError+1
	2247	errorMessage(nError) = 'ERROR: COSP input variable (MODIS simulator): cospIN%asym has not been allocated'
	2248	alloc_status = .false.
	2249	endif
	2250	if (.not. allocated(cospIN%ss_alb)) then
	2251	nError=nError+1
	2252	errorMessage(nError) = 'ERROR: COSP input variable (MODIS simulator): cospIN%ss_alb has not been allocated'
	2253	alloc_status = .false.
	2254	endif
	2255	if (.not. allocated(cospgridIN%sunlit)) then
	2256	nError=nError+1
	2257	errorMessage(nError) = 'ERROR: COSP input variable (MODIS simulator): cospgridIN%sunlit has not been allocated'
	2258	alloc_status = .false.
	2259	endif
	2260	if (.not. alloc_status) then
	2261	Lmodis_subcolumn = .false.
	2262	Lmodis_column = .false.
	2263	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	2264	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	2265	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	2266	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	2267	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	2268	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	2269	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	2270	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	2271	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	2272	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	2273	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	2274	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	2275	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	2276	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	2277	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	2278	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	2279	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	2280	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	2281	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	2282	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	2283	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	2284	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	2285	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	2286	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	2287	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	2288	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	2289	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	2290	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	2291	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	2292	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	2293	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	2294	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	2295	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	2296	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	2297	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	2298	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	2299	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	2300	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	2301	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	2302	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	2303	endif
	2304	endif
	2305
	2306	! RTTOV
	2307	if (Lrttov_column) then
	2308	alloc_status = .true.
	2309	if (.not. allocated(cospgridIN%emis_sfc)) then
	2310	nError=nError+1
	2311	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%emis_sfc has not been allocated'
	2312	alloc_status = .false.
	2313	endif
	2314	if (.not. allocated(cospgridIN%hgt_matrix_half)) then
	2315	nError=nError+1
	2316	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%emis_sfc has not been allocated'
	2317	alloc_status = .false.
	2318	endif
	2319	if (.not. allocated(cospgridIN%u_sfc)) then
	2320	nError=nError+1
	2321	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%u_sfc has not been allocated'
	2322	alloc_status = .false.
	2323	endif
	2324	if (.not. allocated(cospgridIN%v_sfc)) then
	2325	nError=nError+1
	2326	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%v_sfc has not been allocated'
	2327	alloc_status = .false.
	2328	endif
	2329	if (.not. allocated(cospgridIN%skt)) then
	2330	nError=nError+1
	2331	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%skt has not been allocated'
	2332	alloc_status = .false.
	2333	endif
	2334	if (.not. allocated(cospgridIN%phalf)) then
	2335	nError=nError+1
	2336	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%phalf has not been allocated'
	2337	alloc_status = .false.
	2338	endif
	2339	if (.not. allocated(cospgridIN%qv)) then
	2340	nError=nError+1
	2341	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%qv has not been allocated'
	2342	alloc_status = .false.
	2343	endif
	2344	if (.not. allocated(cospgridIN%at)) then
	2345	nError=nError+1
	2346	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%at has not been allocated'
	2347	alloc_status = .false.
	2348	endif
	2349	if (.not. allocated(cospgridIN%land)) then
	2350	nError=nError+1
	2351	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%land has not been allocated'
	2352	alloc_status = .false.
	2353	endif
	2354	if (.not. allocated(cospgridIN%lat)) then
	2355	nError=nError+1
	2356	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%lat has not been allocated'
	2357	alloc_status = .false.
	2358	endif
	2359	if (.not. allocated(cospgridIN%lon)) then
	2360	nError=nError+1
	2361	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%lon has not been allocated'
	2362	alloc_status = .false.
	2363	endif
	2364	if (.not. allocated(cospgridIN%seaice)) then
	2365	nError=nError+1
	2366	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%seaice has not been allocated'
	2367	alloc_status = .false.
	2368	endif
	2369	if (.not. allocated(cospgridIN%pfull)) then
	2370	nError=nError+1
	2371	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%pfull has not been allocated'
	2372	alloc_status = .false.
	2373	endif
	2374	if (.not. allocated(cospgridIN%phalf)) then
	2375	nError=nError+1
	2376	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%phalf has not been allocated'
	2377	alloc_status = .false.
	2378	endif
	2379	if (.not. allocated(cospgridIN%at)) then
	2380	nError=nError+1
	2381	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%at has not been allocated'
	2382	alloc_status = .false.
	2383	endif
	2384	if (.not. allocated(cospgridIN%qv)) then
	2385	nError=nError+1
	2386	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%qv has not been allocated'
	2387	alloc_status = .false.
	2388	endif
	2389	if (.not. allocated(cospgridIN%o3)) then
	2390	nError=nError+1
	2391	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%o3 has not been allocated'
	2392	alloc_status = .false.
	2393	endif
	2394	if (.not. allocated(cospgridIN%tca)) then
	2395	nError=nError+1
	2396	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%tca has not been allocated'
	2397	alloc_status = .false.
	2398	endif
	2399	if (.not. allocated(cospgridIN%cloudIce)) then
	2400	nError=nError+1
	2401	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%cloudIce has not been allocated'
	2402	alloc_status = .false.
	2403	endif
	2404	if (.not. allocated(cospgridIN%cloudLiq)) then
	2405	nError=nError+1
	2406	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%cloudLiq has not been allocated'
	2407	alloc_status = .false.
	2408	endif
	2409	if (.not. allocated(cospgridIN%fl_rain)) then
	2410	nError=nError+1
	2411	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%fl_rain has not been allocated'
	2412	alloc_status = .false.
	2413	endif
	2414	if (.not. allocated(cospgridIN%fl_snow)) then
	2415	nError=nError+1
	2416	errorMessage(nError) = 'ERROR: COSP input variable (RTTOV): cospgridIN%fl_snow has not been allocated'
	2417	alloc_status = .false.
	2418	endif
	2419	if (.not. alloc_status) then
	2420	Lrttov_column = .false.
	2421	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2422	endif
	2423	endif
	2424
	2425	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	2426	! PART 1: Check input array values for out-of-bounds values. When an out-of-bound value
	2427	! is encountered, COSP outputs that are dependent on that input are filled with
	2428	! an undefined value (set in cosp_config.f90) and if necessary, that simulator
	2429	! is turned off.
	2430	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	2431	if (any([Lisccp_subcolumn, Lisccp_column, Lmisr_subcolumn, Lmisr_column, Lmodis_subcolumn, Lmodis_column])) then
	2432	if (any(cospgridIN%sunlit .lt. 0)) then
	2433	nError=nError+1
	2434	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%sunlit contains values out of range (0 or 1)'
	2435	Lisccp_subcolumn = .false.
	2436	Lisccp_column = .false.
	2437	Lmisr_subcolumn = .false.
	2438	Lmisr_column = .false.
	2439	Lmodis_subcolumn = .false.
	2440	Lmodis_column = .false.
	2441	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2442	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2443	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2444	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2445	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2446	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2447	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2448	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2449	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2450	if (associated(cospOUT%misr_fq)) cospOUT%misr_fq(:,:,:) = R_UNDEF
	2451	if (associated(cospOUT%misr_dist_model_layertops)) cospOUT%misr_dist_model_layertops(:,:) = R_UNDEF
	2452	if (associated(cospOUT%misr_meanztop)) cospOUT%misr_meanztop(:) = R_UNDEF
	2453	if (associated(cospOUT%misr_cldarea)) cospOUT%misr_cldarea(:) = R_UNDEF
	2454	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	2455	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	2456	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	2457	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	2458	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	2459	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	2460	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	2461	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	2462	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	2463	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	2464	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	2465	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	2466	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	2467	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	2468	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	2469	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	2470	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	2471	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	2472	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	2473	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	2474	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	2475	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	2476	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	2477	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	2478	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	2479	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	2480	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	2481	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	2482	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	2483	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	2484	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	2485	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	2486	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	2487	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	2488	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	2489	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	2490	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	2491	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	2492	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	2493	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	2494	endif
	2495	endif
	2496
	2497	if (any([Lisccp_subcolumn, Lisccp_column, Lmisr_subcolumn, Lmisr_column, Lrttov_column,&
	2498	Lcalipso_column, Lcloudsat_column, Lradar_lidar_tcc,Llidar_only_freq_cloud, &
	2499	Lcloudsat_tcc, Lcloudsat_tcc2])) then
	2500	if (any(cospgridIN%at .lt. 0)) then
	2501	nError=nError+1
	2502	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%at contains values out of range (at<0), expected units (K)'
	2503	Lisccp_subcolumn = .false.
	2504	Lisccp_column = .false.
	2505	Lmisr_subcolumn = .false.
	2506	Lmisr_column = .false.
	2507	Lrttov_column = .false.
	2508	Lcalipso_column = .false.
	2509	Lcloudsat_column = .false.
	2510	Lradar_lidar_tcc = .false.
	2511	Llidar_only_freq_cloud = .false.
	2512	Lcloudsat_tcc = .false.
	2513	Lcloudsat_tcc2 = .false.
	2514	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2515	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2516	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2517	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2518	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2519	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2520	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2521	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2522	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2523	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2524	if (associated(cospOUT%misr_fq)) cospOUT%misr_fq(:,:,:) = R_UNDEF
	2525	if (associated(cospOUT%misr_dist_model_layertops)) cospOUT%misr_dist_model_layertops(:,:) = R_UNDEF
	2526	if (associated(cospOUT%misr_meanztop)) cospOUT%misr_meanztop(:) = R_UNDEF
	2527	if (associated(cospOUT%misr_cldarea)) cospOUT%misr_cldarea(:) = R_UNDEF
	2528	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2529	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2530	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2531	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2532	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2533	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2534	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	2535	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	2536	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	2537	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	2538	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	2539	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	2540	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	2541	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2542	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2543	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2544	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2545	endif
	2546	endif
	2547	if (any([Lisccp_subcolumn, Lisccp_column, Lrttov_column])) then
	2548	if (any(cospgridIN%pfull .lt. 0)) then
	2549	nError=nError+1
	2550	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%pfull contains values out of range'
	2551	Lisccp_subcolumn = .false.
	2552	Lisccp_column = .false.
	2553	Lrttov_column = .false.
	2554	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2555	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2556	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2557	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2558	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2559	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2560	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2561	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2562	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2563	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2564	endif
	2565	endif
	2566	if (any([Lisccp_subcolumn,Lisccp_column,Lmodis_subcolumn,Lmodis_column,Lcalipso_column,Lrttov_column,&
	2567	LgrLidar532_column,Latlid_column])) then
	2568	if (any(cospgridIN%phalf .lt. 0)) then
	2569	nError=nError+1
	2570	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%phalf contains values out of range'
	2571	Lisccp_subcolumn = .false.
	2572	Lisccp_column = .false.
	2573	Lmodis_subcolumn = .false.
	2574	Lmodis_column = .false.
	2575	Lcalipso_column = .false.
	2576	Lrttov_column = .false.
	2577	Latlid_column = .false.
	2578	LgrLidar532_column = .false.
	2579	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2580	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2581	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2582	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2583	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2584	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2585	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2586	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2587	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2588	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2589	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	2590	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	2591	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	2592	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	2593	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	2594	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	2595	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	2596	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	2597	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	2598	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	2599	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	2600	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	2601	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	2602	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	2603	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	2604	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	2605	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	2606	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	2607	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	2608	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	2609	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	2610	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	2611	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	2612	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	2613	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	2614	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	2615	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	2616	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	2617	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	2618	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	2619	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	2620	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	2621	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	2622	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	2623	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	2624	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	2625	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	2626	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	2627	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	2628	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	2629	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2630	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2631	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2632	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2633	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2634	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2635	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	2636	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	2637	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	2638	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	2639	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	2640	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	2641	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	2642	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	2643	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	2644	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	2645	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	2646	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	2647	endif
	2648	endif
	2649	if (any([Lisccp_subcolumn,Lisccp_column,Lrttov_column])) then
	2650	if (any(cospgridIN%qv .lt. 0)) then
	2651	nError=nError+1
	2652	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%qv contains values out of range'
	2653	Lisccp_subcolumn = .false.
	2654	Lisccp_column = .false.
	2655	Lrttov_column = .false.
	2656	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2657	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2658	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2659	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2660	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2661	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2662	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2663	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2664	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2665	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2666	endif
	2667	endif
	2668	if (any([Lmisr_subcolumn,Lmisr_column,Lcloudsat_subcolumn,Lcloudsat_column,Lcalipso_column,Lradar_lidar_tcc,&
	2669	Llidar_only_freq_cloud,LgrLidar532_column,Latlid_column,Lcloudsat_tcc, Lcloudsat_tcc2])) then
	2670	if (any(cospgridIN%hgt_matrix .lt. -300)) then
	2671	nError=nError+1
	2672	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%hgt_matrix contains values out of range'
	2673	Lmisr_subcolumn = .false.
	2674	Lmisr_column = .false.
	2675	Lcloudsat_subcolumn = .false.
	2676	Lcloudsat_column = .false.
	2677	Lcalipso_column = .false.
	2678	Lradar_lidar_tcc = .false.
	2679	Llidar_only_freq_cloud = .false.
	2680	Lcloudsat_tcc = .false.
	2681	Lcloudsat_tcc2 = .false.
	2682	Latlid_column = .false.
	2683	LgrLidar532_column = .false.
	2684	if (associated(cospOUT%misr_fq)) cospOUT%misr_fq(:,:,:) = R_UNDEF
	2685	if (associated(cospOUT%misr_dist_model_layertops)) cospOUT%misr_dist_model_layertops(:,:) = R_UNDEF
	2686	if (associated(cospOUT%misr_meanztop)) cospOUT%misr_meanztop(:) = R_UNDEF
	2687	if (associated(cospOUT%misr_cldarea)) cospOUT%misr_cldarea(:) = R_UNDEF
	2688	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2689	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2690	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2691	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2692	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2693	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2694	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	2695	if (associated(cospOUT%cloudsat_Ze_tot)) cospOUT%cloudsat_Ze_tot(:,:,:) = R_UNDEF
	2696	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2697	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2698	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2699	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2700	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	2701	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	2702	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	2703	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	2704	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	2705	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	2706	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	2707	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	2708	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	2709	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	2710	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	2711	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	2712	endif
	2713	endif
	2714	if (any([Lrttov_column,Lcloudsat_column,Lcalipso_column,Lradar_lidar_tcc,Llidar_only_freq_cloud, &
	2715	LgrLidar532_column, Latlid_column, Lcloudsat_tcc, Lcloudsat_tcc2])) then
	2716	if (any(cospgridIN%hgt_matrix_half .lt. -300)) then
	2717	nError=nError+1
	2718	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%hgt_matrix_half contains values out of range'
	2719	Lrttov_column = .false.
	2720	Lcloudsat_column = .false.
	2721	Lcalipso_column = .false.
	2722	Lradar_lidar_tcc = .false.
	2723	Llidar_only_freq_cloud = .false.
	2724	Lcloudsat_tcc = .false.
	2725	Lcloudsat_tcc2 = .false.
	2726	Latlid_column = .false.
	2727	LgrLidar532_column = .false.
	2728	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2729	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2730	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2731	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2732	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2733	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2734	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2735	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	2736	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	2737	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	2738	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	2739	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	2740	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	2741	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	2742	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	2743	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	2744	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	2745	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	2746	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	2747	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	2748	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	2749	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	2750	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	2751	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	2752	endif
	2753	endif
	2754	if (any([Lrttov_column,Lcalipso_column,Lparasol_column])) then
	2755	if (any(cospgridIN%land .lt. 0)) then
	2756	nError=nError+1
	2757	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%land contains values out of range'
	2758	Lrttov_column = .false.
	2759	Lcalipso_column = .false.
	2760	Lparasol_column = .false.
	2761	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2762	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	2763	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	2764	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	2765	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	2766	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	2767	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	2768	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	2769	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	2770	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	2771	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	2772	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	2773	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	2774	if (associated(cospOUT%parasolGrid_refl)) cospOUT%parasolGrid_refl(:,:) = R_UNDEF
	2775	endif
	2776	endif
	2777	if (any([Lisccp_subcolumn,Lisccp_column,Lrttov_column])) then
	2778	if (any(cospgridIN%skt .lt. 0)) then
	2779	nError=nError+1
	2780	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%skt contains values out of range'
	2781	Lisccp_subcolumn = .false.
	2782	Lisccp_column = .false.
	2783	Lrttov_column = .false.
	2784	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2785	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2786	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2787	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2788	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2789	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2790	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2791	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2792	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2793	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2794	endif
	2795	endif
	2796
	2797	! RTTOV Inputs
	2798	if (Lrttov_column) then
	2799	if (cospgridIN%zenang .lt. -90. .OR. cospgridIN%zenang .gt. 90) then
	2800	nError=nError+1
	2801	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%zenang contains values out of range'
	2802	Lrttov_column = .false.
	2803	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2804	endif
	2805	if (cospgridIN%co2 .lt. 0) then
	2806	nError=nError+1
	2807	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%co2 contains values out of range'
	2808	Lrttov_column = .false.
	2809	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2810	endif
	2811	if (cospgridIN%ch4 .lt. 0) then
	2812	nError=nError+1
	2813	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%ch4 contains values out of range'
	2814	Lrttov_column = .false.
	2815	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2816	endif
	2817	if (cospgridIN%n2o .lt. 0) then
	2818	nError=nError+1
	2819	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%n2o contains values out of range'
	2820	Lrttov_column = .false.
	2821	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2822	endif
	2823	if (cospgridIN%co.lt. 0) then
	2824	nError=nError+1
	2825	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%co contains values out of range'
	2826	Lrttov_column = .false.
	2827	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2828	endif
	2829	if (any(cospgridIN%o3 .lt. 0)) then
	2830	nError=nError+1
	2831	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%o3 contains values out of range'
	2832	Lrttov_column = .false.
	2833	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2834	endif
	2835	if (any(cospgridIN%emis_sfc .lt. 0. .OR. cospgridIN%emis_sfc .gt. 1)) then
	2836	nError=nError+1
	2837	errorMessage(nError) = 'ERROR: COSP input variable: cospgridIN%emis_sfc contains values out of range'
	2838	Lrttov_column = .false.
	2839	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2840	endif
	2841	if (any(cospgridIN%u_sfc .lt. -100. .OR. cospgridIN%u_sfc .gt. 100.)) then
	2842	nError=nError+1
	2843	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%u_sfc contains values out of range'
	2844	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2845	Lrttov_column = .false.
	2846	endif
	2847	if (any(cospgridIN%v_sfc .lt. -100. .OR. cospgridIN%v_sfc .gt. 100.)) then
	2848	nError=nError+1
	2849	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%v_sfc contains values out of range'
	2850	Lrttov_column = .false.
	2851	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2852	endif
	2853	if (any(cospgridIN%lat .lt. -90 .OR. cospgridIN%lat .gt. 90)) then
	2854	nError=nError+1
	2855	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%lat contains values out of range'
	2856	Lrttov_column = .false.
	2857	if (associated(cospOUT%rttov_tbs)) cospOUT%rttov_tbs(:,:) = R_UNDEF
	2858	endif
	2859	endif
	2860
	2861	! COSP_INPUTS
	2862	if (any([Lisccp_subcolumn,Lisccp_column])) then
	2863	if (cospIN%emsfc_lw .lt. 0. .OR. cospIN%emsfc_lw .gt. 1.) then
	2864	nError=nError+1
	2865	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%emsfc_lw contains values out of range'
	2866	Lisccp_subcolumn = .false.
	2867	Lisccp_column = .false.
	2868	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2869	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2870	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2871	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2872	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2873	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2874	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2875	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2876	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2877	endif
	2878	endif
	2879	if (any([Lisccp_subcolumn,Lisccp_column,Lmisr_subcolumn,Lmisr_column,Lmodis_subcolumn,Lmodis_column])) then
	2880	if (any(cospIN%tau_067 .lt. 0)) then
	2881	nError=nError+1
	2882	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tau_067 contains values out of range'
	2883	Lisccp_subcolumn = .false.
	2884	Lisccp_column = .false.
	2885	Lmisr_subcolumn = .false.
	2886	Lmisr_column = .false.
	2887	Lmodis_subcolumn = .false.
	2888	Lmodis_column = .false.
	2889	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2890	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2891	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2892	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2893	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2894	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2895	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2896	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2897	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2898	if (associated(cospOUT%misr_fq)) cospOUT%misr_fq(:,:,:) = R_UNDEF
	2899	if (associated(cospOUT%misr_dist_model_layertops)) cospOUT%misr_dist_model_layertops(:,:) = R_UNDEF
	2900	if (associated(cospOUT%misr_meanztop)) cospOUT%misr_meanztop(:) = R_UNDEF
	2901	if (associated(cospOUT%misr_cldarea)) cospOUT%misr_cldarea(:) = R_UNDEF
	2902	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	2903	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	2904	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	2905	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	2906	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	2907	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	2908	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	2909	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	2910	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	2911	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	2912	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	2913	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	2914	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	2915	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	2916	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	2917	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	2918	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	2919	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	2920	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	2921	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	2922	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	2923	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	2924	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	2925	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	2926	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	2927	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	2928	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	2929	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	2930	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	2931	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	2932	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	2933	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	2934	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	2935	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	2936	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	2937	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	2938	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	2939	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	2940	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	2941	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	2942	endif
	2943	endif
	2944	if (any([Lisccp_subcolumn,Lisccp_column])) then
	2945	if (any(cospIN%emiss_11 .lt. 0. .OR. cospIN%emiss_11 .gt. 1)) then
	2946	nError=nError+1
	2947	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%emiss_11 contains values out of range'
	2948	Lisccp_subcolumn = .false.
	2949	Lisccp_column = .false.
	2950	if (associated(cospOUT%isccp_totalcldarea)) cospOUT%isccp_totalcldarea(:) = R_UNDEF
	2951	if (associated(cospOUT%isccp_meantb)) cospOUT%isccp_meantb(:) = R_UNDEF
	2952	if (associated(cospOUT%isccp_meantbclr)) cospOUT%isccp_meantbclr(:) = R_UNDEF
	2953	if (associated(cospOUT%isccp_meanptop)) cospOUT%isccp_meanptop(:) = R_UNDEF
	2954	if (associated(cospOUT%isccp_meantaucld)) cospOUT%isccp_meantaucld(:) = R_UNDEF
	2955	if (associated(cospOUT%isccp_meanalbedocld)) cospOUT%isccp_meanalbedocld(:) = R_UNDEF
	2956	if (associated(cospOUT%isccp_boxtau)) cospOUT%isccp_boxtau(:,:) = R_UNDEF
	2957	if (associated(cospOUT%isccp_boxptop)) cospOUT%isccp_boxptop(:,:) = R_UNDEF
	2958	if (associated(cospOUT%isccp_fq)) cospOUT%isccp_fq(:,:,:) = R_UNDEF
	2959	endif
	2960	endif
	2961	if (any([Lmodis_subcolumn,Lmodis_column])) then
	2962	if (any(cospIN%asym .lt. -1. .OR. cospIN%asym .gt. 1)) then
	2963	nError=nError+1
	2964	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%asym contains values out of range'
	2965	Lmodis_subcolumn = .false.
	2966	Lmodis_column = .false.
	2967	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	2968	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	2969	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	2970	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	2971	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	2972	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	2973	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	2974	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	2975	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	2976	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	2977	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	2978	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	2979	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	2980	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	2981	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	2982	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	2983	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	2984	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	2985	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	2986	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	2987	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	2988	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	2989	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	2990	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	2991	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	2992	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	2993	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	2994	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	2995	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	2996	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	2997	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	2998	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	2999	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	3000	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	3001	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	3002	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	3003	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	3004	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	3005	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	3006	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	3007	endif
	3008	if (any(cospIN%ss_alb .lt. 0 .OR. cospIN%ss_alb .gt. 1)) then
	3009	nError=nError+1
	3010	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%ss_alb contains values out of range'
	3011	Lmodis_subcolumn = .false.
	3012	Lmodis_column = .false.
	3013	if (associated(cospOUT%modis_Cloud_Fraction_Total_Mean)) &
	3014	cospOUT%modis_Cloud_Fraction_Total_Mean(:) = R_UNDEF
	3015	if (associated(cospOUT%modis_Cloud_Fraction_Water_Mean)) &
	3016	cospOUT%modis_Cloud_Fraction_Water_Mean(:) = R_UNDEF
	3017	if (associated(cospOUT%modis_Cloud_Fraction_Ice_Mean)) &
	3018	cospOUT%modis_Cloud_Fraction_Ice_Mean(:) = R_UNDEF
	3019	if (associated(cospOUT%modis_Cloud_Fraction_High_Mean)) &
	3020	cospOUT%modis_Cloud_Fraction_High_Mean(:) = R_UNDEF
	3021	if (associated(cospOUT%modis_Cloud_Fraction_Mid_Mean)) &
	3022	cospOUT%modis_Cloud_Fraction_Mid_Mean(:) = R_UNDEF
	3023	if (associated(cospOUT%modis_Cloud_Fraction_Low_Mean)) &
	3024	cospOUT%modis_Cloud_Fraction_Low_Mean(:) = R_UNDEF
	3025	if (associated(cospOUT%modis_Optical_Thickness_Total_Mean)) &
	3026	cospOUT%modis_Optical_Thickness_Total_Mean(:) = R_UNDEF
	3027	if (associated(cospOUT%modis_Optical_Thickness_Water_Mean)) &
	3028	cospOUT%modis_Optical_Thickness_Water_Mean(:) = R_UNDEF
	3029	if (associated(cospOUT%modis_Optical_Thickness_Ice_Mean)) &
	3030	cospOUT%modis_Optical_Thickness_Ice_Mean(:) = R_UNDEF
	3031	if (associated(cospOUT%modis_Optical_Thickness_Total_LogMean)) &
	3032	cospOUT%modis_Optical_Thickness_Total_LogMean(:) = R_UNDEF
	3033	if (associated(cospOUT%modis_Optical_Thickness_Water_LogMean)) &
	3034	cospOUT%modis_Optical_Thickness_Water_LogMean(:) = R_UNDEF
	3035	if (associated(cospOUT%modis_Optical_Thickness_Ice_LogMean)) &
	3036	cospOUT%modis_Optical_Thickness_Ice_LogMean(:) = R_UNDEF
	3037	if (associated(cospOUT%modis_Cloud_Particle_Size_Water_Mean)) &
	3038	cospOUT%modis_Cloud_Particle_Size_Water_Mean(:) = R_UNDEF
	3039	if (associated(cospOUT%modis_Cloud_Particle_Size_Ice_Mean)) &
	3040	cospOUT%modis_Cloud_Particle_Size_Ice_Mean(:) = R_UNDEF
	3041	if (associated(cospOUT%modis_Cloud_Top_Pressure_Total_Mean)) &
	3042	cospOUT%modis_Cloud_Top_Pressure_Total_Mean(:) = R_UNDEF
	3043	if (associated(cospOUT%modis_Liquid_Water_Path_Mean)) &
	3044	cospOUT%modis_Liquid_Water_Path_Mean(:) = R_UNDEF
	3045	if (associated(cospOUT%modis_Ice_Water_Path_Mean)) &
	3046	cospOUT%modis_Ice_Water_Path_Mean(:) = R_UNDEF
	3047	if (associated(cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure)) &
	3048	cospOUT%modis_Optical_Thickness_vs_Cloud_Top_Pressure(:,:,:) = R_UNDEF
	3049	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffICE)) &
	3050	cospOUT%modis_Optical_Thickness_vs_ReffICE(:,:,:) = R_UNDEF
	3051	if (associated(cospOUT%modis_Optical_Thickness_vs_ReffLIQ)) &
	3052	cospOUT%modis_Optical_Thickness_vs_ReffLIQ(:,:,:) = R_UNDEF
	3053	endif
	3054	endif
	3055	if (any([Latlid_subcolumn,Latlid_column])) then
	3056	if (any(cospIN%betatot_atlid .lt. 0)) then
	3057	nError=nError+1
	3058	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%betatot_atlid contains values out of range'
	3059	Latlid_subcolumn = .false.
	3060	Latlid_column = .false.
	3061	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	3062	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	3063	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	3064	if (associated(cospOUT%atlid_beta_tot)) cospOUT%atlid_beta_tot(:,:,:) = R_UNDEF
	3065	if (associated(cospOUT%atlid_beta_mol)) cospOUT%atlid_beta_mol(:,:) = R_UNDEF
	3066	endif
	3067	if (any(cospIN%beta_mol_atlid .lt. 0)) then
	3068	nError=nError+1
	3069	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%beta_mol_atlid contains values out of range'
	3070	Latlid_subcolumn = .false.
	3071	Latlid_column = .false.
	3072	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	3073	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	3074	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	3075	if (associated(cospOUT%atlid_beta_tot)) cospOUT%atlid_beta_tot(:,:,:) = R_UNDEF
	3076	if (associated(cospOUT%atlid_beta_mol)) cospOUT%atlid_beta_mol(:,:) = R_UNDEF
	3077	endif
	3078	if (any(cospIN%tautot_atlid .lt. 0)) then
	3079	nError=nError+1
	3080	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_atlid contains values out of range'
	3081	Latlid_subcolumn = .false.
	3082	Latlid_column = .false.
	3083	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	3084	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	3085	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	3086	if (associated(cospOUT%atlid_beta_tot)) cospOUT%atlid_beta_tot(:,:,:) = R_UNDEF
	3087	if (associated(cospOUT%atlid_beta_mol)) cospOUT%atlid_beta_mol(:,:) = R_UNDEF
	3088	endif
	3089	if (any(cospIN%tau_mol_atlid .lt. 0)) then
	3090	nError=nError+1
	3091	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tau_mol_atlid contains values out of range'
	3092	Latlid_subcolumn = .false.
	3093	Latlid_column = .false.
	3094	if (associated(cospOUT%atlid_cfad_sr)) cospOUT%atlid_cfad_sr(:,:,:) = R_UNDEF
	3095	if (associated(cospOUT%atlid_lidarcld)) cospOUT%atlid_lidarcld(:,:) = R_UNDEF
	3096	if (associated(cospOUT%atlid_cldlayer)) cospOUT%atlid_cldlayer(:,:) = R_UNDEF
	3097	if (associated(cospOUT%atlid_beta_tot)) cospOUT%atlid_beta_tot(:,:,:) = R_UNDEF
	3098	if (associated(cospOUT%atlid_beta_mol)) cospOUT%atlid_beta_mol(:,:) = R_UNDEF
	3099	endif
	3100	endif
	3101
	3102	if (any([LgrLidar532_subcolumn,LgrLidar532_column])) then
	3103	if (any(cospIN%betatot_grLidar532 .lt. 0)) then
	3104	nError=nError+1
	3105	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%betatot_grLidar532 contains values out of range'
	3106	LgrLidar532_subcolumn = .false.
	3107	LgrLidar532_column = .false.
	3108	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	3109	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	3110	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	3111	if (associated(cospOUT%grLidar532_beta_tot)) cospOUT%grLidar532_beta_tot(:,:,:) = R_UNDEF
	3112	if (associated(cospOUT%grLidar532_beta_mol)) cospOUT%grLidar532_beta_mol(:,:) = R_UNDEF
	3113	endif
	3114	if (any(cospIN%beta_mol_grLidar532 .lt. 0)) then
	3115	nError=nError+1
	3116	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%beta_mol_grLidar532 contains values out of range'
	3117	LgrLidar532_subcolumn = .false.
	3118	LgrLidar532_column = .false.
	3119	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	3120	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	3121	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	3122	if (associated(cospOUT%grLidar532_beta_tot)) cospOUT%grLidar532_beta_tot(:,:,:) = R_UNDEF
	3123	if (associated(cospOUT%grLidar532_beta_mol)) cospOUT%grLidar532_beta_mol(:,:) = R_UNDEF
	3124	endif
	3125	if (any(cospIN%tautot_grLidar532 .lt. 0)) then
	3126	nError=nError+1
	3127	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_grLidar532 contains values out of range'
	3128	LgrLidar532_subcolumn = .false.
	3129	LgrLidar532_column = .false.
	3130	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	3131	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	3132	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	3133	if (associated(cospOUT%grLidar532_beta_tot)) cospOUT%grLidar532_beta_tot(:,:,:) = R_UNDEF
	3134	if (associated(cospOUT%grLidar532_beta_mol)) cospOUT%grLidar532_beta_mol(:,:) = R_UNDEF
	3135	endif
	3136	if (any(cospIN%tau_mol_grLidar532 .lt. 0)) then
	3137	nError=nError+1
	3138	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tau_mol_grLidar532 contains values out of range'
	3139	LgrLidar532_subcolumn = .false.
	3140	LgrLidar532_column = .false.
	3141	if (associated(cospOUT%grLidar532_cfad_sr)) cospOUT%grLidar532_cfad_sr(:,:,:) = R_UNDEF
	3142	if (associated(cospOUT%grLidar532_lidarcld)) cospOUT%grLidar532_lidarcld(:,:) = R_UNDEF
	3143	if (associated(cospOUT%grLidar532_cldlayer)) cospOUT%grLidar532_cldlayer(:,:) = R_UNDEF
	3144	if (associated(cospOUT%grLidar532_beta_tot)) cospOUT%grLidar532_beta_tot(:,:,:) = R_UNDEF
	3145	if (associated(cospOUT%grLidar532_beta_mol)) cospOUT%grLidar532_beta_mol(:,:) = R_UNDEF
	3146	endif
	3147	endif
	3148
	3149	if (any([Lcalipso_subcolumn,Lcalipso_column])) then
	3150	if (any(cospIN%betatot_calipso .lt. 0)) then
	3151	nError=nError+1
	3152	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%betatot_calipso contains values out of range'
	3153	Lcalipso_subcolumn = .false.
	3154	Lcalipso_column = .false.
	3155	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3156	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3157	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3158	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3159	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3160	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3161	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3162	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3163	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3164	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3165	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3166	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3167	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3168	endif
	3169	if (any(cospIN%betatot_liq_calipso .lt. 0)) then
	3170	nError=nError+1
	3171	errorMessage(nError) = ('ERROR: COSP input variable: cospIN%betatot_liq_calipso contains values out of range')
	3172	Lcalipso_subcolumn = .false.
	3173	Lcalipso_column = .false.
	3174	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3175	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3176	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3177	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3178	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3179	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3180	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3181	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3182	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3183	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3184	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3185	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3186	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3187	endif
	3188	if (any(cospIN%betatot_ice_calipso .lt. 0)) then
	3189	nError=nError+1
	3190	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%betatot_ice_calipso contains values out of range'
	3191	Lcalipso_subcolumn = .false.
	3192	Lcalipso_column = .false.
	3193	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3194	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3195	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3196	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3197	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3198	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3199	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3200	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3201	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3202	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3203	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3204	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3205	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3206	endif
	3207	if (any(cospIN%tautot_calipso .lt. 0)) then
	3208	nError=nError+1
	3209	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_calipso contains values out of range'
	3210	Lcalipso_subcolumn = .false.
	3211	Lcalipso_column = .false.
	3212	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3213	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3214	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3215	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3216	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3217	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3218	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3219	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3220	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3221	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3222	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3223	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3224	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3225	endif
	3226	if (any(cospIN%tautot_liq_calipso .lt. 0)) then
	3227	nError=nError+1
	3228	errorMessage(nError) = ('ERROR: COSP input variable: cospIN%tautot_liq_calipso contains values out of range')
	3229	Lcalipso_subcolumn = .false.
	3230	Lcalipso_column = .false.
	3231	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3232	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3233	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3234	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3235	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3236	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3237	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3238	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3239	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3240	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3241	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3242	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3243	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3244	endif
	3245	if (any(cospIN%tautot_ice_calipso .lt. 0)) then
	3246	nError=nError+1
	3247	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_ice_calipso contains values out of range'
	3248	Lcalipso_subcolumn = .false.
	3249	Lcalipso_column = .false.
	3250	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3251	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3252	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3253	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3254	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3255	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3256	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3257	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3258	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3259	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3260	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3261	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3262	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3263	endif
	3264	if (any(cospIN%tau_mol_calipso .lt. 0)) then
	3265	nError=nError+1
	3266	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tau_mol_calipso contains values out of range'
	3267	Lcalipso_subcolumn = .false.
	3268	Lcalipso_column = .false.
	3269	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3270	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3271	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3272	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3273	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3274	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3275	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3276	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3277	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3278	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3279	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3280	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3281	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3282	endif
	3283	endif
	3284	if (any([Lcalipso_subcolumn,Lcalipso_column,Lcloudsat_column,Lradar_lidar_tcc, &
	3285	Llidar_only_freq_cloud, Lcloudsat_tcc, Lcloudsat_tcc2])) then
	3286	if (any(cospIN%beta_mol_calipso .lt. 0)) then
	3287	nError=nError+1
	3288	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%beta_mol_calipso contains values out of range'
	3289	Lcalipso_subcolumn = .false.
	3290	Lcalipso_column = .false.
	3291	Lcloudsat_column = .false.
	3292	Lradar_lidar_tcc = .false.
	3293	Llidar_only_freq_cloud = .false.
	3294	Lcloudsat_tcc = .false.
	3295	Lcloudsat_tcc2 = .false.
	3296	if (associated(cospOUT%calipso_cfad_sr)) cospOUT%calipso_cfad_sr(:,:,:) = R_UNDEF
	3297	if (associated(cospOUT%calipso_lidarcld)) cospOUT%calipso_lidarcld(:,:) = R_UNDEF
	3298	if (associated(cospOUT%calipso_lidarcldphase)) cospOUT%calipso_lidarcldphase(:,:,:) = R_UNDEF
	3299	if (associated(cospOUT%calipso_cldlayer)) cospOUT%calipso_cldlayer(:,:) = R_UNDEF
	3300	if (associated(cospOUT%calipso_cldlayerphase)) cospOUT%calipso_cldlayerphase(:,:,:) = R_UNDEF
	3301	if (associated(cospOUT%calipso_lidarcldtmp)) cospOUT%calipso_lidarcldtmp(:,:,:) = R_UNDEF
	3302	if (associated(cospOUT%calipso_srbval)) cospOUT%calipso_srbval(:) = R_UNDEF
	3303	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	3304	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	3305	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	3306	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	3307	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	3308	if (associated(cospOUT%calipso_lidarcldtype)) cospOUT%calipso_lidarcldtype(:,:,:) = R_UNDEF
	3309	if (associated(cospOUT%calipso_cldtype)) cospOUT%calipso_cldtype(:,:) = R_UNDEF
	3310	if (associated(cospOUT%calipso_cldtypetemp)) cospOUT%calipso_cldtypetemp(:,:) = R_UNDEF
	3311	if (associated(cospOUT%calipso_cldtypemeanz)) cospOUT%calipso_cldtypemeanz(:,:) = R_UNDEF
	3312	if (associated(cospOUT%calipso_cldtypemeanzse)) cospOUT%calipso_cldtypemeanzse(:,:) = R_UNDEF
	3313	if (associated(cospOUT%calipso_cldthinemis)) cospOUT%calipso_cldthinemis(:) = R_UNDEF
	3314	endif
	3315	endif
	3316	if (any([Lparasol_subcolumn,Lparasol_column])) then
	3317	if (any(cospIN%tautot_S_liq .lt. 0)) then
	3318	nError=nError+1
	3319	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_S_liq contains values out of range'
	3320	Lparasol_subcolumn = .false.
	3321	Lparasol_column = .false.
	3322	if (associated(cospOUT%parasolPix_refl)) cospOUT%parasolPix_refl(:,:,:) = R_UNDEF
	3323	if (associated(cospOUT%parasolGrid_refl)) cospOUT%parasolGrid_refl(:,:) = R_UNDEF
	3324	endif
	3325	if (any(cospIN%tautot_S_ice .lt. 0)) then
	3326	nError=nError+1
	3327	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%tautot_S_ice contains values out of range'
	3328	Lparasol_subcolumn = .false.
	3329	Lparasol_column = .false.
	3330	if (associated(cospOUT%parasolPix_refl)) cospOUT%parasolPix_refl(:,:,:) = R_UNDEF
	3331	if (associated(cospOUT%parasolGrid_refl)) cospOUT%parasolGrid_refl(:,:) = R_UNDEF
	3332	endif
	3333	endif
	3334	if (any([Lcloudsat_subcolumn,Lcloudsat_column,Lradar_lidar_tcc,Llidar_only_freq_cloud, &
	3335	Lcloudsat_tcc, Lcloudsat_tcc2])) then
	3336	if (any(cospIN%z_vol_cloudsat .lt. 0)) then
	3337	nError=nError+1
	3338	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%z_vol_cloudsat contains values out of range'
	3339	Lcloudsat_subcolumn = .false.
	3340	Lcloudsat_column = .false.
	3341	Lradar_lidar_tcc = .false.
	3342	Llidar_only_freq_cloud = .false.
	3343	Lcloudsat_tcc = .false.
	3344	Lcloudsat_tcc2 = .false.
	3345	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	3346	if (associated(cospOUT%cloudsat_Ze_tot)) cospOUT%cloudsat_Ze_tot(:,:,:) = R_UNDEF
	3347	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	3348	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	3349	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	3350	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	3351	endif
	3352	if (any(cospIN%kr_vol_cloudsat .lt. 0)) then
	3353	nError=nError+1
	3354	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%kr_vol_cloudsat contains values out of range'
	3355	Lcloudsat_subcolumn = .false.
	3356	Lcloudsat_column = .false.
	3357	Lradar_lidar_tcc = .false.
	3358	Llidar_only_freq_cloud = .false.
	3359	Lcloudsat_tcc = .false.
	3360	Lcloudsat_tcc2 = .false.
	3361	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	3362	if (associated(cospOUT%cloudsat_Ze_tot)) cospOUT%cloudsat_Ze_tot(:,:,:) = R_UNDEF
	3363	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	3364	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	3365	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	3366	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	3367	endif
	3368	if (any(cospIN%g_vol_cloudsat .lt. 0)) then
	3369	nError=nError+1
	3370	errorMessage(nError) = 'ERROR: COSP input variable: cospIN%g_vol_cloudsat contains values out of range'
	3371	Lcloudsat_subcolumn = .false.
	3372	Lcloudsat_column = .false.
	3373	Lradar_lidar_tcc = .false.
	3374	Llidar_only_freq_cloud = .false.
	3375	Lcloudsat_tcc = .false.
	3376	Lcloudsat_tcc2 = .false.
	3377	if (associated(cospOUT%cloudsat_cfad_ze)) cospOUT%cloudsat_cfad_ze(:,:,:) = R_UNDEF
	3378	if (associated(cospOUT%cloudsat_Ze_tot)) cospOUT%cloudsat_Ze_tot(:,:,:) = R_UNDEF
	3379	if (associated(cospOUT%lidar_only_freq_cloud)) cospOUT%lidar_only_freq_cloud(:,:) = R_UNDEF
	3380	if (associated(cospOUT%radar_lidar_tcc)) cospOUT%radar_lidar_tcc(:) = R_UNDEF
	3381	if (associated(cospOUT%cloudsat_tcc)) cospOUT%cloudsat_tcc(:) = R_UNDEF
	3382	if (associated(cospOUT%cloudsat_tcc2)) cospOUT%cloudsat_tcc2(:) = R_UNDEF
	3383	endif
	3384	endif
	3385	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	3386	! Part 2: Check input fields array size for consistency. This needs to be done for each
	3387	! simulator
	3388	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	3389	! ISCCP
	3390	if (Lisccp_subcolumn .or. Lisccp_column) then
	3391	if (size(cospIN%frac_out,1) .ne. cospIN%Npoints .OR. &
	3392	size(cospIN%tau_067,1) .ne. cospIN%Npoints .OR. &
	3393	size(cospIN%emiss_11,1) .ne. cospIN%Npoints .OR. &
	3394	size(cospgridIN%skt) .ne. cospIN%Npoints .OR. &
	3395	size(cospgridIN%qv,1) .ne. cospIN%Npoints .OR. &
	3396	size(cospgridIN%at,1) .ne. cospIN%Npoints .OR. &
	3397	size(cospgridIN%phalf,1) .ne. cospIN%Npoints .OR. &
	3398	size(cospgridIN%sunlit) .ne. cospIN%Npoints .OR. &
	3399	size(cospgridIN%pfull,1) .ne. cospIN%Npoints) then
	3400	Lisccp_subcolumn = .false.
	3401	Lisccp_column = .false.
	3402	nError=nError+1
	3403	errorMessage(nError) = 'ERROR(isccp_simulator): The number of points in the input fields are inconsistent'
	3404	endif
	3405	if (size(cospIN%frac_out,2) .ne. cospIN%Ncolumns .OR. &
	3406	size(cospIN%tau_067,2) .ne. cospIN%Ncolumns .OR. &
	3407	size(cospIN%emiss_11,2) .ne. cospIN%Ncolumns) then
	3408	Lisccp_subcolumn = .false.
	3409	Lisccp_column = .false.
	3410	nError=nError+1
	3411	errorMessage(nError) = 'ERROR(isccp_simulator): The number of sub-columns in the input fields are inconsistent'
	3412	endif
	3413	if (size(cospIN%frac_out,3) .ne. cospIN%Nlevels .OR. &
	3414	size(cospIN%tau_067,3) .ne. cospIN%Nlevels .OR. &
	3415	size(cospIN%emiss_11,3) .ne. cospIN%Nlevels .OR. &
	3416	size(cospgridIN%qv,2) .ne. cospIN%Nlevels .OR. &
	3417	size(cospgridIN%at,2) .ne. cospIN%Nlevels .OR. &
	3418	size(cospgridIN%pfull,2) .ne. cospIN%Nlevels .OR. &
	3419	size(cospgridIN%phalf,2) .ne. cospIN%Nlevels+1) then
	3420	Lisccp_subcolumn = .false.
	3421	Lisccp_column = .false.
	3422	nError=nError+1
	3423	errorMessage(nError) = 'ERROR(isccp_simulator): The number of levels in the input fields are inconsistent'
	3424	endif
	3425	endif
	3426
	3427	! MISR
	3428	if (Lmisr_subcolumn .or. Lmisr_column) then
	3429	if (size(cospIN%tau_067,1) .ne. cospIN%Npoints .OR. &
	3430	size(cospgridIN%sunlit) .ne. cospIN%Npoints .OR. &
	3431	size(cospgridIN%hgt_matrix,1) .ne. cospIN%Npoints .OR. &
	3432	size(cospgridIN%at,1) .ne. cospIN%Npoints) then
	3433	Lmisr_subcolumn = .false.
	3434	Lmisr_column = .false.
	3435	nError=nError+1
	3436	errorMessage(nError) = 'ERROR(misr_simulator): The number of points in the input fields are inconsistent'
	3437	endif
	3438	if (size(cospIN%tau_067,2) .ne. cospIN%Ncolumns) then
	3439	Lmisr_subcolumn = .false.
	3440	Lmisr_column = .false.
	3441	nError=nError+1
	3442	errorMessage(nError) = 'ERROR(misr_simulator): The number of sub-columns in the input fields are inconsistent'
	3443	endif
	3444	if (size(cospIN%tau_067,3) .ne. cospIN%Nlevels .OR. &
	3445	size(cospgridIN%hgt_matrix,2) .ne. cospIN%Nlevels .OR. &
	3446	size(cospgridIN%at,2) .ne. cospIN%Nlevels) then
	3447	Lmisr_subcolumn = .false.
	3448	Lmisr_column = .false.
	3449	nError=nError+1
	3450	errorMessage(nError) = 'ERROR(misr_simulator): The number of levels in the input fields are inconsistent'
	3451	endif
	3452	endif
	3453
	3454	! MODIS
	3455	if (Lmodis_subcolumn .or. Lmodis_column) then
	3456	if (size(cospIN%fracLiq,1) .ne. cospIN%Npoints .OR. &
	3457	size(cospIN%tau_067,1) .ne. cospIN%Npoints .OR. &
	3458	size(cospIN%asym,1) .ne. cospIN%Npoints .OR. &
	3459	size(cospIN%ss_alb,1) .ne. cospIN%Npoints) then
	3460	Lmodis_subcolumn = .false.
	3461	Lmodis_column = .false.
	3462	nError=nError+1
	3463	errorMessage(nError) = 'ERROR(modis_simulator): The number of points in the input fields are inconsistent'
	3464	endif
	3465	if (size(cospIN%fracLiq,2) .ne. cospIN%Ncolumns .OR. &
	3466	size(cospIN%tau_067,2) .ne. cospIN%Ncolumns .OR. &
	3467	size(cospIN%asym,2) .ne. cospIN%Ncolumns .OR. &
	3468	size(cospIN%ss_alb,2) .ne. cospIN%Ncolumns) then
	3469	Lmodis_subcolumn = .false.
	3470	Lmodis_column = .false.
	3471	nError=nError+1
	3472	errorMessage(nError) = 'ERROR(modis_simulator): The number of sub-columns in the input fields are inconsistent'
	3473	endif
	3474	if (size(cospIN%fracLiq,3) .ne. cospIN%Nlevels .OR. &
	3475	size(cospIN%tau_067,3) .ne. cospIN%Nlevels .OR. &
	3476	size(cospIN%asym,3) .ne. cospIN%Nlevels .OR. &
	3477	size(cospIN%ss_alb,3) .ne. cospIN%Nlevels) then
	3478	Lmodis_subcolumn = .false.
	3479	Lmodis_column = .false.
	3480	nError=nError+1
	3481	errorMessage(nError) = 'ERROR(modis_simulator): The number of levels in the input fields are inconsistent'
	3482	endif
	3483	endif
	3484
	3485	! CLOUDSAT
	3486	if (Lcloudsat_subcolumn .or. Lcloudsat_column) then
	3487	if (size(cospIN%z_vol_cloudsat,1) .ne. cospIN%Npoints .OR. &
	3488	size(cospIN%kr_vol_cloudsat,1) .ne. cospIN%Npoints .OR. &
	3489	size(cospIN%g_vol_cloudsat,1) .ne. cospIN%Npoints .OR. &
	3490	size(cospgridIN%hgt_matrix,1) .ne. cospIN%Npoints) then
	3491	Lcloudsat_subcolumn = .false.
	3492	Lcloudsat_column = .false.
	3493	nError=nError+1
	3494	errorMessage(nError) = 'ERROR(cloudsat_simulator): The number of points in the input fields are inconsistent'
	3495	endif
	3496	if (size(cospIN%z_vol_cloudsat,2) .ne. cospIN%Ncolumns .OR. &
	3497	size(cospIN%kr_vol_cloudsat,2) .ne. cospIN%Ncolumns .OR. &
	3498	size(cospIN%g_vol_cloudsat,2) .ne. cospIN%Ncolumns) then
	3499	Lcloudsat_subcolumn = .false.
	3500	Lcloudsat_column = .false.
	3501	nError=nError+1
	3502	errorMessage(nError) = 'ERROR(cloudsat_simulator): The number of sub-columns in the input fields are inconsistent'
	3503	endif
	3504	if (size(cospIN%z_vol_cloudsat,3) .ne. cospIN%Nlevels .OR. &
	3505	size(cospIN%kr_vol_cloudsat,3) .ne. cospIN%Nlevels .OR. &
	3506	size(cospIN%g_vol_cloudsat,3) .ne. cospIN%Nlevels .OR. &
	3507	size(cospgridIN%hgt_matrix,2) .ne. cospIN%Nlevels) then
	3508	Lcloudsat_subcolumn = .false.
	3509	Lcloudsat_column = .false.
	3510	nError=nError+1
	3511	errorMessage(nError) = 'ERROR(cloudsat_simulator): The number of levels in the input fields are inconsistent'
	3512	endif
	3513	endif
	3514
	3515	! GROUND LIDAR @ 532nm
	3516	if (LgrLidar532_subcolumn .or. LgrLidar532_column) then
	3517	if (size(cospIN%beta_mol_grLidar532,1) .ne. cospIN%Npoints .OR. &
	3518	size(cospIN%betatot_grLidar532,1) .ne. cospIN%Npoints .OR. &
	3519	size(cospIN%tau_mol_grLidar532,1) .ne. cospIN%Npoints .OR. &
	3520	size(cospIN%tautot_grLidar532,1) .ne. cospIN%Npoints) then
	3521	LgrLidar532_subcolumn = .false.
	3522	LgrLidar532_column = .false.
	3523	nError=nError+1
	3524	errorMessage(nError) = 'ERROR(grLidar532_simulator): The number of points in the input fields are inconsistent'
	3525	endif
	3526	if (size(cospIN%betatot_grLidar532,2) .ne. cospIN%Ncolumns .OR. &
	3527	size(cospIN%tautot_grLidar532,2) .ne. cospIN%Ncolumns) then
	3528	LgrLidar532_subcolumn = .false.
	3529	LgrLidar532_column = .false.
	3530	nError=nError+1
	3531	errorMessage(nError) = 'ERROR(grLidar532_simulator): The number of sub-columns in the input fields are inconsistent'
	3532	endif
	3533	if (size(cospIN%beta_mol_grLidar532,2) .ne. cospIN%Nlevels .OR. &
	3534	size(cospIN%betatot_grLidar532,3) .ne. cospIN%Nlevels .OR. &
	3535	size(cospIN%tau_mol_grLidar532,2) .ne. cospIN%Nlevels .OR. &
	3536	size(cospIN%tautot_grLidar532,3) .ne. cospIN%Nlevels) then
	3537	LgrLidar532_subcolumn = .false.
	3538	LgrLidar532_column = .false.
	3539	nError=nError+1
	3540	errorMessage(nError) = 'ERROR(grLidar532_simulator): The number of levels in the input fields are inconsistent'
	3541	endif
	3542	endif
	3543
	3544	! ATLID
	3545	if (Latlid_subcolumn .or. Latlid_column) then
	3546	if (size(cospIN%beta_mol_atlid,1) .ne. cospIN%Npoints .OR. &
	3547	size(cospIN%betatot_atlid,1) .ne. cospIN%Npoints .OR. &
	3548	size(cospIN%tau_mol_atlid,1) .ne. cospIN%Npoints .OR. &
	3549	size(cospIN%tautot_atlid,1) .ne. cospIN%Npoints) then
	3550	Latlid_subcolumn = .false.
	3551	Latlid_column = .false.
	3552	nError=nError+1
	3553	errorMessage(nError) = 'ERROR(atlid_simulator): The number of points in the input fields are inconsistent'
	3554	endif
	3555	if (size(cospIN%betatot_atlid,2) .ne. cospIN%Ncolumns .OR. &
	3556	size(cospIN%tautot_atlid,2) .ne. cospIN%Ncolumns) then
	3557	Latlid_subcolumn = .false.
	3558	Latlid_column = .false.
	3559	nError=nError+1
	3560	errorMessage(nError) = 'ERROR(atlid_simulator): The number of sub-columns in the input fields are inconsistent'
	3561	endif
	3562	if (size(cospIN%beta_mol_atlid,2) .ne. cospIN%Nlevels .OR. &
	3563	size(cospIN%betatot_atlid,3) .ne. cospIN%Nlevels .OR. &
	3564	size(cospIN%tau_mol_atlid,2) .ne. cospIN%Nlevels .OR. &
	3565	size(cospIN%tautot_atlid,3) .ne. cospIN%Nlevels) then
	3566	Latlid_subcolumn = .false.
	3567	Latlid_column = .false.
	3568	nError=nError+1
	3569	errorMessage(nError) = 'ERROR(atlid_simulator): The number of levels in the input fields are inconsistent'
	3570	endif
	3571	endif
	3572
	3573	! CALIPSO
	3574	if (Lcalipso_subcolumn .or. Lcalipso_column) then
	3575	if (size(cospIN%beta_mol_calipso,1) .ne. cospIN%Npoints .OR. &
	3576	size(cospIN%betatot_calipso,1) .ne. cospIN%Npoints .OR. &
	3577	size(cospIN%betatot_liq_calipso,1) .ne. cospIN%Npoints .OR. &
	3578	size(cospIN%betatot_ice_calipso,1) .ne. cospIN%Npoints .OR. &
	3579	size(cospIN%tau_mol_calipso,1) .ne. cospIN%Npoints .OR. &
	3580	size(cospIN%tautot_calipso,1) .ne. cospIN%Npoints .OR. &
	3581	size(cospIN%tautot_liq_calipso,1) .ne. cospIN%Npoints .OR. &
	3582	size(cospIN%tautot_ice_calipso,1) .ne. cospIN%Npoints) then
	3583	Lcalipso_subcolumn = .false.
	3584	Lcalipso_column = .false.
	3585	nError=nError+1
	3586	errorMessage(nError) = 'ERROR(calipso_simulator): The number of points in the input fields are inconsistent'
	3587	endif
	3588	if (size(cospIN%betatot_calipso,2) .ne. cospIN%Ncolumns .OR. &
	3589	size(cospIN%betatot_liq_calipso,2) .ne. cospIN%Ncolumns .OR. &
	3590	size(cospIN%betatot_ice_calipso,2) .ne. cospIN%Ncolumns .OR. &
	3591	size(cospIN%tautot_calipso,2) .ne. cospIN%Ncolumns .OR. &
	3592	size(cospIN%tautot_liq_calipso,2) .ne. cospIN%Ncolumns .OR. &
	3593	size(cospIN%tautot_ice_calipso,2) .ne. cospIN%Ncolumns) then
	3594	Lcalipso_subcolumn = .false.
	3595	Lcalipso_column = .false.
	3596	nError=nError+1
	3597	errorMessage(nError) = 'ERROR(calipso_simulator): The number of sub-columns in the input fields are inconsistent'
	3598	endif
	3599	if (size(cospIN%beta_mol_calipso,2) .ne. cospIN%Nlevels .OR. &
	3600	size(cospIN%betatot_calipso,3) .ne. cospIN%Nlevels .OR. &
	3601	size(cospIN%betatot_liq_calipso,3) .ne. cospIN%Nlevels .OR. &
	3602	size(cospIN%betatot_ice_calipso,3) .ne. cospIN%Nlevels .OR. &
	3603	size(cospIN%tau_mol_calipso,2) .ne. cospIN%Nlevels .OR. &
	3604	size(cospIN%tautot_calipso,3) .ne. cospIN%Nlevels .OR. &
	3605	size(cospIN%tautot_liq_calipso,3) .ne. cospIN%Nlevels .OR. &
	3606	size(cospIN%tautot_ice_calipso,3) .ne. cospIN%Nlevels) then
	3607	Lcalipso_subcolumn = .false.
	3608	Lcalipso_column = .false.
	3609	nError=nError+1
	3610	errorMessage(nError) = 'ERROR(calipso_simulator): The number of levels in the input fields are inconsistent'
	3611	endif
	3612	endif
	3613
	3614	! PARASOL
	3615	if (Lparasol_subcolumn .or. Lparasol_column) then
	3616	if (size(cospIN%tautot_S_liq,1) .ne. cospIN%Npoints .OR. &
	3617	size(cospIN%tautot_S_ice,1) .ne. cospIN%Npoints) then
	3618	Lparasol_subcolumn = .false.
	3619	Lparasol_column = .false.
	3620	nError=nError+1
	3621	errorMessage(nError) = 'ERROR(parasol_simulator): The number of points in the input fields are inconsistent'
	3622	endif
	3623	if (size(cospIN%tautot_S_liq,2) .ne. cospIN%Ncolumns .OR. &
	3624	size(cospIN%tautot_S_ice,2) .ne. cospIN%Ncolumns) then
	3625	Lparasol_subcolumn = .false.
	3626	Lparasol_column = .false.
	3627	nError=nError+1
	3628	errorMessage(nError) = 'ERROR(parasol_simulator): The number of levels in the input fields are inconsistent'
	3629	endif
	3630	endif
	3631
	3632	! RTTOV
	3633	if (Lrttov_column) then
	3634	if (size(cospgridIN%pfull,1) .ne. cospIN%Npoints .OR. &
	3635	size(cospgridIN%at,1) .ne. cospIN%Npoints .OR. &
	3636	size(cospgridIN%qv,1) .ne. cospIN%Npoints .OR. &
	3637	size(cospgridIN%hgt_matrix_half,1) .ne. cospIN%Npoints .OR. &
	3638	size(cospgridIN%u_sfc) .ne. cospIN%Npoints .OR. &
	3639	size(cospgridIN%v_sfc) .ne. cospIN%Npoints .OR. &
	3640	size(cospgridIN%skt) .ne. cospIN%Npoints .OR. &
	3641	size(cospgridIN%phalf,1) .ne. cospIN%Npoints .OR. &
	3642	size(cospgridIN%qv,1) .ne. cospIN%Npoints .OR. &
	3643	size(cospgridIN%land) .ne. cospIN%Npoints .OR. &
	3644	size(cospgridIN%lat) .ne. cospIN%Npoints) then
	3645	Lrttov_column = .false.
	3646	nError=nError+1
	3647	errorMessage(nError) = 'ERROR(rttov_simulator): The number of points in the input fields are inconsistent'
	3648	endif
	3649	if (size(cospgridIN%pfull,2) .ne. cospIN%Nlevels .OR. &
	3650	size(cospgridIN%at,2) .ne. cospIN%Nlevels .OR. &
	3651	size(cospgridIN%qv,2) .ne. cospIN%Nlevels .OR. &
	3652	size(cospgridIN%hgt_matrix_half,2) .ne. cospIN%Nlevels+1 .OR. &
	3653	size(cospgridIN%phalf,2) .ne. cospIN%Nlevels+1 .OR. &
	3654	size(cospgridIN%qv,2) .ne. cospIN%Nlevels) then
	3655	Lrttov_column = .false.
	3656	nError=nError+1
	3657	errorMessage(nError) = 'ERROR(rttov_simulator): The number of levels in the input fields are inconsistent'
	3658	endif
	3659	endif
	3660	end subroutine cosp_errorCheck
	3661
	3662	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	3663	! END MODULE
	3664	!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	3665	END MODULE MOD_COSP

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

Download in other formats:

Original Format