NOAA-GFDL · MJHarrison-GFDL · Apr 11, 2015
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/.gitignore b/ice_ocean_SIS2/FLOR2/INPUT/.gitignore
@@ -0,0 +1,2 @@
+# Undo *.nc in top-level .gitignore
+!*.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/MOM_channels_FLOR b/ice_ocean_SIS2/FLOR2/INPUT/MOM_channels_FLOR
@@ -0,0 +1,9 @@
+! This file specifies restricted channel widths in MOM.  The order is:
+!  [U|V]_width, min_longitude, max_longitude, min_latitude, max_latitude, width
+
+U_width,  -5.9,  -5.4, 35.8, 36.2,  12000.0  ! Gibraltar
+U_width,  26.4,  26.8, 40.4, 40.6,   5000.0  ! Dardanelles
+
+V_width,  28.8,  29.4, 41.1, 41.2,   2500.0  ! Bosporus - should be 1km wide.
+V_width,  43.25, 43.5, 12.6, 12.8,  10000.0  ! Red Sea, Bab-el-Mendeb
+V_width, 141.5, 141.7, 52.0,  52.15, 2500.0  ! Between Sakhalin & Russia
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/WOA05_pottemp_salt.nc b/ice_ocean_SIS2/FLOR2/INPUT/WOA05_pottemp_salt.nc
@@ -0,0 +1 @@
+.datasets/obs/WOA05_pottemp_salt.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/atmos_mosaic_tile1Xland_mosaic_tile1.nc b/ice_ocean_SIS2/FLOR2/INPUT/atmos_mosaic_tile1Xland_mosaic_tile1.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/atmos_mosaic_tile1Xland_mosaic_tile1.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/atmos_mosaic_tile1Xocean_mosaic_tile1.nc b/ice_ocean_SIS2/FLOR2/INPUT/atmos_mosaic_tile1Xocean_mosaic_tile1.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/atmos_mosaic_tile1Xocean_mosaic_tile1.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/grid_spec.nc b/ice_ocean_SIS2/FLOR2/INPUT/grid_spec.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/grid_spec.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/land_mask.nc b/ice_ocean_SIS2/FLOR2/INPUT/land_mask.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/land_mask.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/land_mosaic_tile1Xocean_mosaic_tile1.nc b/ice_ocean_SIS2/FLOR2/INPUT/land_mosaic_tile1Xocean_mosaic_tile1.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/land_mosaic_tile1Xocean_mosaic_tile1.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/layer_coord.nc b/ice_ocean_SIS2/FLOR2/INPUT/layer_coord.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/ncar_precip.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/ncar_precip.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/ncar_precip.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/ncar_rad.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/ncar_rad.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/ncar_rad.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/ocean_hgrid.nc b/ice_ocean_SIS2/FLOR2/INPUT/ocean_hgrid.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocean_hgrid.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/ocean_mask.nc b/ice_ocean_SIS2/FLOR2/INPUT/ocean_mask.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/ocean_mask.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/ocean_mosaic.nc b/ice_ocean_SIS2/FLOR2/INPUT/ocean_mosaic.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/ocean_mosaic.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/q_10_mod.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/q_10_mod.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/q_10_mod.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/runoff.nc b/ice_ocean_SIS2/FLOR2/INPUT/runoff.nc
@@ -0,0 +1 @@
+.datasets/OM4_025/INPUT/runoff.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/seawifs_1998-2006_smoothed_2X.nc b/ice_ocean_SIS2/FLOR2/INPUT/seawifs_1998-2006_smoothed_2X.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/seawifs_1998-2006_smoothed_2X.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/slp.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/slp.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/slp.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/t_10_mod.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/t_10_mod.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/t_10_mod.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/tidal_amplitude.nc b/ice_ocean_SIS2/FLOR2/INPUT/tidal_amplitude.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/tidal_amplitude.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/topog.nc b/ice_ocean_SIS2/FLOR2/INPUT/topog.nc
@@ -0,0 +1 @@
+/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/topog.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/u_10_mod.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/u_10_mod.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/u_10_mod.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/v_10_mod.clim.nc b/ice_ocean_SIS2/FLOR2/INPUT/v_10_mod.clim.nc
@@ -0,0 +1 @@
+.datasets/CORE/NYF_v2.0/v_10_mod.clim.nc
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/vgrid.py b/ice_ocean_SIS2/FLOR2/INPUT/vgrid.py
@@ -0,0 +1,69 @@
+import numpy
+import math
+
+def dzIter(nk, Htot, dzTop, Huniform, fnPow, prec):
+    """
+    Optimizes the highest ratio dzFn so that the sum(dz)=Htot.
+    """
+    def dzFn(nk, Huniform, dzTop, fnPow, zFac):
+        """
+        Returns dz = dzTop * ( 1 + zFac \int \fn(k) dk ) and sum(dz).
+        """
+        dz = numpy.ones(nk) * dzTop
+        k = math.ceil( Huniform/dzTop )
+        dz = dzTop * numpy.cumprod( 1. + zFac * (fn( numpy.linspace(1, nk, nk).astype(numpy.float64), k, nk )**fnPow) )
+        return dz, numpy.sum(dz)
+
+    def fn(z, z0, z1):
+        """
+        Cosine bell function between z0 and z1 s.t. f(z<z0)=0 and f(z>z1)=0.
+        """
+        zStar = (z-z0)/(z1-z0) # non-dimensional coordinate 0..1
+        zStar = numpy.maximum(0., zStar)
+        zStar = numpy.minimum(1., zStar)
+        return 0.5*(1. - numpy.cos(2.*math.pi*zStar))
+
+    def optimizeZfac(nk, Htot, dzTop, Huniform, fnPow,  prec):
+        """
+        Optimizes the highest ratio dzFn() so that the sum(dz)=Htot.
+        """
+        it = 0
+        zc0 = 0; dz, H0 = dzFn( nk, Huniform, dzTop, fnPow,  zc0 )
+        zc2 = 2; dz, H2 = dzFn( nk, Huniform, dzTop, fnPow,  zc2 )
+        while H2-H0 > prec/8  and it<200: # Binary search
+            zc1 = (zc0 + zc2)/2; dz, H1 = dzFn( nk, Huniform, dzTop, fnPow,  zc1 )
+            if Htot<H1: zc2, H2 = 1.*zc1, 1.*H1
+            else: zc0, H0 = 1.*zc1, 1.*H1
+            it += 1
+        #print 'zFac=',zc1,'max ratio=',(dz[1:]/dz[:-1]).max()
+        return dz
+
+    def roundDz(nk, Htot, dzTop, Huniform, fnPow,  prec):
+        """
+        Returns dz = optimizeZfac() rounded to the precision "prec", and sum(dz).
+        """
+        dz = prec * numpy.rint( optimizeZfac(nk, Htot, dzTop, Huniform, fnPow,  prec)/prec ) # Round to a given precision
+        return dz, numpy.sum(dz)
+
+    dz, H = roundDz( nk, Htot, dzTop, Huniform, fnPow,  prec )
+    return dz
+
+def zFromDz(dz):
+    """
+    Sums dz to return zInterface and zCenter.
+    """
+    zInt = zeros(dz.shape[0]+1)
+    zInt[1:] = -numpy.cumsum(dz)
+    zCenter = 0.5*( zInt[:-1] + zInt[1:] )
+    return zInt, zCenter
+
+nk = 75         # number of levels
+Htot=6500       # deepest ocean point
+Huniform = 0    # upper region of uniform resolution
+dzTop = 2       # thickness of top level
+fnPow = 1.42865 # ???
+prec = .01      # precision to round thicknesses/depths todz = dzIter(nk, Htot, dzTop, Huniform, fnPow, prec)
+dz = dzIter(nk, Htot, dzTop, Huniform, fnPow,  prec)
+print 'dz=',dz
+print 'sum(dz)=',numpy.sum(dz)
+print 'max ratio=',(dz[1:]/dz[:-1]).max()
diff --git a/ice_ocean_SIS2/FLOR2/INPUT/vgrid_75_2m.nc b/ice_ocean_SIS2/FLOR2/INPUT/vgrid_75_2m.nc
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		/lustre/f1/unswept/Matthew.Harrison/MOM6_datasets/FLOR2/ocn_only/atmos_mosaic_tile1Xland_mosaic_tile1.nc