!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
use
grids_example.nc
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
def ax/x=0:151:1 xax1
def ax/y=0:301:1 yax1
!!!!!!!!!!!!!!
!Generating new first and
last cells
let lon_rhoa =
lon_rho[gx=xax1,gy=yax1]
let lon_rhob = if i lt 152 then lon_rhoa else
(lon_rhoa[i=150]+(lon_rhoa[i=150]-lon_rhoa[i=149]))
let lon_rhoc = if j lt 302 then lon_rhob else
(lon_rhob[j=300]+(lon_rhob[j=300]-lon_rhob[j=299]))
let lon_rhod = if i gt 1 then lon_rhoc else
(lon_rhoc[i=2]+(lon_rhoc[i=2]-lon_rhoc[i=3]))
let lon_rhoe = if j gt 1 then lon_rhod else
(lon_rhod[j=2]+(lon_rhod[j=2]-lon_rhod[j=3]))
let lat_rhoa = lat_rho[gx=xax1,gy=yax1]
let lat_rhob = if i lt 152 then lat_rhoa else
(lat_rhoa[i=150]+(lat_rhoa[i=150]-lat_rhoa[i=149]))
let lat_rhoc = if j lt 302 then lat_rhob else
(lat_rhob[j=300]+(lat_rhob[j=300]-lat_rhob[j=299]))
let lat_rhod = if i gt 1 then lat_rhoc else
(lat_rhoc[i=2]+(lat_rhoc[i=2]-lat_rhoc[i=3]))
let lat_rhoe = if j gt 1 then lat_rhod else
(lat_rhod[j=2]+(lat_rhod[j=2]-lat_rhod[j=3]))
!!!!!!!!!!!!
!Calculating limits between
grid cells
let lon_rho2 =
((lon_rhoe[i=@shf:0,j=@shf:0]+lon_rhoe[i=@shf:1,j=@shf:0]+lon_rhoe[i=@shf:0,j=@shf:1]+lon_rhoe[i=@shf:1,j=@shf:1])/4)
let lat_rho2 =
((lat_rhoe[i=@shf:0,j=@shf:0]+lat_rhoe[i=@shf:1,j=@shf:0]+lat_rhoe[i=@shf:0,j=@shf:1]+lat_rhoe[i=@shf:1,j=@shf:1])/4)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
set win/asp=0.5 2
let diff_rho = lon_rho+lat_rho
let diff_rhoe = lon_rhoe+lat_rhoe
def symb hlim -74.8:-74.55
def symb vlim 36.58:36.7
def symb lev (-38.2,-37.9,0.01)
can mode verify
set view left
shade/i=1:150/j=1:300/lev=($lev)/hlim=($hlim)/vlim=($vlim)/nolab
diff_rho,lon_rho,lat_rho
rep/i=1:150 (plot/ov/vs/dash/line/col=1/nolab
lon_rho,lat_rho)
rep/j=1:300 (plot/ov/vs/dash/line/col=1/nolab
lon_rho,lat_rho)
rep/i=148:150 (rep/j=1:3
(plot/size=0.2/col=1/ov/vs/nolab/symb=2
{`lon_rho`},{`lat_rho`}))
label/nouser `($ppl$xlen)/2` `($ppl$ylen)*1.1` 0 0 0.2
Original variables
set view right
shade/i=1:152/j=1:302/lev=($lev)/hlim=($hlim)/vlim=($vlim)/nolab
diff_rhoe,lon_rhoe,lat_rhoe
rep/i=1:151 (plot/ov/vs/dash/line/col=1/nolab
lon_rho2,lat_rho2)
rep/j=1:300 (plot/ov/vs/dash/line/col=1/nolab
lon_rho2,lat_rho2)
rep/i=148:150 (rep/j=1:3
(plot/size=0.2/col=1/ov/vs/nolab/symb=2
{`lon_rho2`},{`lat_rho2`}))
label/nouser `($ppl$xlen)/2` `($ppl$ylen)*1.1` 0 0 0.2
Edited variables
set mode verify