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

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

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

revision 1.14 by cnh, Sun Feb 4 14:38:47 2001 UTC revision 1.15 by cnh, Wed Sep 26 18:09:15 2001 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5  #define USE_FACTORIZED_POLY  #define USE_FACTORIZED_POLY
6    
7    CBOP
8    C     !ROUTINE: FIND_RHO
9    C     !INTERFACE:
10        subroutine FIND_RHO(        subroutine FIND_RHO(
11       I      bi, bj, iMin, iMax, jMin, jMax,  k, kRef, eqn,       I      bi, bj, iMin, iMax, jMin, jMax,  k, kRef, eqn,
12       I      tFld, sFld,       I      tFld, sFld,
13       O      rholoc,       O      rholoc,
14       I      myThid )       I      myThid )
15  C     /==========================================================\  
16  C     | o SUBROUTINE FIND_RHO                                    |  C     !DESCRIPTION: \bv
17  C     |   Calculates [rho(S,T,z)-Rhonil] of a slice              |  C     *==========================================================*
18  C     |==========================================================|  C     | o SUBROUTINE FIND_RHO                                    
19  C     |                                                          |  C     |   Calculates [rho(S,T,z)-Rhonil] of a slice              
20  C     | k - is the Theta/Salt level                              |  C     *==========================================================*
21  C     | kRef - determines pressure reference level               |  C     |                                                          
22  C     |        (not used in 'LINEAR' mode)                       |  C     | k - is the Theta/Salt level                              
23  C     | eqn - determines the eqn. of state: 'LINEAR' or 'POLY3'  |  C     | kRef - determines pressure reference level                
24  C     |                                                          |  C     |        (not used in 'LINEAR' mode)                        
25  C     \==========================================================/  C     | eqn - determines the eqn. of state: 'LINEAR' or 'POLY3'  
26    C     |                                                          
27    C     *==========================================================*
28    C     \ev
29    
30    C     !USES:
31        implicit none        implicit none
32  C     == Global variables ==  C     == Global variables ==
33  #include "SIZE.h"  #include "SIZE.h"
34  #include "EEPARAMS.h"  #include "EEPARAMS.h"
35  #include "PARAMS.h"  #include "PARAMS.h"
36    
37    C     !INPUT/OUTPUT PARAMETERS:
38  C     == Routine arguments ==  C     == Routine arguments ==
39        integer bi,bj,iMin,iMax,jMin,jMax        integer bi,bj,iMin,iMax,jMin,jMax
40        integer k                 ! Level of Theta/Salt slice        integer k                 ! Level of Theta/Salt slice
# Line 36  C     == Routine arguments == Line 45  C     == Routine arguments ==
45        _RL rholoc(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL rholoc(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
46        integer myThid        integer myThid
47    
48    C     !LOCAL VARIABLES:
49  C     == Local variables ==  C     == Local variables ==
50        integer i,j        integer i,j
51        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig
52        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
53    CEOP
54    
55  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
56        DO j=1-OLy,sNy+OLy        DO j=1-OLy,sNy+OLy
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22