202 |
ENDDO |
ENDDO |
203 |
ENDDO |
ENDDO |
204 |
|
|
205 |
|
|
206 |
DO bj=myByLo(myThid),myByHi(myThid) |
DO bj=myByLo(myThid),myByHi(myThid) |
207 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
DO bi=myBxLo(myThid),myBxHi(myThid) |
208 |
|
|
233 |
jMin = 1-OLy+1 |
jMin = 1-OLy+1 |
234 |
jMax = sNy+OLy |
jMax = sNy+OLy |
235 |
|
|
236 |
|
|
237 |
K = 1 |
K = 1 |
238 |
BOTTOM_LAYER = K .EQ. Nr |
BOTTOM_LAYER = K .EQ. Nr |
239 |
|
|
275 |
I myThid ) |
I myThid ) |
276 |
ENDIF |
ENDIF |
277 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
278 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
279 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoKm1, |
280 |
O buoyKm1, |
O buoyKm1, |
281 |
I myThid ) |
I myThid ) |
315 |
I myThid ) |
I myThid ) |
316 |
ENDIF |
ENDIF |
317 |
C-- Calculate buoyancy |
C-- Calculate buoyancy |
318 |
CALL CALC_BUOY( |
CALL CALC_BUOYANCY( |
319 |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
I bi,bj,iMin,iMax,jMin,jMax,K,rhoK, |
320 |
O buoyK, |
O buoyK, |
321 |
I myThid ) |
I myThid ) |
450 |
C & maxval(salt(1:sNx,1:sNy,:,:,:)) |
C & maxval(salt(1:sNx,1:sNy,:,:,:)) |
451 |
C write(0,*) 'dynamics: phiHyd ',minval(phiHyd/(Gravity*Rhonil),mask=phiHyd.NE.0.), |
C write(0,*) 'dynamics: phiHyd ',minval(phiHyd/(Gravity*Rhonil),mask=phiHyd.NE.0.), |
452 |
C & maxval(phiHyd/(Gravity*Rhonil)) |
C & maxval(phiHyd/(Gravity*Rhonil)) |
453 |
|
C CALL PLOT_FIELD_XYZRL( gU, ' GU exiting dyanmics ' , |
454 |
|
C &Nr, 1, myThid ) |
455 |
|
C CALL PLOT_FIELD_XYZRL( gV, ' GV exiting dyanmics ' , |
456 |
|
C &Nr, 1, myThid ) |
457 |
|
C CALL PLOT_FIELD_XYZRL( gS, ' GS exiting dyanmics ' , |
458 |
|
C &Nr, 1, myThid ) |
459 |
|
C CALL PLOT_FIELD_XYZRL( gT, ' GT exiting dyanmics ' , |
460 |
|
C &Nr, 1, myThid ) |
461 |
|
|
462 |
|
|
463 |
RETURN |
RETURN |
464 |
END |
END |