[MITgcm-support] tracer concentration

[Jean-Michel Campin]

Hi,

Matt's formula is not always giving the right integral
(e.g., with CS-grid); better to take directly RAC rather
than DXG.*DYG .

Jean-Michel

On Sat, Feb 26, 2011 at 01:22:47PM -0800, Matthew Mazloff wrote:
> Hi Angela,
>
> I think you need to volume weight your integral
>
> try
>> tr01(:,:,:)=rdmds('TRAC01',currentiter);
>> sum(sum(sum(tr01.*DXG.*DYG.*drf.*hFacC)))
>
>
> -Matt
>
>
> On Feb 26, 2011, at 1:18 PM, Angela Marie Zalucha wrote:
>
>> I checked and trcstat_ptracer01_mean is staying quite constant (less  
>> than
>> 1% difference).  But when I read into matlab the total tracer
>> concentration is not constant using:
>>
>> tr01(:,:,:)=rdmds('TRAC01',currentiter);
>> sum(sum(sum(tr01)))
>>
>>
>>   Angela
>>
>>
>>
>> On Fri, 25 Feb 2011, Jean-Michel Campin wrote:
>>
>>> Hi Angela,
>>>
>>> If you consider the output of ptracer monitor  
>>> (trcstat_ptracer01_mean),
>>> it give the true global average tracer concentration (and account for
>>> hFac, both the topography + r* streching factor).
>>> If not fluid is added (dynstat_eta_mean = constant), and no source- 
>>> sink,
>>> trcstat_ptracer01_mean should stay constant.
>>> Could you check this ?
>>>
>>> And assuming that trcstat_ptracer01_mean is drifting, might want to  
>>> have a look at:
>>> a) all *OPTIONS.h files in your "code" dir
>>> b) all parameter files data* (from your run dir).
>>> If it's OK to put those in a tar file.
>>>
>>> Thanks,
>>> Jean-Michel
>>>
>>> On Fri, Feb 25, 2011 at 12:45:49PM -0500, Angela Marie Zalucha wrote:
>>>> Hi,
>>>>
>>>>  I don't really know what you mean by CD scheme.  I am using r*
>>>> coordinates and my simulation is based on the hs94.cs-32x32x5  
>>>> experiment
>>>> (but with a different radiation scheme and 30 levels).  The ptracers
>>>> package I am using in the default setting except to specify the  
>>>> specific
>>>> initial tracer locations.  My data.ptracers file is
>>>>
>>>> &PTRACERS_PARM01
>>>> PTRACERS_numInUse=1,
>>>> PTRACERS_advScheme=81,
>>>> PTRACERS_diffKh=0.,
>>>> PTRACERS_diffK4=0.E9,
>>>> PTRACERS_diffKr=0.,
>>>> PTRACERS_Iter0=201000,
>>>>
>>>> Is Yuan's suggestion correct (that I need to be multiplying the  
>>>> tracers
>>>> by hfac in the post-processing since I am using r*)?
>>>>
>>>>  Angela
>>>>
>>>> On Fri, 25 Feb 2011, Jean-Michel Campin wrote:
>>>>
>>>>> Hi Angela,
>>>>>
>>>>> MITgcm conserve passive tracer (even with CD-scheme on), i.e.,
>>>>> the total content follows the integrated source-sink, and stay
>>>>> constent if no source-sink.
>>>>> Could you give more details about your problem/set-up (otherwise,
>>>>> will not be able to provide a usefull advise).
>>>>>
>>>>> Thanks,
>>>>> Jean-Michel
>>>>>
>>>>> On Thu, Feb 24, 2011 at 05:44:28PM -0500, Angela Marie Zalucha  
>>>>> wrote:
>>>>>> The tracers are passive.
>>>>>>
>>>>>> There is no change in atmospheric mass with time.
>>>>>>
>>>>>> There are no sources and sinks of tracers.
>>>>>>
>>>>>>
>>>>>>  Angela
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, 24 Feb 2011, Yuan Lian wrote:
>>>>>>
>>>>>>> Are you running the ptracers package with passive tracers? Is any
>>>>>>> change in the atmospheric mass if
>>>>>>> the tracer is passive? Do you have any sink or source of  
>>>>>>> tracers? I
>>>>>>> don't see much problems
>>>>>>> with most of the advection schemes for passive tracers, 
>>>>>>> though some
>>>>>>> deviation exists but should be
>>>>>>> smaller than 0.1% in most cases.
>>>>>>>
>>>>>>> If the total atmospheric mass is changing, make sure the 
>>>>>>> tracer mixing
>>>>>>> ratio is updated properly
>>>>>>> (i.e, scaling the mixing ratio in columns).This should also  
>>>>>>> work if the
>>>>>>> tracer is semi-passive, i.e., the
>>>>>>> mass loading of tracer is considered in the model.
>>>>>>>
>>>>>>> If you have sink and source of tracers, make sure they are
>>>>>>> self-consistent. Or if you are updating tracer
>>>>>>> tendencies in physics grid, make sure the schemes are  
>>>>>>> conservative.
>>>>>>>
>>>>>>> Cheers,
>>>>>>> Yuan
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On 2/24/11 1:17 PM, Angela Marie Zalucha wrote:
>>>>>>>> Hello, I am working with the ptracers package.  I am 
>>>>>>>> finding that
>>>>>>>> after running for about 5 days, the tracers are not 
>>>>>>>> conserving mass
>>>>>>>> at all--the increase in mass is about 50%, which is pretty
>>>>>>>> significant.  I am using advection scheme 81 (I have tried 
>>>>>>>> 7, 30, 33,
>>>>>>>> 77, 80, but with similar results).  Is this typical?  Are  
>>>>>>>> there any
>>>>>>>> ways around this?
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>  Angela
>>>>>>>>
>>>>>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>>>>>> Angela Zalucha, Ph.D.
>>>>>>>> Postdoctoral Associate
>>>>>>>> Department of Earth, Atmospheric,
>>>>>>>>  and Planetary Sciences
>>>>>>>> Massachusetts Institute of Technology
>>>>>>>>
>>>>>>>> Office located at Southwest Research Institute
>>>>>>>> 1050 Walnut Street, Suite 300
>>>>>>>> Boulder, CO 80302
>>>>>>>> USA
>>>>>>>> (720) 208-7211
>>>>>>>>
>>>>>>>>
>>>>>>>> http://web.mit.edu/azalucha/www/index.html
>>>>>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> MITgcm-support mailing list
>>>>>>>> MITgcm-support at mitgcm.org
>>>>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> MITgcm-support mailing list
>>>>>>> MITgcm-support at mitgcm.org
>>>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>>>
>>>>>>
>>>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>>>> Angela Zalucha, Ph.D.
>>>>>> Postdoctoral Associate
>>>>>> Department of Earth, Atmospheric,
>>>>>>  and Planetary Sciences
>>>>>> Massachusetts Institute of Technology
>>>>>>
>>>>>> Office located at Southwest Research Institute
>>>>>> 1050 Walnut Street, Suite 300
>>>>>> Boulder, CO 80302
>>>>>> USA
>>>>>> (720) 208-7211
>>>>>>
>>>>>>
>>>>>> http://web.mit.edu/azalucha/www/index.html
>>>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> _______________________________________________
>>>>>> MITgcm-support mailing list
>>>>>> MITgcm-support at mitgcm.org
>>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>
>>>>> _______________________________________________
>>>>> MITgcm-support mailing list
>>>>> MITgcm-support at mitgcm.org
>>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>>>
>>>>
>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>> Angela Zalucha, Ph.D.
>>>> Postdoctoral Associate
>>>> Department of Earth, Atmospheric,
>>>>  and Planetary Sciences
>>>> Massachusetts Institute of Technology
>>>>
>>>> Office located at Southwest Research Institute
>>>> 1050 Walnut Street, Suite 300
>>>> Boulder, CO 80302
>>>> USA
>>>> (720) 208-7211
>>>>
>>>>
>>>> http://web.mit.edu/azalucha/www/index.html
>>>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>>>
>>>>
>>>>
>>>>
>>>> _______________________________________________
>>>> MITgcm-support mailing list
>>>> MITgcm-support at mitgcm.org
>>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>
>>> _______________________________________________
>>> MITgcm-support mailing list
>>> MITgcm-support at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>>
>>
>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>> Angela Zalucha, Ph.D.
>> Postdoctoral Associate
>> Department of Earth, Atmospheric,
>>   and Planetary Sciences
>> Massachusetts Institute of Technology
>>
>> Office located at Southwest Research Institute
>> 1050 Walnut Street, Suite 300
>> Boulder, CO 80302
>> USA
>> (720) 208-7211
>>
>>
>> http://web.mit.edu/azalucha/www/index.html
>> _-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-
>>
>>
>>
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support