/[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.16 by cnh, Tue Jun 9 16:34:03 1998 UTC
# Line 174  C     uninitialised but inert locations. Line 174  C     uninitialised but inert locations.
174          rhokm1(i,j)  = 0. _d 0          rhokm1(i,j)  = 0. _d 0
175          rhokp1(i,j)  = 0. _d 0          rhokp1(i,j)  = 0. _d 0
176          rhotmp(i,j)  = 0. _d 0          rhotmp(i,j)  = 0. _d 0
177            maskC (i,j)  = 0. _d 0
178         ENDDO         ENDDO
179        ENDDO        ENDDO
180    
# Line 227  C--     Integrate hydrostatic balance fo Line 228  C--     Integrate hydrostatic balance fo
228       I      bi,bj,iMin,iMax,jMin,jMax,1,rhoKm1,rhoKm1,       I      bi,bj,iMin,iMax,jMin,jMax,1,rhoKm1,rhoKm1,
229       U      pH,       U      pH,
230       I      myThid )       I      myThid )
231          DO J=1-Oly,sNy+Oly          DO J=jMin,jMax
232           DO I=1-Olx,sNx+Olx           DO I=iMin,iMax
233            rhoKp1(I,J)=rhoKm1(I,J)            rhoKp1(I,J)=rhoKm1(I,J)
234           ENDDO           ENDDO
235          ENDDO          ENDDO
# Line 243  copt I     bi, bj, iMin, iMax, jMin, jMa Line 244  copt I     bi, bj, iMin, iMax, jMin, jMa
244  copt O     rhoKm1,  copt O     rhoKm1,
245  copt I     myThid )  copt I     myThid )
246  C       rhoKm1=rhoKp1  C       rhoKm1=rhoKp1
247          DO J=1-Oly,sNy+Oly          DO J=jMin,jMax
248           DO I=1-Olx,sNx+Olx           DO I=iMin,iMax
249            rhoKm1(I,J)=rhoKp1(I,J)            rhoKm1(I,J)=rhoKp1(I,J)
250           ENDDO           ENDDO
251          ENDDO          ENDDO
# Line 365  C--     Implicit diffusion Line 366  C--     Implicit diffusion
366         ENDDO         ENDDO
367        ENDDO        ENDDO
368    
369        write(0,*) 'dynamics: pS',minval(cg2d_x),maxval(cg2d_x)        write(0,*) 'dynamics: pS ',minval(cg2d_x(1:sNx,1:sNy,:,:)),
370        write(0,*) 'dynamics: U',minval(uVel(1:sNx,1:sNy,:,:,:)),       &                           maxval(cg2d_x(1:sNx,1:sNy,:,:))
371       &                         maxval(uVel(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: U  ',minval(uVel(1:sNx,1:sNy,:,:,:)),
372        write(0,*) 'dynamics: V',minval(vVel(1:sNx,1:sNy,:,:,:)),       &                           maxval(uVel(1:sNx,1:sNy,:,:,:))
373       &                         maxval(vVel(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: V  ',minval(vVel(1:sNx,1:sNy,:,:,:)),
374        write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),       &                           maxval(vVel(1:sNx,1:sNy,:,:,:))
375       &                         maxval(K13(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K13',minval(K13(1:sNx,1:sNy,:)),
376        write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),  cblk &                           maxval(K13(1:sNx,1:sNy,:))
377       &                         maxval(K23(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K23',minval(K23(1:sNx,1:sNy,:)),
378        write(0,*) 'dynamics: K33',minval(K33(1:sNx,1:sNy,:)),  cblk &                           maxval(K23(1:sNx,1:sNy,:))
379       &                         maxval(K33(1:sNx,1:sNy,:))  cblk  write(0,*) 'dynamics: K33',minval(K33(1:sNx,1:sNy,:)),
380        write(0,*) 'dynamics: gT',minval(gT(1:sNx,1:sNy,:,:,:)),  cblk &                           maxval(K33(1:sNx,1:sNy,:))
381       &                         maxval(gT(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: gT ',minval(gT(1:sNx,1:sNy,:,:,:)),
382        write(0,*) 'dynamics: T',minval(Theta(1:sNx,1:sNy,:,:,:)),       &                           maxval(gT(1:sNx,1:sNy,:,:,:))
383       &                         maxval(Theta(1:sNx,1:sNy,:,:,:))        write(0,*) 'dynamics: T  ',minval(Theta(1:sNx,1:sNy,:,:,:)),
384        write(0,*) 'dynamics: pH',minval(pH/(Gravity*Rhonil)),       &                           maxval(Theta(1:sNx,1:sNy,:,:,:))
385       &                          maxval(pH/(Gravity*Rhonil))  cblk  write(0,*) 'dynamics: pH ',minval(pH/(Gravity*Rhonil)),
386    cblk &                           maxval(pH/(Gravity*Rhonil))
387    
388        RETURN        RETURN
389        END        END
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22