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","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
Received on Mon Jan 31 07:09:31 2011
This archive was generated by hypermail 2.1.8 : Mon Jan 31 2011 - 10:38:24 MST