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

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

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

revision 1.5 by molod, Mon Jun 14 20:34:50 2004 UTC revision 1.7 by molod, Tue Jul 6 19:47:49 2004 UTC
# Line 34  c pe on physics grid refers to bottom ed Line 34  c pe on physics grid refers to bottom ed
34        integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2        integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
35        integer nhms1,nymd1,nhms2,nymd2,imns,ipls        integer nhms1,nymd1,nhms2,nymd2,imns,ipls
36        _RL facm, facp        _RL facm, facp
37          logical alarm
38          external alarm
39    
40        im1 = 1-OLx        im1 = 1-OLx
41        im2 = sNx+OLx        im2 = sNx+OLx
# Line 53  c  Construct the physics grid pressures Line 55  c  Construct the physics grid pressures
55  c                                         (even though dpphy counted bottom up)  c                                         (even though dpphy counted bottom up)
56          do j = 1,sNy          do j = 1,sNy
57          do i = 1,sNx          do i = 1,sNx
58           pephy(i,j,Nrphys+1,bi,bj)=(Ro_surf(i,j,bi,bj)+etaH(i,j,bi,bj))/           pephy(i,j,Nrphys+1,bi,bj)=(Ro_surf(i,j,bi,bj)+etaH(i,j,bi,bj))
      .                       rstarExpC(i,j,bi,bj)  
59           do L = 2,Nrphys+1           do L = 2,Nrphys+1
60           LL = Nrphys+2-L           LL = Nrphys+2-L
61           pephy(i,j,LL,bi,bj)=pephy(i,j,LL+1,bi,bj)-dpphys(i,j,L-1,bi,bj)           pephy(i,j,LL,bi,bj)=pephy(i,j,LL+1,bi,bj)-dpphys(i,j,L-1,bi,bj)
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22