Re: Issues with reading in forecast ensemble data

Mike,

It probably worth a try to change code:

  do hLag = 0, dimsizes(hr_arr)-1
        do x = 0, dimsizes(mem)-1

    fileNames(hLag,x) = "/archive/real_time/gfs/gfs_ens_fcst/"+hr_arr(hLag)+"/ens."+mem(x)

        end do
  end do

   printVarSummary(fileNames)
  ; print(fileNames(0,:))

  print("Adding files together")

   v500in = new((/64, 20,73,144/),"float")

  do y = 0,19
  gfsFiles = addfiles( fileNames(:,y) + ".grib", "r" )

  print(gfsFiles)

  ListSetType( gfsFiles, "join" )

   print("Finished joining files")

  printVarSummary(gfsFiles)

  presLev = 500
   
  print("Constructing forecast variable")

   test = gfsFiles[:]->V_GRD_2_ISBL
   print(test(:,0,0,0))

   printVarSummary(test)

   printVarSummary(v500in(:,0,:,:))

   v500in(:,y,:,:) = gfsFiles[:]->V_GRD_2_ISBL(:,{500},::-1,:)

   printVarSummary(v500in)

end do

   printVarSummary(v500in)

to:

  v500in = new((/64, 20,73,144/),"float")
  do hLag = 0, dimsizes(hr_arr)-1
     do x = 0, dimsizes(mem)-1
        fileName = "/archive/real_time/gfs/gfs_ens_fcst/"+hr_arr(hLag)+"/ens."+mem(x)

        gfsFile = addfile(fileName + ".grib", "r" )

        v500in(hLag,x,:,:) = gfsFile->V_GRD_2_ISBL({500},::-1,:)

        delete(gfsFile)
     end do
  end do

  printVarSummary(v500in)

Wei

On Apr 2, 2014, at 2:54 PM, Michael Ventrice <michael.ventrice@weather.com> wrote:

> Hello,
>
> I am trying to read in GFS ensemble model data in real-time using NCL (Using version 6.0.0). I am coming across an odd issue where I call in all 64 time steps and ensemble members individually, join the files together by looping over each ensemble member, which are dumped into an variable that I created. Now I can see all file paths being correctly called out to forecast hour 384. However after I join files, I am only capable to read 43 time steps where it should read all 64 time steps in.
>
> Oddly enough, I was receiving the same issue when I swap over to the European ensembles, where I can only read in 43 time steps. Is there some limitation that I hitting without knowing?
>
> Thanks for whatever help you may offer. I have attached the snippet of code that I am using and output file from the terminal window.
>
> Best regards,
> Mike V.
>
>
> --
> Michael Ventrice, Ph.D. | Operational Scientist
> w: e: michael.ventrice@weather.com
>
> <script.txt><output.txt>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 3 08:43:03 2014