| 1 | DEFINE SYMBOL MyDir "/ccc/scratch/cont003/gencmip6/p86caub/RUN_DIR/CM5A2-BILAN_NEWOASIS01" |
|---|
| 2 | |
|---|
| 3 | GO use_and_define da ($MyDir)/COLIQRUN_LMDZ_25.nc |
|---|
| 4 | GO use_and_define do ($MyDir)/O_Runoff_oceanx_25.nc |
|---|
| 5 | GO use_and_define dgo /ccc/work/cont003/dsm/p86mart/GRAF/DATA/orca2.3.nc |
|---|
| 6 | GO use_and_define drmp ($MyDir)/rmp_tlmd_to_torc_MOSAIC_rivflu.nc |
|---|
| 7 | |
|---|
| 8 | LET runoff_srf = O_Runoff[D=($do)]*e1t[d=($dgo)]*e2t[D=($dgo)]*maskper[D=($dgo)] |
|---|
| 9 | LET runoff_srf_full = O_Runoff[D=($do)]*e1t[d=($dgo)]*e2t[D=($dgo)] |
|---|
| 10 | |
|---|
| 11 | LET bilo = runoff_srf[I=@SUM,J=@SUM]/1000.0 |
|---|
| 12 | LET bilo_f = runoff_srf_full[I=@SUM,J=@SUM]/1000.0 |
|---|
| 13 | LET bila = coliqrun[D=($da),I=@SUM,J=@SUM] |
|---|
| 14 | |
|---|
| 15 | LET srf_atm = RESHAPE ( src_grid_area[D=($drmp)], coliqrun[D=($da)] ) |
|---|
| 16 | LET srf_oce = RESHAPE ( dst_grid_area[D=($drmp)], O_Runoff[D=($do)] ) * maskper[D=($dgo)] |
|---|
| 17 | |
|---|
| 18 | LIST/K=1:17 bila, bilo, bilo_f, bilo/bila, (bilo-bila)/bila, bilo_f/bila, (bilo_f-bila)/bila |
|---|
| 19 | |
|---|
| 20 | SET WINDOW 1 ; GO page_new -i 1 2 |
|---|
| 21 | SET VIEW 1 ; SHADE reshape ( src_grid_imask[D=($drmp)], coliqrun[D=($da)] ) |
|---|
| 22 | SET VIEW 2 ; SHADE reshape ( dst_grid_imask[D=($drmp)], O_Runoff[D=($do)] ) |
|---|
| 23 | |
|---|
| 24 | |
|---|
| 25 | SET WINDOW 2 ; GO page_new -il 1 1 |
|---|
| 26 | SET VIEW 1 |
|---|
| 27 | LET xx = IF O_Runoff GT 0 THEN O_Runoff |
|---|
| 28 | SHADE xx[D=($do),K=2] |
|---|
