[ncl-talk] Problem with Time dimension after saving to netcdf file in NCL

[Dennis Shea]

[1] It is likely that NCL's "dimension reduction"  is the issue.
Generally, this is a "feature" on NCL. However, sometimes it requires some
extra coding. You must use an 'if' to test. Maybe something like:

   dimr = dimsizes(rain)
   ntim = dimr(0)
   nlat = dimr(1)
   mlon = dimr(2)
   if (ntim.eq.1) then
       rain_aave = new ((/1,nlat,mlon/),
typeof(rain_region),getVarFillVaalue(rain_region))
       rain_aave(0,:,:) =wgt_areaave_Wrap(rain_region,1.0,1.0,0)
   else
       .rain_aave=wgt_areaave_Wrap(rain_region,1.0,1.0,0)
  end if

[2]
undef("add_hyphen")
function add_hyphen(yyyymmdd[*]:integer)
local yyyy, mmdd, mm, dd
begin
   yyyy    = yyyymmdd/10000
   mmdd = yyyymmdd-(yyyymmdd*10000)
   mm     = mmdd/100
   dd       = mmdd-(mm*100)
   return(yyyy+"-"+mm+"-"+dd)
end



On Thu, Aug 29, 2019 at 4:19 AM Lyndz via ncl-talk <ncl-talk at ucar.edu>
wrote:

> Dear NCL-experts,
>
> I trying to run my script (attached in this email). There are two issues
> that I would like to address and ask for help.
>
> *[1] The script runs but the output time dimension is incorrect when only
> one-time step satisfies the conditions. For cases with multiple time steps
> satisfying the conditions, the output is ok. Notice the time is 0 from the
> ncdump below. Any ideas how to solve this issue?*
>
> Here's the ncdump of the output netcdf file when only one time-step
> satisfies the coditions:
> dimensions:
>
> time = UNLIMITED ; // (0 currently)
> latitude = 280 ;
> longitude = 360 ;
> variables:
> float rain(latitude, longitude) ;
> rain:time = 116. ;
> rain:level_description = "Earth surface" ;
> rain:missing_value = -99.9f ;
> rain:_FillValue = -99.9f ;
> rain:units = "mm/day" ;
> rain:long_name = "daily precipitation analysis interpolated onto 0.25deg
> grids [mm/day]" ;
> double latitude(latitude) ;
> latitude:axis = "Y" ;
> latitude:units = "degrees_north" ;
> latitude:long_name = "Latitude" ;
> latitude:standard_name = "latitude" ;
> double longitude(longitude) ;
> longitude:axis = "X" ;
> longitude:units = "degrees_east" ;
> longitude:long_name = "Longitude" ;
> longitude:standard_name = "longitude" ;
>
> // global attributes:
> :creation_date = "Thu Aug 29 18:28:59 JST 2019" ;
> :Conventions = "None" ;
> :source_file = "Daily Aphrodite Rainfall" ;
> :title = "rain_stretch equal rain_threshold=0" ;
> }
>
> *[2] One part of the script save the dates into a text file. Is there a
> way to save the output file as YYYY-MM-DD separated by a dash (-) instead
> of YYYYMMDD?*
>
> I'll appreciate any help about these.
>
> Sincerely,
>
> Lyndz
>
> _______________________________________________
> 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/20190829/83ea8f03/attachment.html>