/[MITgcm]/MITgcm/model/src/do_oceanic_phys.F
ViewVC logotype

Diff of /MITgcm/model/src/do_oceanic_phys.F

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

revision 1.56 by heimbach, Wed Sep 26 04:12:40 2007 UTC revision 1.57 by heimbach, Fri Sep 28 00:52:18 2007 UTC
# Line 142  C     to diagnose surface buoyancy fluxe Line 142  C     to diagnose surface buoyancy fluxe
142  C     term.  Because this call precedes model thermodynamics, temperature  C     term.  Because this call precedes model thermodynamics, temperature
143  C     under sea-ice may not be "exactly" at the freezing point by the time  C     under sea-ice may not be "exactly" at the freezing point by the time
144  C     theta is dumped or time-averaged.  C     theta is dumped or time-averaged.
145    #ifndef ALLOW_AUTODIFF_TAMC
146        IF ( useSEAICE ) THEN        IF ( useSEAICE ) THEN
147  #ifdef ALLOW_AUTODIFF_TAMC  #else
148  CADJ STORE atemp,aqh,precip    = comlev1, key = ikey_dynamics  CADJ STORE atemp,aqh,precip    = comlev1, key = ikey_dynamics
149  CADJ STORE swdown,lwdown       = comlev1, key = ikey_dynamics  CADJ STORE swdown,lwdown       = comlev1, key = ikey_dynamics
150  cph# ifdef EXF_READ_EVAP  cph# ifdef EXF_READ_EVAP
# Line 172  CADJ STORE recip_hfacc         = comlev1 Line 173  CADJ STORE recip_hfacc         = comlev1
173          CALL TIMER_START('SEAICE_MODEL    [DO_OCEANIC_PHYS]', myThid)          CALL TIMER_START('SEAICE_MODEL    [DO_OCEANIC_PHYS]', myThid)
174          CALL SEAICE_MODEL( myTime, myIter, myThid )          CALL SEAICE_MODEL( myTime, myIter, myThid )
175          CALL TIMER_STOP ('SEAICE_MODEL    [DO_OCEANIC_PHYS]', myThid)          CALL TIMER_STOP ('SEAICE_MODEL    [DO_OCEANIC_PHYS]', myThid)
176  #ifdef ALLOW_COST_ICE  #ifdef ALLOW_COST
177          CALL COST_ICE_TEST ( myTime, myIter, myThid )          CALL SEAICE_COST_SENSI ( myTime, myIter, myThid )
178  #endif  #endif
179    #ifndef ALLOW_AUTODIFF_TAMC
180        ENDIF        ENDIF
181    #endif
182  #endif /* ALLOW_SEAICE */  #endif /* ALLOW_SEAICE */
183    
184  #if (defined ALLOW_THSICE) && !(defined ALLOW_ATM2D)  #if (defined ALLOW_THSICE) && !(defined ALLOW_ATM2D)
# Line 396  cph Avoid variable aliasing for adjoint Line 399  cph Avoid variable aliasing for adjoint
399       I             myThid )       I             myThid )
400            ENDIF            ENDIF
401    
 #ifdef ALLOW_AUTODIFF_TAMC  
 ctest# ifndef GM_EXCLUDE_CLIPPING  
 CADJ STORE rhoK   (:,:) = comlev1_bibj_k ,       key=kkey, byte=isbyte  
 ctest# endif  
 CADJ STORE rhoKm1 (:,:) = comlev1_bibj_k ,       key=kkey, byte=isbyte  
 #endif /* ALLOW_AUTODIFF_TAMC */  
402  C--       Implicit Vertical Diffusion for Convection  C--       Implicit Vertical Diffusion for Convection
403  c ==> should use sigmaR !!!  c ==> should use sigmaR !!!
404            IF (k.GT.1 .AND. ivdc_kappa.NE.0.) THEN            IF (k.GT.1 .AND. ivdc_kappa.NE.0.) THEN
# Line 426  c ==> should use sigmaR !!! Line 423  c ==> should use sigmaR !!!
423  C--     end of diagnostic k loop (Nr:1)  C--     end of diagnostic k loop (Nr:1)
424          ENDDO          ENDDO
425    
426    #ifdef ALLOW_AUTODIFF_TAMC
427    CADJ STORE IVDConvCount(:,:,:,bi,bj)
428    CADJ &     = comlev1_bibj, key=itdkey, byte=isbyte
429    #endif
430    
431  C--     Diagnose Mixed Layer Depth:  C--     Diagnose Mixed Layer Depth:
432          IF ( useGMRedi .OR. doDiagsRho.GE.1 ) THEN          IF ( useGMRedi .OR. doDiagsRho.GE.1 ) THEN
433            CALL CALC_OCE_MXLAYER( rhoK, sigmaR,            CALL CALC_OCE_MXLAYER( rhoK, sigmaR,
Legend:
Removed from v.1.56  
changed lines
  Added in v.1.57
  ViewVC Help
Powered by ViewVC 1.1.22