[ecco-support] Fwd: Re: A question about temperature tendency
浩刘
hliu at whoi.edu
Fri Oct 30 21:06:51 EDT 2015
Hi, ecco
Diana helped with me on question 2,how about question 1) and 3) I
have asked?Can you help me with thoes?
Thank you so much.
BTW,1) is If the studied pressure below 100m in the ocean,should I still
use shortwave penetration (derived by oceQsw) to close the temperature
tendency equation?
3)If I want to calculate temperature tendency due to
advection, should I do it in the following step?
(3.1)CellVol = RAC(ix,iy) * DRF(iz) * hFacC(ix,iy,iz);
*(3.2)Adv_tend *= - [ {ADVr_TH(ix,iy,iz) – ADVr_TH(ix,iy,iz+1)}/CellVol +
{ADVx_TH(ix+1,iy,iz) – ADVx_TH(ix,iy,iz)}/CellVol +
{ADVy_TH(ix,iy+1,iz) – ADVy_TH(ix,iy,iz)}/CellVol ]
Best Regards,Hao
-------- Forwarded Message --------
Subject: Re: A question about temperature tendency
Date: Fri, 30 Oct 2015 18:46:44 -0400
From: Diana Lees Spiegel <diana at ocean.mit.edu>
To: XXXXXX <hliu at whoi.edu>
CC: "xliang at mit.eduXXXecco-support"@mit.edu, linxiao <linxiaop at ouc.edu.cn>
Hao,
Re question 2) Yes use calc_UEVNfromUXVY.m As before, I recommend that
you do the calculations on the native grid (using Gael Forget's gcmfaces
matlab package) and then prior to creating a plot use
calc_UEVNfromUXVY.m and regrid.
(If you are only interested in the region where the grid is rectangular
in lat and lon, i.e. Y=[-69.4712 : 56.2021], there are simpler methods.)
Best,
Diana
======================================================================
Diana Spiegel office: 54-1423
Oceanography (EAPS Dept) phone: 617-253-2092 <tel:617-253-2092>
Massachusetts Institute of Technology
======================================================================
On Fri, Oct 30, 2015 at 1:43 PM, 浩刘 <hliu at whoi.edu
<mailto:hliu at whoi.edu>> wrote:
Hi Diana，
Thank you so much on last month support.I mapped wind stress and
Ekman pumping perfectly.
Now I have stuck in three problems on_subsurface_ temperature
tendency .I wonder if you would be so kind to help me .(1) If the
studied pressure below 100m in the ocean,should I still use
shortwave penetration (derived by oceQsw) to close the temperature
tendency equation?(2) ADVx_TH and ADVy_TH are vectors, so if I want
to gridded on my personal map,should I still use calc_UEVNfromUXVY.m
before re-gridding ?(3)If I want to calculate temperature tendency
due to advection, should I do it in the following step:
(3.1)CellVol = RAC(ix,iy) * DRF(iz) * hFacC(ix,iy,iz);
*(3.2)Adv_tend *= - [ {ADVr_TH(ix,iy,iz) –
ADVr_TH(ix,iy,iz+1)}/CellVol +
{ADVx_TH(ix+1,iy,iz) – ADVx_TH(ix,iy,iz)}/CellVol +
{ADVy_TH(ix,iy+1,iz) – ADVy_TH(ix,iy,iz)}/CellVol ]
Thank you again. Looking forward to hearing from you.
Happy Halloween!
Best Regards,Hao
On 2015/9/25 17:45, Diana Lees Spiegel wrote:
> Hello Hao,
>
> Your question raises two issues. Xinfeng has answered your concern
> about the discontinuity by referring you to calc_UEVNfromUXVY.m.
> But I'd like to underscore another issue.
>
> In general, for the best accuracy, it is recommended that you
> perform your calculations on the ECCO V4 original grid. Spatial
> interpolation on the model grid can introduce inaccuracies in
> calculations, especially outside of the region where the model
> uses a rectangular latlon grid (i.e. for j=[62:239] within the
> array format, which corresponds to Y=[-69.4712:56.2021]) Avoiding
> interpolation is particularly crucial if you plan to calculate
> transports.
>
> To calculate on the original model grid you have two choices.
> After loading both oceTAUX and oceTAUY you could use convert2array
> on each, calculate the time means using the array formats of
> oceTAUX and oceTAUY, and then again use convert2array to convert
> each back to gcmfaces format. Or you could simply calculate time
> means directly in the gcmfaces format for both oceTAUX and
> oceTAUY. In either case you would then use calc_UEVNfromUXVY.m to
> calculate oceTAUXE and oceTAUYN from oceTAUX and oceTAUY.
> oceTAUXE would be identical in either case and plots without
> displaying a discontinuity.
>
> For a quick and dirty plot (as in qwckplot.m) you can plot an
> array field in array coordinates using pcolor or imagesc. To plot
> in geographic coordinates see the gcmfaces routine
> plot_one_field.m for examples of convert2pcol.m and other plotting
> methods. Interpolating the array field to a regular grid and
> plotting with your favorite Matlab 2D plotting routine may produce
> results similar to convert2pcol, but I'm not certain the polar cap
> will be correct.
>
> Hope this is clear,
> Diana Spiegel
>
>
> ======================================================================
> Diana Spiegel office: 54-1412
> Oceanography (EAPS Dept) phone: 617-253-2092
> <tel:617-253-2092>
> Massachusetts Institute of Technology
> ======================================================================
>
>
> On Wed, Sep 23, 2015 at 3:57 PM, Xinfeng Liang <xliang at mit.edu
> <mailto:xliang at mit.edu>> wrote:
>
> Hi Hao,
>
> As shown in Figure 3 in the gcmfaces manual
> (http://mitgcm.org/viewvc/*checkout*/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces.pdf),
> the whole domain includes five different faces, three of which
> need to be rotated to get the zonal and meridional components
> of vectors. If I remember correctly, this can be done using
> calc_UEVNfromUXVY.m (included in the gcmfaces package). Please
> let us know whether that works for you.
>
> Cheers,
> Xinfeng
>
> > On Sep 23, 2015, at 3:34 PM, ecco-support-request at MIT.EDU
> <mailto:ecco-support-request at MIT.EDU> wrote:
> >
> > Send ecco-support mailing list submissions to
> > ecco-support at mit.edu <mailto:ecco-support at mit.edu>
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> > http://mailman.mit.edu/mailman/listinfo/ecco-support
> > or, via email, send a message with subject or body 'help' to
> > ecco-support-request at mit.edu
> <mailto:ecco-support-request at mit.edu>
> >
> > You can reach the person managing the list at
> > ecco-support-owner at mit.edu <mailto:ecco-support-owner at mit.edu>
> >
> > When replying, please edit your Subject line so it is more
> specific
> > than "Re: Contents of ecco-support digest..."
> >
> >
> > Today's Topics:
> >
> > 1. A question related to WIND data of ecco v4 (??)
> >
> >
> >
> ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Wed, 23 Sep 2015 11:20:34 -0400
> > From: ?? <hliu at whoi.edu <mailto:hliu at whoi.edu>>
> > Subject: [ecco-support] A question related to WIND data of
> ecco v4
> > To: ecco-support at mit.edu <mailto:ecco-support at mit.edu>
> > Message-ID: <5602C342.4080702 at whoi.edu
> <mailto:5602C342.4080702 at whoi.edu>>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Hi,
> > I mapped a gridded zonal wind stress data ,and
> find a
> > discontinuity around 150E and 40W,attachment is what I did.
> > Why is there a discontinuity?
> > How I do it: First I use gcmfaces to load oceTAUX
> ,then I
> > use convert2array code to transfer to array data ,then I
> griddata them
> > in matlab ,and at last mean the whole period.
> > Looking forward to hearing from you!
> > Sincerely,Hao
> > -------------- next part --------------
> > A non-text attachment was scrubbed...
> > Name: gridded TAUX by ecco version4 from 1993 to 2011.jpg
> > Type: image/jpeg
> > Size: 67296 bytes
> > Desc: not available
> > Url :
> http://mailman.mit.edu/mailman/private/ecco-support/attachments/20150923/0a2d5bd9/attachment.jpg
> >
> > ------------------------------
> >
> > _______________________________________________
> > ecco-support mailing list
> > ecco-support at mit.edu <mailto:ecco-support at mit.edu>
> > http://mailman.mit.edu/mailman/listinfo/ecco-support
> >
> >
> > End of ecco-support Digest, Vol 4, Issue 5
> > ******************************************
>
>
> _______________________________________________
> ecco-support mailing list
> ecco-support at mit.edu <mailto:ecco-support at mit.edu>
> http://mailman.mit.edu/mailman/listinfo/ecco-support
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.mit.edu/mailman/private/ecco-support/attachments/20151030/fe826695/attachment-0001.html
More information about the ecco-support
mailing list