[gradsusr] Computing Precipitation Yearly and Seasonal Means from ecmwf data

Peter Tuju peterenos at ymail.com
Thu Sep 17 09:50:07 EDT 2015


Dear grads users,I am trying out to find the yearly Precipitation Means from the 3 hourly daily ecmwf data spanning from 
1980-01-01 to 2015-06-31 with 12965 time steps but ignoring the last year (2015) which contains missing 
data for the Months from July to December. Below is my code I use;
reinit
sdfopen ecmwf_netcdf.nc
set display color white
set grid off
set gxout shaded
set mpdset hires
set t 1 12
define yr_ave=ave(cp, t+1, t=12965, 12)                            **Starting with the First Year (1980)
d yr_ave

Can you understand what I am trying to do?. What is wrong with my code? Please guide me.

Also I want to plot the Precipitation Seasonal Means from the 3 hourly daily ecmwf data spanning from 
1980-01-01 to 2015-06-31. I want to get the spatial precipitations for the seasons of MAM, JJA, SON, 
OND, DJB and JFM. But I dont know how I can go further. I need your help please>
 _____________
Peter  Enos
Dar es Salaam
T A N Z A N I A
----------------------
      From: "gradsusr-request at gradsusr.org" <gradsusr-request at gradsusr.org>
 To: gradsusr at gradsusr.org 
 Sent: Wednesday, September 16, 2015 8:53 PM
 Subject: gradsusr Digest, Vol 67, Issue 27
   
Send gradsusr mailing list submissions to
    gradsusr at gradsusr.org

To subscribe or unsubscribe via the World Wide Web, visit
    http://gradsusr.org/mailman/listinfo/gradsusr
or, via email, send a message with subject or body 'help' to
    gradsusr-request at gradsusr.org

You can reach the person managing the list at
    gradsusr-owner at gradsusr.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of gradsusr digest..."


Today's Topics:

  1. wind direction difference (?? ????)
  2. mask country by shpfile (Dibas Shrestha)


----------------------------------------------------------------------

Message: 1
Date: Tue, 15 Sep 2015 21:30:51 +0300
From: ?? ???? <shaike07 at gmail.com>
Subject: [gradsusr] wind direction difference
To: gradsusr at gradsusr.org
Message-ID:
    <CAPSLP6_N7TCoFntq+1f41vDZ+BiWWCes1-2__sv792pNzCMxkw at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Hello everyone,
I am trying to plot differences in wind field between two scenarios. I am
using Grads 2.0 a8.
I will appreciate help on how to plot the differences in wind direction? or
better yet how to plot difference in speed and direction simultaneously on
the same map?
thank you
shai

-- 
Shai Kaplan, Ph.D.
ASU-BGU Partnership project manager & Postdoctoral Researcher
Ben-Gurion University of the Negev, Israel
Email: skaplan3 at asu.edu / shaike07 at gmail.com
Phone:  058-5933184
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150915/bbae2e88/attachment-0001.html 

------------------------------

Message: 2
Date: Tue, 15 Sep 2015 06:19:28 +0000 (UTC)
From: Dibas Shrestha <st.dibas at yahoo.com>
Subject: [gradsusr] mask country by shpfile
To: "gradsusr at gradsusr.org" <gradsusr at gradsusr.org>
Message-ID:
    <257905170.102264.1442297968764.JavaMail.yahoo at mail.yahoo.com>
Content-Type: text/plain; charset="utf-8"

Dear grads users, 

I am using TRMM_3B42 daily v7 data for my study. I want to get area average rainfall value for Nepal using country shpfile. Further I want to display only Nepal by using mask-out function. While I am trying to use GDAL and GrADS, following error has been is occurs. I would appreciate your suggestions.? 
? 
Data Request Warning:? Request is completely outside file limits
? Entire grid contents are set to missing data 
? Grid limits of file:???? X = 1 1440? Y = 1 400? Z = 1 1? T = 1 1? E = 1 1 
? Grid limits of request:? X = 320 361? Y = 300 329? Z = 1 1? T = -313919 -313919? E = 1 1 
? Warning issued for variable = r
Define memory allocation size = 10080 bytes

Attached herewith is data file and .gs file is as follows. 
sdfopen nep-map.nc
define ca=band1
sdfopen 3B42_daily.1998.08.07.7.nc
set dfile 2 
set time jan1998
set mpdset hires
set gxout grfill
set lat 25 32
set lon 80 90
define carain=r*ca
d carain

*quit


Thank you in advance. 
Dibas Shrestha 
Nepal Academy of Science and Technology
Nepal
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150915/d5349400/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 3B42_daily.1998.08.07.7.nc
Type: application/x-netcdf
Size: 2321820 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20150915/d5349400/attachment.nc 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nep-map.nc
Type: application/x-netcdf
Size: 4623792 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20150915/d5349400/attachment-0001.nc 

------------------------------

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr


End of gradsusr Digest, Vol 67, Issue 27
****************************************


  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150917/99086ea5/attachment.html 


More information about the gradsusr mailing list