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

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

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

revision 1.6 by molod, Tue May 31 20:14:37 2005 UTC revision 1.7 by molod, Wed Jun 1 18:08:33 2005 UTC
# Line 120  c ----------------------- Line 120  c -----------------------
120        enddo        enddo
121        enddo        enddo
122    
       if(diagnostics_is_on('SDIAG1  ',myid) ) then  
        do j=1,jm  
        do i=1,im  
         tmpdiag(i,j) = float(nthin(i,j))  
        enddo  
        enddo  
        call diagnostics_fill(tmpdiag,'SDIAG1  ',0,1,3,bi,bj,myid)  
       endif  
       if(diagnostics_is_on('SDIAG2  ',myid) ) then  
        do j=1,jm  
        do i=1,im  
         tmpdiag(i,j) = float(nbase(i,j))  
        enddo  
        enddo  
        call diagnostics_fill(tmpdiag,'SDIAG2  ',0,1,3,bi,bj,myid)  
       endif  
   
123  c Compute Topography Sub-Grid Standard Deviation  c Compute Topography Sub-Grid Standard Deviation
124  c        and constrain the Maximum Value  c        and constrain the Maximum Value
125  c ----------------------------------------------  c ----------------------------------------------
# Line 146  c -------------------------------------- Line 129  c --------------------------------------
129        enddo        enddo
130        enddo        enddo
131    
       if(diagnostics_is_on('SDIAG3  ',myid) ) then  
        do j=1,jm  
        do i=1,im  
         tmpdiag(i,j) = phis_std(i,j)  
        enddo  
        enddo  
        call diagnostics_fill(tmpdiag,'SDIAG3  ',0,1,3,bi,bj,myid)  
       endif  
   
132  c Compute Virtual Temperatures  c Compute Virtual Temperatures
133  c ----------------------------  c ----------------------------
134        do L = 1,Lm        do L = 1,Lm
# Line 180  c -------------------------------------- Line 154  c --------------------------------------
154    
155        do n=1,npcs        do n=1,npcs
156    
157        call strip ( phis_std,std,im*jm,istrip,1,n )        call stripit ( phis_std,std,im*jm,im*jm,istrip,1,n )
158    
159        call strip ( pz,ps,im*jm,istrip,1 ,n )        call stripit ( pz,ps,im*jm,im*jm,istrip,1 ,n )
160        call strip ( uz,us,im*jm,istrip,Lm,n )        call stripit ( uz,us,im*jm,im*jm,istrip,Lm,n )
161        call strip ( vz,vs,im*jm,istrip,Lm,n )        call stripit ( vz,vs,im*jm,im*jm,istrip,Lm,n )
162        call strip ( tv,ts,im*jm,istrip,Lm,n )        call stripit ( tv,ts,im*jm,im*jm,istrip,Lm,n )
163        call strip ( pl,plstr,im*jm,istrip,Lm,n )        call stripit ( pl,plstr,im*jm,im*jm,istrip,Lm,n )
164        call strip ( ple,plestr,im*jm,istrip,Lm,n )        call stripit ( ple,plestr,im*jm,im*jm,istrip,Lm,n )
165        call strip ( dpres,dpresstr,im*jm,istrip,Lm,n )        call stripit ( dpres,dpresstr,im*jm,im*jm,istrip,Lm,n )
166        call stripint ( nthin,nthinstr,im*jm,istrip,1,n )        call stripitint ( nthin,nthinstr,im*jm,im*jm,istrip,1,n )
167        call stripint ( nbase,nbasestr,im*jm,istrip,1,n )        call stripitint ( nbase,nbasestr,im*jm,im*jm,istrip,1,n )
168    
169        call GWDD ( ps,us,vs,ts,        call GWDD ( ps,us,vs,ts,
170       .            dragus,dragvs,dragxs,dragys,std,       .            dragus,dragvs,dragxs,dragys,std,
171       .            plstr,plestr,dpresstr,grav,rgas,cp,       .            plstr,plestr,dpresstr,grav,rgas,cp,
172       .            istrip,Lm,nthinstr,nbasestr,lstar )       .            istrip,Lm,nthinstr,nbasestr,lstar )
173    
174        call paste ( dragus,dragu,istrip,im*jm,Lm,n )        call pastit( dragus,dragu,istrip,im*jm,im*jm,Lm,n )
175        call paste ( dragvs,dragv,istrip,im*jm,Lm,n )        call pastit( dragvs,dragv,istrip,im*jm,im*jm,Lm,n )
176        call paste ( dragxs,dragx,istrip,im*jm,1 ,n )        call pastit( dragxs,dragx,istrip,im*jm,im*jm,1 ,n )
177        call paste ( dragys,dragy,istrip,im*jm,1 ,n )        call pastit( dragys,dragy,istrip,im*jm,im*jm,1 ,n )
178    
179        enddo        enddo
180    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22