Hi folks,
I am trying to calculate the rmsd and bias using the following script but i am getting the following errors could you guys please help with the problem which i am facing.::::
yes? set mode metafie
yes? use
pottm2008.ncyes? use
pottmp.2008.ncyes? set memory/size=500
Cached data cleared from memory
yes? set reg/t=10-may-2008:31-may-2008/z=10
yes? let mtemp=temperature[d=1,l=@ave]
yes? let mpot=pottmp[d=2,l=@var]
yes? let sdtemp=temperature[d=1,l=@var]^0.5
yes? let sdpot=pottmp[d=2,l=@var]^0.5
yes? let dU=temperature[d=1]-pottmp[d=2,g=temperature[d=1]]
yes? let bias=dU[l=@ave]
yes? let rmsd=dU[l=@var]^0.5
yes? list mtemp,mpot,sdtemp,sdpot,bias,rmsd
**ERROR: inconsistent sizes of data regions: X axis
Variable 2 is I=001:360
_expression_ is I=0001:2001
yes? go multi_view 3, 1, .8, .1, 0, .2, .1, .1
yes? plo/nolabel/vlim=-100:140/color=black/thick=3 temperature[d=1]
**ERROR: dimensions improperly specified: specified data is not a line
- its a 3D region: "plo/nolabel/vlim=-100:140/color=black/thick=3 temperature[d=1]"
yes? plo/nolabel/color=black/thick=3 temperature[d=1]
**ERROR: dimensions improperly specified: specified data is not a line
- its a 3D region: "plo/nolabel/color=black/thick=3 temperature[d=1]"
yes? plo/ov/nolabel/color=red/thick=3 pottmp[d=2]
**ERROR: dimensions improperly specified: overlay is on a different axis
"plo/ov/nolabel/color=red/thick=3 pottmp[d=2]"
yes?
the grids of the data are as follows:::
yes? use
pottm2008.nc yes? sh g temperature
GRID GOZ1
name axis # pts start end
LONN2551_N551 LONGITUDE 2001mr 130W(-130) 30E(30)
LAT959_2051 LATITUDE 1093 i 39.976S 39.976N
LEV DEPTH (m) 10 i- 0 200
TIME TIME 32 r 08-MAY-2008 00:00 08-JUN-2008 00:00
normal E
normal F
yes? sh d
currently SET data sets:
1> ./
pottm2008.nc (default)
name title I J K L M N
TEMPERATURE
potential temperature (degc) 1:2001 1:1093 1:10 1:32 ... ...
yes? use
pottm2008.nc yes? sh g pottmp
**ERROR: variable unknown or not in data set: POTTMP
yes? use
pottmp.2008.nc yes? sh g pottmp
GRID GCN2
name axis # pts start end
LON LONGITUDE 360mr 0.5E 0.5W
LAT LATITUDE 418 r 74.5S 64.499N
LEVEL DEPTH (m) 40 i- 5 4478
TIME1 TIME 12 i 01-JAN-2008 00:00 01-DEC-2008 00:00
normal E
normal F
yes? sh d pottmp
**ERROR: unknown data set:
yes? sh d
currently SET data sets:
1> ./
pottm2008.nc name title I J K L M N
TEMPERATURE
potential temperature (degc) 1:2001 1:1093 1:10 1:32 ... ...
2> ./
pottmp.2008.nc (default)
name title I J K L M N
DATE ... ... ... 1:12 ... ...
TIMEPLOT ... ... ... 1:12 ... ...
POTTMP Potential temperature 1:360 1:418 1:40 1:12 ... ...
yes?