/[MITgcm]/MITgcm/pkg/seaice/seaice_model.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_model.F

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

revision 1.2 by heimbach, Tue Nov 12 20:47:27 2002 UTC revision 1.3 by dimitri, Thu Dec 5 08:43:03 2002 UTC
# Line 21  C     === Global variables === Line 21  C     === Global variables ===
21  #include "SEAICE_PARAMS.h"  #include "SEAICE_PARAMS.h"
22  #include "SEAICE_EXTERNAL.h"  #include "SEAICE_EXTERNAL.h"
23    
24    #ifdef ALLOW_AUTODIFF_TAMC
25    # include "tamc.h"
26    #endif /* ALLOW_AUTODIFF_TAMC */
27    
28  C     === Routine arguments ===  C     === Routine arguments ===
29  C     myTime - Simulation time  C     myTime - Simulation time
30  C     myIter - Simulation timestep number  C     myIter - Simulation timestep number
# Line 45  C--   Third level model velocity is used Line 49  C--   Third level model velocity is used
49         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
50          DO j=0,sNy+1          DO j=0,sNy+1
51           DO i=0,sNx+1           DO i=0,sNx+1
52            GWATX(I,J,bi,bj)=0.5*(uVel(i+1,j,3,bi,bj)            GWATX(I,J,bi,bj)=HALF*(uVel(i+1,j,3,bi,bj)
53       &                         +uVel(i+1,j+1,3,bi,bj))       &                         +uVel(i+1,j+1,3,bi,bj))
54            GWATY(I,J,bi,bj)=0.5*(vVel(i,j+1,3,bi,bj)            GWATY(I,J,bi,bj)=HALF*(vVel(i,j+1,3,bi,bj)
55       &                         +vVel(i+1,j+1,3,bi,bj))       &                         +vVel(i+1,j+1,3,bi,bj))
56  #ifdef SEAICE_DEBUG  #ifdef SEAICE_DEBUG
57  c          write(*,'(2i4,2i2,f7.1,7f12.3)')  c          write(*,'(2i4,2i2,f7.1,7f12.3)')
# Line 95  C--   Write sea ice restart files Line 99  C--   Write sea ice restart files
99        CALL SEAICE_WRITE_PICKUP ( .FALSE.,        CALL SEAICE_WRITE_PICKUP ( .FALSE.,
100       &     myTime+deltaTClock, myIter+1, myThid )       &     myTime+deltaTClock, myIter+1, myThid )
101    
102    C---------------------------------------------------
103    C OOH NOOOO we need to move the whole stuff
104    C---------------------------------------------------
105    #ifdef ALLOW_AUTODIFF_TAMC
106    CRG CADJ store UICE,VICE,AREA,HEFF,fu,fv,EmPmR,Qnet,Qsw = comlev1_bibj
107    #endif /* ALLOW_AUTODIFF_TAMC */
108    
109    C--   Call sea-ice cost function routine
110    CRG      CALL SEAICE_COST( myTime, myIter, myThid )
111    
112  #endif ALLOW_SEAICE  #endif ALLOW_SEAICE
113    
114        RETURN        RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22