Commit 795a52b9 authored by Carl Schreck's avatar Carl Schreck
Browse files

Switching climo period

parent ff9885f0
......@@ -29,10 +29,10 @@ mkdir -p $LOG_DIR
RUN_DIR=$BASE_DIR/cluster/
EMPTY=\"\"
for VAR_NAME in uwnd chi psi vwnd olr uShear shear pwat
#for VAR_NAME in olr
for VAR_NAME in uwnd chi psi vwnd olr vShear uShear shear pwat
#for VAR_NAME in vShear
do
if [[ $VAR_NAME = "olr" || $VAR_NAME = "pwat" || $VAR_NAME = "uShear" || $VAR_NAME = "shear" ]]; then
if [[ $VAR_NAME = "olr" || $VAR_NAME = "pwat" || $VAR_NAME = "uShear" || $VAR_NAME = "shear" || $VAR_NAME = "vShear" ]]; then
ALL_LEVELS=( $EMPTY )
else
ALL_LEVELS=( 850 200 )
......
......@@ -15,12 +15,13 @@ begin ; main
print_clock( "Here we go!" )
; These are some parameters that could be useful to have up top
varName = "vShear"
outvar = "vwnd850"
invar = str_sub_str(outvar, "wnd", "")
print_clock( "Reading" )
inPath = "~/data/cfsr/clim/" + varName + ".clim.1999.2010.nc"
inPath = "~/data/cfsr/clim/" + invar + ".clim.1991.2020.nc"
inFile = addfile( inPath, "r" )
inData = inFile->$varName$
inData = inFile->$invar$
print_clock( "Reading template" )
......@@ -42,11 +43,11 @@ begin ; main
outData(3::4,:,:) = (/ regridData /)
print_clock("Writing")
outPath = "~/data/monitor/v2/clim/cfs/analysis/" + varName \
+ ".cfsr.mean.clim.daily.1999.2010.nc"
outPath = "~/data/monitor/v2/clim/cfs/analysis/" + outvar \
+ ".cfsr.mean.clim.daily.1991.2020.nc"
system( "rm " + outPath )
outFile = addfile(outPath,"c")
outFile->$varName$ = outData
outFile->$outvar$ = outData
print_clock( "Thank you, come again." )
......
......@@ -15,30 +15,32 @@ begin ; main
print_clock( "Here we go!" )
; These are some parameters that could be useful to have up top
inVar = "v850"
outVar = "vwnd850"
outVar = "vShear"
inVar = str_sub_str(outVar, "wnd", "")
print_clock( "Reading" )
inPath = "~/data/cfsr/clim/" + inVar + ".clim.nc"
inPath = "~/data/cfsr/clim/" + inVar + ".clim.1991.2020.nc"
inFile = addfile( inPath, "r" )
inData = inFile->$inVar$
print_clock( "Reading template" )
tempPath = "~/data/monitor/v2/clim/cfs/analysis/" \
+ "wnd850.cfsr.mean.clim.daily.1999.2010.nc"
; tempPath = "~/data/monitor/v2/clim/cfs/analysis/" \
; + "wnd850.cfsr.mean.clim.daily.1999.2010.nc"
tempPath = "~/data/monitor/v2/clim/change_climo/" + outVar \
+ ".clim.1999.2010.nc"
tempFile = addfile( tempPath, "r" )
tempData = tempFile->VGRD_P0_L100_GLL0
tempData = tempFile->$outVar$
print_clock( "Regridding" )
regridData = area_conserve_remap_Wrap( inData&lon, inData&lat, \
inData, tempData&lon_0, tempData&lat_0(::-1), False )
inData, tempData&lon, tempData&lat, False )
regridData = regridData(:,::-1,:)
print_clock("Writing")
outPath = "~/data/monitor/v2/clim/cfs/analysis/" + outVar \
+ ".clim.nc"
outPath = "~/data/monitor/v2/clim/change_climo/" + outVar \
+ ".clim.1991.2020.nc"
system( "rm " + outPath )
outFile = addfile(outPath,"c")
outFile->$outVar$ = regridData
......
Supports Markdown
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