I use your CAMS data for my fog forecasting project with the AROME model, integrating the effects of CCN and IFN on the forecast. I previously used the below script for directly downloading aerosols in Lambert grid for my domain. However, with the new ADS, the data is now downloaded only in lat/lon projection, which disrupts the continuity of my work.
Could you please help me resolve this issue?
script :
"import cdsapi
import sys
c = cdsapi.Client()
c.retrieve(
‘cams-global-atmospheric-composition-forecasts’,
{
‘variable’: [
‘dust_aerosol_0.03-0.55um_mixing_ratio’, ‘dust_aerosol_0.55-0.9um_mixing_ratio’, ‘dust_aerosol_0.9-20um_mixing_ratio’,
‘hydrophilic_black_carbon_aerosol_mixing_ratio’, ‘hydrophilic_organic_matter_aerosol_mixing_ratio’, ‘hydrophobic_black_carbon_aerosol_mixing_ratio’,
‘hydrophobic_organic_matter_aerosol_mixing_ratio’, ‘sea_salt_aerosol_0.03-0.5um_mixing_ratio’, ‘sea_salt_aerosol_0.5-5um_mixing_ratio’,
‘sea_salt_aerosol_5-20um_mixing_ratio’, ‘sulphate_aerosol_mixing_ratio’,‘specific_humidity’,‘surface_geopotential’,‘surface_pressure’, ‘temperature’,
‘u_component_of_wind’,‘v_component_of_wind’ ,
],
‘model_level’: [
‘1’, ‘2’, ‘3’,
‘4’, ‘5’, ‘6’,
‘7’, ‘8’, ‘9’,
‘10’, ‘11’, ‘12’,
‘13’, ‘14’, ‘15’,
‘16’, ‘17’, ‘18’,
‘19’, ‘20’, ‘21’,
‘22’, ‘23’, ‘24’,
‘25’, ‘26’, ‘27’,
‘28’, ‘29’, ‘30’,
‘31’, ‘32’, ‘33’,
‘34’, ‘35’, ‘36’,
‘37’, ‘38’, ‘39’,
‘40’, ‘41’, ‘42’,
‘43’, ‘44’, ‘45’,
‘46’, ‘47’, ‘48’,
‘49’, ‘50’, ‘51’,
‘52’, ‘53’, ‘54’,
‘55’, ‘56’, ‘57’,
‘58’, ‘59’, ‘60’,
‘61’, ‘62’, ‘63’,
‘64’, ‘65’, ‘66’,
‘67’, ‘68’, ‘69’,
‘70’, ‘71’, ‘72’,
‘73’, ‘74’, ‘75’,
‘76’, ‘77’, ‘78’,
‘79’, ‘80’, ‘81’,
‘82’, ‘83’, ‘84’,
‘85’, ‘86’, ‘87’,
‘88’, ‘89’, ‘90’,
‘91’, ‘92’, ‘93’,
‘94’, ‘95’, ‘96’,
‘97’, ‘98’, ‘99’,
‘100’, ‘101’, ‘102’,
‘103’, ‘104’, ‘105’,
‘106’, ‘107’, ‘108’,
‘109’, ‘110’, ‘111’,
‘112’, ‘113’, ‘114’,
‘115’, ‘116’, ‘117’,
‘118’, ‘119’, ‘120’,
‘121’, ‘122’, ‘123’,
‘124’, ‘125’, ‘126’,
‘127’, ‘128’, ‘129’,
‘130’, ‘131’, ‘132’,
‘133’, ‘134’, ‘135’,
‘136’, ‘137’,
],
‘date’: sys.argv[1]+‘-’+sys.argv[2]+‘-’+sys.argv[3]+‘/’+sys.argv[4]+‘-’+sys.argv[5]+‘-’+sys.argv[6],
‘leadtime_hour’: int(sys.argv[7]),
‘time’: ‘00:00’,
‘type’: ‘forecast’,
‘area’: [
8.3, -7.7,-7.7,
8.3,
],
‘rotation’: [-56.0 , 3.0,
],
‘grid’:[‘0.8’, ‘0.8’,],
‘format’: ‘grib’,
},
‘aer_’+sys.argv[1]+sys.argv[2]+sys.argv[3]+‘_’+sys.argv[7]+‘.grib’)
"