/[MITgcm]/MITgcm/pkg/ecco/cost_trans_zonal.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_trans_zonal.F

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

revision 1.2 by heimbach, Thu Jul 28 13:51:36 2005 UTC revision 1.3 by jmc, Tue Oct 9 00:02:51 2007 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
5    
# Line 9  c     ================================== Line 10  c     ==================================
10  c     SUBROUTINE cost_trans_zonal  c     SUBROUTINE cost_trans_zonal
11  c     ==================================================================  c     ==================================================================
12  c  c
13  c     o Compute zonal transports.  c     o Compute zonal transports.
14  c  c
15  c     ==================================================================  c     ==================================================================
16  c     SUBROUTINE cost_trans_zonal  c     SUBROUTINE cost_trans_zonal
# Line 83  c--   5: P12 - Tasmania 145E Line 84  c--   5: P12 - Tasmania 145E
84        DATA beglat / -67.0, -14.5, -70.0, -67.0, -67.0 /        DATA beglat / -67.0, -14.5, -70.0, -67.0, -67.0 /
85        DATA endlat / -55.0,  -8.5, -31.0, -32.0, -42.0 /        DATA endlat / -55.0,  -8.5, -31.0, -32.0, -42.0 /
86  c  c
87  c      _RL ylat2,beglon2,endlon2        c      _RL ylat2,beglon2,endlon2
88  c      _RL ylat3,beglon3,endlon3        c      _RL ylat3,beglon3,endlon3
89  c      parameter(ylat= 29., beglon=-42., endlon =-2.)  c      parameter(ylat= 29., beglon=-42., endlon =-2.)
90  c      parameter(ylat= 29., beglon=282., endlon =352.)  c      parameter(ylat= 29., beglon=282., endlon =352.)
91  c      parameter(ylat= 29., beglon=-82., endlon =-2.)  c      parameter(ylat= 29., beglon=-82., endlon =-2.)
# Line 231  c Line 232  c
232          musum(isect) = musum(isect) + mu(isect)          musum(isect) = musum(isect) + mu(isect)
233    
234          write(msgbuf,'(a,i3,i5,2i3,5(X,D15.8))')          write(msgbuf,'(a,i3,i5,2i3,5(X,D15.8))')
235       &       'ECCO_TRANS_ZON ', isect, irec, kmin(isect), kmax(isect),       &       'ECCO_TRANS_ZON ', isect, irec, kmin(isect), kmax(isect),
236       &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)       &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)
237          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
238       &     SQUEEZE_RIGHT, mythid )       &     SQUEEZE_RIGHT, mythid )
# Line 241  c Line 242  c
242          call print_message( msgbuf, standardmessageunit,          call print_message( msgbuf, standardmessageunit,
243       &     SQUEEZE_RIGHT, mythid )       &     SQUEEZE_RIGHT, mythid )
244    
245          write(funit,'(i3,i5,2i3,5(X,D22.15))')          write(funit,'(i3,i5,2i3,5(X,D22.15))')
246       &       isect, irec, kmin(isect), kmax(isect),       &       isect, irec, kmin(isect), kmax(isect),
247       &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)       &       tu(isect), su(isect), mu(isect), mumin(isect), mumax(isect)
248    
249  c-- end loop over irec  c-- end loop over irec
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22