/[MITgcm]/MITgcm/pkg/thsice/thsice_solve4temp.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_solve4temp.F

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

revision 1.24 by heimbach, Sat Oct 16 19:22:34 2010 UTC revision 1.25 by heimbach, Thu Oct 21 02:10:34 2010 UTC
# Line 193  c     act4 = ikey_dynamics - 1 Line 193  c     act4 = ikey_dynamics - 1
193    
194  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
195  CADJ STORE flxsw(:,:) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE flxsw(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
196        DO j = jMin, jMax        DO j = 1-OLy, sNy+OLy
197         DO i = iMin, iMax         DO i = 1-OLx, sNx+OLx
198          tic1(i,j) = 0. _d 0          tic1(i,j) = 0. _d 0
199          tic2(i,j) = 0. _d 0          tic2(i,j) = 0. _d 0
200         END DO         END DO
# Line 223  c    &         + sNx*sNy*max1*max2*act3 Line 223  c    &         + sNx*sNy*max1*max2*act3
223  c    &         + sNx*sNy*max1*max2*max3*act4  c    &         + sNx*sNy*max1*max2*max3*act4
224  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
225  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
226  cCADJ STORE  devdt        = comlev1_thsice_1, key=ikey_1  cCADJ STORE  devdt(i,j)        = comlev1_thsice_1, key=ikey_1
227  cCADJ STORE  dFlxdT        = comlev1_thsice_1, key=ikey_1  cCADJ STORE  dFlxdT(i,j)        = comlev1_thsice_1, key=ikey_1
 cCADJ STORE  flxexceptsw  = comlev1_thsice_1, key=ikey_1  
228  cCADJ STORE  flxsw(i,j)   = comlev1_thsice_1, key=ikey_1  cCADJ STORE  flxsw(i,j)   = comlev1_thsice_1, key=ikey_1
229  cCADJ STORE  qic1(i,j)    = comlev1_thsice_1, key=ikey_1  cCADJ STORE  qic1(i,j)    = comlev1_thsice_1, key=ikey_1
230  cCADJ STORE  qic2(i,j)    = comlev1_thsice_1, key=ikey_1  cCADJ STORE  qic2(i,j)    = comlev1_thsice_1, key=ikey_1
# Line 331  C--   Compute coefficients used in quadr Line 330  C--   Compute coefficients used in quadr
330  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|
331    
332  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
333  CADJ STORE devdt    = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE devdt(:,:)    = comlev1_bibj,key=iicekey,byte=isbyte
334  CADJ STORE tsf      = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE tsf(:,:)      = comlev1_bibj,key=iicekey,byte=isbyte
335  #endif  #endif
336    
337  C--   Get surface fluxes over melting surface  C--   Get surface fluxes over melting surface
# Line 380  C     Tsfc converges. Line 379  C     Tsfc converges.
379        ENDIF        ENDIF
380    
381  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
382  CADJ STORE devdt    = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE devdt(:,:)    = comlev1_bibj,key=iicekey,byte=isbyte
383  CADJ STORE dflxdt   = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE dflxdt(:,:)   = comlev1_bibj,key=iicekey,byte=isbyte
384  CADJ STORE flx0exsw = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE flx0exsw(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
385  CADJ STORE flxtexsw = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE flxtexsw(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
386  #endif  #endif
387    
388  C ----- begin iteration  -----  C ----- begin iteration  -----
# Line 423  C-    could add this "ifdef" to hide THS Line 422  C-    could add this "ifdef" to hide THS
422         ENDIF         ENDIF
423    
424  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
425  CADJ STORE devdt(i,j) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE devdt(:,:)    = comlev1_bibj,key=iicekey,byte=isbyte
426  CADJ STORE dflxdt(i,j) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE dflxdt(:,:)   = comlev1_bibj,key=iicekey,byte=isbyte
427  CADJ STORE flxtexsw(i,j) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE flxtexsw(:,:) = comlev1_bibj,key=iicekey,byte=isbyte
428  CADJ STORE iceflag(i,j) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE iceflag(:,:)  = comlev1_bibj,key=iicekey,byte=isbyte
429  CADJ STORE tsf(i,j) = comlev1_bibj,key=iicekey,byte=isbyte  CADJ STORE tsf(:,:)      = comlev1_bibj,key=iicekey,byte=isbyte
430  #endif  #endif
431    
432  C--   Compute new top layer and surface temperatures.  C--   Compute new top layer and surface temperatures.
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22