/[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.5 by heimbach, Thu Sep 20 02:04:45 2012 UTC revision 1.6 by dgoldberg, Thu Sep 27 20:29:00 2012 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "COST_OPTIONS.h"
5  #include "STREAMICE_OPTIONS.h"  #include "STREAMICE_OPTIONS.h"
6    
7  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 65  C     ZERO OUT FLOATING POINT ARRAYS Line 66  C     ZERO OUT FLOATING POINT ARRAYS
66            zv_old_si(i,j,bi,bj) = 0. _d 0            zv_old_si(i,j,bi,bj) = 0. _d 0
67            h_after_uflux_SI(i,j,bi,bj) = 0. _d 0            h_after_uflux_SI(i,j,bi,bj) = 0. _d 0
68  #endif  #endif
69    #ifdef USE_ALT_RLOW
70              R_low_si(i,j,bi,bj) = 0. _d 0
71    #endif
72    
73  #ifdef STREAMICE_HYBRID_STRESS  #ifdef STREAMICE_HYBRID_STRESS
74            do k=1,Nr            do k=1,Nr
75             visc_streamice_full(i,j,k,bi,bj) =             visc_streamice_full(i,j,k,bi,bj) =
# Line 75  C     ZERO OUT FLOATING POINT ARRAYS Line 80  C     ZERO OUT FLOATING POINT ARRAYS
80  #endif  #endif
81           ENDDO           ENDDO
82          ENDDO          ENDDO
83    
84    #ifdef ALLOW_COST_TEST
85            cost_func1_streamice (bi,bj) = 0.0
86    #endif
87    
88         ENDDO         ENDDO
89        ENDDO        ENDDO
90    
# Line 114  C     INIT. INTEGER ARRAYS Line 124  C     INIT. INTEGER ARRAYS
124         ENDDO         ENDDO
125        ENDDO        ENDDO
126    
127    
128    #ifdef USE_ALT_RLOW
129    ! init alternate array for topog
130          IF ( bathyFile .NE. ' ' ) THEN
131            _BARRIER
132    C The 0 is the "iteration" argument. The ' ' is an empty suffix
133           CALL READ_FLD_XY_RS( bathyFile, '',
134         &      R_low_si, 0, myThid )
135          
136          ELSE
137            WRITE(msgBuf,'(A)') 'STREAMICE TOPOG - FILENAME MISSING'
138            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
139         &                    SQUEEZE_RIGHT , 1)
140          ENDIF
141    #endif
142    
143  ! initialize thickness  ! initialize thickness
144            
145        IF ( STREAMICEthickInit.EQ.'PARAM' ) THEN        IF ( STREAMICEthickInit.EQ.'PARAM' ) THEN
# Line 300  C The 0 is the "iteration" argument. The Line 326  C The 0 is the "iteration" argument. The
326        _EXCH_XY_RL(STREAMICE_hmask, myThid )        _EXCH_XY_RL(STREAMICE_hmask, myThid )
327        _EXCH_XY_RL(area_shelf_streamice, myThid )        _EXCH_XY_RL(area_shelf_streamice, myThid )
328        _EXCH_XY_RL(C_basal_friction, myThid )        _EXCH_XY_RL(C_basal_friction, myThid )
329    #ifdef USE_ALT_RLOW
330          _EXCH_XY_RL(R_low_si, myThid )
331    #endif
332    
333  #ifdef STREAMICE_HYBRID_STRESS  #ifdef STREAMICE_HYBRID_STRESS
334    
335        CALL STREAMICE_VISC_BETA (myThid)        CALL STREAMICE_VISC_BETA (myThid)
336    
337  #endif  #endif
338          
339          CALL WRITE_FLD_XY_RL ( "C_basal_fric", "",
340         & C_basal_friction, 0, myThid )
341        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",
342       & H_streamIce, 0, myThid )       & H_streamIce, 0, myThid )
343        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",
# Line 323  C The 0 is the "iteration" argument. The Line 353  C The 0 is the "iteration" argument. The
353       &   C_basal_friction, 0, myThid )       &   C_basal_friction, 0, myThid )
354        CALL WRITE_FLD_XY_RL ( "V_init", "",        CALL WRITE_FLD_XY_RL ( "V_init", "",
355       &   V_streamice, 0, myThid )       &   V_streamice, 0, myThid )
356    #ifdef USE_ALT_RLOW
357          CALL WRITE_FLD_XY_RL ( "R_low_si", "init",
358         & R_low_si, 0, myThid )
359    #endif
360    
361  !       CALL WRITE_FULLARRAY_RL ("H",H_streamice,1,0,0,1,0,myThid)  !       CALL WRITE_FULLARRAY_RL ("H",H_streamice,1,0,0,1,0,myThid)
362  !       CALL WRITE_FULLARRAY_RL ("hmask",STREAMICE_hmask,1,0,0,1,0,myThid)  !       CALL WRITE_FULLARRAY_RL ("hmask",STREAMICE_hmask,1,0,0,1,0,myThid)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22