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

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

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

revision 1.18 by jmc, Tue Dec 5 05:25:08 2006 UTC revision 1.19 by jmc, Fri Sep 5 20:15:28 2008 UTC
# Line 162  C-    Calculate reciprol vertical grid s Line 162  C-    Calculate reciprol vertical grid s
162         recip_drF(k)   = 1. _d 0/drF(k)         recip_drF(k)   = 1. _d 0/drF(k)
163        ENDDO        ENDDO
164    
 C--   Calculate horizontal grid factor for the deep model (<=> dropping the  
 C      shallow atmosphere approximation): only function of the vertical index  
 C-    first: initialise deep-model grid factor:  
       DO k=1,Nr  
         deepFacC(k) = 1. _d 0  
         deepFac2C(k)= 1. _d 0  
         recip_deepFacC(k) = 1. _d 0  
         recip_deepFac2C(k)= 1. _d 0  
       ENDDO  
       DO k=1,Nr+1  
         deepFacF(k) = 1. _d 0  
         deepFac2F(k)= 1. _d 0  
         recip_deepFacF(k) = 1. _d 0  
         recip_deepFac2F(k)= 1. _d 0  
       ENDDO  
       IF ( deepAtmosphere ) THEN  
 C-    set deep-model grid factor:  
        IF ( usingZCoords ) THEN  
         DO k=1,Nr  
           deepFacC(k) = (rSphere+rC(k))*recip_rSphere  
           deepFac2C(k)  = deepFacC(k)*deepFacC(k)  
         ENDDO  
         DO k=1,Nr+1  
           deepFacF(k) = (rSphere+rF(k))*recip_rSphere  
           deepFac2F(k)  = deepFacF(k)*deepFacF(k)  
         ENDDO  
        ELSE  
         STOP 'INI_VERTICAL_GRID: setting deepFac is not coded'  
        ENDIF  
 C-    set reciprocal of deep-model grid factor:  
        DO k=1,Nr  
         recip_deepFacC(k) = 1. _d 0/deepFacC(k)  
         recip_deepFac2C(k)= 1. _d 0/deepFac2C(k)  
        ENDDO  
        DO k=1,Nr+1  
         recip_deepFacF(k) = 1. _d 0/deepFacF(k)  
         recip_deepFac2F(k)= 1. _d 0/deepFac2F(k)  
        ENDDO  
       ENDIF  
   
165        _END_MASTER(myThid)        _END_MASTER(myThid)
166        _BARRIER        _BARRIER
167    
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22