Re: making JJA without using month_to_season function

From: Ping Huang <ph0007_at_nyahnyahspammersnyahnyah>
Date: Wed Oct 05 2011 - 05:20:53 MDT

Hi, Siraj

I'm sorry the former code has an error. These codes should be following:

precip = inFile->PRECT
precip2 = onedtond(precip, (/10, 5, 96, 144/))
precip_JJA = dim_avg_n(precip2(:, 1:3, :, :), 1)

BTW, onedtond can transfer an nd data to another n dimensions data. It is very powerful.

2011-10-53:38 Siraj ul Islam д

> Hi all,
>
> This may be a very easy question. I have a ten year monthly data file and for each year there are only 5 months (May, June, July, August and Sep.) data. In total my 10 year data has (10 x 5) 50 months. Now I want to analyze JJA season each year. As all the built in functions (clmMonLLT, month_to_season) in NCL need 12 months data file (time dimension must be a multiple of 12) and therefore I cannot use them for my data. Is there any other way to average out JJA each year or to calculate climatology for each month?
>
> here is the ncdump out put for my file.
>
>
> netcdf cam4_dsp_may_2000_2009 {
> dimensions:
> lat = 96 ;
> lon = 144 ;
> time = UNLIMITED ; // (50 currently)
> variables:
> double lat(lat) ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> double lon(lon) ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> double time(time) ;
> time:long_name = "time" ;
> time:units = "days since 2000-05-01 00:00:00" ;
> time:calendar = "noleap" ;
> time:bounds = "time_bnds" ;
> float PRECT(time, lat, lon) ;
> PRECT:units = "m/s" ;
> PRECT:long_name = "Total (convective and large-scale) precipitation rate (liq + ice)" ;
> PRECT:cell_methods = "time: mean" ;
> float PS(time, lat, lon) ;
> PS:units = "Pa" ;
> PS:long_name = "Surface pressure" ;
> PS:cell_methods = "time: mean" ;
> float TREFHT(time, lat, lon) ;
> TREFHT:units = "K" ;
> TREFHT:long_name = "Reference height temperature" ;
> TREFHT:cell_methods = "time: mean" ;
> float U200(time, lat, lon) ;
> U200:units = "m/s" ;
> U200:long_name = "Zonal wind at 200 mbar pressure surface" ;
> U200:cell_methods = "time: mean" ;
> float U850(time, lat, lon) ;
> U850:units = "m/s" ;
> U850:long_name = "Zonal wind at 850 mbar pressure surface" ;
> U850:cell_methods = "time: mean" ;
> float V200(time, lat, lon) ;
> V200:units = "m/s" ;
> V200:long_name = "Meridional wind at 200 mbar pressure surface" ;
> V200:cell_methods = "time: mean" ;
> float V850(time, lat, lon) ;
> V850:units = "m/s" ;
> V850:long_name = "Meridional wind at 850 mbar pressure surface" ;
> V850:cell_methods = "time: mean" ;
> float Z500(time, lat, lon) ;
> Z500:units = "m" ;
> Z500:long_name = "Geopotential Z at 500 mbar pressure surface" ;
> Z500:cell_methods = "time: mean" ;
>
> // global attributes:
> :Conventions = "CF-1.0" ;
> :source = "CAM" ;
> :case = "fv1.9x2.5_2000_05" ;
> :title = ;
> :logname = "siraj" ;
> :host = "b158" ;
> :Version = "$Name$" ;
> :revision_Id = "$Id$" ;
> :initial_file = "/home/siraj/ccsm4_working_copy/inputdata/2000/May/cam_initial_day01_00Z.nc" ;
> :topography_file = "/home/siraj/ccsm4_working_copy/inputdata/atm/cam/topo/USGS-gtopo30_1.9x2.5_remap_c050602.nc" ;
> :history = "Wed Sep 21 15:55:24 2011: ncrcat fv1.9x2.5_2000_05.cam2.h0.2000-05.nc fv1.9x2.5_2000_05.cam2.h0.2000-06.nc fv1.9x2.5_2000_05.cam2.h0.2000-07.nc fv1.9x2.5_2000_05.cam2.h0.2000-08.nc fv1.9x2.5_2000_05.cam2.h0.2000-09.nc fv1.9x2.5_2001_05.cam2.h0.2001-05.nc fv1.9x2.5_2001_05.cam2.h0.2001-06.nc fv1.9x2.5_2001_05.cam2.h0.2001-07.nc fv1.9x2.5_2001_05.cam2.h0.2001-08.nc fv1.9x2.5_2001_05.cam2.h0.2001-09.nc fv1.9x2.5_2002_05.cam2.h0.2002-05.nc fv1.9x2.5_2002_05.cam2.h0.2002-06.nc fv1.9x2.5_2002_05.cam2.h0.2002-07.nc fv1.9x2.5_2002_05.cam2.h0.2002-08.nc fv1.9x2.5_2002_05.cam2.h0.2002-09.nc fv1.9x2.5_2003_05.cam2.h0.2003-05.nc fv1.9x2.5_2003_05.cam2.h0.2003-06.nc fv1.9x2.5_2003_05.cam2.h0.2003-07.nc fv1.9x2.5_2003_05.cam2.h0.2003-08.nc fv1.9x2.5_2003_05.cam2.h0.2003-09.nc fv1.9x2.5_2004_05.cam2.h0.2004-05.nc fv1.9x2.5_2004_05.cam2.h0.2004-06.nc fv1.9x2.5_2004_05.cam2.h0.2004-07.nc fv1.9x2.5_2004_05.cam2.h0.2004-08.nc fv1.9x2.5_2004_05.cam2.h0.2004-09.nc fv1.9x2.5_2005_05.cam2.h0.2005-05.nc fv1.9x2.5_2005_05.cam2.h0.2005-06.nc fv1.9x2.5_2005_05.cam2.h0.2005-07.nc fv1.9x2.5_2005_05.cam2.h0.2005-08.nc fv1.9x2.5_2005_05.cam2.h0.2005-09.nc fv1.9x2.5_2006_05.cam2.h0.2006-05.nc fv1.9x2.5_2006_05.cam2.h0.2006-06.nc fv1.9x2.5_2006_05.cam2.h0.2006-07.nc fv1.9x2.5_2006_05.cam2.h0.2006-08.nc fv1.9x2.5_2006_05.cam2.h0.2006-09.nc fv1.9x2.5_2007_05.cam2.h0.2007-05.nc fv1.9x2.5_2007_05.cam2.h0.2007-06.nc fv1.9x2.5_2007_05.cam2.h0.2007-07.nc fv1.9x2.5_2007_05.cam2.h0.2007-08.nc fv1.9x2.5_2007_05.cam2.h0.2007-09.nc fv1.9x2.5_2008_05.cam2.h0.2008-05.nc fv1.9x2.5_2008_05.cam2.h0.2008-06.nc fv1.9x2.5_2008_05.cam2.h0.2008-07.nc fv1.9x2.5_2008_05.cam2.h0.2008-08.nc fv1.9x2.5_2008_05.cam2.h0.2008-09.nc fv1.9x2.5_2009_05.cam2.h0.2009-05.nc fv1.9x2.5_2009_05.cam2.h0.2009-06.nc fv1.9x2.5_2009_05.cam2.h0.2009-07.nc fv1.9x2.5_2009_05.cam2.h0.2009-08.nc fv1.9x2.5_2009_05.cam2.h0.2009-09.nc -v TREFHT,PS,PRECT,V200,U200,V850,U850,Z500 cam4_dsp_may_2000_2009.nc" ;
> :nco_openmp_thread_number = 1 ;
> }
>
>
> --
> Siraj Ul Islam
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Oct 5 05:20:53 2011

This archive was generated by hypermail 2.1.8 : Sun Oct 09 2011 - 13:05:26 MDT