Re: fatal:Div: operator failed, can't continue

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Jan 31 2011 - 08:36:34 MST

Please read the documentation:
http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_rh.shtml

==
"qv" must be kg/kg , you are passing "g/kg"
"P" must be "Pa", you are passing "hPa"

On 1/31/11 7:09 AM, wei huang wrote:
> Dear NCL users,
>
> I got an error to calculating saturated mixing ratio from the wrf water
> vapor mixing ratio. I don’t know my assumption is correct or not. Could
> anyone help me to figure out the problem or any advice about it?
>
> Thanks
>
> wei
> a = addfile("wrfout_02.nc <http://wrfout_02.nc>","r")
>
> P= a->P; perturbation
>
> Pb = a->PB; base state pressure
>
> P= P + Pb; total pressure
>
> T = a->T; perturbation potential temperature
>
> tk = T + 300.
>
> qv = a->QVAPOR/1000; water vapor mixing ratio(g/kg)
>
> RH = wrf_rh (qv, P*.01, tk) ;P*0.01 converts Pa to hPa
>
> printVarSummary(RH)
>
> qvs = (qv/RH)*100; saturated water vapor mixing ratio(g/kg)
>
> printVarSummary(qvs)
>
> end
>
> ;*****************************************
>
> Variable: RH
>
> Type: float
>
> Total Size: 102667500 bytes
>
> 25666875 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes:[24] x [39] x [south_north | 160] x [west_east | 105]
>
> Coordinates:
>
> Number Of Attributes: 2
>
> units :%
>
> description : Relative Humidity
>
> fatal:divide: Division by 0, Can't continue
>
> fatal:Div: operator failed, can't continue
>
> fatal:Execute: Error occurred at or near line 16 in file
>
>
>
> _______________________________________________
> 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
Received on Mon Jan 31 08:36:40 2011

This archive was generated by hypermail 2.1.8 : Mon Jan 31 2011 - 10:38:24 MST