/[MITgcm]/MITgcm/pkg/aim_v23/aim_dyn2aim.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/aim_dyn2aim.F

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

revision 1.3 by jmc, Thu Apr 8 00:14:09 2004 UTC revision 1.4 by jmc, Thu Jul 8 15:51:19 2004 UTC
# Line 101  c        PSA(I2) = rF(K)/atm_po Line 101  c        PSA(I2) = rF(K)/atm_po
101        ENDDO        ENDDO
102    
103  C-    Set cell delta_P fraction (of the full delta.P = drF_k):  C-    Set cell delta_P fraction (of the full delta.P = drF_k):
104        DO K = 1,Nr  #ifdef NONLIN_FRSURF
105         Katm = _KD2KA( K )        IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
106         DO J = 1,sNy         IF ( select_rStar.GT.0 ) THEN
107          DO I = 1,sNx          DO k = 1,Nr
108           I2 = I+(J-1)*sNx           Katm = _KD2KA( k )
109  c        dpFac(I2,Katm) = 1. _d 0           DO j = 1,sNy
110           dpFac(I2,Katm) = hFacC(I,J,K,bi,bj)            DO i = 1,sNx
111  c        IF (hFacC(I,J,K,bi,bj).GT.0. _d 0) THEN             I2 = i+(j-1)*sNx
112  c          dpFac(I2,Katm) = hFacC(I,J,K,bi,bj)             dpFac(I2,Katm) = h0FacC(i,j,k,bi,bj)*rStarFacC(i,j,bi,bj)
 c        ELSE  
113  c          dpFac(I2,Katm) = 1. _d 0  c          dpFac(I2,Katm) = 1. _d 0
114  c        ENDIF            ENDDO
115             ENDDO
116          ENDDO          ENDDO
117         ENDDO         ELSE
118        ENDDO          DO k = 1,Nr
119             Katm = _KD2KA( k )
120             DO j = 1,sNy
121              DO i = 1,sNx
122               I2 = i+(j-1)*sNx
123               IF ( k.EQ.ksurfC(i,j,bi,bj) ) THEN
124                dpFac(I2,Katm) = hFac_surfC(i,j,bi,bj)
125               ELSE
126                dpFac(I2,Katm) = hFacC(i,j,k,bi,bj)
127               ENDIF
128    c          dpFac(I2,Katm) = 1. _d 0
129              ENDDO
130             ENDDO
131            ENDDO
132           ENDIF
133          ELSE
134    #else /* ndef NONLIN_FRSURF */
135          IF (.TRUE.) THEN
136    #endif /* NONLIN_FRSURF */
137            DO k = 1,Nr
138             Katm = _KD2KA( k )
139             DO j = 1,sNy
140              DO i = 1,sNx
141               I2 = i+(j-1)*sNx
142               dpFac(I2,Katm) = hFacC(i,j,k,bi,bj)
143    c          dpFac(I2,Katm) = 1. _d 0
144              ENDDO
145             ENDDO
146            ENDDO
147          ENDIF
148    
149  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
150    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22