RE: [ferret_users] Inconsistent t-axis

["JoshAB1995@xxxxxxxxxxx" <JoshAB1995@xxxxxxxxxxx>]

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