Dear Beate,
I don't believe this issue has been addressed yet by the convention. Clearly
you will need to save the half-level a's and b's, along with ps, but there is no
way currently to indicate that this is a hybrid sigma coordinate, but with full
level p's calculated from the half-level p's using the formula you supplied.
We'll have to give this further thought ... any suggestions out there?
Karl
Beate.Geyer at gkss.de wrote:
>
> Dear member of CF community,
> I got data from GCM-output to convert to netCDF with CF-standard.
> The description of atmosphere_hybrid_sigma_pressure_coordinate in
> model-documentation is the following:
> - definition for half-integer level is (like described at
> http://www.cgd.ucar.edu/cms/eaton/cf-metadata/CF-current.html#sub_app)
> p(n,k-1/2,j,i) = ap(k-1/2) + b(k-1/2)*ps(n,j,i) and
>
> -the definition for the full levels (exponential interpolation instead
> of linear) is
> p(k)=
> exp(1/p(k-1/2)-p(k+1/2))*(p(k-1/2)*ln(p(k-1/2)-p(k+1/2)*ln(p(k+1/2))-1
> for k<k(max)
> = 1/2p(k-1/2)
> for k=k(max)
> where p is time and location depended like above (for half-integer level).
>
> How is the correct way to include the latter information in the
> netCDF-File?
>
> Thanks for your help,
> Beate Geyer
>
> ====================================================
> Dr. Beate Geyer
> GKSS Research Centre
> Institute for Coastal Research
> Max-Planck-Str.
> D-21502 Geesthacht, Germany
>
> E-Mail: Beate.Geyer at gkss.de
> Phone: 04152-872016
> ====================================================
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> CF-metadata mailing list
> CF-metadata at cgd.ucar.edu
> http://www.cgd.ucar.edu/mailman/listinfo/cf-metadata
Received on Wed Sep 21 2005 - 15:51:23 BST
