/[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.3 by molod, Tue Mar 30 18:28:44 2004 UTC revision 1.4 by molod, Tue Mar 30 18:52:39 2004 UTC
# Line 55  c pe on dynamics and physics grid refers Line 55  c pe on dynamics and physics grid refers
55         jdim2 = sNy         jdim2 = sNy
56         dt = 1.         dt = 1.
57    
        CALL TIMER_START('SETUP FIZHI CORR  [STEP_FIZHI_CORR]',mythid)  
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)
60    
# Line 110  c Create a wind magnitude field on the p Line 109  c Create a wind magnitude field on the p
109          enddo          enddo
110         enddo         enddo
111         enddo         enddo
        CALL TIMER_STOP('SETUP FIZHI CORR  [STEP_FIZHI_CORR]',mythid)  
112    
113         CALL TIMER_START('PHYS2DYN          [STEP_FIZHI_CORR]',mythid)         CALL TIMER_START('PHYS2DYN          [STEP_FIZHI_CORR]',mythid)
114         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 168  c    Fifth: Interpolate correction terms Line 166  c    Fifth: Interpolate correction terms
166         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
167    
168          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,
169       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,uphytemp)       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,uphytemp)
170          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,
171       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,vphytemp)       .      1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,1,vphytemp)
172          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,
173       .     1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,thphytemp)       .     1,sNy,bi,bj,windphy,pephy,ksurfC,Nrphys,nlperdyn,0,thphytemp)
174          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,
# Line 180  c    Fifth: Interpolate correction terms Line 178  c    Fifth: Interpolate correction terms
178         CALL TIMER_STOP('DYN2PHYS          [STEP_FIZHI_CORR]',mythid)         CALL TIMER_STOP('DYN2PHYS          [STEP_FIZHI_CORR]',mythid)
179    
180  c    Last: Increment physics state by the correction term  c    Last: Increment physics state by the correction term
        CALL TIMER_START('STEP_PHYSICS      [STEP_FIZHI_CORR]',mythid)  
181         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
182         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
183          call step_physics(uphy,vphy,thphy,sphy,dt,im1,im2,jm1,jm2,          call step_physics(uphy,vphy,thphy,sphy,dt,im1,im2,jm1,jm2,
# Line 189  c    Last: Increment physics state by th Line 186  c    Last: Increment physics state by th
186    
187         enddo         enddo
188         enddo         enddo
        CALL TIMER_STOP('STEP_PHYSICS      [STEP_FIZHI_CORR]',mythid)  
189    
190         return         return
191         end         end
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22