/[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.67 by gforget, Wed Jun 11 18:29:08 2008 UTC revision 1.68 by jmc, Mon Aug 11 22:25:52 2008 UTC
# Line 362  C         slope terms (e.g. GM/Redi tens Line 362  C         slope terms (e.g. GM/Redi tens
362       &         .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN       &         .OR. useSALT_PLUME .OR. doDiagsRho.GE.1 ) THEN
363  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
364              IF ( debugLevel .GE. debLevB )              IF ( debugLevel .GE. debLevB )
365       &       CALL DEBUG_CALL('FIND_RHO',myThid)       &       CALL DEBUG_CALL('FIND_RHO_2D',myThid)
366  #endif  #endif
367  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
368  CADJ STORE theta(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE theta(:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
369  CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE salt (:,:,k,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
370  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
371              CALL FIND_RHO(              CALL FIND_RHO_2D(
372       I        bi, bj, iMin, iMax, jMin, jMax, k, k,       I                 iMin, iMax, jMin, jMax, k,
373       I        theta, salt,       I                 theta(1-OLx,1-OLy,k,bi,bj),
374       O        rhoK,       I                 salt (1-OLx,1-OLy,k,bi,bj),
375       I        myThid )       O                 rhoK,
376         I                 k, bi, bj, myThid )
377    
378              IF (k.GT.1) THEN              IF (k.GT.1) THEN
379  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
380  CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
381  CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte  CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_bibj_k, key=kkey, byte=isbyte
382  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
383               CALL FIND_RHO(               CALL FIND_RHO_2D(
384       I        bi, bj, iMin, iMax, jMin, jMax, k-1, k,       I                 iMin, iMax, jMin, jMax, k,
385       I        theta, salt,       I                 theta(1-OLx,1-OLy,k-1,bi,bj),
386       O        rhoKm1,       I                 salt (1-OLx,1-OLy,k-1,bi,bj),
387       I        myThid )       O                 rhoKm1,
388         I                 k-1, bi, bj, myThid )
389              ENDIF              ENDIF
390  #ifdef ALLOW_DEBUG  #ifdef ALLOW_DEBUG
391              IF ( debugLevel .GE. debLevB )              IF ( debugLevel .GE. debLevB )
Legend:
Removed from v.1.67  
changed lines
  Added in v.1.68
  ViewVC Help
Powered by ViewVC 1.1.22