Hi Satyesh,
Thanks for making me think about the grids! Simple solution:
Let regrid = X1X[gt=etw[d=1]]
Save regrid
Then use the regrid instead of the original salinity dataset and it’s working a treat!
Best,
Josh
Sent from Mail for Windows 10
From: owner-ferret_users@xxxxxxxx <owner-ferret_users@xxxxxxxx> on behalf of JoshAB1995@xxxxxxxxxxx <JoshAB1995@xxxxxxxxxxx>
Sent: Thursday, June 6, 2019 2:37:08 PM
To: Satyesh Ghetiya
Cc: ferret_users@xxxxxxxx
Subject: RE: [ferret_users] Inconsistent t-axis
Hi Satyesh,
Here’s the set data:
yes? use "./RCP45_SBT.nc"
yes? use "./RCP45_SBS.nc"
yes? use "./RCP45_SST.nc"
yes? use "./RCP45_SSS.nc"
yes? sh d
currently SET data sets:
1> ./RCP45_SBT.nc
name title I J K L
ETW sea water potential temparature 1:151 1:151 1:1 1:600
2> ./RCP45_SBS.nc
name title I J K L
X1X sea water salinity at the lowes 1:151 1:151 1:1 1:600
3> ./RCP45_SST.nc
name title I J K L
ETW sea water potential temparature 1:151 1:151 1:1 1:600
4> ./RCP45_SSS.nc (default)
name title I J K L
X1X sea water salinity at the surfa 1:151 1:151 1:1 1:600
Here are the grids:
yes? sh gr etw[d=1]
GRID GGI1
name axis # pts start end subset
LON LONGITUDE 151mi 5W(-5) 10E full
LAT LATITUDE 151 i 48N 63N full
Z Z 1 r 1 1 full
TIME TIME 600 r 31-DEC-2049 16:04 01-DEC-2099 08:35 full
** INTERRUPTED! **
yes? sh grid etw[d=2]
**ERROR: variable unknown or not in data set: ETW[D=2]
yes? sh grid etw[d=3]
GRID GGI1
name axis # pts start end subset
LON LONGITUDE 151mi 5W(-5) 10E full
LAT LATITUDE 151 i 48N 63N full
Z Z 1 r 1 1 full
TIME TIME 600 r 31-DEC-2049 16:04 01-DEC-2099 08:35 full
yes? sh grid X1X[d=2]
GRID GGV1
name axis # pts start end subset
LON LONGITUDE 151mi 5W(-5) 10E full
LAT LATITUDE 151 i 48N 63N full
Z Z 1 r 1 1 full
TIME1 TIME 600 i 01-JAN-2050 00:00 01-DEC-2099 16:30 full
yes? sh grid X1X[d=4]
GRID GGV1
name axis # pts start end subset
LON LONGITUDE 151mi 5W(-5) 10E full
LAT LATITUDE 151 i 48N 63N full
Z Z 1 r 1 1 full
TIME1 TIME 600 i 01-JAN-2050 00:00 01-DEC-2099 16:30 full
yes?
I realise that they’re on different grids, would a simple regrid fix the issue?
I’ll give it a go anyway.
Cheers,
Josh
From: Satyesh Ghetiya
Sent: 06 June 2019 14:26
To: JoshAB1995@xxxxxxxxxxx
Cc: ferret_users@xxxxxxxx
Subject: Re: [ferret_users] Inconsistent t-axis
Can you give results of following commands from your data ?
(1) sh gr etw[d=1]
(2) sh gr X1X[d=2]
(3)sh gr etw[d=3]
(4)sh gr X1X[d=4]
(5) sh da (to show us all the data and variables)
On Thu, Jun 6, 2019 at 6:46 PM JoshAB1995@xxxxxxxxxxx <JoshAB1995@xxxxxxxxxxx> wrote:
Hi everyone,
Running into a strange error when I try to save a file after performing the following script:
use "./RCP85_SST.nc"
use "./RCP85_SSS.nc"
let temp=etw[d=1]*1.00024
let salt=X1X[d=2]
let potemp=THETA_FO(salt,temp,z[g=temp],100)
let surfdens=rho_un(salt,potemp,100)
shade surfdens[l=@ave]
use "./RCP85_SBT.nc"
use "./RCP85_SBS.nc"
let temp2=etw[d=3]*1.00024
let salt2=X1X[d=4]
let potemp2=THETA_FO(salt2,temp2,z[g=temp2],100)
let benthdens=rho_un(salt2,potemp2,100)
shade benthdens[l=@ave]
let strat=surfdens – benthdens
When I use save strat I get an error message saying the following:
yes? save strat
**ERROR: inconsistent sizes of data regions: T axis
TEMP has 72 points (L=01:72)
_expression_ has 600 points (L=001:600)
yes?
Have I written an erroneous script or does anyone know why this might happen? Every dataset has 600 time slices so I’m not sure why it’s saying that ‘temp’ has only 72?
Any help is greatly appreciated, thanks in advance.
Josh
--
G Satyesh