DEFINE SYMBOL MyDir "/ccc/scratch/cont003/gencmip6/p86caub/RUN_DIR/CM5A2-BILAN_NEWOASIS01"

GO use_and_define da   ($MyDir)/COLIQRUN_LMDZ_25.nc
GO use_and_define do   ($MyDir)/O_Runoff_oceanx_25.nc
GO use_and_define dgo  /ccc/work/cont003/dsm/p86mart/GRAF/DATA/orca2.3.nc
GO use_and_define drmp ($MyDir)/rmp_tlmd_to_torc_MOSAIC_rivflu.nc

LET runoff_srf      = O_Runoff[D=($do)]*e1t[d=($dgo)]*e2t[D=($dgo)]*maskper[D=($dgo)]
LET runoff_srf_full = O_Runoff[D=($do)]*e1t[d=($dgo)]*e2t[D=($dgo)]

LET bilo   = runoff_srf[I=@SUM,J=@SUM]/1000.0
LET bilo_f = runoff_srf_full[I=@SUM,J=@SUM]/1000.0
LET bila   = coliqrun[D=($da),I=@SUM,J=@SUM]

LET srf_atm = RESHAPE ( src_grid_area[D=($drmp)], coliqrun[D=($da)] )
LET srf_oce = RESHAPE ( dst_grid_area[D=($drmp)], O_Runoff[D=($do)] ) * maskper[D=($dgo)]

LIST/K=1:17 bila, bilo, bilo_f, bilo/bila, (bilo-bila)/bila, bilo_f/bila, (bilo_f-bila)/bila

SET WINDOW 1 ; GO page_new -i 1 2
SET VIEW 1 ; SHADE reshape ( src_grid_imask[D=($drmp)], coliqrun[D=($da)] )
SET VIEW 2 ; SHADE reshape ( dst_grid_imask[D=($drmp)], O_Runoff[D=($do)] )


SET WINDOW 2 ; GO page_new -il 1 1
SET VIEW 1
LET xx = IF O_Runoff GT 0 THEN O_Runoff
SHADE xx[D=($do),K=2]