/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_readparms.F

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

revision 1.2 by dgoldberg, Mon May 14 16:47:27 2012 UTC revision 1.9 by dgoldberg, Mon Dec 10 02:34:45 2012 UTC
# Line 35  C     iUnit      :: Work variable for IO Line 35  C     iUnit      :: Work variable for IO
35    
36        NAMELIST /STREAMICE_PARM01/        NAMELIST /STREAMICE_PARM01/
37       &     streamice_density, streamice_density_ocean_avg,       &     streamice_density, streamice_density_ocean_avg,
38       &     A_glen_isothermal, n_glen, eps_glen_min,       &     A_glen_isothermal, n_glen, eps_glen_min, eps_u_min,
39       &     C_basal_fric_const, n_basal_friction,       &     C_basal_fric_const, n_basal_friction,
40       &     streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,       &     streamice_vel_update,streamice_cg_tol,streamice_nonlin_tol,
41         &     streamice_nonlin_tol_fp,
42       &     streamice_max_cg_iter, streamice_max_nl_iter,       &     streamice_max_cg_iter, streamice_max_nl_iter,
43       &     STREAMICE_GL_regularize,       &     STREAMICE_GL_regularize,
44       &     STREAMICEthickInit,       &     STREAMICEthickInit,
45         &     STREAMICEsigcoordInit,
46         &     STREAMICEsigcoordFile,
47       &     STREAMICEthickFile,       &     STREAMICEthickFile,
48       &     STREAMICEcalveMaskFile,       &     STREAMICEcalveMaskFile,
49       &     STREAMICEison,       &     STREAMICEison,
# Line 48  C     iUnit      :: Work variable for IO Line 51  C     iUnit      :: Work variable for IO
51       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,
52       &     STREAMICE_GL_regularize, STREAMICE_move_front,       &     STREAMICE_GL_regularize, STREAMICE_move_front,
53       &     STREAMICE_calve_to_mask,       &     STREAMICE_calve_to_mask,
54         &     STREAMICE_diagnostic_only,
55  !     &     STREAMICE_construct_matrix,  !     &     STREAMICE_construct_matrix,
56       &     STREAMICE_lower_cg_tol,       &     STREAMICE_lower_cg_tol,
57       &     streamice_CFL_factor       &     streamice_CFL_factor,
58         &     streamice_adjDump,
59    !      &     STREAMICE_hybrid_stress,
60         &     streamice_bg_surf_slope_x, streamice_bg_surf_slope_y,
61         &     streamice_kx_b_init, streamice_ky_b_init,
62         &     STREAMICEbasalTracConfig,
63         &     STREAMICEbasalTracFile,
64         &     STREAMICEvelOptimFile,
65         &     STREAMICE_ppm_driving_stress,
66         &     STREAMICE_h_ctrl_const_surf,
67         &     streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel,
68         &     streamice_wgt_avthick, streamice_addl_backstress    
69    
70    
71        NAMELIST /STREAMICE_PARM02/        NAMELIST /STREAMICE_PARM02/
72       &     shelf_max_draft,       &     shelf_max_draft,
# Line 82  C     iUnit      :: Work variable for IO Line 98  C     iUnit      :: Work variable for IO
98       &     min_y_CFBC_WEST, max_y_CFBC_WEST,       &     min_y_CFBC_WEST, max_y_CFBC_WEST,
99       &     min_y_CFBC_EAST, max_y_CFBC_EAST,       &     min_y_CFBC_EAST, max_y_CFBC_EAST,
100       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
101       &     flux_bdry_val_WEST, flux_bdry_val_EAST       &     flux_bdry_val_WEST, flux_bdry_val_EAST,
102         &     STREAMICE_NS_periodic, STREAMICE_EW_periodic
103    
104        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
105    
# Line 93  C--   Default values for STREAMICE Line 110  C--   Default values for STREAMICE
110        A_glen_isothermal = 9.461e-18   !  Pa (-1/3) a        A_glen_isothermal = 9.461e-18   !  Pa (-1/3) a
111        n_glen = 3.        n_glen = 3.
112        eps_glen_min = 1.0e-12        eps_glen_min = 1.0e-12
113          eps_u_min = 1.0e-6
114        C_basal_fric_const = 31.71 ! Pa (m/a)-1n        C_basal_fric_const = 31.71 ! Pa (m/a)-1n
115        n_basal_friction = 1.        n_basal_friction = 1.
116        streamice_vel_update = 169200. ! seconds        streamice_vel_update = 169200. ! seconds
117        streamice_cg_tol = 1e-6        streamice_cg_tol = 1e-6
118        streamice_nonlin_tol = 1e-6        streamice_nonlin_tol = 1e-6
119          streamice_nonlin_tol_fp = 1.e-14
120        streamice_max_cg_iter = 2000        streamice_max_cg_iter = 2000
121        streamice_max_nl_iter = 100        streamice_max_nl_iter = 100
122        streamice_n_sub_regularize = 4        streamice_n_sub_regularize = 4
123        streamice_CFL_factor = .5        streamice_CFL_factor = .5
124                streamice_adjDump = 0.
125          streamice_bg_surf_slope_x = .0
126          streamice_bg_surf_slope_y = 0.
127          streamice_kx_b_init = 1.
128          streamice_ky_b_init = 1.
129          streamice_wgt_drift = 0.
130          streamice_wgt_surf = 0.    
131          streamice_wgt_vel = 0.    
132          streamice_wgt_avthick = 0.
133          streamice_addl_backstress = 0.
134    
135        STREAMICEthickInit = 'FILE'        STREAMICEthickInit = 'FILE'
136        STREAMICEthickFile = ' '        STREAMICEthickFile = ' '
137        STREAMICEcalveMaskFile = ' '        STREAMICEcalveMaskFile = ' '
138          STREAMICEsigcoordInit = 'UNIFORM'
139          STREAMICEsigcoordFile = ' '
140          STREAMICEbasalTracConfig = 'UNIFORM'
141          STREAMICEbasalTracFile = ' '
142          STREAMICEvelOptimFile = ''
143    
144        STREAMICEison = .TRUE.        STREAMICEison = .TRUE.
145        STREAMICE_tave_mdsio = .TRUE.        STREAMICE_tave_mdsio = .TRUE.
# Line 117  C--   Default values for STREAMICE Line 151  C--   Default values for STREAMICE
151        STREAMICE_calve_to_mask = .FALSE.        STREAMICE_calve_to_mask = .FALSE.
152  !      STREAMICE_construct_matrix = .TRUE.  !      STREAMICE_construct_matrix = .TRUE.
153        STREAMICE_lower_cg_tol = .FALSE.        STREAMICE_lower_cg_tol = .FALSE.
154          STREAMICE_diagnostic_only = .FALSE.
155          STREAMICE_ppm_driving_stress = .FALSE.
156          STREAMICE_h_ctrl_const_surf = .FALSE.
157    !       STREAMICE_hybrid_stress= .FALSE.
158    
159        min_x_noflow_NORTH = 0.        min_x_noflow_NORTH = 0.
160        max_x_noflow_NORTH = 0.        max_x_noflow_NORTH = 0.
# Line 159  C--   Default values for STREAMICE Line 197  C--   Default values for STREAMICE
197        flux_bdry_val_WEST = 0.        flux_bdry_val_WEST = 0.
198        flux_bdry_val_EAST = 0.        flux_bdry_val_EAST = 0.
199    
200          STREAMICE_NS_periodic = .FALSE.
201          STREAMICE_EW_periodic = .FALSE.
202    
203        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
204        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 174  C     Read parameters from open data fil Line 214  C     Read parameters from open data fil
214       &    'STREAMICE_READPARMS: read first param block'       &    'STREAMICE_READPARMS: read first param block'
215        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
216       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
       READ(UNIT=iUnit,NML=STREAMICE_PARM02)  
217    
218        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
219           READ(UNIT=iUnit,NML=STREAMICE_PARM02)
220          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
221       &    'STREAMICE_READPARMS: read second param block'       &    'STREAMICE_READPARMS: read second param block'
222         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
223       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
        READ(UNIT=iUnit,NML=STREAMICE_PARM03)  
224        ENDIF        ENDIF
225    
226          READ(UNIT=iUnit,NML=STREAMICE_PARM03)
227        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
228       &    'STREAMICE_READPARMS: read third param block'       &    'STREAMICE_READPARMS: read third param block'
229        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22