bringing some ISSM/MITgcm files over from the svn slr repo

Author: menemenlis

test/MITgcm/documents/NOTES

@@ -0,0 +1,54 @@
+The notes below and files in this directory initially originate from the slr svn repo:
+svn --username ????? co http://issm.ess.uci.edu/svn/slr
+
+
+
+
+{{{ Meeting with Tony Lee.
+Leverage model we have: mass ocean + ice. Atmopheric constraints. 
+
+CMIP: projection RCP scenarios. 2D forcings for atmosphere. Ask Peter Gallagher for help.
+Decadal Prediction:  no point of contact/website that has all the decadal prediction archives.  UK Met
+Depresys Doug Smith, NCAR runs (Danabasoglu Gokhan) - includes GFDL among them (Tony Rosatti).
+
+Issue: do we stop there to figure out the enveloppe of regional sl under decadal or CMIP forcings. 
+Or do we do an extra step: implication for inundations. Storm surge wave statistics to give you a bigger enveloppe. 
+We don't have expertise in storm surge or extreme events.
+
+NCAR wave forecast, lead is  Hendrick Tolman NOAA operational forecast system, RFTOOS.
+
+Issue of timeframe? 10 year, 50 years.  Max and minimum number? 
+
+Forcing comes from a coupled system: if the coupled system does not have GRACE errors, ice shelves -> does it make sense to use that forcing?  Mindful about potential limitations: 
+- CMIP6 will cointinuously improve these constraoints, so that's ok. 
+- ongoing effort to modify, alleviate these concerns by MITgcm and CMIP5 constraints (ask Dimitris about this)
+ask Dimitris more about it. 
+
+GMAO uses Poseidon, thinking about changing to MOM3. Since MOM4 got released, nothing has happened.
+
+}}}
+
+
+
+
+From J. Wahr: 
+- wanted to know whether we want to match altimetry or whether we want 
+- spatial resolution issue. Ocean function is 1degx1deb -> need to increase that. 
+- same thing for spherical harmonics that are currently runnning at 1deg. 
+- other issue is that code relies on ocean funcion that does not change in
+  time. This is ok if we are sticking to 10-100 yr projections, but it's not
+  ok if we are targetting 1000 yrs projections where this would become an
+  issue. There is currently no solution to this matter.
+
+Here is the email referring to some of theses points: 
+Isabella, Eric:  I've now modified the sea level code to include rotational
+feedback.  Before sending it to Isabella to be archived at Irvine, I have a
+question.  I'm not sure exactly what JPL wants to do with the code, but right
+now it's set up to output changes in sea level as would be seen by a tide
+gauge.  That is, it outputs: (change in geocentric sea surface height) minus
+(change in geocentric sea floor height).  An altimeter sees, instead, just the
+geocentric sea surface height.  For the JPL applications you have in mind,
+would you rather have the program output the latter (i.e. geocentric sea
+surface height)?
+ 
+-  John

test/MITgcm/documents/README

