<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Gael,<div class=""><br class=""></div><div class="">thanks for the pointers. I will try it out and report on the results.</div><div class=""><br class=""></div><div class="">Cheers,</div><div class=""><br class=""><div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;">Paola<br class=""><br class=""><span style="font-size: 12px;" class="">----------------------------------------------<br class="">Paola Cessi   <br class="">Scripps Institution of Oceanography  <br class="">9500 Gilman Drive #0213               <br class="">La Jolla, CA 92093-0213<br class="">USA               <br class="">Tel: +1 858 534 0622   Fax: +1 858 534 8045<br class="">e-mail: <a href="mailto:pcessi@ucsd.edu" class="">pcessi@ucsd.edu</a>     Web:  <a href="http://www-pord.ucsd.edu/~pcessi" class="">http://www-pord.ucsd.edu/~pcessi</a><br class=""><br class="">Delivery address (Fed-Ex, DHL, UPS, etc.):<br class="">Paola Cessi<br class="">University of California San Diego<br class="">Keck Center - Room 366<br class="">8851 Shellback Way<br class="">La Jolla, CA 92093-0213<br class="">USA</span></div><div style="color: rgb(0, 0, 0); font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-size: 12px;" class=""><br class=""></span></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline">
</div>
<br class=""><div style=""><blockquote type="cite" class=""><div class="">On May 11, 2020, at 6:44 PM, Gael Forget <<a href="mailto:gforget@mit.edu" class="">gforget@mit.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Paola<div class=""><br class=""></div><div class="">The meta file created by this package indeed does not seem to include the list of variables. However, <span style="" class="">I think </span>it should be the same list of 9 variables described in <a href="https://mitgcm.readthedocs.io/en/latest/outp_pkgs/outp_pkgs.html#input-files" class="">https://mitgcm.readthedocs.io/en/latest/outp_pkgs/outp_pkgs.html#input-files</a> . Below is an exert from flt_init_varia.F which is consistent with the documented list. </div><div class=""><br class=""></div><div class="">Its important to note that the first record (i.e. 9 element vector) is special. This is stated in the doc but not excessively clear and kind of tricky imho.</div><div class=""><br class=""></div><div class="">In any event, you should be able to read these files in Matlab using <b class="">rdmds.m</b> as usual (a copy is in <a href="https://github.com/MITgcm/gcmfaces" class="">https://github.com/MITgcm/gcmfaces</a>). Doing it this way may read all tile pickups in one multidimensional array.  </div><div class=""><br class=""></div><div class="">Alternatively if you want more granular control on read and write, you can use <span style="" class=""><b class="">read2memory</b> and <b class="">write2file</b> (also in </span><a href="https://github.com/MITgcm/gcmfaces" class="">gcmfaces</a><span style="" class="">; try e.g. `help </span><span style="" class="">read2memory`). </span></div><div class=""><span style="" class=""><br class=""></span></div><div class=""><span style="" class="">In an example where nrecords is 74 (from the meta file), this </span><span style="" class="">should be</span> something like this :</div><div class=""><br class=""></div><div class="">>> tmp2=read2memory('pickup_flt.ckptA.001.001.data',[9 74],64)</div><div class=""><br class=""></div><div class="">I copy below the first and 2nd records for reference. All of the others look like #2. </div><div class=""><br class=""></div><div class="">And then writing the modified pickup to a new file would look like this :</div><div class=""><br class=""></div><div class="">>> write2file(‘modified/pickup_flt.ckptA.001.001.data',tmp2,64)</div><div class=""><br class=""></div><div class="">Hope this helps. Sorry this stuff is still not documented well enough. Good luck with your experiment.</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Gael</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div class="">>> tmp2(:,1)</div><div class="">          73</div><div class="">       18001</div><div class="">    10800600</div><div class="">           0</div><div class="">           0</div><div class="">         300</div><div class="">           9</div><div class="">           0</div><div class="">           0</div></div><div class=""><br class=""></div><div class=""><div class="">>> tmp2(:,2)</div><div class="">   12.0000</div><div class="">  600.0000</div><div class="">   13.2794</div><div class="">   13.2823</div><div class="">    3.0000</div><div class="">         0</div><div class="">         0</div><div class="">         0</div><div class="">   -1.0000</div></div><div class=""><br class=""></div><div class="">From <span style="" class="">flt_init_varia.F :</span></div><div class=""><br class=""><div class=""><div class="">152               npart( npart_tile(bi,bj),bi,bj) = tmp(1)</div><div class="">153               tstart(npart_tile(bi,bj),bi,bj) = tmp(2)</div><div class="">154               ipart( npart_tile(bi,bj),bi,bj) = ix</div><div class="">155               jpart( npart_tile(bi,bj),bi,bj) = jy</div><div class="">156               kpart( npart_tile(bi,bj),bi,bj) = kz</div><div class="">157               kfloat(npart_tile(bi,bj),bi,bj) = tmp(6)</div><div class="">158               iup(   npart_tile(bi,bj),bi,bj) = tmp(7)</div><div class="">159               itop(  npart_tile(bi,bj),bi,bj) = tmp(8)</div><div class="">160               tend(  npart_tile(bi,bj),bi,bj) = tmp(9)</div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><blockquote type="cite" class=""><div class="">On May 11, 2020, at 8:18 PM, Paola Cessi <<a href="mailto:pcessi@ucsd.edu" class="">pcessi@ucsd.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">Hi Dimitri,</span><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><br class=""></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">I am using matlab. These are pickups files from the flt package, so I am not sure they are the same as the regular MITgcm pickups. They are written and read with the mdsio package.</div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><br class=""></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;">I _think_ that what I need to change is the 9th entry of the “arr" array, i.e. arr(9) written/read by the routines<span class="Apple-converted-space"> </span><span class="final-path" style="background-color: rgb(255, 255, 255); box-sizing: border-box; color: rgb(36, 41, 46);">mdsio_wr_rec_rl.F</span><span class="" style="color: rgb(88, 96, 105); background-color: rgb(255, 255, 255);">/</span><span class="" style="background-color: rgb(255, 255, 255); color: rgb(36, 41, 46);">mdsio_rd_rec_rl.F. That said, I don’t understand how to locate this variable given the .meta file (attached)</span><span class="" style="font-family: -apple-system, system-ui, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji"; font-size: 16px; background-color: rgb(255, 255, 255); color: rgb(36, 41, 46);">.</span></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;"><span class="" style="font-family: -apple-system, system-ui, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji"; font-size: 16px; background-color: rgb(255, 255, 255); color: rgb(36, 41, 46);"><br class=""></span></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;"><font color="#24292e" class=""><span class="" style="background-color: rgb(255, 255, 255);">Many thanks,</span></font></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;"><span class="" style="orphans: 2; widows: 2; background-color: rgb(255, 255, 255);"><font color="#24292e" class=""><br class=""></font></span></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;"></div><span id="cid:135AFDE9-C993-4306-8966-BC78D56FEAC0@fios-router.home" class=""><pickup_flt.ckptA.001.001.meta></span><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; orphans: 2; widows: 2;"></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><br class=""></div><div class="" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><br class=""><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div style="font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Paola<br class=""><br class=""><span class="" style="font-size: 12px;">----------------------------------------------<br class="">Paola Cessi   <br class="">Scripps Institution of Oceanography  <br class="">9500 Gilman Drive #0213               <br class="">La Jolla, CA 92093-0213<br class="">USA               <br class="">Tel: +1 858 534 0622   Fax: +1 858 534 8045<br class="">e-mail: <a href="mailto:pcessi@ucsd.edu" class="">pcessi@ucsd.edu</a>     Web:  <a href="http://www-pord.ucsd.edu/~pcessi" class="">http://www-pord.ucsd.edu/~pcessi</a><br class=""><br class="">Delivery address (Fed-Ex, DHL, UPS, etc.):<br class="">Paola Cessi<br class="">University of California San Diego<br class="">Keck Center - Room 366<br class="">8851 Shellback Way<br class="">La Jolla, CA 92093-0213<br class="">USA</span></div><div class="" style="font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><span class="" style="font-size: 12px;"><br class=""></span></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class=""><div class=""><blockquote type="cite" class=""><div class="">On May 11, 2020, at 4:41 PM, Dimitris Menemenlis <<a href="mailto:dmenemenlis@gmail.com" class="">dmenemenlis@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Paola, do you use matlab or python?<div class=""><br class=""></div><div class="">If matlab and if noone answers with a generic answer, please send your pickup*.meta file and let me know which variable you want to replace and I will send you two lines of matlab code (and two small scripts) that will read and write the desired variable from and to the pickup file.</div><div class=""><br class=""></div><div class="">Cheers, Dimitris<br class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On May 11, 2020, at 3:39 PM, Paola Cessi <<a href="mailto:pcessi@ucsd.edu" class="">pcessi@ucsd.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">I need to change one entry in the pick-up files. Is there a simple sample program to read a pickup file change, one entry and write it again with the entry changed?<div class=""><br class=""></div><div class="">Thanks,</div><div class=""><br class=""><div class=""><div class="" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;"><div class="" style="font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;">Paola<br class=""><br class=""><span class="" style="font-size: 12px;">----------------------------------------------<br class="">Paola Cessi   <br class="">Scripps Institution of Oceanography  <br class="">9500 Gilman Drive #0213               <br class="">La Jolla, CA 92093-0213<br class="">USA               <br class="">Tel: +1 858 534 0622   Fax: +1 858 534 8045<br class="">e-mail: <a href="mailto:pcessi@ucsd.edu" class="">pcessi@ucsd.edu</a>     Web:  <a href="http://www-pord.ucsd.edu/~pcessi" class="">http://www-pord.ucsd.edu/~pcessi</a><br class=""><br class="">Delivery address (Fed-Ex, DHL, UPS, etc.):<br class="">Paola Cessi<br class="">University of California San Diego<br class="">Keck Center - Room 366<br class="">8851 Shellback Way<br class="">La Jolla, CA 92093-0213<br class="">USA</span></div><div class="" style="font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><span class="" style="font-size: 12px;"><br class=""></span></div><br class="Apple-interchange-newline"></div><br class="Apple-interchange-newline"></div><br class=""></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!QoirxRdRyr4vnrPbSnaKY7tFkekGeHWD65KRYp5_8wt09DuCmWRzUfDbDwwgxw$" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class=""></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class=""><a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!QoirxRdRyr4vnrPbSnaKY7tFkekGeHWD65KRYp5_8wt09DuCmWRzUfDbDwwgxw$" class="">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!QoirxRdRyr4vnrPbSnaKY7tFkekGeHWD65KRYp5_8wt09DuCmWRzUfDbDwwgxw$</a><span class="Apple-converted-space"> </span><br class=""></div></blockquote></div><br class=""></div><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">MITgcm-support mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="mailto:MITgcm-support@mitgcm.org" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">MITgcm-support@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""></div></body></html>