<span> </span>Dear Grads Users,<br><br>I am Raju A.<br>I have IMD gridded rainfall.<br>I want to calculate Standard deviation for JJAS.<br>How can i do.<br>Please help me<br>Thanks in advance<br><br>RAJU A
<br>Research Fellow
<br>Indian Institute Of Tropical Meteorology
<br>Dr. Homi Bhabha Road,
<br>PUNE-411008<br><br>On Fri, 23 Jul 2010 15:10:54 +0530  wrote<br>&gt;Send gradsusr mailing list submissions to<br>&gt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;<br>&gt;To subscribe or unsubscribe via the World Wide Web, visit<br>&gt;<a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;or, via email, send a message with subject or body 'help' to<br>&gt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr-request@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr-request@gradsusr.org</a><br>&gt;<br>&gt;You can reach the person managing the list at<br>&gt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr-owner@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr-owner@gradsusr.org</a><br>&gt;<br>&gt;When replying, please edit your Subject line so it is more specific<br>&gt;than "Re: Contents of gradsusr digest..."<br>&gt;<br>&gt;<br>&gt;Today's Topics:<br>&gt;<br>&gt;  1. Re: how to use vint function in GrADS (sushant puranik)<br>&gt;  2. Re: how to use vint function in GrADS (MY)<br>&gt;  3. Re: how to use vint function in GrADS (MY)<br>&gt;  4. Re: Obtaining total values of a dimension on a netCDFfile<br>&gt;   (Becker, Bernd)<br>&gt;<br>&gt;<br>&gt;----------------------------------------------------------------------<br>&gt;<br>&gt;Message: 1<br>&gt;Date: Fri, 23 Jul 2010 11:04:01 +0530<br>&gt;From: sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;<br>&gt;Subject: Re: [gradsusr] how to use vint function in GrADS<br>&gt;To: GrADS Users Forum &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a>&gt;<br>&gt;Message-ID:<br>&gt;&lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=AANLkTingGqD37gATrPBoOtCYi0XXLwcYd8LgPcC+=MjW@mail.gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">AANLkTingGqD37gATrPBoOtCYi0XXLwcYd8LgPcC+=MjW@mail.gmail.com</a>&gt;<br>&gt;Content-Type: text/plain; charset="utf-8"<br>&gt;<br>&gt;Thank you for suggestion. Actually i want to calculate following integral.<br>&gt;<br>&gt;         1000<br>&gt;<br>&gt;Q = (1/g) ? (V q)dp.<br>&gt;<br>&gt;       400<br>&gt;<br>&gt;I used this kind of command in grads<br>&gt;<br>&gt;'define ps=1000'<br>&gt;'define a=vint(ps,vwnd*shum,400)'<br>&gt;'d a/9.81'<br>&gt;with this command i am able to make a plot.<br>&gt;<br>&gt;whether it is correct or not ?<br>&gt;<br>&gt;Thanking you<br>&gt;<br>&gt;<br>&gt;Sushant<br>&gt;<br>&gt;On Fri, Jul 23, 2010 at 10:38 AM, MY &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=funnytyphoon@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">funnytyphoon@gmail.com</a>&gt; wrote:<br>&gt;<br>&gt;&gt; Hi,<br>&gt;&gt;<br>&gt;&gt; please try<br>&gt;&gt;<br>&gt;&gt; 'define ps=1000'<br>&gt;&gt; 'd vint(ps,vwnd*shum,400)'<br>&gt;&gt;<br>&gt;&gt; best wishes,<br>&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt; YING, Ming<br>&gt;&gt; Shanghai Typhoon Institute of China Meteorological Administration<br>&gt;&gt; 166 Puxi Road, Xuhui District, Shanghai, 200030<br>&gt;&gt; PR China<br>&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; 2010/7/23 sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;:<br>&gt;&gt; &gt; Hi<br>&gt;&gt; &gt; I am trying to find vertical integral in GrADS using "vint" command. For<br>&gt;&gt; &gt; that purpose i opened vwnd and specific humidity file as first and second<br>&gt;&gt; &gt; respectively. The q file command for both the variable are as<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; ga-&gt; q file<br>&gt;&gt; &gt; File 1 :mean daily NMC reanalysis (2008)<br>&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt; &gt;   vwnd 8 -999 mean Daily V wind<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; ga-&gt; q file 2<br>&gt;&gt; &gt; File 2 :mean daily NMC reanalysis (2008)<br>&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt; &gt;   shum 8 -999 mean Daily specific humidity<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; after this i used<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; ga-&gt; d vint(1000,vwnd*shum.2,400)<br>&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt; &gt; Error from VINT: Incompatible grids.<br>&gt;&gt; &gt; Contouring: -0 to 2.4e+36 interval 3e+35<br>&gt;&gt; &gt; Segmentation fault<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; after showing this message, the grads window gets closed.<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; Can anyone has solution for this perticular problem.<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; Thanks in advance<br>&gt;&gt; &gt;<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; --<br>&gt;&gt; &gt; Sushant<br>&gt;&gt; &gt;<br>&gt;&gt; &gt;<br>&gt;&gt; &gt; _______________________________________________<br>&gt;&gt; &gt; gradsusr mailing list<br>&gt;&gt; &gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt; &gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt; &gt;<br>&gt;&gt; &gt;<br>&gt;&gt;<br>&gt;&gt; _______________________________________________<br>&gt;&gt; gradsusr mailing list<br>&gt;&gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;<br>&gt;-------------- next part --------------<br>&gt;An HTML attachment was scrubbed...<br>&gt;URL: <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fpipermail%2Fgradsusr%2Fattachments%2F20100723%2F4be8ea92%2Fattachment%2D0001%2Ehtml" target="_new">http://gradsusr.org/pipermail/gradsusr/attachments/20100723/4be8ea92/attachment-0001.html</a> <br>&gt;<br>&gt;------------------------------<br>&gt;<br>&gt;Message: 2<br>&gt;Date: Fri, 23 Jul 2010 15:52:37 +0800<br>&gt;From: MY &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=funnytyphoon@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">funnytyphoon@gmail.com</a>&gt;<br>&gt;Subject: Re: [gradsusr] how to use vint function in GrADS<br>&gt;To: GrADS Users Forum &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a>&gt;<br>&gt;Message-ID:<br>&gt;&lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax@mail.gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax@mail.gmail.com</a>&gt;<br>&gt;Content-Type: text/plain; charset=GB2312<br>&gt;<br>&gt;Hi, Sushant,<br>&gt;<br>&gt;I think 'd a' is enough. Please refer to the document for vint().<br>&gt;<br>&gt;best wishes,<br>&gt;++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;YING, Ming<br>&gt;Shanghai Typhoon Institute of China Meteorological Administration<br>&gt;166 Puxi Road, Xuhui District, Shanghai, 200030<br>&gt;PR China<br>&gt;++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;2010/7/23 sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;:<br>&gt;&gt; Thank you for suggestion. Actually i want to calculate following integral.<br>&gt;&gt;<br>&gt;&gt;         1000<br>&gt;&gt;<br>&gt;&gt; Q = (1/g) ? (V q)dp.<br>&gt;&gt;<br>&gt;&gt;       400<br>&gt;&gt;<br>&gt;&gt; I used this kind of command in grads<br>&gt;&gt;<br>&gt;&gt; 'define ps=1000'<br>&gt;&gt; 'define a=vint(ps,vwnd*shum,400)'<br>&gt;&gt; 'd a/9.81'<br>&gt;&gt; with this command i am able to make a plot.<br>&gt;&gt;<br>&gt;&gt; whether it is correct or not ?<br>&gt;&gt;<br>&gt;&gt; Thanking you<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; Sushant<br>&gt;&gt;<br>&gt;&gt; On Fri, Jul 23, 2010 at 10:38 AM, MY &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=funnytyphoon@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">funnytyphoon@gmail.com</a>&gt; wrote:<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Hi,<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; please try<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; 'define ps=1000'<br>&gt;&gt;&gt; 'd vint(ps,vwnd*shum,400)'<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; best wishes,<br>&gt;&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt;&gt; YING, Ming<br>&gt;&gt;&gt; Shanghai Typhoon Institute of China Meteorological Administration<br>&gt;&gt;&gt; 166 Puxi Road, Xuhui District, Shanghai, 200030<br>&gt;&gt;&gt; PR China<br>&gt;&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; 2010/7/23 sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;:<br>&gt;&gt;&gt; &gt; Hi<br>&gt;&gt;&gt; &gt; I am trying to find vertical integral in GrADS using "vint" command. For<br>&gt;&gt;&gt; &gt; that purpose i opened vwnd and specific humidity file as first and<br>&gt;&gt;&gt; &gt; second<br>&gt;&gt;&gt; &gt; respectively. The q file command for both the variable are as<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; q file<br>&gt;&gt;&gt; &gt; File 1 :mean daily NMC reanalysis (2008)<br>&gt;&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt;&gt; &gt;   vwnd 8 -999 mean Daily V wind<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; q file 2<br>&gt;&gt;&gt; &gt; File 2 :mean daily NMC reanalysis (2008)<br>&gt;&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt;&gt; &gt;   shum 8 -999 mean Daily specific humidity<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; after this i used<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; d vint(1000,vwnd*shum.2,400)<br>&gt;&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt;&gt; &gt; Error from VINT: Incompatible grids.<br>&gt;&gt;&gt; &gt; Contouring: -0 to 2.4e+36 interval 3e+35<br>&gt;&gt;&gt; &gt; Segmentation fault<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; after showing this message, the grads window gets closed.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; Can anyone has solution for this perticular problem.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; Thanks in advance<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; --<br>&gt;&gt;&gt; &gt; Sushant<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; _______________________________________________<br>&gt;&gt;&gt; &gt; gradsusr mailing list<br>&gt;&gt;&gt; &gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt;&gt; &gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; _______________________________________________<br>&gt;&gt;&gt; gradsusr mailing list<br>&gt;&gt;&gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt;&gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; _______________________________________________<br>&gt;&gt; gradsusr mailing list<br>&gt;&gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;<br>&gt;&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;------------------------------<br>&gt;<br>&gt;Message: 3<br>&gt;Date: Fri, 23 Jul 2010 15:52:37 +0800<br>&gt;From: MY &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=funnytyphoon@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">funnytyphoon@gmail.com</a>&gt;<br>&gt;Subject: Re: [gradsusr] how to use vint function in GrADS<br>&gt;To: GrADS Users Forum &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a>&gt;<br>&gt;Message-ID:<br>&gt;&lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax@mail.gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">AANLkTi=_Cf-t2P0qUkjDL3Q0rd2ETxXQvmkuFuf-dkax@mail.gmail.com</a>&gt;<br>&gt;Content-Type: text/plain; charset=GB2312<br>&gt;<br>&gt;Hi, Sushant,<br>&gt;<br>&gt;I think 'd a' is enough. Please refer to the document for vint().<br>&gt;<br>&gt;best wishes,<br>&gt;++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;YING, Ming<br>&gt;Shanghai Typhoon Institute of China Meteorological Administration<br>&gt;166 Puxi Road, Xuhui District, Shanghai, 200030<br>&gt;PR China<br>&gt;++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;2010/7/23 sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;:<br>&gt;&gt; Thank you for suggestion. Actually i want to calculate following integral.<br>&gt;&gt;<br>&gt;&gt;         1000<br>&gt;&gt;<br>&gt;&gt; Q = (1/g) ? (V q)dp.<br>&gt;&gt;<br>&gt;&gt;       400<br>&gt;&gt;<br>&gt;&gt; I used this kind of command in grads<br>&gt;&gt;<br>&gt;&gt; 'define ps=1000'<br>&gt;&gt; 'define a=vint(ps,vwnd*shum,400)'<br>&gt;&gt; 'd a/9.81'<br>&gt;&gt; with this command i am able to make a plot.<br>&gt;&gt;<br>&gt;&gt; whether it is correct or not ?<br>&gt;&gt;<br>&gt;&gt; Thanking you<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; Sushant<br>&gt;&gt;<br>&gt;&gt; On Fri, Jul 23, 2010 at 10:38 AM, MY &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=funnytyphoon@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">funnytyphoon@gmail.com</a>&gt; wrote:<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; Hi,<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; please try<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; 'define ps=1000'<br>&gt;&gt;&gt; 'd vint(ps,vwnd*shum,400)'<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; best wishes,<br>&gt;&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt;&gt; YING, Ming<br>&gt;&gt;&gt; Shanghai Typhoon Institute of China Meteorological Administration<br>&gt;&gt;&gt; 166 Puxi Road, Xuhui District, Shanghai, 200030<br>&gt;&gt;&gt; PR China<br>&gt;&gt;&gt; ++++++++++++++++++++++++++++++++++++++++++++++++++++<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; 2010/7/23 sushant puranik &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=sushantpuranik@gmail.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">sushantpuranik@gmail.com</a>&gt;:<br>&gt;&gt;&gt; &gt; Hi<br>&gt;&gt;&gt; &gt; I am trying to find vertical integral in GrADS using "vint" command. For<br>&gt;&gt;&gt; &gt; that purpose i opened vwnd and specific humidity file as first and<br>&gt;&gt;&gt; &gt; second<br>&gt;&gt;&gt; &gt; respectively. The q file command for both the variable are as<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; q file<br>&gt;&gt;&gt; &gt; File 1 :mean daily NMC reanalysis (2008)<br>&gt;&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08v.nc<br>&gt;&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt;&gt; &gt;   vwnd 8 -999 mean Daily V wind<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; q file 2<br>&gt;&gt;&gt; &gt; File 2 :mean daily NMC reanalysis (2008)<br>&gt;&gt;&gt; &gt;  Descriptor: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt;&gt; &gt;  Binary: /home/Sushant/Desktop/10may08q.nc<br>&gt;&gt;&gt; &gt;  Type = Gridded<br>&gt;&gt;&gt; &gt;  Xsize = 29 Ysize = 25 Zsize = 8 Tsize = 1<br>&gt;&gt;&gt; &gt;  Number of Variables = 1<br>&gt;&gt;&gt; &gt;   shum 8 -999 mean Daily specific humidity<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; after this i used<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; ga-&gt; d vint(1000,vwnd*shum.2,400)<br>&gt;&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt;&gt; &gt; gxstart=0, gxstop=29, gystart=0, gystop=25.<br>&gt;&gt;&gt; &gt; gzstart=0, gzstop=1, gtstart=0, gtstop=1.<br>&gt;&gt;&gt; &gt; Error from VINT: Incompatible grids.<br>&gt;&gt;&gt; &gt; Contouring: -0 to 2.4e+36 interval 3e+35<br>&gt;&gt;&gt; &gt; Segmentation fault<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; after showing this message, the grads window gets closed.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; Can anyone has solution for this perticular problem.<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; Thanks in advance<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; --<br>&gt;&gt;&gt; &gt; Sushant<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt; _______________________________________________<br>&gt;&gt;&gt; &gt; gradsusr mailing list<br>&gt;&gt;&gt; &gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt;&gt; &gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt; &gt;<br>&gt;&gt;&gt;<br>&gt;&gt;&gt; _______________________________________________<br>&gt;&gt;&gt; gradsusr mailing list<br>&gt;&gt;&gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt;&gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt;<br>&gt;&gt; _______________________________________________<br>&gt;&gt; gradsusr mailing list<br>&gt;&gt; <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;&gt; <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;&gt;<br>&gt;&gt;<br>&gt;<br>&gt;<br>&gt;<br>&gt;------------------------------<br>&gt;<br>&gt;Message: 4<br>&gt;Date: Fri, 23 Jul 2010 09:37:41 +0100<br>&gt;From: "Becker, Bernd" <bernd.becker@metoffice.gov.uk><br>&gt;Subject: Re: [gradsusr] Obtaining total values of a dimension on a<br>&gt;netCDFfile<br>&gt;To: "GrADS Users Forum" &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a>&gt;<br>&gt;Message-ID:<br>&gt;&lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=E209E2066317644FB35E5AA0E8C1FB7F079994@EXXMAIL02.desktop.frd.metoffice.com&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">E209E2066317644FB35E5AA0E8C1FB7F079994@EXXMAIL02.desktop.frd.metoffice.com</a>&gt;<br>&gt;<br>&gt;Content-Type: text/plain; charset="us-ascii"<br>&gt;<br>&gt;Ok. <br>&gt;you are getting there.<br>&gt; <br>&gt;try:<br>&gt;set t 1<br>&gt; <br>&gt;d sum(data,t=1,t=32)<br>&gt; <br>&gt;any numbers looking familiar?<br>&gt;Cheers,<br>&gt;Bernd.<br>&gt; <br>&gt;<br>&gt;________________________________<br>&gt;<br>&gt;From: <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr-bounces@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr-bounces@gradsusr.org</a><br>&gt;[<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=mailto:gradsusr-bounces@gradsusr.org" target="_new">mailto:gradsusr-bounces@gradsusr.org</a>] On Behalf Of Thiago Veloso<br>&gt;Sent: 22 July 2010 18:06<br>&gt;To: GrADS Users Forum<br>&gt;Subject: Re: [gradsusr] Obtaining total values of a dimension on a<br>&gt;netCDFfile<br>&gt;<br>&gt;<br>&gt; Dear Bernd, <br>&gt;<br>&gt; Thanks for your reply.<br>&gt;<br>&gt; After trying the procedure you suggested I received lots of<br>&gt;information on display. Each information set corresponding to a specific<br>&gt;day.<br>&gt;<br>&gt; For example, setting t 1 3 generates:<br>&gt;<br>&gt;ga-&gt; set t 1 3<br>&gt;Time values set: 2009:1:1:0 2009:1:3:0 <br>&gt;ga-&gt; d cmorph<br>&gt;Data Type = grid<br>&gt;Dimensions = 0 1<br>&gt;I Dimension = 0 to 4949 Linear -0.0363783 0.0727567<br>&gt;J Dimension = 1 to 1649 Linear -59.9636 0.0727714<br>&gt;Sizes = 4950 1649 8162550<br>&gt;Undef value = -9.99e+08<br>&gt;Undef count = 114789 Valid count = 8047761<br>&gt;Min, Max = 0 241<br>&gt;Cmin, cmax, cint = 0 240 20<br>&gt;Stats[sum,sumsqr,root(sumsqr),n]:   4.20946e+06 8.7492e+07 9353.71<br>&gt;8047761<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/n]:   0.523059 10.8716 3.29721<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.523059 10.8716 3.29721<br>&gt;Stats[(sigma,var)(n)]:   3.25546 10.598<br>&gt;Stats[(sigma,var)(n-1)]:  3.25546 10.598<br>&gt;Data Type = grid<br>&gt;Dimensions = 0 1<br>&gt;I Dimension = 0 to 4949 Linear -0.0363783 0.0727567<br>&gt;J Dimension = 1 to 1649 Linear -59.9636 0.0727714<br>&gt;Sizes = 4950 1649 8162550<br>&gt;Undef value = -9.99e+08<br>&gt;Undef count = 89326 Valid count = 8073224<br>&gt;Min, Max = 0 217<br>&gt;Cmin, cmax, cint = 0 200 20<br>&gt;Stats[sum,sumsqr,root(sumsqr),n]:   3.78218e+06 7.38863e+07 8595.71<br>&gt;8073224<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/n]:   0.468485 9.15202 3.02523<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.468485 9.15202 3.02523<br>&gt;Stats[(sigma,var)(n)]:   2.98874 8.93254<br>&gt;Stats[(sigma,var)(n-1)]:  2.98874 8.93254<br>&gt;Data Type = grid<br>&gt;Dimensions = 0 1<br>&gt;I Dimension = 0 to 4949 Linear -0.0363783 0.0727567<br>&gt;J Dimension = 1 to 1649 Linear -59.9636 0.0727714<br>&gt;Sizes = 4950 1649 8162550<br>&gt;Undef value = -9.99e+08<br>&gt;Undef count = 177477 Valid count = 7985073<br>&gt;Min, Max = 0 200<br>&gt;Cmin, cmax, cint = 0 200 20<br>&gt;Stats[sum,sumsqr,root(sumsqr),n]:   3.8412e+06 9.12816e+07 9554.14<br>&gt;7985073<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/n]:   0.481048 11.4315 3.38105<br>&gt;Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 0.481048 11.4315 3.38106<br>&gt;Stats[(sigma,var)(n)]:   3.34666 11.2001<br>&gt;Stats[(sigma,var)(n-1)]:  3.34666 11.2001<br>&gt;ga-&gt; <br>&gt;<br>&gt; Is this information enough?<br>&gt;<br>&gt; Best wishes,<br>&gt;<br>&gt; Thiago.<br>&gt;<br>&gt;<br>&gt;--- On Thu, 22/7/10, Becker, Bernd <bernd.becker@metoffice.gov.uk><br>&gt;wrote:<br>&gt;<br>&gt;<br>&gt;<br>&gt;From: Becker, Bernd <bernd.becker@metoffice.gov.uk><br>&gt;Subject: Re: [gradsusr] Obtaining total values of a dimension on<br>&gt;a netCDF file<br>&gt;To: "GrADS Users Forum" &lt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a>&gt;<br>&gt;Date: Thursday, 22 July, 2010, 12:20<br>&gt;<br>&gt;<br>&gt;sdfopen file.nc<br>&gt;set gxout stats<br>&gt; <br>&gt;set t 1 32<br>&gt;d var<br>&gt; <br>&gt;what does that produce?<br>&gt; <br>&gt;Bernd.<br>&gt;<br>&gt;________________________________<br>&gt;<br>&gt;From: <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr-bounces@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr-bounces@gradsusr.org</a><br>&gt;[<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=mailto:gradsusr-bounces@gradsusr.org" target="_new">mailto:gradsusr-bounces@gradsusr.org</a>] On Behalf Of Thiago Veloso<br>&gt;Sent: 22 July 2010 15:18<br>&gt;To: <a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;Subject: [gradsusr] Obtaining total values of a dimension on a<br>&gt;netCDF file<br>&gt;<br>&gt;<br>&gt; Dear GrADS colleagues, <br>&gt;<br>&gt; I am working with a netCDF file which contains precipitation values<br>&gt;for 31 days. What I am trying to do is to use NCO to obtain the sum of<br>&gt;the "time" dimension (thus, the 31 days), but with no success until the<br>&gt;present moment.<br>&gt;<br>&gt; Could anybody give some directions?<br>&gt;<br>&gt; Best regards,<br>&gt;<br>&gt; Thiago.<br>&gt;<br>&gt;<br>&gt;-----Inline Attachment Follows-----<br>&gt;<br>&gt;<br>&gt;_______________________________________________<br>&gt;gradsusr mailing list<br>&gt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;<a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;<br>&gt;<br>&gt;<br>&gt;-------------- next part --------------<br>&gt;An HTML attachment was scrubbed...<br>&gt;URL: <a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fpipermail%2Fgradsusr%2Fattachments%2F20100723%2Fed4915ee%2Fattachment%2Ehtml" target="_new">http://gradsusr.org/pipermail/gradsusr/attachments/20100723/ed4915ee/attachment.html</a> <br>&gt;<br>&gt;------------------------------<br>&gt;<br>&gt;_______________________________________________<br>&gt;gradsusr mailing list<br>&gt;<a href="/prism/writemail?&amp;mode=mail_to_individual&amp;email=gradsusr@gradsusr.org&amp;output=web&amp;els=91de4925ef95f4d64786fcbaeb67dd7d" target="_new">gradsusr@gradsusr.org</a><br>&gt;<a href="http://www.rediffmail.com/cgi-bin/red.cgi?red=http%3A%2F%2Fgradsusr%2Eorg%2Fmailman%2Flistinfo%2Fgradsusr" target="_new">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>&gt;<br>&gt;<br>&gt;End of gradsusr Digest, Vol 5, Issue 50<br>&gt;***************************************<br>&gt;</bernd.becker@metoffice.gov.uk></bernd.becker@metoffice.gov.uk></bernd.becker@metoffice.gov.uk><br><A HREF="http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?" target="_blank"><IMG SRC="http://sigads.rediff.com/RealMedia/ads/adstream_nx.ads/www.rediffmail.com/signatureline.htm@Middle"></A>