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

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

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

revision 1.1 by heimbach, Thu Mar 29 15:59:21 2012 UTC revision 1.7 by dgoldberg, Thu Sep 27 20:29:01 2012 UTC
# Line 67  c************************************** Line 67  c**************************************
67            zu_old_si(i,j,bi,bj) = 0. _d 0            zu_old_si(i,j,bi,bj) = 0. _d 0
68            zv_old_si(i,j,bi,bj) = 0. _d 0            zv_old_si(i,j,bi,bj) = 0. _d 0
69            h_after_uflux_si(i,j,bi,bj) = 0. _d 0            h_after_uflux_si(i,j,bi,bj) = 0. _d 0
70    #ifdef STREAMICE_HYBRID_STRESS
71              streamice_taubx (i,j,bi,bj) = 0. _d 0
72              streamice_tauby (i,j,bi,bj) = 0. _d 0
73    #endif
74           ENDDO           ENDDO
75          ENDDO          ENDDO
76         ENDDO         ENDDO
# Line 74  c************************************** Line 78  c**************************************
78    
79  #endif  #endif
80    
 !       PRINT *, "ICE STREAM TIMESTEP ", iLoop, "; ",  
 !      &  Time/86400/365, " years"  
81    
82  !       do while (time_step_remain .gt. 0. _d 0)  !       do while (time_step_remain .gt. 0. _d 0)
83  !  !
# Line 133  c************************************** Line 135  c**************************************
135    
136                
137    
       CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT )  
138    
139          
140    
141  !       if (time_step_int .gt. 1000) then  !       if (time_step_int .gt. 1000) then
142    
143        WRITE(msgBuf,'(A,I10.10,E9.2,A)')        WRITE(msgBuf,'(A,I10.10,E9.2,A)')
144       &  'streamice solo_time_step: nIter',       &  'streamice solo_time_step: nIter',
145       &  myIter, myTime/86400.0/365.0, 'seconds'       &  myIter, myTime/86400.0/365.0, 'seconds'
   
146        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
147       &                   SQUEEZE_RIGHT , 1)       &                   SQUEEZE_RIGHT , 1)
148    
149          CALL STREAMICE_DUMP( mytime, myiter, myThid )
150    
151    
152    
153        CALL STREAMICE_VELMASK_UPD (myThid)        CALL STREAMICE_VELMASK_UPD (myThid)
154        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
155        CALL STREAMICE_VEL_SOLVE( myThid )        CALL STREAMICE_VEL_SOLVE( myThid )
156  !       endif  !       endif
157    
158          
159          if(.not.STREAMICE_diagnostic_only) THEN
160    
161           CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT )  
162    
163          endif
164          
165    
166  !       enddo  !       enddo
167    
168    !       CALL STREAMICE_VELMASK_UPD (myThid)
169    !       CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
170    !       CALL STREAMICE_VEL_SOLVE( myThid )
171    
172  !       modelEnd = myTime.EQ.endTime .OR. myIter.EQ.nEndIter  !       modelEnd = myTime.EQ.endTime .OR. myIter.EQ.nEndIter
173    
# Line 173  c************************************** Line 186  c**************************************
186  !  !
187  !  !
188  !       IF ( useDiagnostics ) THEN  !       IF ( useDiagnostics ) THEN
 !       PRINT *, "CALLING FILL DIAGNOSTICS"  
189  !        CALL DIAGNOSTICS_FILL(U_streamice,'SI_Uvel ',  !        CALL DIAGNOSTICS_FILL(U_streamice,'SI_Uvel ',
190  !      &      0,1,0,1,1,myThid)  !      &      0,1,0,1,1,myThid)
191  !        CALL DIAGNOSTICS_FILL(V_streamice,'SI_Vvel ',  !        CALL DIAGNOSTICS_FILL(V_streamice,'SI_Vvel ',
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22