[Thread Prev][Thread Next][Index]

RE: [ferret_users] seasonal anomaly timeseries for 30yrs



Hello Julliet,

this sounds more like your time axis is not a calendar ferret recognizes and hence not know how to regrid on the new axis. 

Cheers,

Joerg

 


From: owner-ferret_users@xxxxxxxx [owner-ferret_users@xxxxxxxx] on behalf of afwande juliet [afwandej965@xxxxxxxxx]
Sent: 01 December 2015 01:34
To: Ulrich Jacques Diasso
Cc: ferret_users
Subject: Re: [ferret_users] seasonal anomaly timeseries for 30yrs

hi thanks
I am wondering if space was the problem coz I get the regridding error. I am I supposed to use seasonal regrid or monthly. My data is already seasonal average
yes? use precip_MAM_mean.nc
yes? show data
     currently SET data sets:
    1> ./precip_MAM_mean.nc  (default)
 name     title                             I         J         K         L         M         N
 PRECIP_MAM
          MAM Seasonal Mean precip         1:144     1:72      ...       1:30      ...       ...
       (X=0E:0E(360), Y=90S:90N)
 
yes? use climatological_axes
 *** NOTE: regarding /opt/Ferret/V693/go/climatological_axes.cdf ...
 *** NOTE: Climatological axes SEASONAL_REG, MONTH_REG, MONTH_IRREG, MONTH_GREGORIAN, MONTH_NOLEAP, MONTH_360_DAY, MONTH_ALL_LEAP  defined
yes? cancel data climatological_axes
yes? define grid/like=precip_MAM mygrid
yes? let climatology = precip_MAM[gt=seasonal_reg@mod]
yes? let anomaly = precip_MAM - climatology[gt=mygrid]
yes? plot anomaly[x=22:55@ave,y=-12:22@ave]
 **ERROR: regridding: only @ASN regridding between calendar and
          non-calendar axes: CLIMATOLOGY


On Mon, Nov 30, 2015 at 9:21 AM, Ulrich Jacques Diasso <udiasso@xxxxxxxxx> wrote:
Hi Juliet,
Since anomaly is the variable you want to plot, There should be space between plot and anomaly
like this:
plot anomaly[x=...,y=...]

Cheers


On Mon, Nov 30, 2015 at 4:33 PM, afwande juliet <afwandej965@xxxxxxxxx> wrote:
Dear users,
I have GPCP.nc data for 44years. I have used seasonal mask to create data precip_MAM_mean.nc that has only 30yrs data for MAM season only.

Now I want to plot the timeseries anomaly for this data as below  but get the error. Can someone help point out difference . I have tried to follow this example http://www.pmel.noaa.gov/maillists/tmap/ferret_users/fu_99/msg00291.html
but not able.
[anyuola@login01 validation]$ ferret
  NOAA/PMEL TMAP
  FERRET v6.93 
  Linux 2.6.32-504.el6.x86_64 64-bit - 11/13/14
   1-Dec-15 00:14    

yes? use precip_MAM_mean.nc
yes? show data grids
yes? show data
     currently SET data sets:
    1> ./precip_MAM_mean.nc  (default)
 name     title                             I         J         K         L         M         N
 PRECIP_MAM
          MAM Seasonal Mean precip         1:144     1:72      ...       1:30      ...       ...
       (X=0E:0E(360), Y=90S:90N)
 
yes? show grids precip_mam
    GRID GPN1
 name       axis              # pts   start                end
 LON       LONGITUDE          144mr   1.25E                1.25W
 LAT       LATITUDE            72 r   88.75S               88.75N
 normal    Z
 TYEAR     T (years)           30 r   1980                 2009
 normal    E
 normal    F
yes? use climatological_axes
 *** NOTE: regarding /opt/Ferret/V693/go/climatological_axes.cdf ...
 *** NOTE: Climatological axes SEASONAL_REG, MONTH_REG, MONTH_IRREG, MONTH_GREGORIAN, MONTH_NOLEAP, MONTH_360_DAY, MONTH_ALL_LEAP  defined
yes? cancel data climatological_axes
yes? define grid/like=precip_MAM mygrid
yes? let climatology = precip_MAM[gt=seasonal_reg@mod]
yes? let anomaly = precip_MAM - climatology[gt=mygrid]
yes? plot/anomaly [x=22:55@ave,y=-12:22@ave]
 **ERROR: unknown command qualifier: anomaly
yes? plot anomaly [x=22:55@ave,y=-12:22@ave]
 **ERROR: command syntax: anomaly [x=22:55@ave,y=-12:22@ave]
          omaly [x=




--
Ulrich Jacques DIASSO
Climate Change assessment Expert
Montoring for Environment and Security in Africa (MESA)
African Centre of Meteorological Applications for Development (ACMAD)
85, Avenue des Ministeres
Niamey -  NIGER
Cell Phone: +227 91 63 80 91
                       +226 76 17 76 00



[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement