Commit 84be2433 authored by Carl Schreck's avatar Carl Schreck

Automated Nightly Commit - Wed Jun 13 00:00:21 EDT 2018

parent 005c7375
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; subtract_global.ncl
; Carl Schreck (cjschrec@ncsu.edu)
; March 2017
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Description: Subtract the global mean for each time step from each point
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
load "$CJS_NCL_LIB/lib.cjs_graphics.ncl"
load "$CJS_NCL_LIB/lib.array.ncl"
load "$CJS_NCL_LIB/print_clock.ncl"
begin ; main
print_clock( "Here we go!" )
; These are some parameters that could be useful to have up top
print_clock( "Reading" )
inPath = "~/data/cfsr/7dayAnom/t850.7dayAnom.nc"
inFile = addfile( inPath, "r" )
timeUnits = inFile->time@units
minTime = cd_inv_calendar( 1979, 01, 01, 00, 0, 0, timeUnits, 0 )
maxTime = cd_inv_calendar( 1989, 12, 31, 00, 0, 0, timeUnits, 0 )
data = inFile->t850({minTime:maxTime},:,:)
print_clock( "Writing" )
outPath = "~/data/cfsr/7dayAnom/t850.7dayAnom.1979.1989.nc"
system( "rm " + outPath )
outFile = addfile( outPath, "c" )
outFile->t850 = data
print_clock( "Thank you, come again." )
end; main
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment