Re: problem with dpres_hybrid_ccm

Hybrid coefficients come in 2 forms.

NCL assumes the coefficients are of the form:
    p = a(k)*p0 + b(k)*ps
This is the form 'our' climate model uses.

The 'other' form is
    p = A(k) + B(k)*ps

===
  B(k) = b(k)
  A(k) = a(k)*p0
===

Hence, to convert to the form expected by dpres_hybrid_ccm
and press_hybrid_ccm ...

hyai=f2->hyai
hyam=f2->hyam
     :
p0 = 100000
hyam= hyam/p0
hyai= hyai/p0

----
I always recommend that people use 'printVarSummary(...)'
when debugging a code.
A simple script is attached.
Regards
D
On 1/15/12 10:33 AM, laura.niederdrenk wrote:
> Hello,
>
> I have a problem with the function dpres_hybrid_ccm. It seems to produce
> nonsense values in my case. I attach a simple example ncl-script with
> data from one grid point only (and one time step). Although all input
> values seem to be correct the output of dpres_hybrid_ccm is definitely
> wrong. Could you please tell me where my mistake is!?
>
> Thank you very much in advance.
>
> Laura
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk