Re: Average initial CLM dataset

From: Lei <dream916_at_nyahnyahspammersnyahnyah>
Date: Fri, 19 Jun 2009 15:38:46 -0500

Hi Dennis,
   Thanks for your suggestions. One more question before I try to test it,
if the initial datasets contain many variables, how can I make sure the new
dataset file contains the same contents as the old dataset files. Thanks,

Lei

A sample of variables in the initial dataset is provided as follows,
*****************
variables:
      integer mcdate ( ncl_scalar )
         long_name : current date as 8 digit integer (YYYYMMDD)

      integer mcsec ( ncl_scalar )
         long_name : current seconds of current date
         units : s

      double grid1d_lon ( gridcell )
         long_name : gridcell longitude
         units : degrees_east

      double grid1d_lat ( gridcell )
         long_name : gridcell latitude
         units : degrees_north

      integer grid1d_ixy ( gridcell )
         long_name : 2d longitude index of corresponding gridcell

      integer grid1d_jxy ( gridcell )
         long_name : 2d latitude index of corresponding gridcell

      double land1d_lon ( landunit )
         long_name : landunit longitude
         units : degrees_east

      double land1d_lat ( landunit )
         long_name : landunit latitude
         units : degrees_north

      integer land1d_ixy ( landunit )
         long_name : 2d longitude index of corresponding landunit

      integer land1d_jxy ( landunit )
         long_name : 2d latitude index of corresponding landunit

      integer land1d_gi ( landunit )
         long_name : 1d grid index of corresponding landunit

      double land1d_wtxy ( landunit )
         long_name : landunit weight relative to corresponding gridcell

      integer land1d_ityplun ( landunit )
         long_name : landunit type (vegetated,urban,lake,wetland or
glacier)

      double cols1d_lon ( column )
         long_name : column longitude
         units : degrees_east

      double cols1d_lat ( column )
         long_name : column latitude
         units : degrees_north

      integer cols1d_ixy ( column )
         long_name : 2d longitude index of corresponding column

      integer cols1d_jxy ( column )
         long_name : 2d latitude index of corresponding column

      integer cols1d_gi ( column )
         long_name : 1d grid index of corresponding column

      integer cols1d_li ( column )
         long_name : 1d landunit index of corresponding column

      double cols1d_wtxy ( column )
         long_name : column weight relative to corresponding gridcell

      double cols1d_wtlnd ( column )
         long_name : column weight relative to corresponding landunit

      integer cols1d_ityplun ( column )
         long_name : column landunit type (vegetated,urban,lake,wetland
or glacier)

      double pfts1d_lon ( pft )
         long_name : pft longitude
         units : degrees_east

      double pfts1d_lat ( pft )
         long_name : pft latitude
         units : degrees_north

      integer pfts1d_ixy ( pft )
         long_name : 2d longitude index of corresponding pft

      integer pfts1d_jxy ( pft )
         long_name : 2d latitude index of corresponding pft

      integer pfts1d_gi ( pft )
         long_name : 1d grid index of corresponding pft

      integer pfts1d_li ( pft )
         long_name : 1d landunit index of corresponding pft

      integer pfts1d_ci ( pft )
         long_name : 1d column index of corresponding pft

      double pfts1d_wtxy ( pft )
         long_name : pft weight relative to corresponding gridcell

      double pfts1d_wtlnd ( pft )
         long_name : pft weight relative to corresponding landunit

      double pfts1d_wtcol ( pft )
         long_name : pft weight relative to corresponding column

      integer pfts1d_itypveg ( pft )
         long_name : pft vegetation type

      integer pfts1d_ityplun ( pft )
         long_name : pft landunit type (vegetated,urban,lake,wetland or
glacier)

      double T_VEG ( pft )
         long_name : vegetation temperature
         units : K

      double T_GRND ( column )
         long_name : ground temperature
         units : K

      double EFLX_LWRAD_OUT ( pft )
         long_name : emitted infrared (longwave) radiation
         units : watt/m^2

      double H2OCAN ( pft )
         long_name : canopy water
         units : kg/m2

      double H2OSNO ( column )
         long_name : snow water
         units : kg/m2

On Thu, Jun 18, 2009 at 6:14 PM, Dennis Shea <shea_at_ucar.edu> wrote:

> [1]
> Use addfiles with "join"
> See Example 2
> http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml
>
> [2]
> If you have: x(column, levtot) on each file and you have (say) 4 files
>
> you will have
> [ncl_join | 4] x [column | ?] x [levtot | ?]
>
> [3]
> avex = dim_avg_Wrap(x(column|:,levtot|:,ncl_join|:)) ; v5.1.0
>
> or, better, with 5.1.1
>
> avex = dim_avg_n_Wrap(x, 0)
>
> In eother case,
>
> avex(column,levtot)
>
> Possibly, you may also uecatse the netCDF operators [ncecat]
>
> Good luck
>
>
>
> Lei wrote:
>
>> Hello,
>> I am working on simulation of the impact of soil moisture on
>> precipitation. I have used CAM-CLM3.0 to create a number of initial datasets
>> ( such as camrun1.clm2.i.***.nc). I am trying to average variable values in
>> the different initial datasets to a new initial dataset and have no idea on
>> how to achieve this. The difficulties are: these initial datasets contain
>> many variables and variable dimensions are like (column, levtot), not like
>> (levtot, lat, lon). I would appreciate any suggestions or comments. Thank
>> you!
>>
>> Lei
>> Texas A&M University
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> 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 Fri Jun 19 2009 - 14:38:46 MDT

This archive was generated by hypermail 2.2.0 : Mon Jun 22 2009 - 15:36:10 MDT