/[MITgcm]/MITgcm/pkg/aim/phy_vdifsc.F
ViewVC logotype

Diff of /MITgcm/pkg/aim/phy_vdifsc.F

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

revision 1.3 by cnh, Tue May 29 19:28:53 2001 UTC revision 1.4 by cnh, Mon Jun 18 17:39:58 2001 UTC
# Line 84  C Line 84  C
84            QTENVD(J,K) = 0.            QTENVD(J,K) = 0.
85          ENDDO          ENDDO
86        ENDDO        ENDDO
87    
88  c  c
89  C  C
90  C *****************************************  C *****************************************
91  C *****************************************  C *****************************************
92  Cchdbg  Cchdbg
93        if(ifirst) then  C     if(ifirst) then
94          xindconv=0.  C       xindconv=0.
95          xindconv1=0.  C       xindconv1=0.
96          npas=0  C       npas=0
97          ifirst=.FALSE.  C       ifirst=.FALSE.
98        endif  C     endif
99        npas = npas +1  C     npas = npas +1
100  Cchdbg  Cchdbg
101  C ******************************************  C ******************************************
102  C *****************************************  C *****************************************
# Line 115  C Line 116  C
116          RNL1(J) = -DSIG(NLEVxy(J,myThid))/DSIG(NL1(J))          RNL1(J) = -DSIG(NLEVxy(J,myThid))/DSIG(NL1(J))
117         ENDIF         ENDIF
118        ENDDO        ENDDO
119    
120  C  C
121  C  C
122  C New writing of the Conditional stability  C New writing of the Conditional stability
# Line 128  C -------------------------------------- Line 130  C --------------------------------------
130        ENDDO        ENDDO
131  C  C
132        DO J=1,NGP        DO J=1,NGP
133           stab(J)=0.
134         IF ( NLEVxy(J,myThid) .GT. 0 ) THEN         IF ( NLEVxy(J,myThid) .GT. 0 ) THEN
135          dThdp=(Th(J,NL1(J))-Th(J,NLEVxy(J,myThid)))          dThdp=(Th(J,NL1(J))-Th(J,NLEVxy(J,myThid)))
136       &              *((Prefw(NLEVxy(J,myThid))/Pground)**(RD/CP))*CP       &              *((Prefw(NLEVxy(J,myThid))/Pground)**(RD/CP))*CP
# Line 151  C Line 154  C
154  C ***************************************************  C ***************************************************
155  C ***************************************************  C ***************************************************
156  C chdbg  C chdbg
157            if(J.ge.6336 .and. J.eq.6348) then  C         if(J.ge.6336 .and. J.eq.6348) then
158               xindconv=xindconv+1./13.  C            xindconv=xindconv+1./13.
159            endif  C         endif
160            if(J.ge.4160 .and. J.eq.4172) then  C         if(J.ge.4160 .and. J.eq.4172) then
161               xindconv1=xindconv1+1./13.  C            xindconv1=xindconv1+1./13.
162            endif  C         endif
163            if(npas.eq.960 .and. J.eq.1) then  C         if(npas.eq.960 .and. J.eq.1) then
164              write(0,*) 'xindconv=',xindconv  C           write(0,*) 'xindconv=',xindconv
165              write(0,*) 'xindconv1=',xindconv1  C           write(0,*) 'xindconv1=',xindconv1
166            endif  C         endif
167  Cchdbg  Cchdbg
168  C ****************************************************  C ****************************************************
169  C ****************************************************  C ****************************************************
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22