Re: ENS_PROB question

• To: gembud@xxxxxxxxxxxxxxxx
• Subject: Re: ENS_PROB question
• From: david.novak@xxxxxxxx
• Date: Tue, 22 May 2007 19:35:25 -0400

Nevermind - I had the cint=10. It should be  cint=0.1

Dave

----- Original Message -----
From: david.novak@xxxxxxxx
Date: Tuesday, May 22, 2007 7:25 pm
Subject: ENS_PROB question

Hello,

I'm trying to calculate the probability of exceedance for a 12 member ensemble. My script is appended. When running this script I get: [DE 4] Number of MEMBERS = 12 Grid file: {namm,gfsm,srefn1m,srefn2m,srefp1m,srefp2m,namw,gfsw,srefn1w,srefn2w,srefp1w,srefp2w} GRID IDENTIFIER:

   TIME1             TIME2         LEVL1 LEVL2   VCORD PARM
060211/1200F036                         0          NONE EXX__1

GAREA: 37;-80;45;-68 SCALE: 0

MINIMUM AND MAXIMUM VALUES     0.00     1.00
LINE CONTOURS:
LEVELS:         0.00

BUT nothing but a blank map

So looking closer at the ENS_PROB, I saw it is actaully a multivariate function. So to make it work, I thought I needed two conditions - in this case PTOT>25 but <150 or GFUNC = ENS_PROB(AND(GT(PTOT,25),LT(PTOT,150))Upon runing this I get.

[DE -5]  No ensemble file referenced in ensemble diagnostic.
[DG -70]  Ensemble field cannot be computed.

What am I mising?
Thanks!
Dave
________________

#!/bin/tcsh

# probability of exceedance model precip
# 1 panel


set g1 = 'namm'
set g2 = 'gfsm'
set g3 = 'srefn1m'
set g4 = 'srefn2m'
set g5 = 'srefp1m'
set g6 = 'srefp2m'
set g7 = 'namw'
set g8 = 'gfsw'
set g9 = 'srefn1w'
set g10 = 'srefn2w'
set g11 = 'srefp1w'
set g12 = 'srefp2w'

set device = xw
#set device = gf
#set device = "ps||8.5;11|c"
set time = '060211/1200F036'

nice gdcntr << endgdcntr
GDATTIM  =  $time
GLEVEL   =  0
GVCORD   =  none
GFUNC    = ENS_PROB(GT(PTOT,25))
GDFILE   = {$g1,$g2,$g3,$g4,$g5,$g6,$g7,$g8,$g9,$g10,$g11,$g12}
CINT     =  1
LINE     =  3
TITLE    =  1// ^
DEVICE   =  $device
PROJ     =  mer
garea    = 37;-80;45;-68
CLEAR    =  y
clrbar = 1/V///.75;.02/
SCALE    =  0

PANEL TEXT = 1

# FINT     = 25;50;75;100;125;150;175;200;225
FINT     = 5;10;15;20;25;30;35;40;45;50;55;60;65;70
# FINT     = 4;8;12;16;20;24;28;32;36;40;44;48
FLINE    = 0;26;4;24;21;22;23;20;17;13;14;11;10;9
CTYPE    =  f
LATLON   = 0
MAP      = 1/1/3
r

endgdcntr

gpend

mv gempak.gif '/home/dnovak/research/Ensembles/12Feb_2006/plots/12Feb_probprecip.gif'

• References:
  • ENS_PROB question
    • From: david . novak