Read 4 times daily observation data

Dear all; I would like to confirm the following matter. I want to cut the 4 times daily observation data (http://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBSearch.pl?&Dataset=NCEP/DOE+AMIP-II+Reanalysis+%28Reanalysis-2%29&Dataset=NCEP/DOE+AMIP-II+Reanalysis+%28Reanalysis-2%29+&Variable=U-wind). I want to take only the data between 80E-160E and 10N-10S. The purpose in only to reduce the file size for next processing. Thus, I want confirm whether the following code is correct or not in term of NCL code: ------------------ ; Read and cut the 4 times daily data ; take only data between 10N-10S and 80E-160E load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" begin   print( "Here we go! " + systemfunc( "date" ) )  ; mis = -999  ; mis@_FillValue = -999   ; open the file     fileName = "uwnd.2001.nc"   inFile = addfile( fileName, "r" )   obsPerDay = 4       latN = 10   latS = -latN     lonW = 80   lonE = 160   ; read the data   print( "Reading the data..." )   dataShort = inFile->uwnd  (:,:,{latS:latN},{lonW:lonE})     data = short2flt( dataShort )     ; open the file to write to     fileName = "uwndnew.2001.nc"   system( "rm " + fileName )   setfileoption( "nc", "Format", "LargeFile" )   outFile = addfile( fileName, "c" )   print( "Writing data... " + systemfunc( "date" ) )   outFile->uwnd = data     print( "Thank you, come again." + systemfunc( "date" ) ) end ------------------- Thanks for best help.. JukY

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Feb 26 22:43:53 2012