Hi,
In stead of giving single scalar value this _expression_ gives me multiple areas for each latitude. 2ndly I doubts that integral for x and y counts grids along the boundaries twice. Any suggestion? Thanks, Dogar Date: Mon, 28 Oct 2013 09:03:01 -0700 From: j_karmacharya@xxxxxxxxx Subject: Re: [ferret_users] How to divide a variable by area used in calculation To: couragetoknowgc@xxxxxxxxxxx; ferret_users@xxxxxxxx Hi Dogar, If your vaiable has lat-long information, I think this should work: Create a variable with unitary value from the desired variable and use @din to compute area: e.g.: let mask=if heat_anm_ctr1[l=1] then 1 list mask[x=@din.y=@din] ! give area in square m Jagadish On Monday, October 28, 2013 1:02 PM, Star Physicist <couragetoknowgc@xxxxxxxxxxx> wrote: Dear Ferret users, I am calculating Ocean Heat Content (OHC) as follows and wanted to divide it by area used in calculation. How can I do that: set reg/x=-180:180/y=60S:60N !Defining rho x temp let heat=rho[d=1] *temp[d=1] !Integrating over the entire volume let OHC=heat_anm_ctr1[x=@din,y=@din,z=@din]*3990 Now I want to divide this OHC by corresponding area to get the unit as Joule/m^2. How to get the area used in above calculation. Thanks, Dogar Date: Mon, 28 Oct 2013 06:23:11 +0530 Subject: [ferret_users] netcdf point datasets in ferret From: ghanshamsangar@xxxxxxxxx To: ferret_users@xxxxxxxx Is it possible to open point datasets in netcdf format following CF conventions (see discrete sampling geometry in cf metedata conventions1.6 document available online). I could not find online documentation regarding same. I am usinfg ferret 6.84. The dataset is satellite derived wind vectors having relevant fields (ucomp, vcomp, wind speed, Windows direction, latitude, longitude, date/time of each observation).
Ghansham
|