/[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.7 by dgoldberg, Thu Sep 27 20:29:01 2012 UTC revision 1.8 by dgoldberg, Mon Oct 8 02:12:47 2012 UTC
# Line 79  c************************************** Line 79  c**************************************
79  #endif  #endif
80    
81    
 !       do while (time_step_remain .gt. 0. _d 0)  
 !  
 !        n_interm = n_interm + 1  
 !        Imin_ratio = 1.0e-16  
 !  
 !        DO bj = myByLo(myThid), myByHi(myThid)  
 !         DO bi = myBxLo(myThid), myBxHi(myThid)  
 !          DO j=1,sNy  
 !           DO i=1,sNx  
 !              
 !            local_u_max = 0  
 !            local_v_max = 0  
 !              
 !            if (STREAMICE_hmask (i,j,bi,bj) .eq. 1.0) then  
 !             do ki = 0,1  
 !              do kj = 0,1  
 !               local_u_max =  
 !      &         max (local_u_max, abs(U_streamice(i+ki,j+kj,bi,bj)))  
 !               local_v_max =  
 !      &         max (local_v_max, abs(V_streamice(i+ki,j+kj,bi,bj)))  
 !              enddo  
 !             enddo  
 !  
 !             Iratio = max (local_u_max/dxF(i,j,bi,bj) ,  
 !      &       local_v_max/dyF(i,j,bi,bj))  
 !             Imin_ratio = max (Imin_ratio, Iratio)  
 !            endif  
 !    
 !           enddo  
 !          enddo  
 !         enddo  
 !        enddo  
 !  
 !        CALL GLOBAL_MAX_R8 (Imin_ratio, myThid)  
 !        min_ratio = 1./(Imin_ratio+1.e-12)  
 !        time_step_int =  
 !      &  min(streamice_CFL_factor*min_ratio*(365*86400),  
 !      &   deltaT)  
 !          
 !        if (time_step_int .lt. min_time_step) then  
 !         WRITE(msgBuf,'(A,A,E9.2)') 'streamice solo_time_step:',  
 !      &                       ' abnormally small timestep ',  
 !      &                       streamice_CFL_factor*min_ratio*(365*86400)  
 !         CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
 !      &                    SQUEEZE_RIGHT , 1)  
 !        endif  
 !          
 !        if (time_step_int .ge. time_step_remain) then  
 !         time_step_int = time_step_remain  
 !         time_step_remain = 0.0  
 !        else  
 !         time_step_remain = time_step_remain - time_step_int  
 !        endif  
   
         
   
   
   
   
 !       if (time_step_int .gt. 1000) then  
82    
83        WRITE(msgBuf,'(A,I10.10,E9.2,A)')        WRITE(msgBuf,'(A,I10.10,E9.2,A)')
84       &  'streamice solo_time_step: nIter',       &  'streamice solo_time_step: nIter',
# Line 148  c************************************** Line 88  c**************************************
88    
89        CALL STREAMICE_DUMP( mytime, myiter, myThid )        CALL STREAMICE_DUMP( mytime, myiter, myThid )
90    
91    ! NEW DIRECTIVES - DNG
92    #ifdef ALLOW_AUTODIFF_TAMC
93    CADJ STORE float_frac_streamice  = comlev1, key = ikey_dynamics,
94    CADJ &     kind = isbyte
95    CADJ STORE surf_el_streamice  = comlev1, key = ikey_dynamics,
96    CADJ &     kind = isbyte
97    CADJ STORE base_el_streamice  = comlev1, key = ikey_dynamics,
98    CADJ &     kind = isbyte
99    #endif
100    ! NEW DIRECTIVES - DNG
101    
102          if (myIter.eq.0) then
103           CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
104          endif
105    
       CALL STREAMICE_VELMASK_UPD (myThid)  
106        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )        CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
107          CALL STREAMICE_VELMASK_UPD (myThid)
108        CALL STREAMICE_VEL_SOLVE( myThid )        CALL STREAMICE_VEL_SOLVE( myThid )
109  !       endif        
   
   
110        if(.not.STREAMICE_diagnostic_only) THEN        if(.not.STREAMICE_diagnostic_only) THEN
111    
112         CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT )           CALL STREAMICE_ADVECT_THICKNESS ( myThid, deltaT )  
113    
114        endif        endif
115                
116    ! CALL AT END INSTEAD OF BEGINNING - DNG
117          CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )
118    
 !       enddo  
   
 !       CALL STREAMICE_VELMASK_UPD (myThid)  
 !       CALL STREAMICE_UPD_FFRAC_UNCOUPLED ( myThid )  
 !       CALL STREAMICE_VEL_SOLVE( myThid )  
   
 !       modelEnd = myTime.EQ.endTime .OR. myIter.EQ.nEndIter  
   
 ! #ifdef ALLOW_DIAGNOSTICS  
 !  
 ! ! this stuff was to replace what is in forward_step, but not  
 ! ! forward_step will be called  
 !    
 ! !       IF ( useDiagnostics ) THEN  
 ! !         CALL DIAGNOSTICS_SWITCH_ONOFF( myTime, myiter, myThid )  
 ! ! ! C--   State-variables diagnostics  
 ! ! !         CALL TIMER_START('DO_STATEVARS_DIAGS  [FORWARD_STEP]',myThid)  
 ! ! !         CALL DO_STATEVARS_DIAGS( myTime, 0, myIter, myThid )  
 ! ! !         CALL TIMER_STOP ('DO_STATEVARS_DIAGS  [FORWARD_STEP]',myThid)  
 ! !       ENDIF  
 !  
 !  
 !       IF ( useDiagnostics ) THEN  
 !        CALL DIAGNOSTICS_FILL(U_streamice,'SI_Uvel ',  
 !      &      0,1,0,1,1,myThid)  
 !        CALL DIAGNOSTICS_FILL(V_streamice,'SI_Vvel ',  
 !      &      0,1,0,1,1,myThid)  
 !        CALL DIAGNOSTICS_FILL(H_streamice,'SI_Thick',  
 !      &      0,1,0,1,1,myThid)  
 !        CALL DIAGNOSTICS_FILL(area_shelf_streamice,'SI_area ',  
 !      &      0,1,0,1,1,myThid)  
 !        CALL DIAGNOSTICS_FILL(float_frac_streamice,'SI_float',  
 !      &      0,1,0,1,1,myThid)  
 !  
 !       ENDIF  
 !  
 ! !       IF ( useDiagnostics ) THEN  
 ! !           CALL TIMER_START('DO_STATEVARS_DIAGS  [FORWARD_STEP]',myThid)  
 ! !           CALL DO_STATEVARS_DIAGS( myTime, 2, myIter, myThid )  
 ! !           CALL TIMER_STOP ('DO_STATEVARS_DIAGS  [FORWARD_STEP]',myThid)  
 ! !       ENDIF  
 !  
 ! #endif /* ALLOW_DIAGNOSTICS */  
   
 !       myIter = nIter0 + iLoop  
 !       myTime = startTime + deltaTClock * float(iLoop)  
 !  
 !       CALL TIMER_START('DO_THE_MODEL_IO     [FORWARD_STEP]',myThid)  
 !       CALL DO_THE_MODEL_IO( modelEnd, myTime, myIter, myThid )  
 !       CALL TIMER_STOP ('DO_THE_MODEL_IO     [FORWARD_STEP]',myThid)  
119    
120    
121  #endif  #endif
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22