/[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.11 by dgoldberg, Sun Jan 27 21:05:50 2013 UTC revision 1.12 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC
# Line 45  CEOP Line 45  CEOP
45    
46  C     ZERO OUT FLOATING POINT ARRAYS  C     ZERO OUT FLOATING POINT ARRAYS
47    
48    
49        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
50         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
51          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
# Line 61  C     ZERO OUT FLOATING POINT ARRAYS Line 62  C     ZERO OUT FLOATING POINT ARRAYS
62            mass_ice_streamice(i,j,bi,bj) = 0. _d 0            mass_ice_streamice(i,j,bi,bj) = 0. _d 0
63            BDOT_streamice(i,j,bi,bj) = 0. _d 0            BDOT_streamice(i,j,bi,bj) = 0. _d 0
64            ADOT_streamice(i,j,bi,bj) = 0. _d 0            ADOT_streamice(i,j,bi,bj) = 0. _d 0
65  !           C_basal_friction(i,j,bi,bj) = C_basal_fric_const            C_basal_friction(i,j,bi,bj) = C_basal_fric_const
66  !           A_glen(i,j,bi,bj) = A_glen_isothermal            B_glen(i,j,bi,bj) = B_glen_isothermal
67            H_streamice_prev(i,j,bi,bj) = 0. _d 0            H_streamice_prev(i,j,bi,bj) = 0. _d 0
68  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
69            ru_old_si(i,j,bi,bj) = 0. _d 0            ru_old_si(i,j,bi,bj) = 0. _d 0
# Line 132  C     INIT. INTEGER ARRAYS Line 133  C     INIT. INTEGER ARRAYS
133    
134  #ifdef USE_ALT_RLOW  #ifdef USE_ALT_RLOW
135  ! init alternate array for topog  ! init alternate array for topog
136        IF ( bathyFile .NE. ' ' ) THEN        IF ( STREAMICEtopogFile .NE. ' ' ) THEN
137          _BARRIER          _BARRIER
138  C The 0 is the "iteration" argument. The ' ' is an empty suffix  C The 0 is the "iteration" argument. The ' ' is an empty suffix
139         CALL READ_FLD_XY_RS( bathyFile, '',         CALL READ_FLD_XY_RS( STREAMICEtopogFile, '',
140       &      R_low_si, 0, myThid )       &      R_low_si, 0, myThid )
141                
142        ELSE        ELSE
# Line 232  C The 0 is the "iteration" argument. The Line 233  C The 0 is the "iteration" argument. The
233                STREAMICE_hmask(i,j,bi,bj) = 0. _d 0                STREAMICE_hmask(i,j,bi,bj) = 0. _d 0
234               ENDIF               ENDIF
235               Do k=1,Nr               Do k=1,Nr
236               STREAMICE_ctrl_mask(i,j,bi,bj,k) = 1. _d 0               STREAMICE_ctrl_mask(i,j,k,bi,bj) = 1. _d 0
237               enddo               enddo
238              ENDIF              ENDIF
239             ENDDO             ENDDO
# Line 289  C The 0 is the "iteration" argument. The Line 290  C The 0 is the "iteration" argument. The
290          _BARRIER          _BARRIER
291    
292          CALL READ_FLD_XY_RL( STREAMICEGlenConstFile, ' ',          CALL READ_FLD_XY_RL( STREAMICEGlenConstFile, ' ',
293       &      A_glen, 0, myThid )       &      B_glen, 0, myThid )
294    
295         ELSE         ELSE
296          WRITE(msgBuf,'(A)') 'INIT GLEN - FILENAME MISSING'          WRITE(msgBuf,'(A)') 'INIT GLEN - FILENAME MISSING'
# Line 303  C The 0 is the "iteration" argument. The Line 304  C The 0 is the "iteration" argument. The
304           DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
305            DO j=1,sNy            DO j=1,sNy
306             DO i=1,sNx             DO i=1,sNx
307              A_glen(i,j,bi,bj) = A_glen_isothermal              B_glen(i,j,bi,bj) = B_glen_isothermal
308             ENDDO             ENDDO
309            ENDDO            ENDDO
310           ENDDO           ENDDO
# Line 400  C The 0 is the "iteration" argument. The Line 401  C The 0 is the "iteration" argument. The
401        _EXCH_XY_RL(STREAMICE_hmask, myThid )        _EXCH_XY_RL(STREAMICE_hmask, myThid )
402        _EXCH_XY_RL(area_shelf_streamice, myThid )        _EXCH_XY_RL(area_shelf_streamice, myThid )
403        _EXCH_XY_RL(C_basal_friction, myThid )        _EXCH_XY_RL(C_basal_friction, myThid )
404          _EXCH_XY_RL(B_glen, myThid )
405  #ifdef USE_ALT_RLOW  #ifdef USE_ALT_RLOW
406        _EXCH_XY_RL(R_low_si, myThid )        _EXCH_XY_RL(R_low_si, myThid )
407  #endif  #endif
# Line 424  C The 0 is the "iteration" argument. The Line 426  C The 0 is the "iteration" argument. The
426                
427        CALL WRITE_FLD_XY_RL ( "C_basal_fric", "",        CALL WRITE_FLD_XY_RL ( "C_basal_fric", "",
428       & C_basal_friction, 0, myThid )       & C_basal_friction, 0, myThid )
429          CALL WRITE_FLD_XY_RL ( "B_glen_sqrt", "",
430         & B_glen, 0, myThid )
431        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",        CALL WRITE_FLD_XY_RL ( "H_streamice", "init",
432       & H_streamIce, 0, myThid )       & H_streamIce, 0, myThid )
433        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",        CALL WRITE_FLD_XY_RL ( "area_shelf_streamice", "init",
# Line 452  C The 0 is the "iteration" argument. The Line 456  C The 0 is the "iteration" argument. The
456  !       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)
457  !       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)
458  !       CALL WRITE_FULLARRAY_RL ("umask",STREAMICE_umask,1,0,0,1,0,myThid)  !       CALL WRITE_FULLARRAY_RL ("umask",STREAMICE_umask,1,0,0,1,0,myThid)
459    
460    
461                
462  #endif /* ALLOW_STREAMICE */  #endif /* ALLOW_STREAMICE */
463    
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22