1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
|
retrieve,
source=compute.req.inpdata,
class=od,
expver=0001,
type=an,
date=20150322,
time=12,
levtyp=sfc,
repres=gg,
level=off,
param=2t,
fieldset=2t
retrieve,
source=compute.req.inpdata,
repres = sh,
levtype = pl,
levelist = 1000,
grid = 640,
gaussian = reduced,
param = 157,
fieldset = rh_unclip
compute,
fieldset=rh,
formula="0.01*min(max(rh_unclip,2),100)"
compute,
fieldset=logx,
formula="log(rh) + 17.502*(2t-273.16)/(2t-32.19)"
compute,
fieldset=2d,
formula="(17.502*273.16-32.19*logx)/(17.502-logx)"
write,
fieldset=2d,
target = compute.req.data
|