Skip to content

access_ACCESS_OM2_atm_filtered

Jump to bottom
ars599 edited this page Oct 13, 2018 · 5 revisions

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

  2. break into nx * ny blocks

    cdo selindexbox,$xS,$xE,$yS,$yE $ifile ${ofile}.nc

  3. 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

Clone this wiki locally