319 |
|
|
320 |
! finish initialize glen constant |
! finish initialize glen constant |
321 |
|
|
322 |
|
! initialize BDOT |
323 |
|
|
324 |
|
IF ( STREAMICEBdotConfig.EQ.'FILE' ) THEN |
325 |
|
|
326 |
|
IF ( STREAMICEBdotFile .NE. ' ' ) THEN |
327 |
|
_BARRIER |
328 |
|
|
329 |
|
CALL READ_FLD_XY_RL( STREAMICEBdotFile, ' ', |
330 |
|
& BDOT_STREAMICE, 0, myThid ) |
331 |
|
|
332 |
|
ELSE |
333 |
|
WRITE(msgBuf,'(A)') 'INIT GLEN - FILENAME MISSING' |
334 |
|
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit, |
335 |
|
& SQUEEZE_RIGHT , 1) |
336 |
|
ENDIF |
337 |
|
|
338 |
|
! ELSE IF (STREAMICEBdotConfig.EQ.'UNIFORM' ) THEN |
339 |
|
! |
340 |
|
! DO bj = myByLo(myThid), myByHi(myThid) |
341 |
|
! DO bi = myBxLo(myThid), myBxHi(myThid) |
342 |
|
! DO j=1,sNy |
343 |
|
! DO i=1,sNx |
344 |
|
! BDOT_streamice(i,j,bi,bj) = B_glen_isothermal |
345 |
|
! ENDDO |
346 |
|
! ENDDO |
347 |
|
! ENDDO |
348 |
|
! ENDDO |
349 |
|
|
350 |
|
ELSEIF (STREAMICEBdotConfig.EQ.' ' ) THEN |
351 |
|
|
352 |
|
DO bj = myByLo(myThid), myByHi(myThid) |
353 |
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
354 |
|
DO j=1,sNy |
355 |
|
DO i=1,sNx |
356 |
|
BDOT_streamice(i,j,bi,bj) = 0. _d 0 |
357 |
|
ENDDO |
358 |
|
ENDDO |
359 |
|
ENDDO |
360 |
|
ENDDO |
361 |
|
|
362 |
|
ELSE |
363 |
|
|
364 |
|
WRITE(msgBuf,'(A)') 'INIT MELT RATE - NOT IMPLENTED' |
365 |
|
CALL PRINT_ERROR( msgBuf, myThid) |
366 |
|
STOP 'ABNORMAL END: S/R STREAMICE_INIT_VAR' |
367 |
|
ENDIF |
368 |
|
|
369 |
|
! finish initialize BDOT |
370 |
|
|
371 |
|
|
372 |
! initialize basal traction |
! initialize basal traction |
373 |
|
|
374 |
IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN |
IF ( STREAMICEbasalTracConfig.EQ.'FILE' ) THEN |