Diff of /MITgcm_contrib/dgoldberg/streamice/STREAMICE.h

Parent Directory | Revision Log | Revision Graph | Patch

--- MITgcm_contrib/dgoldberg/streamice/STREAMICE.h	2012/09/04 21:11:44	1.6
+++ MITgcm_contrib/dgoldberg/streamice/STREAMICE.h	2012/09/18 17:06:48	1.7
@@ -10,7 +10,9 @@
      & C_basal_fric_const, n_basal_friction, streamice_input_flux_unif,
      & streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol,
      & streamice_nonlin_tol_fp,
-     & streamice_CFL_factor, streamice_adjDump
+     & streamice_CFL_factor, streamice_adjDump, 
+     & streamice_bg_surf_slope_x, streamice_bg_surf_slope_y,
+     & streamice_kx_b_init, streamice_ky_b_init
       _RL streamice_density, streamice_density_ocean_avg
       _RL A_glen_isothermal, n_glen, eps_glen_min
       _RL C_basal_fric_const
@@ -21,6 +23,8 @@
       _RL streamice_nonlin_tol_fp
       _RL streamice_CFL_factor
       _RL streamice_adjDump
+      _RL streamice_bg_surf_slope_x, streamice_bg_surf_slope_y
+      _RL streamice_kx_b_init, streamice_ky_b_init
       
 C     parms for parameterized initial thickness
 C     SHELF_MAX_DRAFT: max thickness of ice in m
@@ -63,10 +67,22 @@
       CHARACTER*(MAX_LEN_FNAM) STREAMICEthickFile
       CHARACTER*(MAX_LEN_FNAM) STREAMICEthickInit
       CHARACTER*(MAX_LEN_FNAM) STREAMICEcalveMaskFile
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordInit
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEsigcoordFile
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEdelsigFile
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracConfig
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEbasalTracFile
+      CHARACTER*(MAX_LEN_FNAM) STREAMICEvelOptimFile
       COMMON /STREAMICE_PARM_C/ 
      &     STREAMICEthickInit,
      &     STREAMICEthickFile,
-     &     STREAMICEcalveMaskFile
+     &     STREAMICEcalveMaskFile,
+     &     STREAMICEsigcoordInit,
+     &     STREAMICEsigcoordFile,
+     &     STREAMICEdelsigFile,
+     &     STREAMICEbasalTracConfig,
+     &     STREAMICEbasalTracFile,
+     &     STREAMICEvelOptimFile
      
 C     -------------------------- LOGICAL PARAMS ---------------------------------------------------
 
@@ -79,6 +95,17 @@
       LOGICAL STREAMICE_calve_to_mask
       LOGICAL STREAMICE_construct_matrix
       LOGICAL STREAMICE_lower_cg_tol
+      LOGICAL STREAMICE_diagnostic_only
+      
+C     The following parameters specify periodic boundary conditions. 
+C     For now this will completely override all other boundary conditions
+C     and apply to the entire boundary
+
+      LOGICAL STREAMICE_NS_periodic
+      LOGICAL STREAMICE_EW_periodic
+      
+C      LOGICAL STREAMICE_hybrid_stress
+      
       COMMON /STREAMICE_PARM_L/
      & STREAMICEison,
      & STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
@@ -86,7 +113,9 @@
      & STREAMICE_GL_regularize, STREAMICE_move_front,
      & STREAMICE_calve_to_mask,
      & STREAMICE_construct_matrix,
-     & STREAMICE_lower_cg_tol
+     & STREAMICE_lower_cg_tol,
+     & STREAMICE_NS_periodic, STREAMICE_EW_periodic,
+     & STREAMICE_diagnostic_only
 
 C     -------------------------- AND NOW ARRAYS ---------------------------------------------------
 
@@ -177,7 +206,12 @@
      &     STREAMICE_dummy_array,
      &     C_basal_friction,
      &     A_glen,
-     &     BDOT_streamice
+     &     BDOT_streamice,
+     &     streamice_taubx, streamice_tauby,
+     &     streamice_u_surf, streamice_v_surf,
+     &     visc_streamice_full, streamice_omega, streamice_basal_geom,
+     &     streamice_vert_shear_uz, streamice_vert_shear_vz,
+     &     streamice_sigma_coord, streamice_delsigma
       _RL H_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL U_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL V_streamice           (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
@@ -196,6 +230,22 @@
       _RL C_basal_friction    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL A_glen    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       
+C      The following arrays are used for the hybrid stress balance      
+      
+      _RL streamice_taubx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL streamice_tauby (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL streamice_u_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL streamice_v_surf (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL streamice_omega (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL streamice_basal_geom 
+     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+      _RL visc_streamice_full 
+     & (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
+      _RL streamice_vert_shear_uz (Nr)
+      _RL streamice_vert_shear_vz (Nr)
+      _RL streamice_sigma_coord (Nr)
+      _RL streamice_delsigma (Nr)
+      
 !! IMPORTANT: MELT RATE IN METERS PER YEAR
 !! POSITIVE WHERE MELTING
       _RL BDOT_streamice (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)