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

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

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

revision 1.35 by adcroft, Thu Jun 7 15:53:21 2001 UTC revision 1.36 by jmc, Fri Jul 6 21:39:37 2001 UTC
# Line 46  C     sumArea - Work variable used to co Line 46  C     sumArea - Work variable used to co
46        _RL     sumArea        _RL     sumArea
47        _RL     aC, aCw, aCs        _RL     aC, aCw, aCs
48    
 C-- Initialise -Boyancy at surface level : Bo_surf  
 C   Bo_surf is defined as d/dr(Phi_surf) and set to g/rtoz (linear free surface)  
 C     with rtoz = conversion factor from r-unit to z-unit  (=horiVertRatio)  
 C   an acurate formulation will include P_surf and T,S_surf effects on rho_surf:  
 C    z-ocean (rtoz=1) : Bo_surf = - Boyancy = gravity * rho_surf/rho_0  
 C    p-atmos (rtoz=rho_c*g) : Bo_surf = (1/rho)_surf  
 C--  
       IF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN  
         DO bj=myByLo(myThid),myByHi(myThid)  
          DO bi=myBxLo(myThid),myBxHi(myThid)  
           DO J=1-Oly,sNy+Oly  
            DO I=1-Olx,sNx+Olx  
              Bo_surf(I,J,bi,bj) = recip_rhoConst  
              recip_Bo(I,J,bi,bj) = rhoConst  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
       ELSE  
 C-  gBaro = gravity (except for External mode test with reduced gravity)  
         DO bj=myByLo(myThid),myByHi(myThid)  
          DO bi=myBxLo(myThid),myBxHi(myThid)  
           DO J=1-Oly,sNy+Oly  
            DO I=1-Olx,sNx+Olx  
              Bo_surf(I,J,bi,bj) = gBaro  
              recip_Bo(I,J,bi,bj) = 1. _d 0 / gBaro  
            ENDDO  
           ENDDO  
          ENDDO  
         ENDDO  
       ENDIF  
   
 C--   Update overlap regions  
       _EXCH_XY_R8(Bo_surf, myThid)  
       _EXCH_XY_R8(recip_Bo, myThid)  
   
49  C--   Initialise laplace operator  C--   Initialise laplace operator
50  C     aW2d: integral in Z Ax/dX  C     aW2d: integral in Z Ax/dX
51  C     aS2d: integral in Z Ay/dY  C     aS2d: integral in Z Ay/dY
Legend:
Removed from v.1.35  
changed lines
  Added in v.1.36
  ViewVC Help
Powered by ViewVC 1.1.22