[ferret_users] More Seasonal Climatology

Hi guys

I am trying to plot seasonal climatology for rainfall from IPCC AR4 . I am having trouble running the models which don’t have a noleap calendar axis.

Some models have a time axis under PR (for rainfall) and details are as follows:

A/ TIME : Time range 1:2640; no calendar axis specified

If I use:

Use “my_data’ ! 4 runs

use climatological_axes
cancel data climatological_axes
define axis/units=days/t0=01-jan-0000/t=1:2640:30/modulo/calendar=360_day seas_ax={330,420,510,600}
list t[gt=seas_ax]
!Shows the following:

!variable: T

! axis seas_ax

!subset : 4 points (time)

!calendar :360_day

!01-Dec / 1: 330.0

!01-Mar / 2: 420.0

!01-Jun / 3: 510.0

!01-Sep / 4: 600.0

let climate1=PR1[gt=seas_ax@asn]
let climate2=PR2[gt=seas_ax@asn]
let climate3=PR3[gt=seas_ax@asn]
let climate4=PR4[gt=seas_ax@asn]
let Sclimate=(climate1+climate2+climate3+climate4)/4

fill/line/key/level=(0,2000,100)/title="Summer Rain"/l=1 Sclimate*3 !etc

When run in ferret comes up with :

**ERROR : invalid command: attempt to regrid on a compressed axis

CONTOUR/FILL/ line/ key/level=(0,2000,100)/title=”summer rain”/l=1 Sclimate*3

Command file, command group, or repeat execution aborted

Show grid PR gives:

GRID GLY5

name axis # pts start end

LON LONGITUDE 128mr 0E 3.75W

LAT LATITUDE 64 i 87.863S 87.863N

Normal Z

TIME T 2640 r 1 2640

B/ TIME 360_day calendar with time range 15:50745

If I use:

Use “my_data’ ! one run

use climatological_axes
cancel data climatological_axes
define axis/units=days/t0=1-jan-1901/t=14415:50745:30/calendar=360_day seas_ax = {330,420,510,600}
let Sclimate=PR1[gt=seas_ax@asn]
list t[gt=seas_ax]
!Shows the following:

!variable: T

! axis seas_ax

!subset : 4 points (time)

!calendar :360_day

!01-Dec-1901 / 1: 330.0

!01-Mar-1902 / 2: 420.0

!01-Jun-1902 / 3: 510.0

!01-Sep-1902 / 4: 600.0

fill/line/key/level=(0,2000,100)/title="Summer Rain"/l=1 Sclimate*3
fill/line/key/level=(0,2000,100)/title="Autumn Rain"/l=2 Sclimate*3
fill/line/key/level=(0,2000,100)/title="Winter Rain"/l=3 Sclimate*3
fill/line/key/level=(0,2000,100)/title="Spring Rain"/l=4 Sclimate*3

Runs the files ok but shows “no valid data” on the plots .

Show grid PR gives:

GRID GLP4

name axis # pts start end

LON LONGITUDE 96mr 0E 3.75W

LAT LATITUDE 72 i 90S 90N

Normal Z

TIME T (360_day) 1692 r 15 50745

C/. TIME 360_day calendar with time range 16-jan-1860 00:00 16-dec-2000 00:00

If I use:

Use “my_data’

use climatological_axes
cancel data climatological_axes
define axis/t/units=months/t0=01-jan-0000/t="16-jan-1860:12:00:00":"16-dec-2000:12:00:00"/edges/modulo/calendar=360_day seas_ax={11,14,17,20,23}
let Sclimate=PR1[gt=seas_ax@mod]
list t[gt=seas_ax]
!Shows the following:

!variable: T

! axis seas_ax

!subset : 4 points (time)

!calendar :360_day

!16-JAN / 1: 12.50

!16-APR / 2: 15.50

!16-JUL / 3: 18.50

!16-OCT / 4: 21.50

! Comes up with error:

**ERROR : invalid command: attempt to regrid on a compressed axis

CONTOUR/FILL/ line/ key/level=(0,2000,100)/title=”summer rain”/l=1 Sclimate*3

Command file, command group, or repeat execution aborted

Show grid PR gives:

GRID GPH4

name axis # pts start end

LON LONGITUDE 96mr 0E 3.75W

LAT 1 LATITUDE 48 i 87.159S 87.159N

Normal Z

TIME TIME 1692 r 16-Jan-1860 00:00 16-Dec-2000 00:00

Much thanks –hope you can help.

Andy Cottrill

Post Graduate Student

Department of Biological & Physical Sciences

Faculty of Sciences

University of Southern Queensland

Toowoomba Qld 4350