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

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

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

revision 1.4 by dgoldberg, Mon Feb 11 00:57:31 2013 UTC revision 1.5 by dgoldberg, Thu Mar 7 15:23:19 2013 UTC
# Line 89  C       The 0 is the "iteration" argumen Line 89  C       The 0 is the "iteration" argumen
89       &                    SQUEEZE_RIGHT , 1)       &                    SQUEEZE_RIGHT , 1)
90        ENDIF        ENDIF
91    
92          
93    #ifdef ALLOW_COST
94          IF ( STREAMICEcostMaskFile .NE. ' ') THEN
95            _BARRIER
96    C       The 0 is the "iteration" argument. The ' ' is an empty suffix
97            CALL READ_FLD_XY_RL( STREAMICEcostMaskFile, ' ',
98         &   STREAMICE_cost_mask, 0, myThid )
99           ELSE
100            WRITE(msgBuf,'(A)') 'COST MASK - NOT IMPLENTED'
101            CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
102         &                    SQUEEZE_RIGHT , 1)
103    
104           DO bj = myByLo(myThid), myByHi(myThid)
105            DO bi = myBxLo(myThid), myBxHi(myThid)
106             DO j=1,sNy
107              DO i=1,sNx
108               STREAMICE_cost_mask (i,j,bi,bj) = 1.0
109              ENDDO
110             ENDDO
111            ENDDO
112           ENDDO
113    
114          ENDIF
115    #endif
116          
117    
118  !    READ IN FILES FOR BOUNDARY MASKS AND DIRICH VALUES  !    READ IN FILES FOR BOUNDARY MASKS AND DIRICH VALUES
119  #ifdef STREAMICE_GEOM_FILE_SETUP  #ifdef STREAMICE_GEOM_FILE_SETUP
120    
# Line 408  C         INITIALIZE BOUNDARY CONDS AT N Line 434  C         INITIALIZE BOUNDARY CONDS AT N
434       &  'XY', Nr, 1, .FALSE., 0, mythid, dummyRS )       &  'XY', Nr, 1, .FALSE., 0, mythid, dummyRS )
435  #endif  #endif
436    
437    #ifdef ALLOW_COST
438    !      _EXCH_XY_RL(STREAMICE_ctrl_mask, myThid )
439          CALL WRITE_FLD_XY_RS ( 'maskCost', '',
440         & STREAMICE_cost_mask, 0, myThid )
441    #endif
442    
443    
444        _EXCH_XY_RL(k1AtC_str, myThid )        _EXCH_XY_RL(k1AtC_str, myThid )
445        _EXCH_XY_RL(k2AtC_str, myThid )        _EXCH_XY_RL(k2AtC_str, myThid )
# Line 422  C         INITIALIZE BOUNDARY CONDS AT N Line 454  C         INITIALIZE BOUNDARY CONDS AT N
454        Xquad (2) = .5 * (1.+1./sqrt(3.))        Xquad (2) = .5 * (1.+1./sqrt(3.))
455    
456        CALL STREAMICE_INIT_PHI( myThid )        CALL STREAMICE_INIT_PHI( myThid )
       PRINT *, "GOT HERE END INIT FIXED"  
457                
458    
459  #endif  #endif
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22