/[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.6 by heimbach, Thu Sep 20 02:04:45 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        CALL STREAMICE_VELMASK_UPD (myThid)        CALL STREAMICE_VELMASK_UPD (myThid)
152        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
153        CALL STREAMICE_VEL_SOLVE( myThid )        CALL STREAMICE_VEL_SOLVE( myThid )
154  !       endif  !       endif
155    
156          
157          if(.not.STREAMICE_diagnostic_only) THEN
158    
159           CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT )  
160    
161          endif
162          
163    
164  !       enddo  !       enddo
165    
166    !       CALL STREAMICE_VELMASK_UPD (myThid)
167    !       CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
168    !       CALL STREAMICE_VEL_SOLVE( myThid )
169    
170  !       modelEnd = myTime.EQ.endTime .OR. myIter.EQ.nEndIter  !       modelEnd = myTime.EQ.endTime .OR. myIter.EQ.nEndIter
171    
# Line 173  c************************************** Line 184  c**************************************
184  !  !
185  !  !
186  !       IF ( useDiagnostics ) THEN  !       IF ( useDiagnostics ) THEN
 !       PRINT *, "CALLING FILL DIAGNOSTICS"  
187  !        CALL DIAGNOSTICS_FILL(U_streamice,'SI_Uvel ',  !        CALL DIAGNOSTICS_FILL(U_streamice,'SI_Uvel ',
188  !      &      0,1,0,1,1,myThid)  !      &      0,1,0,1,1,myThid)
189  !        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.6
  ViewVC Help
Powered by ViewVC 1.1.22