/[MITgcm]/MITgcm/pkg/fizhi/fizhi_rayleigh.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_rayleigh.F

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

revision 1.7 by molod, Tue Dec 14 19:56:45 2004 UTC revision 1.8 by molod, Sat May 21 23:50:13 2005 UTC
# Line 35  C ************************************** Line 35  C **************************************
35    
36        implicit none        implicit none
37    
 #ifdef ALLOW_DIAGNOSTICS  
 #include "SIZE.h"  
 #include "DIAGNOSTICS_SIZE.h"  
 #include "DIAGNOSTICS.h"  
 #endif  
   
38        integer myid,im,jm,lm,bi,bj        integer myid,im,jm,lm,bi,bj
39        _RL zsurf(im,jm),pres(im,jm,lm),pkap(im,jm,lm)        _RL zsurf(im,jm),pres(im,jm,lm),pkap(im,jm,lm)
40        _RL pekap(im,jm,lm+1)        _RL pekap(im,jm,lm+1)
# Line 52  C ************************************** Line 46  C **************************************
46        _RL z(im,jm,lm)        _RL z(im,jm,lm)
47        _RL dz(im,jm,lm)        _RL dz(im,jm,lm)
48        _RL cpog, cpinv, virtcon, getcon, dampcoef        _RL cpog, cpinv, virtcon, getcon, dampcoef
49    #ifdef ALLOW_DIAGNOSTICS
50          logical  diagnostics_is_on
51          external diagnostics_is_on
52          _RL tmpdiag(im,jm)
53    #endif
54    
55  C **********************************************************************  C **********************************************************************
56  C ****   APPLY RAYLEIGH FRICTION TO WIND (INCLUDE HEATING)           ***  C ****   APPLY RAYLEIGH FRICTION TO WIND (INCLUDE HEATING)           ***
# Line 99  C ************************************** Line 98  C **************************************
98    
99  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
100        do L=1,lm        do L=1,lm
101         if( irfu.ne.0 ) then  
102          if(diagnostics_is_on('RFU     ',myid) ) then
103         do j=1,jm         do j=1,jm
104         do i=1,im         do i=1,im
105          qdiag(i,j,irfu+L-1,bi,bj) = qdiag(i,j,irfu+L-1,bi,bj) +          tmpdiag(i,j) = rfu(i,j,L)*86400
      .                  rfu(i,j,L)*86400  
106         enddo         enddo
107         enddo         enddo
108         endif         call diagnostics_fill(tmpdiag,'RFU     ',L,1,3,bi,bj,myid)
109          endif
110    
111         if( irfv.ne.0 ) then        if(diagnostics_is_on('RFV     ',myid) ) then
112         do j=1,jm         do j=1,jm
113         do i=1,im         do i=1,im
114          qdiag(i,j,irfv+L-1,bi,bj) = qdiag(i,j,irfv+L-1,bi,bj) +          tmpdiag(i,j) = rfv(i,j,L)*86400
      .                  rfv(i,j,L)*86400  
115         enddo         enddo
116         enddo         enddo
117         endif         call diagnostics_fill(tmpdiag,'RFV     ',L,1,3,bi,bj,myid)
118          endif
119    
120         if( irft.ne.0 ) then        if(diagnostics_is_on('RFT     ',myid) ) then
121         do j=1,jm         do j=1,jm
122         do i=1,im         do i=1,im
123          qdiag(i,j,irft+L-1,bi,bj) = qdiag(i,j,irft+L-1,bi,bj) +          tmpdiag(i,j) = rft(i,j,L)*86400
      .                  rft(i,j,L)*86400  
124         enddo         enddo
125         enddo         enddo
126         endif         call diagnostics_fill(tmpdiag,'RFT     ',L,1,3,bi,bj,myid)
       enddo  
 #endif  
   
 #ifdef ALLOW_DIAGNOSTICS  
       if( (bi.eq.1) .and. (bj.eq.1) ) then  
        nrfu = nrfu + 1  
        nrfv = nrfv + 1  
        nrft = nrft + 1  
127        endif        endif
128    
129          enddo
130  #endif  #endif
131    
132        return        return
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22