--- MDTF.py Starting POD MJO_teleconnection Found files: /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/CM4_historical_LONG.pr.day.nc /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/CM4_historical_LONG.rlut.day.nc /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/CM4_historical_LONG.ua.day.nc /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/CM4_historical_LONG.ua.day.nc /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/CM4_historical_LONG.zg.day.nc Env vars: CASENAME: CM4_historical_LONG CLDTOT_var: CLDTOT CODE_ROOT: /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics DATADIR: /vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG FIRSTYR: 1980 FSNTOA_var: FSNTOA ICEFRAC_var: ICEFRAC LANDFRAC_var: LANDFRAC LASTYR: 2014 MODEL_DATA_ROOT: /vftmp/Wenhao.Dong/pid14579/inputdata/model/ OBS_DATA: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection OBS_DATA_ROOT: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data OBS_DATA_SOURCE: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data OUTPUT_DIR: /archive/wnd/MDTF-GFDL/TEST/OUTPUT/ POD_HOME: /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection RGB: /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/src/rgb WK_DIR: /vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection WORKING_DIR: /vftmp/Wenhao.Dong/pid14579/wkdir conda_env_root: /home/mdteam/anaconda/envs conda_root: /home/mdteam/anaconda convert_flags: -crop 0x0+5+5 convert_output_fmt: png data_manager: GFDL_auto dry_run: False environment_manager: GFDL_conda evspsbl_var: evspsbl file_transfer_timeout: 3000 frepp: False hfls_var: hfls hfss_var: hfss ignore-component: False lat_coord: lat lat_var: lat lev_coord: plev19 lev_var: plev19 lon_coord: lon lon_var: lon make_variab_tar: False model: CMIP_GFDL mrsos_var: mrsos omega500_var: wap pr_conversion_factor: 1 pr_var: pr prc_conversion_factor: 1 prc_var: prc prls_conversion_factor: 1 prls_var: PRECL prw_var: prw ps_var: ps psl_var: PSL qa_var: hus qsat_int_var: qsat_int r_lib_root: ./envs/r_libs rlds_var: rlds rlus_var: rlus rlut_var: rlut rsds_var: rsds rsdt_var: rsdt rsus_var: rsus rsut_var: rsut save_nc: True save_ps: False ta_var: ta tas_var: ta tauu_var: TAUX tave_var: tave test_mode: False time_coord: time time_var: time ts_var: TS u200_var: ua u250_var: ua u850_var: ua u_var: ua ua_var: ua v200_var: va v250_var: va v850_var: va v_var: va va_var: va variab_dir: /vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014 venv_root: ./envs/venv verbose: 1 wap_var: wap z250_var: zg zg_var: zg --- MDTF.py calling POD MJO_teleconnection Found conda at /home/mdteam/anaconda/bin/conda Found program python. Found program ncl. Found contributed.ncl in /home/mdteam/anaconda/envs/_MDTF-diagnostics-NCL/lib/ncarg/nclscripts. Found gsn_code.ncl in /home/mdteam/anaconda/envs/_MDTF-diagnostics-NCL/lib/ncarg/nclscripts. Found gsn_csm.ncl in /home/mdteam/anaconda/envs/_MDTF-diagnostics-NCL/lib/ncarg/nclscripts. Found shea_util.ncl in /home/mdteam/anaconda/envs/_MDTF-diagnostics-NCL/lib/ncarg/nclscripts. Found diagnostics_cam.ncl in /home/mdteam/anaconda/envs/_MDTF-diagnostics-NCL/lib/ncarg/nclscripts. ======================================================================= Execution of MJO Teleconnection Diagnotics is started from here ----------------------------------------------------------------------- Following input data file are found CM4_historical_LONG.pr.day.nc CM4_historical_LONG.rlut.day.nc CM4_historical_LONG.ua.day.nc CM4_historical_LONG.ua.day.nc CM4_historical_LONG.zg.day.nc ----------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_RMM_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_RMM_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) --------------------------------------------------------------------------------- (0) Model: CM4_historical_LONG (0) Input data path(Model) :/vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/ (0) Input data path(obs) :/home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection/ (0) Output data path :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Input data file for OLR :CM4_historical_LONG.rlut.day.nc (0) Input data file for U850:CM4_historical_LONG.ua.day.nc (0) Input data file for U250:CM4_historical_LONG.ua.day.nc (0) --------------------------------------------------------------------------------- (0) Date range from year 1980 to 2014 (0) Starting date: Year:1980 Month:1 Day:1 (0) End date: Year:2014 Month:12 Day:31 (0) First doy: 1980001 last doy: 2014365 (0) Total days: 12775 (0) ------------------------------------------- (0) Reading data for 250mb wind (0) U250 data has been read (0) Reading data for 850mb wind (0) U850 data has been read (0) ------------------------------------------- (0) Now reading OLR data (0) ++++++++++++++++++++ Variable: latt2 Type: double Total Size: 720 bytes 90 values Number of Dimensions: 1 Dimensions and sizes: [lat | 90] Coordinates: lat: [ -89.. 89] Number Of Attributes: 4 long_name : latitude units : degrees_N axis : Y bounds : lat_bnds Variable: latt Type: double Total Size: 720 bytes 90 values Number of Dimensions: 1 Dimensions and sizes: [lat | 90] Coordinates: lat: [ -89.. 89] Number Of Attributes: 4 long_name : latitude units : degrees_N axis : Y bounds : lat_bnds (0) ++++++++++++++++++++ (0) OLR data has been read (0) ------------------------------------------- (0) remove interannual variability (0) Interannual variability has been removed (0) ------------------------------------------- (0) remove seasonal cycle (0) calculate raw climatology then smooth it (0) no missing values (0) no missing values (0) no missing values (0) seasonal cycle has been removed (0) ------------------------------------------- (0) averaging the data over 15S-15N Variable: u200 Type: float Total Size: 7358400 bytes 1839600 values Number of Dimensions: 2 Dimensions and sizes: [time | 12775] x [lon | 144] Coordinates: Number Of Attributes: 3 _FillValue : -999 long_name : Anomalies from Daily Climatology average_op_ncl : dim_avg_n over dimension(s): lat Variable: olr Type: float Total Size: 7358400 bytes 1839600 values Number of Dimensions: 2 Dimensions and sizes: [time | 12775] x [lon | 144] Coordinates: Number Of Attributes: 3 _FillValue : -999 long_name : Anomalies from Daily Climatology average_op_ncl : dim_avg_n over dimension(s): lat (0) Done with latitudinal averaging (0) ------------------------------------------- (0) interpolate data to be of size 144 to match observed RMM EOF data Variable: u200 Type: float Total Size: 7358400 bytes 1839600 values Number of Dimensions: 2 Dimensions and sizes: [time | 12775] x [lon | 144] Coordinates: lon: [1.25..358.75] Number Of Attributes: 3 _FillValue : -999 long_name : Anomalies from Daily Climatology average_op_ncl : dim_avg_n over dimension(s): lat Variable: olr Type: float Total Size: 7358400 bytes 1839600 values Number of Dimensions: 2 Dimensions and sizes: [time | 12775] x [lon | 144] Coordinates: lon: [1.25..358.75] Number Of Attributes: 3 _FillValue : -999 long_name : Anomalies from Daily Climatology average_op_ncl : dim_avg_n over dimension(s): lat (0) Interpolation completed (0) ------------------------------------------- (0) read in EOF data (0) compute RMMs (0) calculate MJO phase (0) calculate amplitude (0) Data is saved into txt format at:/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) In following file: (0) CM4_historical_LONG_RMMs.txt (0) Process Completed (0) ================================================== ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_geop_hgt_comp_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_geop_hgt_comp_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) -------------------------------------------------------------------------- (0) Sig: True (0) Model: CM4_historical_LONG (0) Input data path :/vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/ (0) Output data path :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Input data file for Z250 :CM4_historical_LONG.zg.day.nc (0) t-statistics data table is read from:/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/t_table_data.txt (0) -------------------------------------------------------------------------- (0) data table for t stat is read in (0) Total days: 12775 (0) First doy: 1980001 last doy: 2014365 (0) remove seasonal cycle (0) Total days: 3150 (0) First doy: 1980001 last doy: 2014365 (0) -------------------------------------------------------------------------- (0) Reading in RMM indices data (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 45.4 (1) 43.6 (2) 55.6 (3) 51.4 (4) 45.4 (5) 50.6 (6) 56.6 (7) 57 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 45 (1) 43.6 (2) 56.2 (3) 52 (4) 45.8 (5) 50.6 (6) 57.2 (7) 56.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44.6 (1) 44.2 (2) 56.8 (3) 52.6 (4) 46.2 (5) 51 (6) 57.2 (7) 56.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44 (1) 44.6 (2) 57.4 (3) 53.6 (4) 46.2 (5) 51.6 (6) 57 (7) 56.6 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44 (1) 44.6 (2) 58.2 (3) 54.4 (4) 46.8 (5) 51.2 (6) 57.2 (7) 55.2 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44.2 (1) 44.4 (2) 58.8 (3) 55 (4) 47 (5) 50.8 (6) 57.6 (7) 53.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44.2 (1) 44.4 (2) 58.8 (3) 55.4 (4) 47.2 (5) 50.8 (6) 57.6 (7) 53.4 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44.4 (1) 44.2 (2) 59 (3) 56 (4) 47.4 (5) 51.2 (6) 57.8 (7) 52.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44.2 (1) 44.4 (2) 59.4 (3) 56 (4) 47.4 (5) 51.8 (6) 57.6 (7) 52.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44 (1) 44.6 (2) 59.8 (3) 55.6 (4) 47.4 (5) 51.8 (6) 57.4 (7) 52.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44 (1) 45 (2) 60.4 (3) 55.2 (4) 47.2 (5) 51.8 (6) 57.2 (7) 53 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.8 (1) 45.4 (2) 61 (3) 54.6 (4) 47 (5) 51.8 (6) 57.4 (7) 53.2 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.4 (1) 46 (2) 61.2 (3) 54.2 (4) 46.6 (5) 52 (6) 57.4 (7) 53.2 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.6 (1) 46.2 (2) 61.4 (3) 54.2 (4) 46.2 (5) 52.4 (6) 57.4 (7) 53 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.6 (1) 46.4 (2) 61.4 (3) 54.8 (4) 45.2 (5) 52.6 (6) 57.6 (7) 53 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.4 (1) 47 (2) 61.4 (3) 55.2 (4) 44.4 (5) 53 (6) 57.2 (7) 53.2 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.6 (1) 47.4 (2) 62 (3) 54.8 (4) 43.6 (5) 53.2 (6) 57.2 (7) 52.6 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.6 (1) 48 (2) 62.2 (3) 54.2 (4) 43 (5) 53.8 (6) 57.2 (7) 51.8 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 43.8 (1) 48.8 (2) 61.8 (3) 53.4 (4) 42.4 (5) 54.2 (6) 57.2 (7) 51.2 (0) significance testing (0) MJO average # of days: 5 Variable: n (subsection) Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) 44 (1) 49.6 (2) 61 (3) 52.8 (4) 42.4 (5) 54.4 (6) 57 (7) 50.6 (0) (0) Writing NetCDF files (0) Writing NetCDF file for composite of geopotential height for each phase (0) Writing NetCDF file for t-statistic values (0) Composite of geopotential height at 250-hPa are genreted now (0) Data is saved as netcdf format at:/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) In following files: (0) geop_compositesP_hgt250_CM4_historical_LONG.nc (0) tstatP_hgt250_CM4_historical_LONG.nc (0) -------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_prec_comp_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_prec_comp_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) -------------------------------------------------------------------------- (0) Model: CM4_historical_LONG (0) Input data path :/vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/ (0) Output data path :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Input data file for Precip :CM4_historical_LONG.pr.day.nc (0) -------------------------------------------------------------------------- (0) Precipitation (kg m-2 s-1) : min=0 max=0.00371879 Variable: pr Type: float Total Size: 662256000 bytes 165564000 values Number of Dimensions: 3 Dimensions and sizes: [time | 12775] x [lat | 90] x [lon | 144] Coordinates: time: [47450.5..60224.5] lat: [ -89.. 89] lon: [1.25..358.75] Number Of Attributes: 9 long_name : Precipitation units : kg m-2 s-1 missing_value : 1e+20 _FillValue : -999 cell_methods : time: mean cell_measures : area: area time_avg_info : average_T1,average_T2,average_DT standard_name : precipitation_flux interp_method : conserve_order1 (0) 0 (0) Precipitation (kg m-2 s-1) : min=0 max=0.00371879 (0) =============== (0) convert to mm/day (0) after convert Variable: pr Type: float Total Size: 662256000 bytes 165564000 values Number of Dimensions: 3 Dimensions and sizes: [time | 12775] x [lat | 90] x [lon | 144] Coordinates: time: [47450.5..60224.5] lat: [ -89.. 89] lon: [1.25..358.75] Number Of Attributes: 9 long_name : Precipitation units : kg m-2 s-1 missing_value : 1e+20 _FillValue : -999 cell_methods : time: mean cell_measures : area: area time_avg_info : average_T1,average_T2,average_DT standard_name : precipitation_flux interp_method : conserve_order1 (0) 0 (0) Precipitation (kg m-2 s-1) : min=0 max=321.303 (0) =============== (0) Total days: 12775 (0) First doy: 1980001 last doy: 2014365 (0) time indices 0 12774 19800101 20141231 (0) remove seasonal cycle (0) AAAAAAAA (0) 0 (0) BBBBBBBB (0) 0 (0) Total days: 3150 (0) First doy: 1980001 last doy: 2014365 (0) Reading RMM indices data file in (0) write output of daily composited precipitation data (0) Composites of precipitation anomalies are genreted now (0) Data is saved as netcdf format at:/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) In following file: (0) PR_composites_CM4_historical_LONG.nc (0) -------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_U250_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_U250_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) -------------------------------------------------------------------------- (0) ** Calculating 250-hPa Zonal Wind RMS Error ** (0) -------------------------------------------------------------------------- (0) starting year1980 (0) Last year2014 (0) Reading DJF mean ERA-I(1979-2016) zonal wind data ..... (0) Input data path :/vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/ (0) Output data path :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Input data file for Observed U250 :ERAi.1979_2016.DJF.mean.nc (0) Input data file forCM4_historical_LONGU250 :CM4_historical_LONG.ua.day.nc (0) -------------------------------------------------------------------------- (0) calculating DJF mean (0) Done (0) -------------------------------------------------------------------------- (0) Reading Model 250hPa zonal wind data (0) Isolating DJF days -- Excluding 29 feb, Jan, Feb of first year and Dec. of last year (0) Total DJF days: 3060 (0) latitude dimension: 90 (0) longitude dimension: 144 (0) Calculating DJF mean of Model data (0) Done (0) -------------------------------------------------------------------------- (0) Interpoling to reference resolution -- (0) Done (0) Calculating Diffrenerce between model and observations (0) Done (0) Calculating RMS error (0) U250 zonal wind RMSE over (15N-60N,110E-120W): CM4_historical_LONG 4.49516 (0) RMS error (0) U250 Zonal wind RMSE aroud Jet maxima at Model between (15N-60N,110E-120W) 1.81523 (0) RMSE calculation for model: CM4_historical_LONG is completed (0) -------------------------------------------------------------------------- (0) Reading zonal wind RMS error data for CMIP5 Models (0) Data Path: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection/ (0) RMS error data file for CMIP5 Models: U250_RMS_jetext and U250_RMS (0) Adding New model RMSE values to CMIP5 RMSE array (0) Ascii file for U250-hPa RMSE are updated with new model's U250-hPa RMSE values (0) Updated RMSE files are now saved at: /vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Updated RMSE ascii files are: U250_RMS_jetext and U250_RMS_updated (0) Process of U250-hPa RMSE caclulation is completed now (0) -------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_daig_Corr_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_daig_Corr_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) ------------------------------------------------------------------------------------ (0) ** Calculating teleconnection patttern correlation of all MJO phase for Z250-hPa ** (0) ------------------------------------------------------------------------------------ (0) Domain: Pacific (0) Dataset: CM4_historical_LONG (0) Reading Geopotential height data for model : CM4_historical_LONG (0) Input data path (rean.) :/home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection/ (0) Input data path (Model) :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) Input data file observation geop composite :geop_compositesP_hgt250_CM4_historical_LONG.nc (0) Input data file model geop composite :geop_compositesP_hgt250_CM4_historical_LONG.nc (0) Reading Geopotential height data for reanalysis (0) data reading is completed (0) ------------------------------------------------------------------------------------ (0) interpolate (0) interpolation done (0) ------------------------------------------------------------------------------------ (0) isolating Pacific domain (0) done (0) ------------------------------------------------------------------------------------ (0) calculate the cross-correlation between the patterns (0) done (0) writing correlation values into an ascii file (0) Output path: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection/ (0) output file : corrMJO_Pac_250mb_CM4_historical_LONG (0) ---------------------------------------------------------------------------------- (0) Reading phase pattern correlation data for CMIP5 Models (0) Adding New model pattern Corr. values to CMIP5 Pattern Corr. array (0) Ascii file for pattern correlation is updated with new model's Pat. Corr. values (0) saving upadated pat. corr. ascii file at :/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/ (0) saving upadated pat. corr. file name :ccr_Z250comp_CMIP5_updated (0) Process of Pattern correlation calculation is completed now (0) ---------------------------------------------------------------------------------- (0) ---------------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_EWR_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_EWR_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) ---------------------------------------------------------------------------------- (0) ** Calculating Ratio of the eastward and westward moving power of precipitation ** (0) ---------------------------------------------------------------------------------- (0) Dataset: CM4_historical_LONG (0) Reading daily precipitation data (0) Input data path :/vftmp/Wenhao.Dong/pid14579/inputdata/model/CM4_historical_LONG/day/ (0) Input data file observation geop composite :CM4_historical_LONG.pr.day.nc (0) Input data path (Observations EWR) :/home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/MJO_teleconnection/ (0) Input data file CMIP5 models EWR :EWratio_wf_CMIP5 (0) Done (0) ---------------------------------------------------------------------------------- (0) original file (0) Precipitation (kg m-2 s-1) : min=0 max=0.00371879 (0) 0 (0) ===================== (0) convert to mm/day (0) after conversion (0) Precipitation (kg m-2 s-1) : min=0 max=321.303 (0) 0 (0) removing seasonal cycle (0) AAAAAAAA (0) 0 (0) BBBBBBBB (0) 0 (0) Precipitation (kg m-2 s-1) : min=0 max=321.303 (0) calculating raw and smooth climatology .... Variable: P Type: float Total Size: 73584000 bytes 18396000 values Number of Dimensions: 3 Dimensions and sizes: [time_p | 12775] x [lat | 10] x [lon | 144] Coordinates: time_p: [47450.5..60224.5] lat: [ -9.. 9] lon: [1.25..358.75] Number Of Attributes: 4 calendar : noleap units : kg m-2 s-1 long_name : Anomalies: Precipitation _FillValue : -999 (0) (0) Anomalies: Precipitation (kg m-2 s-1) : min=-21.4312 max=317.912 (0) Done (0) ---------------------------------------------------------------------------------- (0) ================= (0) SUM_WEST=0.120577 (0) ================= (0) CM4_historical_LONG:(Eastward power/Westward power) and Eastward power: (0) E/W: 2.77554 east power: 0.334666 (0) ---------------------------------------------------------------------------------- (0) Reading E/W ratio data for CMIP5 Models (0) Adding New model E/W ratio values to CMIP5 E/W ratio array (0) Ascii file for E/W ratio is updated with new model's precip E/W ratio values (0) Now saved at: /vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/netCDF/EWratio_wf_CMIP5_updated (0) Process of E/W ratio calculation is completed now (0) ---------------------------------------------------------------------------------- (0) ---------------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_fig1_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_fig1_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) ----------------------------------------------------------------------------------- (0) Plotting phase composites for precipitation and geopotential height all phases (0) ----------------------------------------------------------------------------------- (0) precipitation and geopotential height composites are read in (0) Figures for precipitation and geopotential height composites are now saved (0) at:/vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection/model/PS/ (0) ----------------------------------------------------------------------------------- ('Calling ', '/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_fig2_MDTF.ncl') NCL routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/MJO_teleconnection/mjo_diag_fig2_MDTF.ncl Copyright (C) 1995-2018 - All Rights Reserved University Corporation for Atmospheric Research NCAR Command Language Version 6.5.0 The use of this software is governed by a License Agreement. See http://www.ncl.ucar.edu/ for more details. (0) -------------------------------------------------------------------- (0) ploting of MJO teleconnection diagnosis figures (0) -------------------------------------------------------------------- (0) read in pattern correlation data (0) read in E/W ratio data (0) read in zonal wind RMS error data (0) 250hPa zonal mean wind error plot (0) -------------------------------------------------------------------- Variable: avdataT Type: float Total Size: 64 bytes 16 values Number of Dimensions: 2 Dimensions and sizes: [2] x [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 Variable: avdataP Type: float Total Size: 32 bytes 8 values Number of Dimensions: 1 Dimensions and sizes: [8] Coordinates: Number Of Attributes: 1 _FillValue : 9.96921e+36 (0) ploting of figures are done ----------------------------------------------------------------------------- |----Execution of MJO Teleconnections diagnostics module is completed now----| ============================================================================= Check: /vftmp/Wenhao.Dong/pid14579/wkdir/MDTF_CM4_historical_LONG_1980_2014/MJO_teleconnection now you can open index.html in browser to see the results -----------------------------------------------------------------------------