COSP v1.2.2 User`s manual - CFMIP
Contents
1. hydrometeor density kgm If used then set HCLASS_APM and HCLASS_BPM to 1 e HCLASS P1 HCLASS P2 HCLASS P3 these parameters depend on the type of dis tribution For the modified gamma distribution used in the UM P1 is the total particle number concentration P2 is the particle mean diameter um and P3 is the distribution width a 1 One of the parameters P1 P2 must be specified and the other one should be set to 1 P3 must be specified The settings for DMIN and DMAX are ignored in the current version for all distributions except for power law Except when the power law distribution is used particle size is fixed to vary from zero to infinity Since COSP v0 2 a capability of Quickbeam to pass the effective radius as input parameter is used In that case the settings in HCLASS_P 1 3 are defaults If the input R is zero at any spatial or hydrometeor coordinate at which there is condensate then the HCALSS default is used Hence if the effective radius is not zero when there is hydrometeor present the values in HCLASS P2 are not used The default values in the COSP HCLASS table reflect those used by Roj Marchand to run the simulator for the MMF Marchand et al 2009 4 Configuration for CFMIP 2 experiments The directory cfmip2 contains the namelists with the configuration for the CFMIP 2 experi ments These files are also available on the CFMIP web site There are two different configu rations 11 e Long time2. Climate and Forecast CF Metadata Convention and fulfill the require ments of the climate community s standard model experiments The Coupled Model Intercom parison Project Phase 5 CMIP5 has requested COSP outputs to be included into a subset of CMIP5 experiments COSP is open source software and can be downloaded from the CFMIP website without charge The document is organised as follows Section 2 provides information on the namelists that are used to configure COSP Section 3 discusses how to set up the microphysical settings Section 4 gives some details on the configuration of COSP for CFMIP 2 experiments Appendix A shows the structure of the NetCDF input data files This document is still under development and therefore is not complete although hope it will still be useful in its current form It is encouraged to read the README txt file that is included with COSP along with this user s manual 2 Configuration setting the COSP namelists The user interaction with COSP is done via namelists This section provides information on the namelists that are used to configure COSP 2 1 COSP INPUT namelist This namelist is located in file cosp_input_nl txt and it contains the input arguments for COSP and all the simulators Table 1 contains a description of the variables in this namelist For details on RTTOV variables please refer to RTTOV documentation Thttp cmip pemdi linl gov emip5 experiment_design html http Awww cfmi
3. Optical depth of stratiform cloud at 0 67 micron dtau_s _FillValue 1 e 30f dtau_s units 1 float dtau_c level point dtau_c long_name Optical depth of convective cloud at 0 67 micron dtau_c _FillValue 1 e 30f dtau_c units 1 float dem_s level point dem_s long_name Longwave emissivity of stratiform cloud at 10 5 micron dem_s _FillValue 1 e 30f dem_s units 1 float dem_c level point dem_c long_name Longwave emissivity of convective cloud at 10 5 micron dem_c _FillValue 1 e 30f dem_c units 1 float skt point skt long_name Skin temperature skt _FillValue 1 e 30f skt units K float sunlit point sunlit long_name Day points sunlit _FillValue 1 e 30f sunlit units 1 float u_wind point u_wind long_name eastward_wind u_wind _FillValue 1 e 30f u_wind units m s 1 float v_wind point v_wind long_name northward_wind v_wind _FillValue 1 e 30f v_wind units m s 1 float mr_ozone level point mr_ozone long_name mass_fraction_of_ozone_in_air mr_ozone _FillValue 1 e 30f mr_ozone units kg kg float emsfc_lw emsfc_lw long_name Surface emissivity at 10 5 micron emsfc_lw _FillValue 1 e 30f emsfc_lw units 1 global attributes 16 fraction title COSP inputs UKMO N320L50 Conventions CF 1 0 description This is the CDL structure of the COSP
4. series long inline txt This is the configuration for the 30 yr monthly and daily means from ISCCP and CALIPSO PARASOL These are global gridded data computed from model gridded inputs with the simulators run inline The production version for these experiments is COSP v1 2 2 e Short time series short_offline txt This is the configuration for the 1 yr time series both for the curtain outputs and global gridded monthly means from curtain outputs Outputs from CloudSat and CALIPSO PARASOL are requested It is hoped that v1 3 will be the production version for these experiments It will contain the final version of the MIP tables released by PCMDI 5 Using your own cloud generator Currently COSP only includes treatment for cloud precipitation overlap but not subgrid vari ability Please see Section 6 5 of the README txt file if you require this extra capability Acknowledgements COSP is a collaborative effort and many people have been involved in the development of the software Thanks to M J Webb S Bony H Chepfer J L Dufresne S A Klein Y Zhang R Marchand J M Haynes R Pincus and V O John Appendix A Structure of the NetCDF input data files The structure of the input data NetCDF files are listed below Examples of these files are distributed with COSP namely cosp_input_um nc for 1D mode and cosp_input_um_2d nc for 2D mode The 1D mode represents data along a trajectory like the orbit track The 2D mo
5. 1 Ry 9 Qu2 Hlas 3 9 3 2 Mixing ratios from precipitation fluxes The radar reflectivities are computed from the hydormeteor mixing ratios However as large scale models typically diagnose precipitation fluxes there exists the possibility of passing pre cipitation fluxes and let COSP convert them into mixing ratios before calling the radar simula tor The variable use_precipitation_fluxes in the COSP_INPUT namelist controls whether the COSP should do this conversion or use the input mixing ratios instead Expanding and integrating Eq 3 1 the expression for the precipitation flux as a function of the mixing ratio and the parameters that define the PSD is given by lz T 1 ax a T FE paz z ax be de PQx b ba 1 i 10 az bz 1 Nazar Az by 1 Solving for the mixing ratio gives 1 p Ga 1 Po where nu Ty ax br dx 1 T2 T az be 1 gE F MarGel 2 Currently the subroutine 11n1 pf_to_mr is used to convert precipitation fluxes to mixing ratios It uses the method from Khairoutdinov and Randall 2003 which is a little bit more restrictive than the method presented here because it assumes constant densities for all hy drometeors and the PSD is a Marshall Palmer distribution a special case of the one presented here Table 4 gives the relationship between the formulation used here and the method used in Khairoutdinov and Randall 2003 Using this correponde
6. COSP user s manual Version 1 2 2 A Bodas Salcedo Met Office Hadley Centre FitzRoy Rd Exeter EX1 3PB United Kingdom March 30 2010 British Crown Copyright 2010 Contents 1 Introduction 1 2 Configuration setting the COSP namelists 2 2 1 GOSP INPUT MAMGISt gt w 0 2 2 02 osent vais ete Ak ELE wl Dog oe a 2 227 CMOR MAMEIST ics So ore tte Vw at Rad Oe Wea eal ae ruse 5 23 COSP OUTPUT namelist vhs este ge do eye oe Bad ahaa ag ae Le 6 3 Microphysical settings 7 9 1 Effective radius eu Aan BM Ark APS a a eA Ae es Ah mn ER 8 3 2 Mixing ratios from precipitation fluxes 9 3 3 Setting the HCLASS table 10 4 Configuration for CFMIP 2 experiments 11 5 Using your own cloud generator 12 1 Introduction The Cloud Feedback Model Intercomparison Project CFMIP Observation Simulator Package COSP is a modular piece of software whose main aim is to enable the simulation of data from several satellite borne sensors from model variables It is written almost entirely in Fortran 90 and it is conceptually divided into three steps First the gridbox mean profiles are broken into subcolumns Then the vertical profiles of individual subcolumns are passed to individual instrument simulators e g lidar forward model ISCCP similator Finally a statistical module gathers the outputs from all the instruments and builds statistics that can be compared to similar statistics
7. Mon Weather Rev 127 10 2514 2531 1999 Marchand R J Haynes G G Mace T Ackerman and G Stephens A comparison of simulated cloud radar output from the multiscale modeling framework global climate model with CloudSat cloud radar observations J Geophys Res 114 DOOA20 doi 10 1029 2008JD009790 2009 Pincus R C Hannay S A Klein K M Xu and R Hemler Overlap assumptions for assumed probability distribution function cloud schemes in large scale models J Geophys Res 110 D15S09 doi doi 10 1029 2004JD005100 2005 Saunders R M Matricardi and P Brunel An improved fast radiative transfer model for assim ilation of satellite radiance observations Q J R Meteorol Soc 125 1407 1425 1999 Webb M C Senior S Bony and J J Morcrette Combining ERBE and ISCCP data to assess clouds in the Hadley Centre ECMWF and LMD atmospheric climate models Clim Dyn 17 905 922 2001 Zhang Y S A Klein J Boyle and G G Mace Evaluation of tropical cloud and pre cipitation statistics of CAM3 using CloudSat and CALIPSO data J Geophys Res doi 10 1029 2009JD012006 2010 21
8. This namelist is located in file cmor cosp_cmor_nl txt and Table 2 details its variables It is expected that this namelist will be expanded in COSPv1 3 to include all the attributes that are required by the CMIP5 Table 2 CMOR namelist INPATH Directory where the MIP table is located OUTPATH Directory where the outputs will be written START DATE Experiment start date MODEL ID String with your model name or id EXPERIMENT_ID Type of experiment This has to be one of those listed in the variable expt_id_ok in the MIP table INSTITUTION Your institution SOURCE Data source e g model version id of your model run CALENDAR Calendar type used by the model REALIZATION Realisation within an ensemble of runs for a given experi ment CONTACT Contact details HISTORY What CMOR has done to the user supplied data e g transforming its units or rearranging its order to be consis tent with the MIP requirements You can live this blank COMMENT Extra comments that may help the interpretation of the data REFERENCES Papers or other references describing the model TABLE Name of the MIP table This has to be consistent with the mode used to run COSP which is defined by the input file Different tables are needed for 1D and 2D models The current list of table distributed with COSP are CMIP5_cfShr MIP table for 1D mode This is a modified version with extra variables of the CMIP5 cf3hr distributed with the CMOR2 library
9. alue float mr_ccliq level mr_ccliq units mr_ccliq long_name mr_ccliq FillValue float mr_ccice level mr_ccice units mr_ccice long_name mr_ccice FillValue float fl_lsrain level kg fl_lsrain long_name fl_lsrain FillValue float fl_lssnow level kg fl_lssnow long_name fl_lssnow FillValue float fl_lsgrpl level kg fl_lsrain units fl_lssnow units fl_lsgrpl units lon fic_humidity Of lon ive_humidity Of lon 3 total_cloud_amount 30f lon ective_cloud_amount 30f lat lon kg mixing ratio_large_scale_cloud_liquid 1 e 30f lat lon 3 kg kg mixing_ratio_large_scale_cloud_ice 1 e 30f lat lon kg kg mixing_ratio_convective_cloud_liquid 1 e 30f lat lon 3 kg kg mixing_ratio_convective_cloud_ice 1 e 30f lat lon m 2 8 1 flux_large_scale_cloud_rain 1 e 30f lat lon m 2 s 1 flux_large_scale_cloud_snow 1 e 30f lat lon m 2 s 1 18 fl_lsgrpl long_name flux_large_scale_cloud_graupel fl_lsgrpl FillValue 1 e 30f float fl_ccrain level lat lon fl_ccrain units kg m 2 s 1 fl_ccrain long_ name flux_convective_cloud_rain fl_ccrain FillValue 1 e 30f float fl_ccsnow level lat lon fl_ccsnow units kg m 2 s 1 fl_ccsnow long_name flux_convective_cloud_snow fl_ccsnow FillValue 1 e 30f float orography lat lon
10. as calcu lated by the ISCCP Simulator Ltauisccp Mean Optical Depth as Calculated by the ISCCP Simulator Lpctisccp ISSCP Mean Cloud Top Pressure Flags for CALIPSO simulator outputs Latb532 Lidar Attenuated Total Backscatter 532 nm LcfadLidarsr532 CALIPSO Scattering Ratio CFAD Lclcalipso2 CALIPSO Cloud Fraction Undetected by CloudSat Lclcalipso CALIPSO Cloud Area Fraction Lclhcalipso CALIPSO High Level Cloud Fraction Lcllcalipso CALIPSO Low Level Cloud Fraction Lclmcalipso CALIPSO Mid Level Cloud Fraction Lcltcalipso CALIPSO Total Cloud Fraction 6 LparasolRefl PARASOL Reflectance LlidarBetaMol532 Lidar Molecular Backscatter 532 nm Flags for CloudSat simulator outputs Lcfaddbze94 CloudSat Radar Reflectivity CFAD Ldbze94 CloudSat Radar Reflectivity Flags for CALIPSO CloudSat combined outputs Lcltlidarradar Lidar and Radar Total Cloud Fraction Flags for other outputs Lfracout Subcolumn output from SCOPS Flags for MODIS simulator outputs Lclhmodis MODIS High Level Cloud Fraction Lclimodis MODIS Ice Cloud Fraction Lcllmodis MODIS Low Level Cloud Fraction Lelmmodis MODIS Mid Level Cloud Fraction Lelmodis MODIS Cloud Area Fraction Lcltmodis MODIS Total Cloud Fraction Lclwmodis MODIS Liquid Cloud Fraction Liwpmodis MODIS Cloud Ice Water Path Llwpmodis MODIS Cloud Liquid Water Path Lpctmodis MODIS Cloud Top Pressure Lreffclimodis MODIS Ice Cloud Particle Size Lreffclwmodis MODIS Liquid Cloud Particle Siz
11. at float float float float float float orography point orography long_name orography orography _FillValue 1 e 30f orography units m psfc point psfc long_name surface_pressure psfc _FillValue 1 e 30f psfc units Pa height level point height long_name height_in_full_levels height _FillValue 1 e 30f height units m height_half level point height_half long_name height_in_half_levels height_half _FillValue 1 e 30f height_half units m T_abs level point T_abs long_name air_temperature T_abs _FillValue 1 e 30f T_abs units K qv level point qv long_name specific_humidity qv _FillValue 1 e 30f qv units 4 rh level point rh long_name relative _humidity_liquid_water rh _FillValue 1 e 30f rh units pfull level point pfull long_name p_in_full_ levels pfull _FillValue 1 e 30f pfull units Pa phalf level point phalf long_name p_in_half_levels phalf _FillValue 1 e 30f phalf units Pa mr_lsliq level point mr_lsliq long_name mixing_ratio_large_scale_cloud_liquid mr_lsliq _FillValue 1 e 30f mr_lsliq units kg kg mr_lsice level point mr_lsice long_name mixing_ratio_large_scale_cloud_ice 14 3 mr_lsice _FillValue 1 e 30f mr_lsice units kg kg float mr_ccliq level point mr_ccliq long_name mixing_ratio_convect
12. at emsfc_lw emsfc_lw units 1 emsfc_lw long_name Surface emissivity at 10 5 micron fraction emsfc_lw FillValue 1 e 30f float mr_ozone level lat lon mr_ozone units kg kg mr_ozone long_name mass_fraction_of_ozone_in_air mr_ozone FillValue 1 e 30f float u_wind lat lon u_wind units m s 1 u_wind long_name eastward_wind u_wind FillValue 1 e 30f float v_wind lat lon v_wind units m s 1 v_wind long_name northward_wind v_wind FillValue 1 e 30f References Chepfer H S Bony D Winker M Chiriaco J L Dufresne and G S ze Use of CALIPSO lidar observations to evaluate the cloudiness simulated by a climate model Geophys Res Lett 35 L15 704 doi 10 1029 2008GL034207 2008 Haynes J M QuickBeam radar simulation software Colorado State University Fort Collins CO USA v1 1 ed 2007 Haynes J M R T Marchand Z Luo A Bodas Salcedo and G L Stephens A multi purpose radar simulation package Quickbeam Bull Am Meteorol Soc 88 11 1723 1727 doi 10 1175 BAMS 88 11 1723 2007 20 Khairoutdinov M F and D A Randall Cloud resolving modeling of the ARM summer 1997 IOP Model formulation results uncertainties and sensitivities J Atmos Sci 60 4 607 625 doi 10 1175 1520 0469 2003 060 0607 CRMOTA 2 0 CO 2 2003 Klein S A and C Jakob Validation and sensitivities of frontal clouds simulated by the ECMWF model
13. de is a gridded lat lon input suitable for model outputs This is the Common Data Language CDL structure of the COSP input NetCDF file in 1D mode netcdf cosp_input_um dimensions point 1236 level 50 hydro 9 variables short year point year long_name year year _FillValue 32767s year units yr byte month point 12 month long_name month month _FillValue 127b byte day point day long_name day day _FillValue 127b day units day byte hour point hour long_name hour hour _FillValue 127b hour units hr byte minute point minute long_name minute minute _FillValue 127b minute units min float second point second long_name second second _FillValue 1 e 30f second units s float t point t long_name t t _FillValue 1 e 30f t units min float tUM point tUM long_name tUM tUM _FillValue 1 e 30f tUM units min float lst point lst long_name lst lst _FillValue 1 e 30f lst units h float lon point lon long_name longitude lon _FillValue 1 e 30f lon units degree_east float lat point lat long_name latitude lat _FillValue 1 e 30f lat units degree_north float landmask point landmask long_name landmask landmask _FillValue 1 e 30f landmask units i 13 float float float float flo
14. e Ltauilogmodis MODIS Ice Cloud Optical Thickness Log10 Mean Ltauimodis MODIS Ice Cloud Optical Thickness Ltautlogmodis MODIS Total Cloud Optical Thickness Log10 Mean Ltautmodis MODIS Total Cloud Optical Thickness Ltauwlogmodis MODIS Liquid Cloud Optical Thickness Log10 Mean Ltauwmodis MODIS Liquid Cloud Optical Thickness Flags for RTTOV outputs Ltbrttov Mean clear sky brightness temperature as calculated by RTTOV 3 Microphysical settings This section discusses how to set up the COSP microphysical settings This is is particularly important for the computation of the radar reflectivities as they are strongly dependent on the paricle size This section should be read in conjunction with Section 4 of the QuickBeam User s Guide In the following discussion lets assume that the particle size distribution PSD n D Shttp reef atmos colostate edu haynes radarsim userguide pdf for a particle of diameter D is defined as a gamma function DD Nog D e 2P 1 where no is the intercept parameter is the slope parameter a is the constant shape parameter x can be either R for rain a for aggregates c for ice crystals or g for graupel For a single moment scheme the intercept parameter is assumed constant or a simple function of Ax Nor gga 2 where naz and np are constants The terminal fall velocity of a precipitating particle V D can be expressed as a function of diameter V D D 2 3 where cr d
15. e a more general version of this conversion similar to the one presented in this document 3 3 Setting the HCLASS table The microphysical assumptions for the radar simulation in COSP are stored in the HCLASS table in cosp_constants f90 The meaning of the HCLASS constants are given in the Quick beam Users guide Haynes 2007 For the sake of completeness here we also give an overview and the settings The HCLASS table consists of several lines each one stored in a different variable These variables are vectors with as many elements as number of hydrom eteors so that the settings for each hydrometeor can be set up independently These variables are e HCLASS_TYPE Set to 1 for modified gamma distribution 2 for exponential distribution 10 3 for power law distribution 4 for monodisperse distribution 5 for lognormal distribution Set to a negative number to ignore the hydrometeor class defined in that position HCLASS COL Reserved for future use value is ignored HCLASS PHASE Set to 0 for liquid 1 for ice HCLASS CP Not used in COSP HCLASS_DMIN The minimum drop size for this class um ignored for monodisperse HCLASS_DMAX The maximum drop size for this class um ignored for monodisperse HCLASS_APM The a coefficient in in the mass diameter relationship If used then set HCLASS_RHO to 1 e HCLASS_BPM The b coefficient in in the mass diameter relationship If used then set HCLASS_RHO to 1 e HCLASS_RHO
16. for the off line CFMIP2 experiments CMIP5_cf3hr cmor1 the same as CMIP5_cf8hr but to be used when linking COSP with CMOR1 3 COSP table 2D table to be used in 2D mode COSP_table_2D cmor1 same as COSP_table_2D but to be used when linking COSP with CMOR1 3 MAXTSTEPS Maximum number of records that can be recorded to the output files CMOR will issue an error and stop if you try to write more records 2 3 COSP_OUTPUT namelist This is the namelist that sets up output related variables see Table 3 It controls the instrument simulators that are run and the list of variables to be written to file If a simulator is switched off then none of its outputs are written out independently of the status of the logical flags of the output variables associated with that particular simulator Table 3 COSP_OUTPUT namelist Logical flags that control which simulators are run Lradar_sim Llidar_sim Lisccp_sim Lmisr_sim Lmodis sim Lrttov_sim Flags for ISCCP simulator outputs Lalbisccp Lboxptopisccp Lboxtauisccp Lclisccp Lcltisccp Lmeantbclrisccp Lmeantbisccp ISSCP Mean Cloud Albedo Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator Optical Depth in Each Column as Calculated by the ISCCP Simulator ISSCP Cloud Area Fraction ISSCP Total Cloud Fraction Mean clear sky 10 5 micron brightness temperature as cal culated by the ISCCP Simulator Mean all sky 10 5 micron brightness temperature
17. from observations The scheme that we use to break the grid box mean profiles of cloud water contents is the Subgrid Cloud Overlap Profile Sampler SCOPS a technique developed for the International Satellite Cloud Climatology Project ISCCP simulator Klein and Jakob 1999 Webb et al 2001 SCOPS uses a pseudo random sampling process fully consistent with the maximum random and maximum random cloud overlap assumptions used in many models e g Pincus et al 2005 Maximum overlap is applied to the convective cloud and maximum random is used for large scale cloud Zhang et al 2010 have developed a simple algorithm that pro vides sub grid distribution of precipitation fluxes compatible with the cloud distribution output by SCOPS and the gridbox mean precipitation fluxes simulated by the model The current version of COSP includes simulators for the following instruments CloudSat radar Haynes et al 2007 CALIPSO lidar Chepfer et al 2008 ISCCP Klein and Jakob 1999 Webb et al 2001 the Multiangle Imaging SpectroRadiometer MISR and the Moder ate Resolution Imaging Spectroradiometer MODIS The fast radiative transfer code RTTOV Saunders et al 1999 can also be linked to COSP to produce clear sky brightness tempera tures for many different channels of past and current infrared and passive microwave radiome ters The Climate Model Output Rewriter CMOR library is used to write the ouputs to NetCDF files that comply with the
18. input NetCDF file in 2D mode netcdf cosp_input_um_2d dimensions lon 17 lat 9 level 38 bnds 2 hydro 9 variables float lon lon lon axis X lon units degrees_east lon long_name longitude lon bounds lon_bnds float lat lat lat axis Y lat units degrees_north lat long_name latitude lat bounds lat_bnds float lon_bnds lon bnds float lat_bnds lat bnds float height level lat lon height units m height long_name height_in_full_levels height FillValue 1 e 30f float pfull level lat lon pfull units Pa pfull long_ name p_in_full_levels pfull FillValue 1 e 30f float phalf level lat lon phalf units Pa phalf long_name p_in_half_levels phalf FillValue 1 e 30f float T_abs level lat lon T_abs units K T_abs long_name air_temperature T_abs FillValue 1 e 30f 17 float qv level lat qv units kg kg qv long_name speci qv FillValue 1 e 3 float rh level lat nan rh long_name relat rh FillValue 1 e 3 float tca level lat rh units tca units 1 tca long_name tca FillValue 1 e float cca level lat cca units 1 cca long_name conv cca FillValue 1 e float mr_lsliq level kg mr_lsliq long_name mr_lsliq units mr_lsliq FillValue float mr_lsice level mr_lsice units mr_lsice long_name mr_lsice FillV
19. ive_cloud_liquid mr_ccliq _FillValue kg kg 1 e 30f gt mr_ccliq units mr_ccice level point mr_ccice long_name mixing _ratio_convective_cloud_ice mr_ccice _FillValue kg kg float 1 e 30f mr_ccice units float fl_lsrain level point fl_lsrain long_name flux_large_scale_cloud_rain fl_lsrain _FillValue kg m 2 s 1 1 e 30f fl_lsrain units fl_lssnow level point flux_large_scale_cloud_snow float fl_lssnow long_name fl_lssnow _FillValue 1 e 30f fl_lssnow units kg m 2 s 1 float fl_lsgrpl level point fl_lsgrpl long_name flux_large_scale_cloud_graupel fl_lsgrpl _FillValue 1 e 30f fl_lsgrpl units kg m 2 s 1 float fl_ccrain level point fl_ccrain long_name flux_convective_cloud_rain fl_ccrain _FillValue 1 e 30f fl_ccrain units kg m 2 s 1 float fl_ccsnow level point fl_ccsnow long_name flux_convective_cloud_snow fl_ccsnow _FillValue 1 e 30f fl_ccsnow units kg m 2 s 1 float tca level point tca long_name total_cloud_amount tca _FillValue 1 e 30f tca units 0 1 float cca level point convective_cloud_amount cca long_name cca _FillValue 0 1 1 e 30f cca units Reff hydro level point Reff long_name hydrometeor_effective_radius Reff _FillValue Reff units m float 1 e 30f 15 float dtau_s level point dtau_s long_name
20. nce and Eq 10 Eq A19 from Khairoutdinov and Randall 2003 can be easily reproduced If your model does not follow the Khairoutdinov and Randall 2003 microphysics then you will have to set use_precipitation_fluxes false and fill in the arrays gbx mr_hydro i with the precipitation mixing ratios in cosp_test i is the index of each precipitation class Model COSP manual Khairoutdinov and Randall 2003 nz D Nax Nom Nba 0 Qx 0 Vex Nri V D CH Am dz bm Gz 0 5 mz D Ay Pm by 3 Table 4 Correspondence between parameters in the formulation used in this manual and those used in Khairoutdinov and Randall 2003 They are classified in three groups according to their relation to the PSD terminal fall speed or mass size relationship _LSRAIN LSSNOW I CVRAIN _ CVSNOW LSGRPL The standard list of hydrometeors is defined in cosp_constants f90 integer parameter I_LSCLIQ integer parameter I_LSCICE integer parameter I_LSRAIN integer parameter I_LSSNOW integer parameter I_CVCLIQ integer parameter I_CVCICE integer parameter I_CVRAIN integer parameter I_CVSNOW AN DOF U NB integer parameter I_LSGRPL Another possibility involves modifying 11n1 pf_to_mr so that the flux to mixing ratio con version is appropriate for your model Depending on the model s microphysics this option may be simpler than passing mixing ratios as inputs COSP v1 3 will includ
21. orography units m orography long_name orography orography FillValue 1 e 30f float landmask lat lon landmask units 1 landmask long_name land_mask landmask FillValue 1 e 30f float height_half level lat lon height_half units m height_half long_name height_in_half_levels height_half FillValue 1 e 30f float psfc lat lon psfc units Pa psfc long_name surface_pressure psfc FillValue 1 e 30f float Reff hydro level lat lon Reff units m Reff long_name hydrometeor_effective_radius Reff FillValue 1 e 30f float dtau_s level lat lon dtau_s units 1 dtau_s long_name Optical depth of stratiform cloud at 0 67 micron dtau_s FillValue 1 e 30f float dtau_c level lat lon dtau_c units 1 dtau_c long_name Optical depth of convective cloud at 0 67 micro dtau_c FillValue 1 e 30f float dem_s level lat lon dem_s units i dem_s long_name Longwave emissivity of stratiform cloud at 10 5 micron dem_s FillValue 1 e 30f 19 float dem_c level lat lon dem_c units 1 dem_c long_name Longwave emissivity of convective cloud at 10 5 micron dem_c FillValue 1 e 30f float skt lat lon skt units K skt long_name Skin temperature skt FillValue 1 e 30f float sunlit lat lon sunlit units 1 sunlit long_name Day points sunlit FillValue 1 e 30f flo
22. p net Table 1 COSP_INPUT namelist General configuration variables CMOR_NL NPOINTS NPOINTS_IT NCOLUMNS NLEVELS USE_VGRID NLR CSAT_VGRID FINPUT Name of CMOR namelist Section 2 2 Number of gridpoints to be processed This has to coincide with the number of points of the NetCDF input file in 2D mode lat lon For 1D curtain mode there is no restric tion Maximum number of gridpoints to be processed in one iter ation This helps to reduce the amount of memory used by COSP If you find memory faults reduce this number Number of subcolumns used for each profile Number of levels This must be the same number as in the input NetCDF file If false the outputs are written on model levels If this is set to true then a vertical grid evenly spaced in altitude is used If true then you need to define number of levels with Nir Number of levels in statistical outputs only used if USE_VGRID true Set to true for CloudSat vertical grid This is just a stan dard grid of 40 levels evenly spaced at CloudSat vertical resolution 480 m This only applies if USE _VGRID true Input NetCDF file This is the input file with all the input vari ables to that your COSP executable will read and process Inputs related to radar simulations RADAR_FREQ SURFACE_RADAR use mie tables use_gas_abs do_ray melt_lay k2 use _reff use_precipitation_fluxes Frequency GHz used in the radar sim
23. propriate to compare to ISCCP IR only algortihm i e you can compare to nighttime ISCCP data with this option ISCCP_TOPHEIGHT_DIRECTIONdirection for finding atmosphere pressure level with interpolated temperature equal to the radiance determined cloud top temperature 1 find the lowest altitude highest pres sure level with interpolated temperature equal to the radiance determined cloud top temperature 2 find the highest altitude lowest pressure level with interpolated temperature equal to the radiance de termined cloud top temperature This is the default value since V4 0 of the ISCCP simulator ONLY APPLICABLE IF top_height EQUALS 1 or 3 Inputs related to RTTOV simulations Platform Satellite Instrument Nchannels Channels Surfem ZenAng CO2 CH4 N20 CO Satellite platform number Satellite Instrument Number of channels to be computed Channel numbers please be sure that you supply Nchan nels Surface emissivity please be sure that you supply Nchan nels Satellite Zenith Angle degrees Mixing ratio of C O2 Mixing ratio of CH4 Mixing ratio of N2O Mixing ratio of CO 2 2 CMOR namelist The CMOR2 library is used to write the ouputs to NetCDF files that comply with the CF Meta data Convention and fulfill the requirements of the climate community s standard model exper iments for CMIP5 The namelist CMOR is used to passed all the metadata that the calls to the CMOR library require
24. ulations Radar position surface 1 spaceborne 0 Use a precomputed lookup table yes 1 no 0 Include gaseous absorption yes 1 no 0 Calculate output Rayleigh refl 1 not 0 This should be set to 0 as the Rayleigh reflectivity is not output by COSP Melting layer model off 0 on 1 Dielectric factor of water 1 use frequency dependent default True if you want effective radius to be used by radar simu lator always used by lidar true True if precipitation fluxes are input to the algorithm Inputs related to lidar simulations Nprmts_ max_ hydro Naero Max number of parameters for hydrometeor size distribu tions Number of aerosol species Not used Nprmts_max_aero lidar ice type OVERLAP Max number of parameters for aerosol size distributions Not used Ice particle shape in lidar calculations 0 ice spheres 1 ice non spherical Overlap type 1 max 2 rand 3 max rand Inputs related to ISCCP simulations ISCCP_TOPHEIGHT 1 adjust top height using both a computed infrared bright ness temperature and the visible optical depth to adjust cloud top pressure Note that this calculation is most ap propriate to compare to ISCCP data during sunlit hours 2 do not adjust top height that is cloud top pres sure is the actual cloud top pressure in the model 3 adjust top height using only the computed infrared brightness temperature Note that this calculation is most ap
25. z hy and Gz are constants p is the air density kg m and po is a reference density of 1 29 We assume a power law relating the mass of the particle to the diameter M D a D 4 The mass diameter relation for rain simply assumes a spherical drop with a density equal to that for liquid water 1000 kg m 3 1 Effective radius COSP requires effective radius as input for CALIPSO and CloudSat Default values can be used although it is recommended to use values that are consistent with the model s micro physics You can use the default values by setting to zero the input array of effective radii The defaults are 30 um for the lidar and the values defined in HCLASS_P1 for CloudSat see details below In order to compute the effective radius it is necessary to be able to infer the particle size distribution This requires to being able to obtain the parameter Az from the model variables specific humidities or precipitation fluxes The ith moment of the PSD is given by i aoe ax i 1 Hz f D Nz D dD Mor al 5 When the hydrometeor mixing ratio is available the value of A is given by fr 1 2 Re me j Pr For precipitation fluxes the flux can be related to the PSD by RF T MD Dina D dD 7 0 Using Eqs 1 2 4 and solving this integral for A gives 1 Ga az tbr tde npp tI ax Cy 2 naal x by de 1 SiS 8 7 F x The effective radius is then given by 3 Hx T ax 4
Download Pdf Manuals
Related Search