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

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

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

revision 1.31 by molod, Thu May 18 23:04:46 2006 UTC revision 1.32 by molod, Mon Jun 19 23:37:34 2006 UTC
# Line 48  c pe on dynamics and physics grid refers Line 48  c pe on dynamics and physics grid refers
48         _RL tempLdiag(sNx,sNy,Nrphys,Nsx,Nsy)         _RL tempLdiag(sNx,sNy,Nrphys,Nsx,Nsy)
49         _RL tempLdiag2(sNx,sNy,Nrphys,Nsx,Nsy)         _RL tempLdiag2(sNx,sNy,Nrphys,Nsx,Nsy)
50         _RL tempdiag(sNx,sNy,Nsx,Nsy)         _RL tempdiag(sNx,sNy,Nsx,Nsy)
51         _RL slp(sNx,sNy)         _RL slp(sNx,sNy),meantemp(sNx,sNy)
52    
53         integer i, j, L, Lbotij, bi, bj         integer i, j, L, Lbotij, bi, bj
54         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
55         _RL tempij(1-olx:sNx+olx,1-oly:sNy+oly)         _RL tempij(1-olx:sNx+olx,1-oly:sNy+oly)
56           _RL grav, getcon
57    
58           grav = getcon('GRAVITY')
59         idim1 = 1-OLx         idim1 = 1-OLx
60         idim2 = sNx+OLx         idim2 = sNx+OLx
61         jdim1 = 1-OLy         jdim1 = 1-OLy
# Line 233  c Do not use a zero field as the top edg Line 235  c Do not use a zero field as the top edg
235         if(diagnostics_is_on('SLP     ',myThid) ) then         if(diagnostics_is_on('SLP     ',myThid) ) then
236          do j = 1,sNy          do j = 1,sNy
237          do i = 1,sNx          do i = 1,sNx
238           tempdiag(i,j,bi,bj) = Phihydlow(i,j,bi,bj)           tempdiag(i,j,bi,bj) = topoZ(i,j,bi,bj)*grav
239          do L = 1,Nrphys+1          do L = 1,Nrphys
240           tempLdiag(i,j,L,bi,bj) = pephy4fiz(i,j,L,bi,bj)           tempLdiag(i,j,L,bi,bj) = pephy4fiz(i,j,L,bi,bj)/100.
241           tempLdiag2(i,j,L,bi,bj) = thphy(i,j,L,bi,bj)           tempLdiag2(i,j,L,bi,bj) = thphy(i,j,L,bi,bj) *
242         .        (1.+0.609*sphy(i,j,L,bi,bj))
243          enddo          enddo
244            tempLdiag(i,j,Nrphys+1,bi,bj)=pephy4fiz(i,j,Nrphys+1,bi,bj)/100.
245          enddo          enddo
246          enddo          enddo
247          call slprs(tempdiag,tempLdiag(1,1,1,bi,bj),          call slprs(tempdiag(1,1,bi,bj),tempLdiag(1,1,1,bi,bj),
248       .  tempLdiag2(1,1,1,bi,bj),fracland(1,1,bi,bj),sNx,sNy,Nrphys,slp)       .  tempLdiag2(1,1,1,bi,bj),fracland(1,1,bi,bj),sNx,sNy,Nrphys,slp)
249          do j = 1,sNy          do j = 1,sNy
250          do i = 1,sNx          do i = 1,sNx
Legend:
Removed from v.1.31  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.22