[ncl-talk] Era-Interim daily precipitation

Ehsan Taghizadeh ehsantaghizadeh at yahoo.com
Thu May 31 08:42:16 MDT 2018


 Thank you so much for your nice replies.With you comment, it seems work correctly. However I've added some changes to the script to have an output file same as input file, I mean about dimensions, attributes, and so on.However I have some difficulties about it and I'll be thankful if I could have your help, again.For example aboutfout->tp24_1!1 = tp!1fout->tp24_1!2 = tp!2I encountered error. It seems they have conflict with following lines:fout->latitude = fin->latitudefout->longitude = fin->longitude
Beside that I want to have "tp24_1" with dimension "Geo2D" to plot.I attached new version of script, and I'll be thankful if I hear from you.
SincerleyEhsan    On Thursday, May 31, 2018, 6:44:48 PM GMT+4:30, Dennis Shea <shea at ucar.edu> wrote:  
 
 The golden rule of  data processing is look at your data

You have:

tp=f->tp     ; <=== type 'short' [ packed data]

However, your file has:

      short tp ( time, latitude, longitude )
         scale_factor : 5.505649013956001e-06
         add_offset :   0.1803980955369768
         _FillValue :   -32767
         missing_value :        -32767
         units :        m
         long_name :    Total precipitation

You *MUST* unpack the variable before computing the sums.

tp = short2flt(f->tp)

printVarSummary(tp)
printMinMax(tp,0)

===
https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml

 




On Thu, May 31, 2018 at 7:19 AM, Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com> wrote:

 Dears,
I've tried sum 12 hour precipitation to obtain 24 hour precipitation using attached script.However it seems doesn't work correctly!I had to sum tp(0,:,:)+tp(1,:,:), also tp(2,:,;)+tp(3,;,:), and so on.I uploaded "2_2_ERAI_Prec_2015_Jan_00_12_12.nc" in ftp.cgd.ucar.edu.I'll be thankful if I hear from you.
SincerelyEhsan
    On Wednesday, May 16, 2018, 8:23:32 AM GMT+4:30, Dennis Shea <shea at ucar.edu> wrote:  
 
 This is book keeping. Something like:

loop  each year
   loop each month
      loop each day
   
        sum all values occurring between 00:00 and 11:59 inclusive ==>PRC(time,lat,lon,0)
        sum all values occurring between 12:00 and 23:59 inclusive ==>PRC(time,lat,lon,1)

     end day
   end month
 end year


On Sun, May 13, 2018 at 9:07 AM, Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com> wrote:

Hi every body;Could I have any help to get daily precipitation from Era-Interim, please?I mean, a way to sum "accumulated precipitation from 00:00 to 12:00" with "accumulated precipitation from 12:00 to 24:00"?
SincerelyEhsan
______________________________ _________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/ mailman/listinfo/ncl-talk



  
______________________________ _________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/ mailman/listinfo/ncl-talk



  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180531/33efba61/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: scripte1.ncl
Type: application/octet-stream
Size: 1132 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180531/33efba61/attachment.obj>


More information about the ncl-talk mailing list