@@ -0,0 +1,85 @@
+------------------------------------------------------------------------------
+
+SLR
+Sea Level Rise  model: a framework for modeling sea level rise using 
+JPL/UCI software in ice sheet modeling (ISSM), post-glacial rebound (PGR 
+from Erik Ivins), ocean circulation (MITgcm) and sea-level models (UCI 
+model from J. Wahr and I. Velicogna and JPL Marsal model, Felix Landerer)
+
+------------------------------------------------------------------------------
+
+Jet Propulsion Laboratory
+
+Eric    Larour (ISSM project manager and SLR lead developer) 
+Erik    Ivins  (PGR lead developer)
+Dimitris Menemenlis (ECCO-ICE PI and MITgcm developer)
+Felix   Landerer (MARSAL scientist/developer)
+Isabella Velicogna (UCI SLM model, scientist/developer)
+
+Email: [email protected]
+
+------------------------------------------------------------------------------
+
+Website: http://slr.jpl.nasa.gov/ for additional information.
+
+------------------------------------------------------------------------------
+
+svn --username ????? co http://issm.ess.uci.edu/svn/slr
+
+will download this directory and subdirectories under a version system called svn.
+
+Once you have downloaded the slr code, you should add the following to your bashrc, or shrc or kshrc: 
+export SLR_DIR=$HOME/slr/  (update to your own directory)
+source $SLR_DIR/etc/environment.sh
+
+if you are under tcsh or csh, do this to your .cshrc: 
+setenv  SLR_DIR '/Users/larour/slr'  (update to your own directory)
+source $SLR_DIR/Etc/environment.sh
+
+Download and install ISSM as per instructions in
+slr/components/issm/README
+
+What we have in the code right now: 
+README: this file
+etc/ where we have our variable settings
+components: where we have all our software ISSM, PGR, MITgcm, GEOS-5.
+external packages: all our external libraries that we will be working with now. 
+
+Web-based tracked svn changes are here: https://slr-svn.jpl.nasa.gov/trac/slr
+
+------------------------------------------------------------------------------
+
+License
+
+SLR is California Institute of Technology Copyright
+
+
+Copyright (c) 2002-2011, California Institute of Technology.
+All rights reserved.  Based on Government Sponsored Research under contracts
+NAS7-1407 and/or NAS7-03001.
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+1. Redistributions of source code must retain the above copyright notice,
+   this list of conditions and the following disclaimer.
+2. Redistributions in binary form must reproduce the above copyright
+   notice, this list of conditions and the following disclaimer in the
+   documentation and/or other materials provided with the distribution.
+3. Neither the name of the California Institute of Technology (Caltech),
+   its operating division the Jet Propulsion Laboratory (JPL), the National
+   Aeronautics and Space Administration (NASA), nor the names of its
+   contributors may be used to endorse or promote products derived from
+   this software without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
+ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
+WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+DISCLAIMED. IN NO EVENT SHALL THE CALIFORNIA INSTITUTE OF TECHNOLOGY BE LIABLE
+FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+------------------------------------------------------------------------------

test/MITgcm/documents/issm_helene_notes.docx

@@ -0,0 +1,45 @@
+% read test4003 MItgcm output
+ISSM_DIR=getenv('ISSM_DIR');
+eval (['cd ' ISSM_DIR '/test/NightlyRun/run'])
+fnm=dir('surfDiag.*.data');
+fld={'SHIfwFlx (kg/m^2/s) Ice shelf fresh water flux (positive upward)', ...
+     'SHIhtFlx (W/m^2   ) Ice shelf heat flux  (positive upward)', ...
+     'SHIgammT (m/s     ) Ice shelf exchange coefficient for theta', ...
+     'SHIgammS (m/s     ) Ice shelf exchange coefficient for salt', ...
+     'SHI_mass (kg/m^2  ) dynamic ice shelf mass for surface load anomaly', ...
+     'SHIRshel (m       ) depth of shelfice', ...
+     'SI_Uvel  (m/a     ) Ice stream x-velocity', ...
+     'SI_Vvel  (m/a     ) Ice stream y-velocity', ...
+     'SI_Thick (m       ) Ice stream thickness', ...
+     'SI_hmask (none    ) Ice stream thickness mask', ...
+     'SI_float (none    ) Ice stream grounding ind', ...
+     'SHIuStar (m/s     ) Friction velocity at bottom of ice shelf'}
+nx=3; ny=200; nf=length(fld); nt=length(fnm);
+Diag=zeros(nx,ny,nf,nt);
+for t=1:nt
+    Diag(:,:,:,t)=readbin(fnm(t).name,[nx ny nf]);
+end
+
+% plot output
+mkdir figs
+orient tall
+wysiwyg
+colormap(cmap)
+for t=1:nt
+    clf
+    for f=1:nf
+        subplot(nf,1,f)
+        if t==1
+            mypcolor(Diag(:,:,f,t));
+        else
+            mypcolor(Diag(:,:,f,t)-Diag(:,:,f,1));
+        end
+        colorbar
+        title(fld{f},'Interpreter','none')
+        set(gca,'YTickLabel',[])
+        if f<nf
+            set(gca,'XTickLabel',[])
+        end
+    end
+    eval(['print -djpeg figs/surfDiag' myint2str(t)])
+end