Main menu

Navigation

Adding new molecule (OCS) to the initial state file

2 posts / 0 new
Last post
myous005@...
Adding new molecule (OCS) to the initial state file

Greeting,

I'm trying to run the WACCM for OCS molecule and two of its isotopologues. however, in the input data file for the initial state (cami_2000-01-01_10x15_L66_c041121.nc), doesn't include (lev='hybrid level at midpoints (1000*(A+B)) matrix) for OCS molecule. my question would be how can I add it for a new molecule such as OCS to the input Netcdf file? or is it necessary to have lev matrix for OCS?

below is the content of the Netcdf file.

Format:

           classic

Dimensions:

           time = 1

           lat  = 19

           lon  = 24

           lev  = 66

           slon = 24

           slat = 18

           ilev = 67

Variables:

    LANDM_COSLAT

           Size:       24x19x1

           Dimensions: lon,lat,time

           Datatype:   double

           Attributes:

                       long_name  = 'land ocean transition mask: ocean (0), continent (1), transition (0-1)'

                       units      = 'none'

                       _FillValue = 1e+36

                       from_hires = 'true'

    time        

           Size:       1x1

           Dimensions: time

           Datatype:   single

           Attributes:

                       long_name = 'time'

                       units     = 'days since 0000-01-01 00:00:00'

                       calendar  = 'noleap'

    lat         

           Size:       19x1

           Dimensions: lat

           Datatype:   single

           Attributes:

                       long_name = 'latitude'

                       units     = 'degrees_north'

    lon         

           Size:       24x1

           Dimensions: lon

           Datatype:   single

           Attributes:

                       long_name = 'longitude'

                       units     = 'degrees_east'

    Q           

           Size:       24x66x19x1

           Dimensions: lon,lev,lat,time

           Datatype:   double

           Attributes:

                       long_name = 'Q constituent tracer'

                       units     = 'kg/kg'

    lev         

           Size:       66x1

           Dimensions: lev

           Datatype:   single

           Attributes:

                       long_name = 'hybrid level at midpoints (1000*(A+B))'

    H2O         

           Size:       24x66x19x1

           Dimensions: lon,lev,lat,time

           Datatype:   double

           Attributes:

                       long_name = 'H2O'

                       units     = 'kg/kg'

    CO2         

           Size:       24x66x19x1

           Dimensions: lon,lev,lat,time

           Datatype:   double

           Attributes:

                       long_name = 'CO2'

                       units     = 'kg/kg'

.

.

.

.

   

 

Regards,

 

Mehdi

mmills

Please post your questions about WACCM in one of the WACCM forums so that I will see them:

https://bb.cgd.ucar.edu/forums/science-issues-0

https://bb.cgd.ucar.edu/forums/software-and-run-time-issues-0

If your IC file does not include OCS, you can run the model for a decade or more to wait for OCS to be advected from the surface (lower boundary condition) to the stratosphere. Alternately, you could add OCS to your IC file as a uniform value equal to your lower boundary condition mixing ratio. It would quickly oxidize in the stratosphere to produce a reasonable profile. Or you can provide a more realistic initial condition based on observations or other information.

Mike Mills
WACCM Liaison
Atmospheric Chemistry Division
NCAR Foothills Lab
Boulder, Colorado USA

Log in or register to post comments

Who's new

  • wurjicnip@...
  • 479266557@...
  • fmoyano@...
  • asayeed@...
  • nbren12@...