/[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.13 by dgoldberg, Sat Apr 6 17:43:41 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_geom_file_setup,
56         &     STREAMICE_diagnostic_only,
57  !     &     STREAMICE_construct_matrix,  !     &     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         &     PETSC_PRECOND_TYPE, PETSC_SOLVER_TYPE    
81    
82    
83        NAMELIST /STREAMICE_PARM02/        NAMELIST /STREAMICE_PARM02/
84       &     shelf_max_draft,       &     shelf_max_draft,
# Line 82  C     iUnit      :: Work variable for IO Line 110  C     iUnit      :: Work variable for IO
110       &     min_y_CFBC_WEST, max_y_CFBC_WEST,       &     min_y_CFBC_WEST, max_y_CFBC_WEST,
111       &     min_y_CFBC_EAST, max_y_CFBC_EAST,       &     min_y_CFBC_EAST, max_y_CFBC_EAST,
112       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,       &     flux_bdry_val_SOUTH, flux_bdry_val_NORTH,
113       &     flux_bdry_val_WEST, flux_bdry_val_EAST       &     flux_bdry_val_WEST, flux_bdry_val_EAST,
114         &     STREAMICE_NS_periodic, STREAMICE_EW_periodic
115    
116        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
117    
# Line 90  C--   Default values for STREAMICE Line 119  C--   Default values for STREAMICE
119    
120        streamice_density = 917.        streamice_density = 917.
121        streamice_density_ocean_avg = 1024.        streamice_density_ocean_avg = 1024.
122        A_glen_isothermal = 9.461e-18   !  Pa (-1/3) a        B_glen_isothermal = 9.461e-18   !  Pa (-1/3) a
123        n_glen = 3.        n_glen = 3.
124        eps_glen_min = 1.0e-12        eps_glen_min = 1.0e-12
125          eps_u_min = 1.0e-6
126        C_basal_fric_const = 31.71 ! Pa (m/a)-1n        C_basal_fric_const = 31.71 ! Pa (m/a)-1n
127        n_basal_friction = 1.        n_basal_friction = 1.
128        streamice_vel_update = 169200. ! seconds        streamice_vel_update = 169200. ! seconds
129        streamice_cg_tol = 1e-6        streamice_cg_tol = 1e-6
130        streamice_nonlin_tol = 1e-6        streamice_nonlin_tol = 1e-6
131          streamice_nonlin_tol_fp = 1.e-14
132        streamice_max_cg_iter = 2000        streamice_max_cg_iter = 2000
133        streamice_max_nl_iter = 100        streamice_max_nl_iter = 100
134        streamice_n_sub_regularize = 4        streamice_n_sub_regularize = 4
135        streamice_CFL_factor = .5        streamice_CFL_factor = .5
136                streamice_adjDump = 0.
137          streamice_bg_surf_slope_x = .0
138          streamice_bg_surf_slope_y = 0.
139          streamice_kx_b_init = 1.
140          streamice_ky_b_init = 1.
141          streamice_wgt_drift = 0.
142          streamice_wgt_tikh = 0.
143          streamice_wgt_surf = 0.    
144          streamice_wgt_vel = 0.    
145          streamice_wgt_avthick = 0.
146          streamice_addl_backstress = 0.0
147          streamice_smooth_gl_width = 0.0
148    
149        STREAMICEthickInit = 'FILE'        STREAMICEthickInit = 'FILE'
150        STREAMICEthickFile = ' '        STREAMICEthickFile = ' '
151        STREAMICEcalveMaskFile = ' '        STREAMICEcalveMaskFile = ' '
152          STREAMICEsigcoordInit = 'UNIFORM'
153          STREAMICEsigcoordFile = ' '
154          STREAMICEbasalTracConfig = 'UNIFORM'
155          STREAMICEbasalTracFile = ' '
156          STREAMICEvelOptimFile = ''
157          STREAMICEtopogFile = ''
158          STREAMICEhmaskFile = ''
159          STREAMICEuFaceBdryFile = ''
160          STREAMICEvFaceBdryFile = ''
161          STREAMICEuDirichValsFile = ''
162          STREAMICEvDirichValsFile = ''
163          STREAMICEGlenConstFile = ''
164          STREAMICEcostMaskFile = ''
165          STREAMICEGlenConstConfig = 'UNIFORM'
166          PETSC_PRECOND_TYPE = 'PCBJACOBI'
167          PETSC_SOLVER_TYPE = 'KSPCG'
168    
169        STREAMICEison = .TRUE.        STREAMICEison = .TRUE.
170        STREAMICE_tave_mdsio = .TRUE.        STREAMICE_tave_mdsio = .TRUE.
# Line 115  C--   Default values for STREAMICE Line 174  C--   Default values for STREAMICE
174        STREAMICE_GL_regularize = .FALSE.        STREAMICE_GL_regularize = .FALSE.
175        STREAMICE_move_front = .FALSE.        STREAMICE_move_front = .FALSE.
176        STREAMICE_calve_to_mask = .FALSE.        STREAMICE_calve_to_mask = .FALSE.
177    !       STREAMICE_geom_file_setup = .FALSE.
178  !      STREAMICE_construct_matrix = .TRUE.  !      STREAMICE_construct_matrix = .TRUE.
179        STREAMICE_lower_cg_tol = .FALSE.        STREAMICE_lower_cg_tol = .FALSE.
180          STREAMICE_diagnostic_only = .FALSE.
181          STREAMICE_ppm_driving_stress = .FALSE.
182          STREAMICE_h_ctrl_const_surf = .FALSE.
183    !       STREAMICE_hybrid_stress= .FALSE.
184    
185        min_x_noflow_NORTH = 0.        min_x_noflow_NORTH = 0.
186        max_x_noflow_NORTH = 0.        max_x_noflow_NORTH = 0.
# Line 159  C--   Default values for STREAMICE Line 223  C--   Default values for STREAMICE
223        flux_bdry_val_WEST = 0.        flux_bdry_val_WEST = 0.
224        flux_bdry_val_EAST = 0.        flux_bdry_val_EAST = 0.
225    
226          STREAMICE_NS_periodic = .FALSE.
227          STREAMICE_EW_periodic = .FALSE.
228    
229        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'        WRITE(msgBuf,'(A)') 'STREAMICE_READPARMS: opening data.streamice'
230        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 174  C     Read parameters from open data fil Line 240  C     Read parameters from open data fil
240       &    'STREAMICE_READPARMS: read first param block'       &    'STREAMICE_READPARMS: read first param block'
241        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
242       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
       READ(UNIT=iUnit,NML=STREAMICE_PARM02)  
243    
244        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN        IF (TRIM(STREAMICEthickInit) .eq. "PARAM") THEN
245           READ(UNIT=iUnit,NML=STREAMICE_PARM02)
246          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
247       &    'STREAMICE_READPARMS: read second param block'       &    'STREAMICE_READPARMS: read second param block'
248         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
249       &    SQUEEZE_RIGHT , 1)       &    SQUEEZE_RIGHT , 1)
        READ(UNIT=iUnit,NML=STREAMICE_PARM03)  
250        ENDIF        ENDIF
251    
252          READ(UNIT=iUnit,NML=STREAMICE_PARM03)
253        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
254       &    'STREAMICE_READPARMS: read third param block'       &    'STREAMICE_READPARMS: read third param block'
255        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22