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

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

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

revision 1.1 by molod, Thu Jan 29 14:22:24 2004 UTC revision 1.2 by molod, Mon Feb 23 20:34:38 2004 UTC
# Line 43  c pe on dynamics and physics grid refers Line 43  c pe on dynamics and physics grid refers
43    
44         integer i, j, L, Lbotij, bi, bj         integer i, j, L, Lbotij, bi, bj
45         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
46           _RL dt
47    
48         im1 = 1-OLx         im1 = 1-OLx
49         im2 = sNx+OLx         im2 = sNx+OLx
# Line 52  c pe on dynamics and physics grid refers Line 53  c pe on dynamics and physics grid refers
53         idim2 = sNx         idim2 = sNx
54         jdim1 = 1         jdim1 = 1
55         jdim2 = sNy         jdim2 = sNy
56           dt = 1.
57    
58         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
59         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
# Line 152  c    Fifth: Interpolate correction terms Line 154  c    Fifth: Interpolate correction terms
154         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
155    
156          call dyn2phys(udyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,          call dyn2phys(udyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,
157       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,uphytemp)       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,uphytemp)
158          call dyn2phys(vdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,          call dyn2phys(vdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,
159       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,vphytemp)       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,vphytemp)
160          call dyn2phys(thdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,          call dyn2phys(thdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,
161       .     1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,thphytemp)       .     1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,thphytemp)
162          call dyn2phys(sdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,          call dyn2phys(sdyntemp,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,1,sNx,
163       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,sphytemp)       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,sphytemp)
 c    Last: Increment physics state by the correction term  
164    
165          call step_physics(uphy,vphy,thphy,sphy,im1,im2,jm1,jm2,Nrphys,  c    Last: Increment physics state by the correction term
166       .   Nsx,Nsy,1,sNx,1,sNy,bi,bj,uphytemp,vphytemp,thphytemp,sphytemp)          call step_physics(uphy,vphy,thphy,sphy,dt,im1,im2,jm1,jm2,
167         .   Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj,
168         .                            uphytemp,vphytemp,thphytemp,sphytemp)
169    
170         enddo         enddo
171         enddo         enddo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22