/[MITgcm]/MITgcm/model/src/dynamics.F
ViewVC logotype

Diff of /MITgcm/model/src/dynamics.F

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

revision 1.14 by cnh, Mon Jun 8 21:43:01 1998 UTC revision 1.15 by adcroft, Tue Jun 9 15:58:36 1998 UTC
# Line 227  C--     Integrate hydrostatic balance fo Line 227  C--     Integrate hydrostatic balance fo
227       I      bi,bj,iMin,iMax,jMin,jMax,1,rhoKm1,rhoKm1,       I      bi,bj,iMin,iMax,jMin,jMax,1,rhoKm1,rhoKm1,
228       U      pH,       U      pH,
229       I      myThid )       I      myThid )
230          DO J=1-Oly,sNy+Oly          DO J=jMin,jMax
231           DO I=1-Olx,sNx+Olx           DO I=iMin,iMax
232            rhoKp1(I,J)=rhoKm1(I,J)            rhoKp1(I,J)=rhoKm1(I,J)
233           ENDDO           ENDDO
234          ENDDO          ENDDO
# Line 243  copt I     bi, bj, iMin, iMax, jMin, jMa Line 243  copt I     bi, bj, iMin, iMax, jMin, jMa
243  copt O     rhoKm1,  copt O     rhoKm1,
244  copt I     myThid )  copt I     myThid )
245  C       rhoKm1=rhoKp1  C       rhoKm1=rhoKp1
246          DO J=1-Oly,sNy+Oly          DO J=jMin,jMax
247           DO I=1-Olx,sNx+Olx           DO I=iMin,iMax
248            rhoKm1(I,J)=rhoKp1(I,J)            rhoKm1(I,J)=rhoKp1(I,J)
249           ENDDO           ENDDO
250          ENDDO          ENDDO
# Line 365  C--     Implicit diffusion Line 365  C--     Implicit diffusion
365         ENDDO         ENDDO
366        ENDDO        ENDDO
367    
368        write(0,*) 'dynamics: pS',minval(cg2d_x),maxval(cg2d_x)        write(0,*) 'dynamics: pS ',minval(cg2d_x(1:sNx,1:sNy,:,:)),
369        write(0,*) 'dynamics: U',minval(uVel(1:sNx,1:sNy,:,:,:)),       &                           maxval(cg2d_x(1:sNx,1:sNy,:,:))
370       &                         maxval(uVel(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: U  ',minval(uVel(1:sNx,1:sNy,:,:,:)),
371        write(0,*) 'dynamics: V',minval(vVel(1:sNx,1:sNy,:,:,:)),       &                           maxval(uVel(1:sNx,1:sNy,:,:,:))
372       &                         maxval(vVel(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: V  ',minval(vVel(1:sNx,1:sNy,:,:,:)),
373        write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),       &                           maxval(vVel(1:sNx,1:sNy,:,:,:))
374       &                         maxval(K13(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),
375        write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),  cblk &                           maxval(K13(1:sNx,1:sNy,:))
376       &                         maxval(K23(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),
377        write(0,*) 'dynamics: K33',minval(K33(1:sNx,1:sNy,:)),  cblk &                           maxval(K23(1:sNx,1:sNy,:))
378       &                         maxval(K33(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K33',minval(K33(1:sNx,1:sNy,:)),
379        write(0,*) 'dynamics: gT',minval(gT(1:sNx,1:sNy,:,:,:)),  cblk &                           maxval(K33(1:sNx,1:sNy,:))
380       &                         maxval(gT(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: gT ',minval(gT(1:sNx,1:sNy,:,:,:)),
381        write(0,*) 'dynamics: T',minval(Theta(1:sNx,1:sNy,:,:,:)),       &                           maxval(gT(1:sNx,1:sNy,:,:,:))
382       &                         maxval(Theta(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: T  ',minval(Theta(1:sNx,1:sNy,:,:,:)),
383        write(0,*) 'dynamics: pH',minval(pH/(Gravity*Rhonil)),       &                           maxval(Theta(1:sNx,1:sNy,:,:,:))
384       &                          maxval(pH/(Gravity*Rhonil))  cblk  write(0,*) 'dynamics: pH ',minval(pH/(Gravity*Rhonil)),
385    cblk &                           maxval(pH/(Gravity*Rhonil))
386    
387        RETURN        RETURN
388        END        END
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22