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

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

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

revision 1.2 by cnh, Sun Sep 6 14:47:43 1998 UTC revision 1.5 by cnh, Fri Nov 6 22:44:43 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5        SUBROUTINE CALC_BUOYANCY(        SUBROUTINE CALC_BUOYANCY(
6       I      bi, bj, iMin, iMax, jMin, jMax,  k, rhoLoc,       I      bi, bj, iMin, iMax, jMin, jMax,  k, rhoLoc,
# Line 36  C--   == Routine arguments == Line 36  C--   == Routine arguments ==
36  C--   == Local variables ==  C--   == Local variables ==
37        INTEGER i,j        INTEGER i,j
38        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig        _RL refTemp,refSalt,sigRef,tP,sP,deltaSig
 ! ------------------------------------------------------------------------------  
39    
40        IF    ( buoyancyRelation .eq. 'ATMOSPHERIC'  ) THEN        IF    ( buoyancyRelation .eq. 'ATMOSPHERIC'  ) THEN
41    
42         DO j=jMin,jMax         DO j=jMin,jMax
43          DO i=iMin,iMax          DO i=iMin,iMax
44           buoy(i,j)=-specVol_S(K)/theta_S(K)*(theta(i,j,k,bi,bj)-tRef(K))           buoy(i,j)=(theta(i,j,k,bi,bj)-tRef(K))/tRef(K)
45          ENDDO          ENDDO
46         ENDDO         ENDDO
47                
# Line 56  C--   == Local variables == Line 55  C--   == Local variables ==
55    
56        ELSE        ELSE
57    
58         STOP 'CALC_BUOANCY: variable "buoyancyRelation" has an illegal value'         STOP
59         &'CALC_BUOANCY: variable "buoyancyRelation" has an illegal value'
60    
61        ENDIF        ENDIF
62    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22