<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Dafydd,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">This is not the correct way to apply the perturbation to ECCO V4r4 for a few reasons:<o:p></o:p></p>
<p class="MsoNormal">1) The xx_<VAR>.effective is the dimensional adjusted field, including the first guess, so using xx_<VAR> would be double counting. Note <VAR> here are control variables for initial conditions (theta, salt, uvel, vvel, and etan) and mixing
parameters (diffkr, kapgm, and kapredi). <o:p></o:p></p>
<p class="MsoNormal">2) The current V4r4 setup reads in non-dimensional control adjustment (xx_<VAR> files).<o:p></o:p></p>
<p class="MsoNormal">3) The model then applies a smoothing operator to the non-dimensional xx, and then rescale it to dimensional xx.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">So one way to achieve what you wanted to do is to perturb the dimensional control adjustment and then let the model read in the dimensional adjustment without applying rescaling and smoothing . The steps are as follows.
<o:p></o:p></p>
<p class="MsoNormal">1) Perturb the dimensional control adjustment for initial theta, xx_theta.dim.0000000129.data that can be found in
<a href="https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/other/adjustments">
https://ecco.jpl.nasa.gov/drive/files/Version4/Release4/other/adjustments</a> along with other dimensional control adjustments.<o:p></o:p></p>
<p class="MsoNormal">2) Copy the perturbed field to xx_theta.0000000129.data (You may want to keep a copy of the original xx_theta.0000000129.data for later use)<o:p></o:p></p>
<p class="MsoNormal">3) Make changes to data.ctrl by replacing <o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc(1,1)='WC01',<o:p></o:p></p>
<p class="MsoNormal">with <o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc(1,1)=' noscaling' ,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The part in data.ctrl related to xx_theta is shown below for before and after the change;<o:p></o:p></p>
<p class="MsoNormal">Before:<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_weight(1) = 'Theta_weights_nonseasonal_rmsv2_areascaled.bin',<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_file(1)='xx_theta',<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_bounds(1:5,1)=-2.0,-1.9,39.,40.,0.,<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc(1,1)='WC01',<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc_i(1,1)=1,<o:p></o:p></p>
<p class="MsoNormal">mult_genarr3d(1) = 1.,<o:p></o:p></p>
<p class="MsoNormal">After:<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_weight(1) = 'Theta_weights_nonseasonal_rmsv2_areascaled.bin',<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_file(1)='xx_theta',<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_bounds(1:5,1)=-2.0,-1.9,39.,40.,0.,<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc(1,1)=' noscaling' ,<o:p></o:p></p>
<p class="MsoNormal">xx_genarr3d_preproc_i(1,1)=1,<o:p></o:p></p>
<p class="MsoNormal">mult_genarr3d(1) = 1.,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Then you can rerun the model with the changes above. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Regards,<o:p></o:p></p>
<p class="MsoNormal">Ou<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><ecco-support-bounces@mit.edu> on behalf of Dafydd Stephenson <dafydd@ucar.edu><br>
<b>Reply-To: </b>"ECCO support list, wider membership" <ecco-support@mit.edu><br>
<b>Date: </b>Thursday, July 22, 2021 at 2:45 PM<br>
<b>To: </b>"ecco-support@mit.edu" <ecco-support@mit.edu><br>
<b>Subject: </b>[EXTERNAL] [ecco-support] Perturbing ECCOv4r4 run/relationship between xx_<VAR> and xx_<VAR>.effective<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal"><span lang="EN-GB">Hi all,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">I have managed to re-run the ECCOv4r4 simulation and am now interested in applying perturbations to the run. I have so far modified xx_theta.0000000129.data to increase the temperature everywhere by 1K as a test, but
the outcome in xx_theta.effective.0000000129.data and T.0000000001.data is in many places 30+K higher than in the unperturbed run.
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">My understanding (hopefully correct) is that a rescaling is applied to xx_<VAR> to generate xx_<VAR>.effective, so I should be able to invert this rescaling to apply an “effective” perturbation of my choosing? However,
I can’t find a weighting file amongst the input files which in any way resembles the transformation applied to my 1K perturbation.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">I’m interested to know (i) whether this is the correct way to apply to a perturbation to the ECCO run and (ii) if so, how to modify xx_<VAR> such that xx_<VAR>.effective corresponds to the perturbation I want.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Thanks!<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-GB">Dafydd<o:p></o:p></span></p>
</div>
</body>
</html>