[MITgcm-support] Sample pTracer using flt package
钱钰坤
qianyk at mail3.sysu.edu.cn
Fri Jun 1 06:19:06 EDT 2018
Hi all,
I've inserted an passive tracer (using ptracer package) and deployed
Lagrangian particles (using flt package) simultaneously into a dynamic flow.
I wanted to sample the tracer value along each Lagrangian float but I only
found the following five lines of code in flt_traj.F:
CALL FLT_BILINEAR (ix,jy,uu,uVel, kp,1,bi,bj,myThid)
CALL FLT_BILINEAR (ix,jy,vv,vVel, kp,2,bi,bj,myThid)
CALL FLT_BILINEAR2D(ix,jy,pp,etaN, 0,bi,bj,myThid)
CALL FLT_BILINEAR (ix,jy,tt,theta, kp,0,bi,bj,myThid)
CALL FLT_BILINEAR (ix,jy,ss,salt, kp,0,bi,bj,myThid)
which means the flt can only sample the basic flow states (u, v, eta, theta, and salt).
I wonder if there is an easy way to achieve this by modifying the code slightly as:
CALL FLT_BILINEAR (ix,jy,ss, pTracer(, , , , , iTracer), kp,0,bi,bj,myThid)
and using the result to replace the salt which I do not care about?
If yes, how to modify the #include or variable declarations at the beginning of flt_traj.F?
Many thanks.
------------------
Best regards
Yu-Kun Qian (钱钰坤)
Center for Monsoon and Environment Research
Department of Atmospheric Sciences
School of Environmental Science and Engineering
Sun Yat-sen University
No. 135 Xingang West Road, Haizhu District
Guangzhou, 510275, P.R. China
Tel; 020-84115227
Email: qianyk at mail3.sysu.edu.cn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180601/28218de6/attachment.html>
More information about the MITgcm-support
mailing list