[Thread Prev][Thread Next][Index]

Re: plotting IPCC stfmmc properly



Hi Mark,
The two datasets seem to disagree as to what the direction of the dimensions are: the first has J=1:189, K=1:31 and the second has I=1:189, J=1:131.  Ferret won't change the direction of axes when you do a regridding operation. I think you should be able to use the SET DATA/ORDER= to change the way one of the datasets is read, and then use the  regridding using @ASN as you suggest.

Ansley

Mark.Collier@csiro.au wrote:

Hi

 

The IPCC variable stfmmc has definitions for j and k which ferret cannot understand.

 

I am trying to to a shade of stfmmc which has properly drawn depth and latitude axis.

 

I thought it might be possible to use @asn – just assign to a grid which works.

 

My data is:

 

yes? show data

     currently SET data sets:

    1> /cs/datastore/csdar/cec006/ocean/ipcc/data/stfmmc/stfmmc_O1b_1991-2000_20

C3M_run1.nc

 name     title                             I         J         K         L

 STFMMC_ATLANTIC

          V_ATLANTIC[Z=@IIN,X=@DIN]        ...       1:189     1:31      1:120

       (X=0E:0E(360), Z=0:5000)

 STFMMC_INDIAN

          V_INDIAN[Z=@IIN,X=@DIN]          ...       1:189     1:31      1:120

       (X=0E:0E(360), Z=0:5000)

 STFMMC_PACIFIC

          V_PACIFIC[Z=@IIN,X=@DIN]         ...       1:189     1:31      1:120

       (X=0E:0E(360), Z=0:5000)

 STFMMC_GLOBAL

          V_GLOBAL[Z=@IIN,X=@DIN]          ...       1:189     1:31      1:120

       (X=0E:0E(360), Z=0:5000)

 

    2> data/20C3M/O1b/run1/stfmmc_O1bc_1961-1990_20C3M_run1.nc  (default)

 name     title                             I         J         K         L

 STFMMC   Meridional Overturning Streamfu  1:189     1:31      1:4       1:12

 

The first file was formed in ferret, the second is concocted, so that output from ferret it modified to be submitted to AR4.

 

If I shade stfmmc_atlantic[l=@ave] it looks perfect.

 

If I shade stfmmc[k=1,l=@ave] it uses a simple X/Y axis (and other features like labelling is lost).

 

I thought I could assign the array from stfmmc onto stfmmc_atlantic,

 

define grid/like=stfmmc_atlantic[d=1] ngrid

 

and then

shade stfmmc[l=1,d=2,k=1,g=ngrid@asn] – it plots something, but it is LATITUDE on the Y axis versus the simple latitude range on the X axis.

 

Is there a simple way of getting the shaded plot right?

Thanks,

Mark.


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement