load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" ;interim_hr_to_day.ncl ;Mike Ventrice ;MVentrice@albany.edu ;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ; Takes a large era-interim.nc file (6-hrly) and coverts to daily averaged steps ;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ;>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> begin ; read in the data print( "Reading the data..." ) inFile = addfile( "/ct13/ventrice/data/interim/raw/g.200.anom.nc", "r" ) timeUnits = "hours since 1800-01-01 00:00:00" inStart = ut_inv_calendar( 1989, 01, 01, 00, 0, 0, timeUnits, 0 ) inEnd = ut_inv_calendar( 2010, 06, 30, 00, 0, 0, timeUnits, 0 ) times = inFile->time lat = inFile->lat lon = inFile->lon var = inFile->g(0:1,:,:) ;ut = ut_calendar(times,0) print( "Writing the anomalies..." ) fileName = "/ct13/ventrice/data/interim/raw/g.200.daily.anom.nc" system( "rm " + fileName ) setfileoption( "nc", "Format", "LargeFile" ) outFile = addfile( fileName, "c" ) time_24 = times(::4) dimNames = (/ "time", "lat", "lon" /) dimSizes = (/ dimsizes(time_24), dimsizes(lat), dimsizes(lon) /) dimUnlim = (/ False, False, False /) filedimdef( outFile, dimNames, dimSizes, dimUnlim ) filevardef( outFile, "time", typeof(time_24), "time" ) filevardef( outFile, "lat", typeof(lat), "lat" ) filevardef( outFile, "lon", typeof(lon), "lon" ) filevardef( outFile, "g", typeof(var), (/ "time", "lat", "lon" /) ) filevarattdef( outFile, "time", time_24 ) filevarattdef( outFile, "lat", lat ) filevarattdef( outFile, "lon", lon ) outFile->time = (/ time_24 /) outFile->lat = (/ lat /) outFile->lon = (/ lon /) print( "Averaging..." ) out = new((/dimsizes(time_24),dimsizes(lat),dimsizes(lon)/), "float") do n = 0, dimsizes(times)-1, 4 print( n+" of "+dimsizes((times)-1) ) no = n/4 data1 = inFile->g(n:n+3,:,:) data = dim_avg_n_Wrap(data1, 0) out(no, :,: ) = data end do delete(data) outFile->g = (/out/) print("Thank you, come again..") end ;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<