/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_interp_vert.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_interp_vert.F

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

revision 1.2 by molod, Tue Oct 25 16:11:20 2005 UTC revision 1.5 by molod, Fri Nov 18 20:25:18 2005 UTC
# Line 93  C Use the biggest of nplevs 1-3 (or any Line 93  C Use the biggest of nplevs 1-3 (or any
93        _RL pksrf(sNx,sNy)        _RL pksrf(sNx,sNy)
94        _RL p        _RL p
95        _RL kappa        _RL kappa
96          _RL oneRL
97  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
98        INTEGER jpoint1,ipoint1        INTEGER jpoint1,ipoint1
99        INTEGER jpoint2,ipoint2        INTEGER jpoint2,ipoint2
# Line 105  C---+----1----+----2----+----3----+----4 Line 106  C---+----1----+----2----+----3----+----4
106    
107          IF(fflags(listId)(2:2).eq.'P') then          IF(fflags(listId)(2:2).eq.'P') then
108            kappa = getcon('KAPPA')            kappa = getcon('KAPPA')
109              oneRL = 1. _d 0
110    
111  c If nonlinear free surf is active, need averaged pressures  c If nonlinear free surf is active, need averaged pressures
112  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
# Line 127  C if fizhi is being  used, may need to g Line 129  C if fizhi is being  used, may need to g
129             endif             endif
130    
131             if(.not. foundp) then             if(.not. foundp) then
132              WRITE(msgBuf,'(3A)') 'DIAGNOSTICS_INTERP_VERT: ',              WRITE(msgBuf,'(4A)') 'DIAGNOSTICS_INTERP_VERT: ',
133       .    ' Have asked for pressure interpolation but have not ',       .    ' Have asked for pressure interpolation but have not ',
134       .    ' Activated surface and 3D pressure diagnostic, ',       .    ' Activated surface and 3D pressure diagnostic, ',
135       .    ' RSURF and PRESSURE'       .    ' RSURF and PRESSURE'
# Line 137  C if fizhi is being  used, may need to g Line 139  C if fizhi is being  used, may need to g
139    
140             DO bj = myByLo(myThid), myByHi(myThid)             DO bj = myByLo(myThid), myByHi(myThid)
141              DO bi = myBxLo(myThid), myBxHi(myThid)              DO bi = myBxLo(myThid), myBxHi(myThid)
142               call getdiag(1.,undef,qtmpsrf(1-OLx,1-OLy,bi,bj),               call getdiag(oneRL,undef,qtmpsrf(1-OLx,1-OLy,bi,bj),
143       .                       jpoint1,0,ipoint1,0,bi,bj,myThid)       .                       jpoint1,0,ipoint1,0,bi,bj,myThid)
144              ENDDO              ENDDO
145             ENDDO             ENDDO
# Line 197  C Load p to the kappa into a temporary a Line 199  C Load p to the kappa into a temporary a
199               DO k = 1,nplevs1               DO k = 1,nplevs1
200                p = plevs1(k)                p = plevs1(k)
201                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,
202       .                                                 nlevels(listId) )       .                                     nlevels(listId),myThid )
203               ENDDO               ENDDO
204              elseif(fflags(listId)(3:3).eq.'2')then              elseif(fflags(listId)(3:3).eq.'2')then
205               nlevsout = nplevs2               nlevsout = nplevs2
206               DO k = 1,nplevs2               DO k = 1,nplevs2
207                p = plevs2(k)                p = plevs2(k)
208                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,
209       .                                                 nlevels(listId) )       .                                     nlevels(listId),myThid )
210               ENDDO               ENDDO
211              elseif(fflags(listId)(3:3).eq.'3')then              elseif(fflags(listId)(3:3).eq.'3')then
212               nlevsout = nplevs3               nlevsout = nplevs3
213               DO k = 1,nplevs3               DO k = 1,nplevs3
214                p = plevs3(k)                p = plevs3(k)
215                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,                call prestopres(qprs(1,1,k),qinp,pkz,pksrf,0.,p,sNx,sNy,
216       .                                                 nlevels(listId) )       .                                     nlevels(listId),myThid )
217               ENDDO               ENDDO
218              endif              endif
219    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22