Re: 3 hrly MM5-CCSM monthly means

That is correct, Dennis. I still used these functions, but only at the preparation stage of CCSM3 data for MM5 initializations.

-----Original Message-----
From: Dennis Shea [mailto:shea@ucar.edu]
Sent: Wednesday, December 01, 2010 4:31 PM
To: Gurer, Kemal@ARB
Cc: Heather McIntyre; ncl-talk@ucar.edu
Subject: Re: 3 hrly MM5-CCSM monthly means

Also the Climate Data Operators (CDO)
   https://code.zmaw.de/projects/cdo

could possibly be used

%> cdo monmean mm5_ccsm.nc ofile

Some examples of using the CDO, NCO, wget, wgrib from within NCL
are at:
      http://www.ncl.ucar.edu/Applications/system.shtml

On 12/1/10 3:17 PM, Gurer, Kemal@ARB wrote:
> Hello Heather,
>
> I do a similar work that you are describing using the readv3.f file that comes with MM5 distribution, rather than using any ncl defined functions to calculate any averages. I attached a gzipped tar file to this e-mail to show how I calculate monthly averaged surface temperature. If you wish to calculate and plot the averages using that method, unzip and untar the attached file, modify the fortran file for the size of your MM5's array dimensions, compile and run it with your MM5 file as input. Then plot the average surface temperature using the attached ncl script. If you wish to calculate the averages of other variables, then add the appropriate variables to the attached fortran file, calculate the averages of those variables in a similar way, write the output to separate files, create another ncl script for your other variable(s), and plot them in the same way that I did for temperature. Let me know if you have any questions. Good luck.
>
> Kemal.
>
> -----Original Message-----
> From: ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu] On Behalf Of Dennis Shea
> Sent: Wednesday, December 01, 2010 1:15 PM
> To: Heather McIntyre
> Cc: ncl-talk@ucar.edu
> Subject: Re: 3 hrly MM5-CCSM monthly means
>
> Not sure what "MM5-CCSM" files look like.
>
> There is an undocumented function in contributed.ncl
>
> undef ("calculate_monthly_values")
> function calculate_monthly_values (x:numeric, arith:string,
> nDim[1]:integer, opt[1]:logical)
>
> ; calculate monthly values [avg, sum, min, max]
> ; x: numeric array of 5D or less [eg: time,lev,lat,lon]
> ; *must* have time coordinate recognized by ut_calendar
> ; if 5D [case,time,lev,lat,lon]
> ;
> ; arith: "avg" [also, "ave"], "sum","min","max" others may be added later
> ;
> ; nDim : scalar integer that specifies the 'time' dimension [generally, 0]
> ;
> ; opt : option ... not used here
> ;
> ; Sample usage: x(time,lat,lon) where time are n-hrly or daily values.
> ; xMonthAvg = calculate_monthly_values(x, "avg", 0, False)
> ; xMonthSum = calculate_monthly_values(x, "sum", 0, False)
> ; xMonthMin = calculate_monthly_values(x, "min", 0, False)
> ; xMonthMax = calculate_monthly_values(x, "max", 0, False)
> ;
>
>
> On 12/1/10 10:46 AM, Heather McIntyre wrote:
>> Hello,
>> I am trying to calculate monthly means from 3hrly MM5-CCSM output data.
>> I've used the ut_calendar function but can't find a function that
>> converts daily data into monthly means.
>> Any suggestions?
>> cheers,
>>
>>
>>
>> _______________________________________________
>> 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
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Dec 1 17:36:05 2010