/[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.4 by dgoldberg, Tue Sep 25 03:05:53 2012 UTC revision 1.5 by dgoldberg, Thu Sep 27 20:29:01 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          print *, "RLOWSI", R_low_si(20,1,1,1)
37    
38        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
39         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
40          DO j=1-Oly,sNy+Oly          DO j=1-Oly,sNy+Oly
41           DO i=1-Olx,sNx+Olx           DO i=1-Olx,sNx+Olx
42    #ifdef USE_ALT_RLOW
43              OD = -1.0 * R_low_si (i,j,bi,bj) -
44         &     H_streamice(i,j,bi,bj) * rhoi/rhow
45    #else
46            OD = -1.0 * R_low (i,j,bi,bj) -            OD = -1.0 * R_low (i,j,bi,bj) -
47       &     H_streamice(i,j,bi,bj) * rhoi/rhow       &     H_streamice(i,j,bi,bj) * rhoi/rhow
48    #endif
 !           if (i.eq.1) then  
 !             print *, R_low(1,1,1,1)  
 !           endif  
49    
50    
51            IF (OD .ge. 0. _d 0) THEN                      IF (OD .ge. 0. _d 0) THEN          
52    
53  c         ice thickness does not take up whole ocean column -> floating  c         ice thickness does not take up whole ocean column -> floating
54             float_frac_streamice(i,j,bi,bj) = 0.0             float_frac_streamice(i,j,bi,bj) = 0.0
55    #ifdef USE_ALT_RLOW
56               base_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)+OD
57    #else
58             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
59    #endif
60             surf_el_streamice(i,j,bi,bj) =             surf_el_streamice(i,j,bi,bj) =
61       &      (1-rhoi/rhow)*H_streamice(i,j,bi,bj)       &      (1-rhoi/rhow)*H_streamice(i,j,bi,bj)
62            ELSE            ELSE
63    
64    
65             float_frac_streamice(i,j,bi,bj) = 1.0             float_frac_streamice(i,j,bi,bj) = 1.0
66    #ifdef USE_ALT_RLOW
67               base_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)
68               surf_el_streamice(i,j,bi,bj) = R_low_si(i,j,bi,bj)
69    #else
70             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)
71             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)
72    #endif
73       &      + H_streamice(i,j,bi,bj)       &      + H_streamice(i,j,bi,bj)
74            ENDIF            ENDIF
75           ENDDO           ENDDO

Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5

  ViewVC Help
Powered by ViewVC 1.1.22