/[MITgcm]/MITgcm/pkg/kpp/kpp_routines.F
ViewVC logotype

Diff of /MITgcm/pkg/kpp/kpp_routines.F

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

revision 1.37 by atn, Fri Dec 21 02:54:34 2007 UTC revision 1.38 by jmc, Mon Aug 11 22:28:06 2008 UTC
# Line 1502  c calculate density, alpha, beta in surf Line 1502  c calculate density, alpha, beta in surf
1502  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
1503  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
1504  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */
1505        call FIND_RHO(        CALL FIND_RHO_2D(
1506       I     bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, 1, 1,       I     1-OLx, sNx+OLx, 1-OLy, sNy+OLy, 1,
1507       I     theta, salt,       I     theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
1508       O     WORK1,       O     WORK1,
1509       I     myThid )       I     1, bi, bj, myThid )
1510  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE
1511  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
1512  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
# Line 1552  CHPF$  INDEPENDENT, NEW (RHOK,RHOKM1,RHO Line 1552  CHPF$  INDEPENDENT, NEW (RHOK,RHOKM1,RHO
1552  CADJ STORE theta(:,:,k,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE theta(:,:,k,bi,bj) = comlev1_kpp_k, key=kkppkey
1553  CADJ STORE salt (:,:,k,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE salt (:,:,k,bi,bj) = comlev1_kpp_k, key=kkppkey
1554  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */
1555           call FIND_RHO(           CALL FIND_RHO_2D(
1556       I        bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, K, K,       I        1-OLx, sNx+OLx, 1-OLy, sNy+OLy, k,
1557       I        theta, salt,       I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
1558       O        RHOK,       O        RHOK,
1559       I        myThid )       I        k, bi, bj, myThid )
1560    
1561  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE
1562  CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE theta(:,:,k-1,bi,bj) = comlev1_kpp_k, key=kkppkey
1563  CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE salt (:,:,k-1,bi,bj) = comlev1_kpp_k, key=kkppkey
1564  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */
1565           call FIND_RHO(           CALL FIND_RHO_2D(
1566       I        bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, K-1, K,       I        1-OLx, sNx+OLx, 1-OLy, sNy+OLy, k,
1567       I        theta, salt,       I        theta(1-OLx,1-OLy,k-1,bi,bj),salt(1-OLx,1-OLy,k-1,bi,bj),
1568       O        RHOKM1,       O        RHOKM1,
1569       I        myThid )       I        k-1, bi, bj, myThid )
1570    
1571  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE
1572  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE theta(:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
1573  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey  CADJ STORE salt (:,:,1,bi,bj) = comlev1_kpp_k, key=kkppkey
1574  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */  #endif /* KPP_AUTODIFF_EXCESSIVE_STORE */
1575           call FIND_RHO(           CALL FIND_RHO_2D(
1576       I        bi, bj, 1-OLx, sNx+OLx, 1-OLy, sNy+OLy, 1, K,       I        1-OLx, sNx+OLx, 1-OLy, sNy+OLy, k,
1577       I        theta, salt,       I        theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
1578       O        RHO1K,       O        RHO1K,
1579       I        myThid )       I        1, bi, bj, myThid )
1580    
1581  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE  #ifdef KPP_AUTODIFF_EXCESSIVE_STORE
1582  CADJ STORE rhok  (:,:)          = comlev1_kpp_k, key=kkppkey  CADJ STORE rhok  (:,:)          = comlev1_kpp_k, key=kkppkey
# Line 1608  CADJ STORE rho1k (:,:)          = comlev Line 1608  CADJ STORE rho1k (:,:)          = comlev
1608  c     work1 - density of t(k-1)  & s(k-1) at depth 1  c     work1 - density of t(k-1)  & s(k-1) at depth 1
1609  c     work2 - density of t(k  )  & s(k  ) at depth 1  c     work2 - density of t(k  )  & s(k  ) at depth 1
1610  c     work3 - density of t(1)-.8 & s(1  ) at depth 1  c     work3 - density of t(1)-.8 & s(1  ) at depth 1
1611              call FIND_RHO(              CALL FIND_RHO_2D(
1612       I           bi, bj, 1, sNx, 1, sNy, K-1, 1, theta, salt,       I           1, sNx, 1, sNy, 1,
1613         I           theta(1-OLx,1-OLy,k-1,bi,bj),
1614         I           salt (1-OLx,1-OLy,k-1,bi,bj),
1615       O           WORK1,       O           WORK1,
1616       I           myThid )       I           k-1, bi, bj, myThid )
1617              call FIND_RHO(              CALL FIND_RHO_2D(
1618       I           bi, bj, 1, sNx, 1, sNy, K  , 1, theta, salt,       I           1, sNx, 1, sNy, 1,
1619         I           theta(1-OLx,1-OLy,k,bi,bj),
1620         I           salt (1-OLx,1-OLy,k,bi,bj),
1621       O           WORK2,       O           WORK2,
1622       I           myThid )       I           k, bi, bj, myThid )
1623              DO J = 1, sNy              DO J = 1, sNy
1624                 DO I = 1, sNx                 DO I = 1, sNx
1625                    IF ( k .LE. klowC(I,J,bi,bj) .AND.                    IF ( k .LE. klowC(I,J,bi,bj) .AND.
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.22