281 |
|
|
282 |
! finish initialize thickness |
! finish initialize thickness |
283 |
|
|
284 |
! initialize basal traction |
! initialize glen constant |
285 |
|
|
286 |
IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN |
IF ( STREAMICEGlenConstConfig.EQ.'FILE' ) THEN |
287 |
|
|
288 |
IF ( STREAMICEbasalTracFile .NE. ' ' ) THEN |
IF ( STREAMICEGlenConstFile .NE. ' ' ) THEN |
289 |
_BARRIER |
_BARRIER |
290 |
C The 0 is the "iteration" argument. The ' ' is an empty suffix |
|
291 |
CALL READ_FLD_XY_RL( STREAMICEbasalTracFile, ' ', |
CALL READ_FLD_XY_RL( STREAMICEGlenConstFile, ' ', |
292 |
& C_basal_friction, 0, myThid ) |
& A_glen, 0, myThid ) |
293 |
|
|
294 |
ELSE |
ELSE |
295 |
WRITE(msgBuf,'(A)') 'INIT THICKNESS - FILENAME MISSING' |
WRITE(msgBuf,'(A)') 'INIT GLEN - FILENAME MISSING' |
296 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
297 |
& SQUEEZE_RIGHT , 1) |
& SQUEEZE_RIGHT , 1) |
298 |
ENDIF |
ENDIF |
299 |
|
|
300 |
ELSE IF (STREAMICEbasalTracConfig.EQ.'UNIFORM' ) THEN |
ELSE IF (STREAMICEGlenConstConfig.EQ.'UNIFORM' ) THEN |
301 |
|
|
302 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
303 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
316 |
STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR' |
STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR' |
317 |
ENDIF |
ENDIF |
318 |
|
|
319 |
! initialize glen constant |
! finish initialize glen constant |
320 |
|
|
321 |
IF ( STREAMICEGlenConstConfig.EQ.'FILE' ) THEN |
! initialize basal traction |
322 |
|
|
323 |
IF ( STREAMICEGlenConstFile .NE. ' ' ) THEN |
IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN |
324 |
|
|
325 |
|
IF ( STREAMICEbasalTracFile .NE. ' ' ) THEN |
326 |
_BARRIER |
_BARRIER |
327 |
C The 0 is the "iteration" argument. The ' ' is an empty suffix |
|
328 |
CALL READ_FLD_XY_RL( STREAMICEGlenConstFile, ' ', |
CALL READ_FLD_XY_RL( STREAMICEbasalTracFile, ' ', |
329 |
& A_glen, 0, myThid ) |
& C_basal_friction, 0, myThid ) |
330 |
|
|
331 |
ELSE |
ELSE |
332 |
WRITE(msgBuf,'(A)') 'INIT THICKNESS - FILENAME MISSING' |
WRITE(msgBuf,'(A)') 'INIT C_BASAL - FILENAME MISSING' |
333 |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
334 |
& SQUEEZE_RIGHT , 1) |
& SQUEEZE_RIGHT , 1) |
335 |
ENDIF |
ENDIF |
392 |
STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR' |
STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR' |
393 |
ENDIF |
ENDIF |
394 |
|
|
395 |
! finish initialize glen |
! finish initialize basal trac |
396 |
|
|
397 |
CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) |
CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid ) |
398 |
|
|