/[MITgcm]/MITgcm/pkg/mom_common/mom_quasihydrostatic.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_quasihydrostatic.F

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

revision 1.4 by jmc, Mon Feb 5 03:16:23 2007 UTC revision 1.5 by jmc, Tue Dec 18 01:18:01 2012 UTC
# Line 41  C  myThid        :: my Thread Id number Line 41  C  myThid        :: my Thread Id number
41        INTEGER myThid        INTEGER myThid
42    
43  C !OUTPUT/MODIFIED VARIABLES: ==========================================  C !OUTPUT/MODIFIED VARIABLES: ==========================================
44  C  effectiveBuoy :: Density anomaly  C  effectiveBuoy :: Density (z-coord) / specific volume (p-coord) anomaly
45        _RL effectiveBuoy(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL effectiveBuoy(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46    
47  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
48  C  i,j           :: loop indices  C  i,j           :: loop indices
# Line 54  C  scalingFactor :: scaling factor (from Line 54  C  scalingFactor :: scaling factor (from
54        _RL scalingFactor        _RL scalingFactor
55  CEOP  CEOP
56    
57        scalingFactor=rhoConst*gravitySign*recip_gravity        IF ( usingZCoords ) THEN
58    C--   Z-coordinate case: Input is density anomaly
59    
60            scalingFactor=rhoConst*gravitySign*recip_gravity
61  C-    to reproduce old (wrong) results:  C-    to reproduce old (wrong) results:
62  c     scalingFactor=gravitySign*recip_gravity  c       scalingFactor=gravitySign*recip_gravity
63    
64          ELSEIF ( fluidIsWater ) THEN
65    C--   P-coordinate, oceanic case: Input is specific-volume anomaly
66    
67            scalingFactor = recip_rhoConst*recip_gravity
68    c       scalingFactor = rVel2wUnit(k) <-- @ interface = wrong location
69    C-    should use expression below (consistent with omega <-> w-velocity
70    C      conversion) but rhoRef(k) = rho(tRef,sRef,p) is computed
71    C     in S/R SET_REF_STATE but is not stored:
72    c       scalingFactor = ( oneRL / rhoRef(k) )*recip_gravity
73    
74          ELSE
75    C--   P-coord., Ideal-Gas case: Input is virtual potential temp. anomaly
76    C     (see White & Bromley, QJRMS 1995)
77            scalingFactor = tRef(k)*recip_gravity
78    
79          ENDIF
80    
81        IF ( use3dCoriolis ) THEN        IF ( use3dCoriolis ) THEN
82         DO j=jMin,jMax         DO j=jMin,jMax
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22