[MITgcm-support] ptracers help/manual
Christopher L. Wolfe
clwolfe at ucsd.edu
Mon May 21 17:12:26 EDT 2007
Manfredi,
Thanks. This is a big help.
Christopher
-----------------------------------------------------------
Dr. Christopher L. Wolfe 858-534-4560
Physical Oceanography Research Division OAR 357
Scripps Institution of Oceanography, UCSD clwolfe at ucsd.edu
-----------------------------------------------------------
On May 21, 2007, at 1:14 PM, Manfredi Manizza wrote:
>
> Hi Christopher,
>
> I hope I rightly understood what you need.
>
> I guess you want to use an input file of a assive tracer where you set
> a
> value you pick and you want to keep it fixed in time for the whole
> simulation as a constant source of tracer that you want to follow as
> advected and diffused by the circulation fields.
>
> IF this IS the case you need to have :
>
> File # 1 - An input file corresponding to your tracer source X by Y
> (your domain) with values you set that could be the same everywhere or
> what is more appropriate for your purpose like a source only in a
> fraction of the domain. Let's call it FILE_TRC.bin
>
> File # 2 - A correpsoning file with same dimension where in the
> position (i,j) of your domain you put 1 if you want to have a constant
> source or 0 if you do not.
> Let's call it RESTORE_FILE.bin
>
> File # 3 - IF you run a regional domain you would need a "sponge
> layer" at the boundaries in order to prevent that what goes outside
> enters back from the other side. Let's call RELAX_FILE.bin
>
> To generate this third file you need to set the value all around the
> boundary equal to 1 and apply an exponential fuction that sets the
> value of decay of the tracer that is maximum at the boundary (1)
> and then decreases as soon as you move inside your domain.
>
> If this third point is not clear, I can send you the matlab script
> that does that so you can adjust it to your domain.
>
> So then, once you have the files ready, you follow this below :
>
> In the "data.ptracers" (in the run direcory) file, you MUST have :
>
> &PTRACERS_PARM01
> PTRACERS_numInUse=1,
> PTRACERS_Iter0= 0,
> #
> # The following option is now obsolete:
> # PTRACERS_write_mnc=.TRUE.,
> #
> # tracer 1 -
> PTRACERS_names(1)='trc1',
> PTRACERS_long_names(1)='trc1',
> PTRACERS_units(1)='mol/m^3',
> PTRACERS_advScheme(1)=77,
> PTRACERS_diffKh(1)=0.E3,
> PTRACERS_diffKr(1)=1.E-5,
> PTRACERS_useGMRedi(1)=.FALSE. ,
> PTRACERS_useKPP(1)=.FALSE. ,
> PTRACERS_initialFile(1)='FILE_TRC.BIN',
> &
>
> In your "data.rbcs" (in the run directory) you MUST have :
> # Relax Boundary Conditions parameters
> &RBCS_PARM01
> tauRelaxT = 0.,
> tauRelaxS = 0.,
> relaxMaskFile(1) = ' ',
> relaxMaskFile(2) = ' ',
> relaxMaskFile(3) = 'RELAX_FILE.bin',
> relaxTFile = ' ',
> relaxSFile = ' ',
> rbcsIniter = 0,
> rbcsForcingPeriod = 15552000.,
> rbcsForcingCycle = 0.,
> useRBCtemp = .FALSE.,
> useRBCsalt = .FALSE.,
> useRBCptracers = .TRUE.,
> &end
>
> &RBCS_PARM02
> useRBCptrnum(1)=.TRUE.
> tauRelaxPTR(1) = 1200.,
> relaxPtracerFile(1) = 'RESTORE_FILE.bin',
> &end
>
> In this way you should be able to not to touch the code
> for your experiments just changing those input files
> should be enough.
>
> I guess this should work as it has been working for my regional Artcic
> Ocean regional set-up if I understood what you needed.
>
> Cheers
> Manfredi
>
>
>
>
>
> On Mon, 21 May 2007, Christopher L. Wolfe wrote:
>
>>
>> Hello modelers,
>>
>> Does anyone know if there is any documentation for the ptracers
>> package? All I can find in the manual is the offline tutorial and a
>> few comments in the package source, but that's not totally helpful.
>> The tutorial shows how to set up ptracers as an initial value
>> problem, but I want to use it as a boundary value problem.
>> Specifically, I want to fix the surface distribution of the tracers
>> and see how the surface values are transported through the model. Is
>> this possible with the ptracers package? If so, how does one do it?
>>
>> Thanks,
>> Christopher
>>
>> -----------------------------------------------------------
>> Dr. Christopher L. Wolfe 858-534-4560
>> Physical Oceanography Research Division OAR 357
>> Scripps Institution of Oceanography, UCSD clwolfe at ucsd.edu
>> -----------------------------------------------------------
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>
> --
>
>
> Manfredi Manizza
> Postdoctoral Associate
> Room 54-1423
> Dept. of Earth, Atmospheric
> and Planetary Sciences
> Massachusetts Institute of Technology
> 77 Massachusetts Avenue
> Cambridge, MA 02139-4307, USA
> phone: 617-452-2977
> fax : 617-253-4464
> email: mmanizza at ocean.mit.edu
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
More information about the MITgcm-support
mailing list