[RegCNET] ICBC data 2005
Sara Rauscher
srausche at ictp.it
Wed Mar 1 13:57:06 CET 2006
Hi Shekhar,
NNRP1 is available up to 2006. If you ftp to the CDC ftp site
(ftp.cdc.noaa.gov - use anonymous login and your email as a password and
then go to /Datasets/ncep.reanalysis/) you will see that data is
available through 2006.
At the same ftp site NNRP2 data is available through 2004 (same ftp
site, /Datasets/ncep.reanalysis2/).
cheers,
sara
M.Sudhanshu Shekhar wrote:
> Hi Bi / all,
>
> I want to do some experiments with 2005 July initial/BC data.
> NNRP1/NNRP2 is present upto 2004/2002. How to get that.
>
> Shekhar
>
>
>
>
> #################################
> *Dr. M. S. Shekhar*
> Centre for Atmospheric Sciences
> Indian Institute of Technology Delhi
> Hauz Khas, New Delhi-110 016
> India
> Phone: (011) 26596052
> email: shekar at cas.iitd.ernet.in <mailto:shekar at cas.iitd.ernet.in>
> hotmail: gmraja75 at hotmail.com <mailto:gmraja75 at hotmail.com>
> #################################
>
> ------------------------------------------------------------------------
> From: /Sara Rauscher <srausche at ictp.it>/
> To: /regcnet at lists.ictp.it, ?????? <dolgikh at meteo.kz>/
> Subject: /Re: [RegCNET] Re: Clarification about 'postproc.in'/
> Date: /Tue, 28 Feb 2006 18:27:51 +0100/
> >Hi,
> >
> >Regarding reading netcdf files produced by postproc in Grads --
> >("unfortunately, the Netcdf files, which you get from postproc,
> are not
> >readable in Grads") -- I have never had problems reading netCDF files
> >created with postproc with Grads. Svetlana's problem with the "no
> >discernable X coordinate" may be specific to her files -- Svetlana,
> >would you send me the netCDF file and/or additional information?
> Maybe I
> >can help.
> >
> >Grads should read these netCDF files in, but the problem is that they
> >don't contain the correct lat/lon information. You can create a data
> >descriptor file that will supply the correct lat/lon information for
> >your netcdf file to Grads.
> >
> >For example, I created the following file using postproc:
> ATM200208AVG.nc
> >
> >If I use the ncdump utility (ncdump -h ATM200208AVG.nc) I can
> view the
> >file information:
> >
> >netcdf ATM200208AVG {
> >dimensions:
> > lon = 49 ;
> > lat = 32 ;
> > level = 18 ;
> > time = UNLIMITED ; // (1 currently)
> >variables:
> > float lon(lon) ;
> > lon:long_name = "Longitude" ;
> > lon:units = "degrees_east" ;
> > lon:actual_range = -4.888423f, 32.60124f ;
> > float lat(lat) ;
> > lat:long_name = "Latitude" ;
> > lat:units = "degrees_north" ;
> > lat:actual_range = 37.03973f, 54.29553f ;
> > float level(level) ;
> > level:long_name = "Height_Index" ;
> > level:units = "level" ;
> > level:actual_range = 1050.f, 0.f ;
> > double time(time) ;
> > time:long_name = "Time" ;
> > time:units = "hours since 1900-1-1 00:00:0.0" ;
> > time:actual_range = 899544., 899544. ;
> > float U(time, level, lat, lon) ;
> > U:long_name = "Zonal Wind" ;
> > U:units = "m/s" ;
> > U:missing_value = -1.e+30f ;
> > float V(time, level, lat, lon) ;
> > V:long_name = "Meridional Wind" ;
> > V:units = "m/s" ;
> > V:missing_value = -1.e+30f ;
> > float TK(time, level, lat, lon) ;
> > TK:long_name = "Temperature" ;
> > TK:units = "K" ;
> > TK:missing_value = -1.e+30f ;
> > float QD(time, level, lat, lon) ;
> > QD:long_name = "Mixing Ratio" ;
> > QD:units = "kg/kg" ;
> > QD:missing_value = -1.e+30f ;
> > float QC(time, level, lat, lon) ;
> > QC:long_name = "Cloud Mixing Ratio" ;
> > QC:units = "kg/kg" ;
> > QC:missing_value = -1.e+30f ;
> > float MSE(time, level, lat, lon) ;
> > MSE:long_name = "Moist Static Energy" ;
> > MSE:units = "m2/s2" ;
> > MSE:missing_value = -1.e+30f ;
> > float RH(time, level, lat, lon) ;
> > RH:long_name = "Relative Humidity" ;
> > RH:units = "fraction" ;
> > RH:missing_value = -1.e+30f ;
> > float HGT(time, level, lat, lon) ;
> > HGT:long_name = "Geopotential Height" ;
> > HGT:units = "m" ;
> > HGT:missing_value = -1.e+30f ;
> > float PS(time, lat, lon) ;
> > PS:long_name = "Surface Pressure" ;
> > PS:units = "hPa" ;
> > PS:missing_value = -1.e+30f ;
> > float RT(time, lat, lon) ;
> > RT:long_name = "Total Precip" ;
> > RT:units = "mm/day" ;
> > RT:missing_value = -1.e+30f ;
> > float TGRND(time, lat, lon) ;
> > TGRND:long_name = "Ground Temperature" ;
> > TGRND:units = "K" ;
> > TGRND:missing_value = -1.e+30f ;
> > float SMT(time, lat, lon) ;
> > SMT:long_name = "Total Soil Water" ;
> > SMT:units = "mm" ;
> > SMT:missing_value = -1.e+30f ;
> > float RB(time, lat, lon) ;
> > RB:long_name = "Base Flow" ;
> > RB:units = "mm/day" ;
> > RB:missing_value = -1.e+30f ;
> >
> >// global attributes:
> > :domxmin = -4.888423f ;
> > :domxmax = 32.60124f ;
> > :domymin = 37.03973f ;
> > :domymax = 54.29553f ;
> > :domzmin = 1050.f ;
> > :domzmax = 0.f ;
> >
> >With information about the model grid (from OUT_HEAD.CTL), I can take
> >this information and make a data descriptor file that I can give any
> >name I like (say, atmavg). Once I write the descriptor file, I
> can open
> >the file in grads using "open atmavg", as indicated below:
> >
> >DSET ^ATM200208AVG.nc
> >dtype netcdf
> >undef -1.e+30
> >pdef 49 32 lcc 45.39 13.48 24.50 16.00 30.00 60.00 13.48 60000.
> 60000.
> >xdef 175 linear -9.61 0.2703
> >ydef 78 linear 34.70 0.2703
> >zdef 18 levels 74.08 122.24 175.22 237.83 310.08 387.14 464.20
> 541.26 618.32 695.38 767.62 830.23 883.21 926.56 960.27 984.35
> 998.80 1008.43
> >1008.43
> >tdef 1 linear 0z01aug2002 30dy
> >vars 13
> >U=>u 18 t,z,y,x westerly wind
> >V=>v 18 t,z,y,x southerly wind
> >TK=>tk 18 t,z,y,x air temperature
> >QD=>qd 18 t,z,y,x air specific humidity
> >QC=>qc 18 t,z,y,x cloud water mixing ratio
> >MSE=>mse 18 t,z,y,x moist static energy
> >RH=>rh 18 t,z,y,x relative humidity
> >HGT=>hgt 18 t,z,y,x geopotential height
> >PSA=>psa 0 t,y,x surface pressure
> >RT=>rt 0 t,y,x total precipitation
> >TGRND=>tgrnd 0 t,y,x groud temperature in BATS
> >SMT=>smt 0 t,y,x total soil water in mm H2O
> >RB=>rb 0 t,y,x accumulated infiltration
> >endvars
> >
> >
> >Note that descriptor files (for netcdf!) work with Grads version 1.9
> >and above ONLY - you can see the documentation at
> >http://grads.iges.org/grads/gadoc/descriptorfile.html.
> >
> >I hope this is helpful, please let me know if you encounter
> problems or
> >it doesn't work.
> >
> >cheers,
> >sara
> >
> >
> >Moetasim wrote:
> >
> > >Hi,
> > >As far as I know...unfortunately, the Netcdf files, which you
> get from
> > >postproc, are not readable in Grads. You have two options:
> > >1) Use ferret/NCL instead of Grads
> > >2) Write CTL file for files which you get after postprocessing
> in Grads format.
> > >Writing CTL file should not be very tricky. YOu can use ncdump
> to check the
> > >contents/list of variables of Netcdf postproc files..The
> sequence of variables
> > >in Grads formatted files will be same. You should use any
> ATM???.CTL file as
> > >template for ATM postproc file and similarly RAD and SRF CTL
> files as template
> > >for writing CTL files for your postprocessed Grads formatted
> files. Remember
> > >that you just need to change variable info and time info in CTL
> files...Rest
> > >should almost remain same...
> > >Try writing one.. If you are successful then do post your CTL
> files at
> > >RegCNET.
> > >And if you are not successful then................. you can ask
> again..or may
> > >be someone else have already tried for such CTL files so he can
> post for you...
> > >I hope you will give a try to any of the options I have
> mentioned above.
> > >
> > >
> > >Quoting ?????? <dolgikh at meteo.kz>:
> > >
> > >
> > >
> > >>Dear Moet.
> > >>
> > >>Thanks a lot for clarification about cru.param.
> > >>
> > >>But we face one more problem and can not go further.
> > >>
> > >>We have a problem with opening NetCDF files (ATM199703AVG.nc,
> etc) in GRads.
> > >>
> > >>The error is:
> > >>ga-> sdfopen ATM199703AVG.nc
> > >>Scanning self-describing file: ATM199703AVG.nc
> > >>SDF file has no discernable X coordinate.
> > >>
> > >>But when we try to open PRE199703.nc everything is ok.
> > >>May be it is possible to change parameters in postproc.in to
> get coordinates
> > >>
> > >>for output .nc file?
> > >>Is there description of postproc.in file anywhere?
> > >>it is possible to get output files in GRADS format (.dat), but
> we can't
> > >>recieve .ctl file with .dat file.
> > >>
> > >>I would be very grateful if you agree to help us again.
> > >>Thank you in advance
> > >>Best regards
> > >>Svetlana
> > >>
> > >>
> > >>
> > >>
> > >
> > >_______________________________________________
> > >RegCNET mailing list
> > >RegCNET at lists.ictp.it
> > >https://lists.ictp.it/mailman/listinfo/regcnet
> > >
> > >
> >
> >--
> >
> >Sara A. Rauscher
> >Physics of Weather and Climate Section
> >The Abdus Salam International Centre for Theoretical Physics
> >Strada Costiera 11
> >34014 Trieste ITALY
> >ph: +39 040 2240 225 (please note change) fax: +39 040 2240 449
> >email: srausche at ictp.it
> >
> >
> >_______________________________________________
> >RegCNET mailing list
> >RegCNET at lists.ictp.it
> >https://lists.ictp.it/mailman/listinfo/regcnet
>
>------------------------------------------------------------------------
>
>_______________________________________________
>RegCNET mailing list
>RegCNET at lists.ictp.it
>https://lists.ictp.it/mailman/listinfo/regcnet
>
>
--
Sara A. Rauscher
Physics of Weather and Climate Section
The Abdus Salam International Centre for Theoretical Physics
Strada Costiera 11
34014 Trieste ITALY
ph: +39 040 2240 225 (please note change) fax: +39 040 2240 449
email: srausche at ictp.it
More information about the RegCNET
mailing list