/[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.1 by heimbach, Thu Mar 29 15:59:21 2012 UTC revision 1.12 by dgoldberg, Thu Mar 7 15:23:19 2013 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,       &     B_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         &     STREAMICEcostMaskFile,
50       &     STREAMICEison,       &     STREAMICEison,
51       &     STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,       &     STREAMICE_dump_mdsio, STREAMICE_tave_mdsio,
52       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,       &     STREAMICE_dump_mnc, STREAMICE_tave_mnc,
53       &     STREAMICE_GL_regularize, STREAMICE_move_front,       &     STREAMICE_GL_regularize, STREAMICE_move_front,
54       &     STREAMICE_calve_to_mask,       &     STREAMICE_calve_to_mask,
55       &     STREAMICE_construct_matrix,  !      &     STREAMICE_geom_file_setup,
56         &     STREAMICE_diagnostic_only,
57    !     &     STREAMICE_construct_matrix,
58       &     STREAMICE_lower_cg_tol,       &     STREAMICE_lower_cg_tol,
59       &     streamice_CFL_factor       &     streamice_CFL_factor,
60         &     streamice_adjDump,
61    !      &     STREAMICE_hybrid_stress,
62         &     streamice_bg_surf_slope_x, streamice_bg_surf_slope_y,
63         &     streamice_kx_b_init, streamice_ky_b_init,
64         &     STREAMICEbasalTracConfig,
65         &     STREAMICEbasalTracFile,
66         &     STREAMICEvelOptimFile,
67         &     STREAMICEtopogFile,
68         &     STREAMICEhmaskFile,
69         &     STREAMICEuFaceBdryFile,
70         &     STREAMICEvFaceBdryFile,
71         &     STREAMICEuDirichValsFile,
72         &     STREAMICEvDirichValsFile,
73         &     STREAMICEGlenConstFile, STREAMICEGlenConstConfig,
74         &     STREAMICE_ppm_driving_stress,
75         &     STREAMICE_h_ctrl_const_surf,
76         &     streamice_wgt_drift,streamice_wgt_surf,streamice_wgt_vel,
77         &     streamice_wgt_avthick, streamice_wgt_tikh,
78         &     streamice_addl_backstress,
79         &     streamice_smooth_gl_width    
80    
81    
82        NAMELIST /STREAMICE_PARM02/        NAMELIST /STREAMICE_PARM02/
83       &     shelf_max_draft,       &     shelf_max_draft,
# Line 82  C     iUnit      :: Work variable for IO Line 109  C     iUnit      :: Work variable for IO
109       &     min_y_CFBC_WEST, max_y_CFBC_WEST,       &     min_y_CFBC_WEST, max_y_CFBC_WEST,
110       &     min_y_CFBC_EAST, max_y_CFBC_EAST,       &     min_y_CFBC_EAST, max_y_CFBC_EAST,
111       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
112       &     flux_bdry_val_WEST, flux_bdry_val_EAST       &     flux_bdry_val_WEST, flux_bdry_val_EAST,
113         &     STREAMICE_NS_periodic, STREAMICE_EW_periodic
114    
115        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
116    
# Line 90  C--   Default values for STREAMICE Line 118  C--   Default values for STREAMICE
118    
119        streamice_density = 917.        streamice_density = 917.
120        streamice_density_ocean_avg = 1024.        streamice_density_ocean_avg = 1024.
121        A_glen_isothermal = 9.461e-18   !  Pa (-1/3) a        B_glen_isothermal = 9.461e-18   !  Pa (-1/3) a
122        n_glen = 3.        n_glen = 3.
123        eps_glen_min = 1.0e-12        eps_glen_min = 1.0e-12
124          eps_u_min = 1.0e-6
125        C_basal_fric_const = 31.71 ! Pa (m/a)-1n        C_basal_fric_const = 31.71 ! Pa (m/a)-1n
126        n_basal_friction = 1.        n_basal_friction = 1.
127        streamice_vel_update = 169200. ! seconds        streamice_vel_update = 169200. ! seconds
128        streamice_cg_tol = 1e-6        streamice_cg_tol = 1e-6
129        streamice_nonlin_tol = 1e-6        streamice_nonlin_tol = 1e-6
130          streamice_nonlin_tol_fp = 1.e-14
131        streamice_max_cg_iter = 2000        streamice_max_cg_iter = 2000
132        streamice_max_nl_iter = 100        streamice_max_nl_iter = 100
133        streamice_n_sub_regularize = 4        streamice_n_sub_regularize = 4
134        streamice_CFL_factor = .5        streamice_CFL_factor = .5
135                streamice_adjDump = 0.
136          streamice_bg_surf_slope_x = .0
137          streamice_bg_surf_slope_y = 0.
138          streamice_kx_b_init = 1.
139          streamice_ky_b_init = 1.
140          streamice_wgt_drift = 0.
141          streamice_wgt_tikh = 0.
142          streamice_wgt_surf = 0.    
143          streamice_wgt_vel = 0.    
144          streamice_wgt_avthick = 0.
145          streamice_addl_backstress = 0.0
146          streamice_smooth_gl_width = 0.0
147    
148        STREAMICEthickInit = 'FILE'        STREAMICEthickInit = 'FILE'
149        STREAMICEthickFile = ' '        STREAMICEthickFile = ' '
150        STREAMICEcalveMaskFile = ' '        STREAMICEcalveMaskFile = ' '
151          STREAMICEsigcoordInit = 'UNIFORM'
152          STREAMICEsigcoordFile = ' '
153          STREAMICEbasalTracConfig = 'UNIFORM'
154          STREAMICEbasalTracFile = ' '
155          STREAMICEvelOptimFile = ''
156          STREAMICEtopogFile = ''
157          STREAMICEhmaskFile = ''
158          STREAMICEuFaceBdryFile = ''
159          STREAMICEvFaceBdryFile = ''
160          STREAMICEuDirichValsFile = ''
161          STREAMICEvDirichValsFile = ''
162          STREAMICEGlenConstFile = ''
163          STREAMICEcostMaskFile = ''
164          STREAMICEGlenConstConfig = 'UNIFORM'
165    
166        STREAMICEison = .TRUE.        STREAMICEison = .TRUE.
167        STREAMICE_tave_mdsio = .TRUE.        STREAMICE_tave_mdsio = .TRUE.
# Line 115  C--   Default values for STREAMICE Line 171  C--   Default values for STREAMICE
171        STREAMICE_GL_regularize = .FALSE.        STREAMICE_GL_regularize = .FALSE.
172        STREAMICE_move_front = .FALSE.        STREAMICE_move_front = .FALSE.
173        STREAMICE_calve_to_mask = .FALSE.        STREAMICE_calve_to_mask = .FALSE.
174        STREAMICE_construct_matrix = .TRUE.  !       STREAMICE_geom_file_setup = .FALSE.
175    !      STREAMICE_construct_matrix = .TRUE.
176        STREAMICE_lower_cg_tol = .FALSE.        STREAMICE_lower_cg_tol = .FALSE.
177          STREAMICE_diagnostic_only = .FALSE.
178          STREAMICE_ppm_driving_stress = .FALSE.
179          STREAMICE_h_ctrl_const_surf = .FALSE.
180    !       STREAMICE_hybrid_stress= .FALSE.
181    
182        min_x_noflow_NORTH = 0.        min_x_noflow_NORTH = 0.
183        max_x_noflow_NORTH = 0.        max_x_noflow_NORTH = 0.
# Line 159  C--   Default values for STREAMICE Line 220  C--   Default values for STREAMICE
220        flux_bdry_val_WEST = 0.        flux_bdry_val_WEST = 0.
221        flux_bdry_val_EAST = 0.        flux_bdry_val_EAST = 0.
222    
223          STREAMICE_NS_periodic = .FALSE.
224          STREAMICE_EW_periodic = .FALSE.
225    
226        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
227        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 174  C     Read parameters from open data fil Line 237  C     Read parameters from open data fil
237       &    'STREAMICE_READPARMS: read first param block'       &    'STREAMICE_READPARMS: read first param block'
238        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
239       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
       READ(UNIT=iUnit,NML=STREAMICE_PARM02)  
240    
241        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
242           READ(UNIT=iUnit,NML=STREAMICE_PARM02)
243          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
244       &    'STREAMICE_READPARMS: read second param block'       &    'STREAMICE_READPARMS: read second param block'
245         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
246       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
        READ(UNIT=iUnit,NML=STREAMICE_PARM03)  
247        ENDIF        ENDIF
248    
249          READ(UNIT=iUnit,NML=STREAMICE_PARM03)
250        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
251       &    'STREAMICE_READPARMS: read third param block'       &    'STREAMICE_READPARMS: read third param block'
252        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22