<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt">Hi Daniele<br><br>your method works: the "set time" command is the crucial point.<br>By defining a new variable for the same field of each new opened file (or by using the .#of-ctl-file at the end of the variable name), one can create the original ensemble set.<br><div><br>Grazie!<br><br>D.<br></div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br><div style="font-family: arial,helvetica,sans-serif; font-size: 13px;"><font size="2" face="Tahoma"><hr size="1"><b><span style="font-weight: bold;">Da:</span></b> Daniele Gandini <daniele.gandini@arpa.piemonte.it><br><b><span style="font-weight: bold;">A:</span></b> GrADS Users Forum <gradsusr@gradsusr.org><br><b><span style="font-weight: bold;">Inviato:</span></b> Mar 24 agosto 2010, 09:33:53<br><b><span
style="font-weight: bold;">Oggetto:</span></b> Re: [gradsusr] ensemble reforecast runs<br></font><br>Dear Daniele Mastrangelo,<br><br>you can do so: put each ensemble member in a different<br>ctl file and then:<br><br>1) open file 1<br>2) display var.1<br>3) open file 2<br>4) set time 12z29Aug1990<br>5) display var.2<br>6) open file 3<br>7) set time 12z29Aug1991<br>8) display var.3<br><br>and so on.<br><br>It looks boring but in a script file you can use<br>a while cycle on ensemble files, i.e.<br>display var.i (i=1,21)<br><br>Buona fortuna.<br><br>Daniele Gandini<br><br><br><br>----- Messaggio originale -----<br>Da: "Daniele Mastrangelo" <<a ymailto="mailto:daniele_mastrangelo@yahoo.it" href="mailto:daniele_mastrangelo@yahoo.it">daniele_mastrangelo@yahoo.it</a>><br>A: "GrADS Users Forum" <<a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Inviato: Luned́, 23 agosto 2010
18:02:59<br>Oggetto: Re: [gradsusr] ensemble reforecast runs<br><br><br><br><br>Hi Wesley, <br><br>thanks for your reply. <br><br>My previous email was not totally clear. Actually, I have 21 ensemble members with the same forecast length; these members are produced initializing the model at the same hour, day, and month (but) of 21 different years (12 UTC 29 aug 1989, 12 UTC 29 aug 1990, ... , 12 UTC 29 aug 2009). <br><br>I tried, however, to apply something similar to your point 2): I aggregated 2 members (e1, e2) in a single ctl file and changed by hand the year of the second member in the EDEF section, making it the same as the first member <br>... <br>tdef 67 linear 12Z29aug1989 12hr <br>zdef 1 linear 1 1 <br>vars 8 <br>... <br>ENDVARS <br>EDEF 2 <br>e1 67 12Z29aug1989 <br>e2 67 12Z29aug1989 ---1989 instead of 1990--- <br>ENDEDEF <br><br>Gribmap ran with no error messages but then, in GrADS, I could not use the e2 fields: "Cannot contour grid - all
undefined values" <br><br>Is there some other way I could try? <br><br>Best regards, <br>Daniele <br><br><br><br><br>Da: Wesley Ebisuzaki <<a ymailto="mailto:Wesley.Ebisuzaki@noaa.gov" href="mailto:Wesley.Ebisuzaki@noaa.gov">Wesley.Ebisuzaki@noaa.gov</a>> <br>A: GrADS Users Forum <<a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>> <br>Inviato: Lun 23 agosto 2010, 16:43:08 <br>Oggetto: Re: [gradsusr] ensemble reforecast runs <br><br>Daniele, <br><br>g2ctl assumes that the length of the forecast is the same for all <br>ensemble members. <br>If not, you can do either <br><br>(1) You could try using the longest ensemble as the 1st ensemble. <br>(It may work.) <br><br>(2) You have to set the TDEF by hand. <br><br>Wesley Ebisuzaki <br><br><br>Daniele Mastrangelo wrote: <br>> Dear GrADS users, <br>> <br>> I'd like to use GrADS v 2.0.a8 to evaluate a set of 21 reforecast <br>>
ensemble runs. These runs are obtained running the same global model <br>> initialized at the same hour and date of 21 different years. <br>> <br>> I obtained a single ctl file through "g2ctl.pl -ens" (g2ctl.pl v <br>> 0.0.4o) but, when I used gribmap to produce the corresponding idx <br>> file, I got the following error message: <br>> <br>> Open Error: ensemble 1 (e2) extends beyond the time axis limits <br>> The data file was not opened. <br>> <br>> I see that the problem is that these runs cover different periods; <br>> moreover, I get the same error when I try to open two of the 21 <br>> different reforecast runs in a normal GrADS session (as file 1 and <br>> file 2). Therefore, my question is: can GrADS treat a set reforecast <br>> runs used as ensemble runs? <br>> If it cannot directly do that, is there a way to bypass the problem? <br>> <br>> Thank you, <br>> Daniele <br>> <br>>
------------------------------------------------------------------------ <br>> <br>> _______________________________________________ <br>> gradsusr mailing list <br>> <a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <br>> <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> <br>> <br><br>_______________________________________________ <br>gradsusr mailing list <br><a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a> <br><br><br>_______________________________________________<br>gradsusr mailing list<br><a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a
href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br><br>_______________________________________________<br>gradsusr mailing list<br><a ymailto="mailto:gradsusr@gradsusr.org" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br></div></div>
</div><br>
</body></html>