--- MDTF.py Starting POD SM_ET_coupling
Found files: 
/vftmp/Wenhao.Dong/pid6243/inputdata/model/ESM4_historical_D1_LONG/mon/ESM4_historical_D1_LONG.mrsos.mon.nc
/vftmp/Wenhao.Dong/pid6243/inputdata/model/ESM4_historical_D1_LONG/mon/ESM4_historical_D1_LONG.evspsbl.mon.nc
/vftmp/Wenhao.Dong/pid6243/inputdata/model/ESM4_historical_D1_LONG/mon/ESM4_historical_D1_LONG.pr.mon.nc
 Env vars: 
CASENAME: ESM4_historical_D1_LONG
CLDTOT_var: CLDTOT
CODE_ROOT: /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics
DATADIR: /vftmp/Wenhao.Dong/pid6243/inputdata/model/ESM4_historical_D1_LONG
FIRSTYR: 1980
FSNTOA_var: FSNTOA
ICEFRAC_var: ICEFRAC
LANDFRAC_var: LANDFRAC
LASTYR: 2014
MODEL_DATA_ROOT: /vftmp/Wenhao.Dong/pid6243/inputdata/model/
OBS_DATA: /home/Oar.Gfdl.Mdteam/DET/analysis/mdtf/obs_data/SM_ET_coupling
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/SM_ET_coupling
RGB: /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/src/rgb
WK_DIR: /vftmp/Wenhao.Dong/pid6243/wkdir/MDTF_ESM4_historical_D1_LONG_1980_2014/SM_ET_coupling
WORKING_DIR: /vftmp/Wenhao.Dong/pid6243/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/pid6243/wkdir/MDTF_ESM4_historical_D1_LONG_1980_2014
venv_root: ./envs/venv
verbose: 1
wap_var: wap
z250_var: zg
zg_var: zg
 --- MDTF.py calling POD SM_ET_coupling

Found conda at /home/mdteam/anaconda/bin/conda
Found program python.
Found program Rscript.
Found R package colorRamps.
Found R package maps.
Found R package fields.
Found R package akima.
Found R package ncdf4.
monthly soil moisture file found
computing SM-ET coupling
--------- Starting SM_ET coupling generate figures (using R)----------------------------
running
  '/home/mdteam/anaconda/envs/_MDTF-diagnostics-R/lib/R/bin/R --slave --no-restore --vanilla --file=/arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/SM_ET_coupling/SM_ET_coupling.R'

Loading required package: spam
Loading required package: dotCall64
Loading required package: grid
Spam version 2.3-0 (2019-09-13) is loaded.
Type 'help( Spam)' or 'demo( spam)' for a short introduction 
and overview of this package.
Help for individual functions is also obtained by adding the
suffix '.spam' to the function name, e.g. 'help( chol.spam)'.

Attaching package: 'spam'

The following objects are masked from 'package:base':

    backsolve, forwardsolve

See https://github.com/NCAR/Fields for
 an extensive vignette, other supplements and source code 
There were 50 or more warnings (use warnings() to see the first 50)
There were 50 or more warnings (use warnings() to see the first 50)
Error in ncvar_put(nc, varncsave, diff) : 
  ncvar_put: error: you asked to write 12851 values, but the passed data array only has 12780 entries!
Execution halted
R routine /arch2/wnd/MDTF-GFDL/NEW20191125/mdtf/MDTF-diagnostics/diagnostics/SM_ET_coupling/SM_ET_coupling.R 
 [1] "Taking lon/lat from model"
[1] "/vftmp/Wenhao.Dong/pid6243/inputdata/model/ESM4_historical_D1_LONG/mon/ESM4_historical_D1_LONG.mrsos.mon.nc"
[1] "Taking soil moisture from model, top-10cm"
[1] "Taking lon/lat from model"
[1] "Getting the number of years in the model file... should be 35 (i.e., 1980-2014)"
[1] "Defining model mask"
[1] "Taking ET from model"
[1] "Correlating SM and ET"
[1] "Plotting model correlation"
[1] "plotting model"
null device 
          1 
[1] "Regridding GLEAM to same as model"
[1] "Plotting GLEAM correlation"
[1] "plotting obs"
null device 
          1 
[1] "Plotting MODEL - GLEAM correlation"
null device 
          1 
[1] "Take precipitation from model"
[1] "Plotting model mean precipitation"
[1] "plotting model"
null device 
          1 
[1] "Regridding GLEAM Precip"
[1] "Plotting GLEAM Pr"
null device 
          1 
[1] "Plotting Pr diff, Model - GLEAM"
null device 
          1 
[1] "Correcting model coupling for precipitation biases"
[1] "loading the CMIP5, 2x2 data"
[1] "regridding cor(SM,ET) from the model to 2x2"
[1] "regridding JJA SM and ET"
[1] "regridding DJF SM and ET"
[1] "computing 2x2 correlation"
[1] "Regridding model P to 2x2"
[1] "Regridding GLEAM coupling to 2x2"
[1] "Regridding GLEAM P to 2x2"
[1] " CORRECTING for Pr DIFFERENCE"
[1] "veryfing CMIP5 models lists match"
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[31] TRUE TRUE TRUE TRUE TRUE TRUE TRUE
[1] "Plotting model coupling corrected for P, on 2x2"
[1] "plotting model"
null device 
          1 
[1] "Plotting GLEAM coupling on 2x2"
[1] "plotting model"
null device 
          1 
[1] "Plotting model - GLEAM coupling on 2x2"
[1] "plotting model"

--------- Finished SM_ET coupling generate figures----------------------------