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] |
---|