Commit 3201e46c authored by Carl Schreck's avatar Carl Schreck

Automated Nightly Commit - Sat Jan 26 00:00:20 EST 2019

parent 5c59a982
......@@ -21,23 +21,23 @@ begin
print_clock( "Here we go! " )
; These are some parameters that could be useful to have up top
varName = "sst"
varName = "skt"
newFile = True
doCalcAnom = True
; basePath = "~/data/cfsr/monthly/" + varName
; basePath = "~/data/ersst/v5/sst.mnmean"
basePath = "~/data/oisst/oisst.month"
inPath = basePath + ".total.nc"
climPath = basePath + ".clim.1981.2010.nc"
anomPath = basePath + ".anom.nc"
stdPath = basePath + ".std.nc"
; basePath = "~/data/ncep1/monthly/"
; inPath = basePath + varName + ".mon" + ".mean.nc"
; climPath = basePath + varName + ".mon" + ".clim.nc"
; anomPath = basePath + varName + ".mon" + ".anom.nc"
; stdPath = basePath + varName + ".mon" + ".std.nc"
; basePath = "~/data/oisst/oisst.month"
; inPath = basePath + ".total.nc"
; climPath = basePath + ".clim.1981.2010.nc"
; anomPath = basePath + ".anom.nc"
; stdPath = basePath + ".std.nc"
basePath = "~/data/ncep1/monthly/"
inPath = basePath + varName + ".mon" + ".mean.nc"
climPath = basePath + varName + ".mon" + ".clim.nc"
anomPath = basePath + varName + ".mon" + ".anom.nc"
stdPath = basePath + varName + ".mon" + ".std.nc"
climStart = 1981001
climEnd = 2010366
......
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; convert_iri.ncl
; Carl Schreck (cjschrec@ncsu.edu)
; January 2019
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Description: Make an IRI ncep1 monthly file look like an esrl one
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
load "$CJS_NCL_LIB/lib.cjs_graphics.ncl"
load "$CJS_NCL_LIB/lib.array.ncl"
begin ; main
print_clock( "Here we go!" )
; These are some parameters that could be useful to have up top
outVar = "skt"
inVar = "temp"
print_clock( "Reading" )
inPath = "~/data/ncep1/iri/" + outVar + ".nc"
inFile = addfile( inPath, "r" )
data = inFile->$inVar$
print_clock( "Fixing coords" )
data@_FillValue = default_fillvalue(typeof(data))
data!0 = "time"
data!1 = "lat"
data!2 = "lon"
newTime = cd_convert( data&time, "days since 1800-01-01" )
add_time_atts( newTime )
delete(data&time)
data&time = newTime
print_clock( "Writing" )
outPath = "~/data/ncep1/monthly/" + outVar + ".mon.mean.nc"
system( "rm -f " + outPath )
outFile = addfile( outPath, "c" )
outFile->$outVar$ = data
print_clock( "Closing" )
delete( outFile )
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