Re: Issues with reading in forecast ensemble data

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Thu Apr 03 2014 - 09:49:01 MDT

Also ..

1] For the record: There is no 43 time step limit in NCL.

Some minor code comments

[2]
hr_arr = (/"00","06","12",....,"372", "378", "384"/)
     can be replaced with
hr_arr = sprinti("%0.2i", ispan(0,384,6)) ; "00",..,"12",..,"96","102",.

[3] mem = (/"01","02","03",...,"18","19","20"/)
     can be replaced with
     mem = sprinti("%0.2i", ispan(0,20,1))

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

On 4/3/14, 8:42 AM, Wei Huang wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 3 09:49:12 2014

This archive was generated by hypermail 2.1.8 : Thu Apr 03 2014 - 13:36:27 MDT