[ferret_users] More Climatology

[Andrew Cottrill <cottrill@xxxxxxxxxx>]

Hi guys

 

I have had success today with most of the noleap calendar axes, a blip in the file must have stopped the file running yesterday.

However, I am still trying to complete seasonal analyses on files with 360_day (numbered and calendar)

and non-calendar axis…

 

A/  Calendar =360_day    : time range 15:50745

For A/ 

If I use:

Use “my_ data

 

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

!Shows the following:

!variable T:

!           axis seas_ax

!subset     : 4 points (time)

!calendar   :360_day

!01-Dec-1900      /     1:  330.0

!01-Mar-1901      /     2:  420.0

!01-Jun-1901      /     3:  510.0  

!01-Sep-1901      /     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 .

 

 

B/  Calendar : =360_day    : time range 16-jan-1860 00:00 to 16-dec-1999 00:00

 

Use “my data”

use climatological_axes
cancel data climatological_axes
define axis/units=days/t="01-jan-1900":"31-dec-1999":12/calendar=360_day seas_ax

let PR1_new=PR1[gt=seas_ax@asn]
let PR2_new=PR2[gt=seas_ax@asn]

use "../../../../data/IPCC/20c3m/precip/monthly/pcmdi.ipcc4.ukmo_hadcm3.20c3m.run1.monthly.pr_A1.nc"
use "../../../../data/IPCC/20c3m/precip/monthly/pcmdi.ipcc4.ukmo_hadcm3.20c3m.run2.monthly.pr_A1.nc"
let PR1_noleap=PR1_new[gt=`PR,return=taxis`]
let PR2_noleap=PR2_new[gt=`PR,return=taxis`]

let climate1=PR1_noleap[gt=seas_ax@mod]
let climate2=PR2_noleap[gt=seas_ax@mod]
let Sclimate=(climate1+climate2)/2
 

etc
 

Runs and complains about modulo axis!  Regridding to a non modulo axis………….
No luck here either with suggestions from b/4

 

 

C/  Non calendar axes with time range: 1:1812 :

Show data gives:

 

Name                title                   I           J           K          L         

LON_BNDS                               1:2        1:90      …         …

LAT_BNDS                                1:2        1:60      …         …

TIME                 time                  …         …         …         1:1812             

                        (invalid coordinate axis)

TIME_BNDS

                                                1:2        …         …         1:1812

PR                    Precipitation      1:90      1:60      …         1:1812

 

Time range : 1:1812

 

Hope you can help.

Cheers….

 

Andy Cottrill

Post Graduate Student

Department of Biological & Physical Sciences

Faculty of Sciences

University of Southern Queensland

Toowoomba         Qld   4350

 

Phone 0419 0912 673