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

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

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

revision 1.2 by dgoldberg, Tue Sep 4 21:11:44 2012 UTC revision 1.6 by heimbach, Fri Sep 28 12:31:31 2012 UTC
# Line 33  CEOP Line 33  CEOP
33        rhoi = streamice_density        rhoi = streamice_density
34        rhow = streamice_density_ocean_avg        rhow = streamice_density_ocean_avg
35    
36    #ifdef USE_ALT_RLOW
37          print *, "RLOWSI", R_low_si(20,1,1,1)
38    #endif
39    
40        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
41         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
42          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
43           DO i=1-Olx,sNx+Olx           DO i=1-Olx,sNx+Olx
44    #ifdef USE_ALT_RLOW
45              OD = -1.0 * R_low_si (i,j,bi,bj) -
46         &     H_streamice(i,j,bi,bj) * rhoi/rhow
47    #else
48            OD = -1.0 * R_low (i,j,bi,bj) -            OD = -1.0 * R_low (i,j,bi,bj) -
49       &     H_streamice(i,j,bi,bj) * rhoi/rhow       &     H_streamice(i,j,bi,bj) * rhoi/rhow
50    #endif
 !           if (i.eq.1) then  
 !             print *, R_low(1,1,1,1)  
 !           endif  
51    
52    
53            IF (OD .ge. 0. _d 0) THEN                      IF (OD .ge. 0. _d 0) THEN          
54    
55  c         ice thickness does not take up whole ocean column -> floating  c         ice thickness does not take up whole ocean column -> floating
56             float_frac_streamice(i,j,bi,bj) = 0.0             float_frac_streamice(i,j,bi,bj) = 0.0
57    #ifdef USE_ALT_RLOW
58               base_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)+OD
59    #else
60             base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)+OD             base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)+OD
61    #endif
62             surf_el_streamice(i,j,bi,bj) =             surf_el_streamice(i,j,bi,bj) =
63       &      (1-rhoi/rhow)*H_streamice(i,j,bi,bj)       &      (1-rhoi/rhow)*H_streamice(i,j,bi,bj)
64            ELSE            ELSE
65    
66    
67             float_frac_streamice(i,j,bi,bj) = 1.0             float_frac_streamice(i,j,bi,bj) = 1.0
68    #ifdef USE_ALT_RLOW
69               base_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)
70               surf_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)
71    #else
72             base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)             base_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)
73             surf_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)             surf_el_streamice(i,j,bi,bj) = R_low(i,j,bi,bj)
74    #endif
75       &      + H_streamice(i,j,bi,bj)       &      + H_streamice(i,j,bi,bj)
76            ENDIF            ENDIF
77           ENDDO           ENDDO
# Line 65  c         ice thickness does not take up Line 79  c         ice thickness does not take up
79         ENDDO         ENDDO
80        ENDDO        ENDDO
81    
82        _EXCH_XY_RL(float_frac_streamice, myThid )         _EXCH_XY_RL(float_frac_streamice, myThid )
83        _EXCH_XY_RL(base_el_streamice, myThid )         _EXCH_XY_RL(base_el_streamice, myThid )
84        _EXCH_XY_RL(surf_el_streamice, myThid )         _EXCH_XY_RL(surf_el_streamice, myThid )
85    
86                
87    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22