Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

A quetion about Aerosol Mass data?

njujyq@163_com

New Member
The file:

AerosolMass_V_64x128_clim_c031022.nc

have the varable of Aerosol mass with cordinate

[time | 12] x [lev | 28] x [lat | 64] x [lon | 128]

and the vertical cordinate is a missing data

it has a vertical cordinate virable

Variable: hybi (file variable)
Type: double
Total Size: 232 bytes
29 values
Number of Dimensions: 1
Dimensions and sizes: [ilev | 29]
Coordinates:
Number Of Attributes: 0
(0) 0
(1) 0.006570000201463699
(2) 0.01386000029742718
(3) 0.02308999933302402
(4) 0.0346900001168251
(5) 0.04919999837875366
(6) 0.06723000109195709
(7) 0.08945000171661377
(8) 0.1165399998426437
(9) 0.1491599977016449
(10) 0.187830001115799
(11) 0.2328599989414215
(12) 0.2842099964618683
(13) 0.3413699865341187
(14) 0.4033400118350983
(15) 0.4686000049114227
(16) 0.535290002822876
(17) 0.6013500094413757
(18) 0.6648200154304504
(19) 0.7240099906921387
(20) 0.7777299880981445
(21) 0.8252699971199036
(22) 0.8664199709892273
(23) 0.9013500213623047
(24) 0.9305400252342224
(25) 0.9545900225639343
(26) 0.9741799831390381
(27) 0.9900000095367432
(28) 1

which is 29 levels

what is the really vetical cordinate of the aerosol data?

With the understanding ,I can generate aerosl data myself

Many thanks!
 

strandwg

Moderator
Staff member
superjyq said:
The file:

AerosolMass_V_64x128_clim_c031022.nc

have the varable of Aerosol mass with cordinate

[time | 12] x [lev | 28] x [lat | 64] x [lon | 128]

and the vertical cordinate is a missing data [...]

I looked at the original file from which it was derived (as noted in the file's header), namely:

AerosolMass_192x94_L28_clim_c030418.nc

and it shows (via 'ncdump -c'):

lev = 2.7, 10.1, 18.3, 28.8, 41.8, 58, 78.2, 102.8, 132.6, 168.2, 210.1,
258.2, 312.5, 372, 435.7, 501.7, 568.1, 632.9, 694.3, 750.8, 801.4,
845.8, 883.8, 915.9, 942.5, 964.4, 982.1, 995 ;

This may or may not be useful to you.
 

njujyq@163_com

New Member
But I do not find the information in the data

So the data is in pressure cordinate

And with surface pressure and hybi

then the CAM3 will load it and change it to the sigma levels

Am I right?
 

eaton

CSEG and Liaisons
The aerosol mass file in question is read by code in the file src/physics/cam1/prescribed_aerosols.F90. That is the place to find details of how the pressure interpolation is being done.

It looks to me that the aerosol data is on a sigma coordinate, hence only hybi part of CAM's hybrid coordinate is present. hybi is the sigma part of the hybrid coordinate (i.e., it multiplies the surface pressure), and it is the values at the layer interfaces rather than the layer midpoints. So the file contains 28 levels with interface pressures given by the values of hybi and PS in the file.

In the prescribed_aerosols module the data on the sigma levels are being interpolated to the actual pressures in the CAM physics columns which is contained in the pint variable.
 

njujyq@163_com

New Member
eaton said:
The aerosol mass file in question is read by code in the file src/physics/cam1/prescribed_aerosols.F90. That is the place to find details of how the pressure interpolation is being done.

It looks to me that the aerosol data is on a sigma coordinate, hence only hybi part of CAM's hybrid coordinate is present. hybi is the sigma part of the hybrid coordinate (i.e., it multiplies the surface pressure), and it is the values at the layer interfaces rather than the layer midpoints. So the file contains 28 levels with interface pressures given by the values of hybi and PS in the file.

In the prescribed_aerosols module the data on the sigma levels are being interpolated to the actual pressures in the CAM physics columns which is contained in the pint variable.

Many thanks
 
I'm trying to figure out the atmospheric heating from that aerosol dataset AerosolMass_V_64x128_clim_c031022.ncI have a CAM3.5 that would read the cyclical aerosol forcing on this grid: 1.9x2.5_L26. Is there anything ready for that? I could not find it.
 

eaton

CSEG and Liaisons
There is a 1.9x2.5 dataset available in the data repository here:https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/inputdata/atm/cam/rad/AerosolMass_V_1.9x2.5_clim_c040811.nc 
 
Top