wiki:access/ACCESS-OM2/atm_filtered

Purpose: try to apply high pass filter to the atmospheric variables

  1. on raijin: concatenate

cdo cat /g/data1/ua8/JRA55-do/v1-3/${var}*.nc ${var}.1958-2018.07092018.nc

  1. break into nx * ny blocks
    cdo selindexbox,$xS,$xE,$yS,$yE $ifile ${ofile}.nc
  1. move outputs to ruby and apply filtering then break into each year

(make_run_ncl_bp_concat_u_10.sh call get_nc_example5.ncl)

;=====high pass=====

opt = 0
ihp = 1
sigma = 1.
nWgt = 2 * 100 * pNum + 1
fca = 1./(pNum*365*7) ; start freq 7 years 3 hourly data
fcb = 1./(pNum*365*5/12) ; last freq 5 months 3 hourly data
wgt = filwgts_lanczos (nWgt, ihp, fca, -999, sigma)

;=====band pass=====

  1. add on climatology data. (make_run_u10.sh)

create climatology extended file
cdo yhourmean ${varArray}.${prdArray}.nc ${varArray}.${prdArray}.hrclim.nc
cdo yhoursub ${varArray}.${prdArray}.nc ${varArray}.${prdArray}.hrclim.nc ${varArray}.${prdArray}.hranom.nc
cdo sub ${varArray}.${prdArray}.nc ${varArray}.${prdArray}.hranom.nc ${varArray}.${prdArray}.hrrepeat.nc

Loop throughout years cdo selyear,$i/$i ${varArray}.${prdArray}.hrrepeat.nc ${varArray}_clim/${varArray}.${i}.hrrepeat.nc
cdo add ${varArray}_clim/${varArray}.${i}.hrrepeat.nc ${varArray}/${varArray}.bandpass.${i}.nc ${varArray}/${varArray}.bandpass.${i}.09102018.nc

Change time calendar from standard to gregorian
ncatted -a calendar,time,o,c,gregorian ${varArray}/${varArray}.bandpass.${i}.09102018.nc

Last modified 2 years ago Last modified on Oct 13, 2018 11:57:58 PM