/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/STREAMICE.h

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

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.4 by heimbach,
Fri Jul 27 21:07:13 2012 UTC
+revision 1.13 by dgoldberg,
Thu Oct  4 15:40:16 2012 UTC
 Line 1
- C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
+ C---+----1--+-+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
  #ifdef ALLOW_STREAMICE
 Line 9 
 C     -------------------------- REAL PA
       & A_glen_isothermal, n_glen, eps_glen_min,
       & C_basal_fric_const, n_basal_friction, streamice_input_flux_unif,
       & streamice_vel_update, streamice_cg_tol, streamice_nonlin_tol,
-      & streamice_CFL_factor, streamice_adjDump
+      & streamice_nonlin_tol_fp,
+      & 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
-Line 17 
 C     -------------------------- REAL PA
+Line 20 
 C     -------------------------- REAL PA
        _RL streamice_input_flux_unif
        _RL streamice_vel_update
        _RL streamice_cg_tol, streamice_nonlin_tol
+       _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
-Line 61 
 C     -------------------------- CHAR PA
+Line 67 
 C     -------------------------- CHAR PA
        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 ---------------------------------------------------
-Line 77 
 C     -------------------------- LOGICAL
+Line 95 
 C     -------------------------- LOGICAL
        LOGICAL STREAMICE_calve_to_mask
        LOGICAL STREAMICE_construct_matrix
        LOGICAL STREAMICE_lower_cg_tol
+       LOGICAL STREAMICE_diagnostic_only
+       LOGICAL STREAMICE_ppm_driving_stress
+       LOGICAL STREAMICE_h_ctrl_const_surf
+ 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,
-Line 84 
 C     -------------------------- LOGICAL
+Line 115 
 C     -------------------------- LOGICAL
       & 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,
+      & STREAMICE_ppm_driving_stress,
+      & STREAMICE_h_ctrl_const_surf
  C     -------------------------- AND NOW ARRAYS ---------------------------------------------------
-Line 173 
 C    REAL ARRAYS
+Line 208 
 C    REAL ARRAYS
       &     u_bdry_values_SI,
       &     v_bdry_values_SI,
       &     STREAMICE_dummy_array,
-      &     C_basal_friction
+      &     C_basal_friction,
+      &     A_glen,
+      &     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,
+ #ifdef USE_ALT_RLOW
+      &     R_low_si,
+ #endif
+      &     streamice_sigma_coord, streamice_delsigma,
+      &     H_streamice_prev
        _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)
-Line 190 
 C    REAL ARRAYS
+Line 238 
 C    REAL ARRAYS
        _RL u_bdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
        _RL v_bdry_values_SI    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
        _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)
+ #ifdef USE_ALT_RLOW
+       _RL R_low_si    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+ #endif
+ 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)
+       _RL H_streamice_prev (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
        _RL STREAMICE_dummy_array (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+       COMMON /STREAMICE_COST_RL/
+      &       cost_func1_streamice
+       _RL cost_func1_streamice(nSx,nSy)
  C    NOTES :
  C      REAL ARRAYS THAT COMPRISE "STATE":
  C       H_streamice,

 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.13
 Legend:



Removed from v.1.4
 


changed lines


 
Added in v.1.13
-Removed from v.1.4
+Added in v.1.13

	ViewVC Help
Powered by ViewVC 1.1.22