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

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

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

revision 1.8 by dgoldberg, Wed Oct 10 15:03:10 2012 UTC revision 1.9 by dgoldberg, Wed Nov 21 01:04:21 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    #ifdef ALLOW_COST
4  #include "COST_OPTIONS.h"  # include "COST_OPTIONS.h"
5    #endif
6  #include "STREAMICE_OPTIONS.h"  #include "STREAMICE_OPTIONS.h"
7    
8  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 38  C     I,J,bi,bj - Loop counters Line 39  C     I,J,bi,bj - Loop counters
39        INTEGER col_y, col_x        INTEGER col_y, col_x
40        _RL slope_pos, c1, x, y, lenx, leny        _RL slope_pos, c1, x, y, lenx, leny
41        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
42          _RS     dummyRS
43    
44  CEOP  CEOP
45    
46  C     ZERO OUT FLOATING POINT ARRAYS  C     ZERO OUT FLOATING POINT ARRAYS
# Line 225  C The 0 is the "iteration" argument. The Line 228  C The 0 is the "iteration" argument. The
228                area_shelf_streamice(i,j,bi,bj) = 0. _d 0                area_shelf_streamice(i,j,bi,bj) = 0. _d 0
229                STREAMICE_hmask(i,j,bi,bj) = 0. _d 0                STREAMICE_hmask(i,j,bi,bj) = 0. _d 0
230               ENDIF               ENDIF
231                 Do k=1,Nr
232                 STREAMICE_ctrl_mask(i,j,bi,bj,k) = 1. _d 0
233                 enddo
234              ENDIF              ENDIF
235             ENDDO             ENDDO
236            ENDDO            ENDDO
# Line 331  C The 0 is the "iteration" argument. The Line 337  C The 0 is the "iteration" argument. The
337        _EXCH_XY_RL(R_low_si, myThid )        _EXCH_XY_RL(R_low_si, myThid )
338  #endif  #endif
339    
340  #ifdef STREAMICE_HYBRID_STRESS  !#ifdef STREAMICE_HYBRID_STRESS
341    
342    !      CALL STREAMICE_VISC_BETA (myThid)
343    
344        CALL STREAMICE_VISC_BETA (myThid)  ! DNG THIS CALL IS TO INITIALISE VISCOSITY
345    !     TO AVOID POSSIBLE ADJOINT INSTABILITIES
346    !     IT IS WRITTEN OVER IN FIRST TIMESTEP
347    
348    #ifdef ALLOW_AUTODIFF
349    
350           CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
351           CALL STREAMICE_VELMASK_UPD (myThid)
352           CALL STREAMICE_VEL_SOLVE( myThid )
353    
354  #endif  #endif
355    
356    !#endif
357                
358        CALL WRITE_FLD_XY_RL ( "C_basal_fric", "",        CALL WRITE_FLD_XY_RL ( "C_basal_fric", "",
359       & C_basal_friction, 0, myThid )       & C_basal_friction, 0, myThid )
# Line 345  C The 0 is the "iteration" argument. The Line 363  C The 0 is the "iteration" argument. The
363       & area_shelf_streamice, 0, myThid )       & area_shelf_streamice, 0, myThid )
364        CALL WRITE_FLD_XY_RL ( "STREAMICE_hmask", "init",        CALL WRITE_FLD_XY_RL ( "STREAMICE_hmask", "init",
365       & STREAMICE_hmask, 0, myThid )       & STREAMICE_hmask, 0, myThid )
366    #ifdef ALLOW_CTRL
367          CALL ACTIVE_WRITE_GEN_RS( 'maskCtrlst', STREAMICE_ctrl_mask,
368         &  'XY', Nr, 1, .FALSE., 0, mythid, dummyRS )
369    #endif
370    !      call active_write_xyz( 'maskCtrlS', STREAMICE_ctrl_mask, 1, 0,
371    !     & mythid, dummy)
372  !       CALL STREAMICE_VELMASK_UPD (myThid)  !       CALL STREAMICE_VELMASK_UPD (myThid)
373  !       CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )  !       CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
374  !       CALL STREAMICE_VEL_SOLVE( myThid )  !       CALL STREAMICE_VEL_SOLVE( myThid )
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22