/[MITgcm]/MITgcm/pkg/gridalt/make_phys_grid.F
ViewVC logotype

Diff of /MITgcm/pkg/gridalt/make_phys_grid.F

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

revision 1.1 by molod, Thu Jan 29 14:23:09 2004 UTC revision 1.2 by molod, Mon Feb 23 20:33:59 2004 UTC
# Line 58  c Require 12 bottom levels (300 mb worth Line 58  c Require 12 bottom levels (300 mb worth
58       .            2500.000,2500.000,2500.000,5000.000,5000.000,5000.000/       .            2500.000,2500.000,2500.000,5000.000,5000.000,5000.000/
59         _RL deltap, dpstar_accum         _RL deltap, dpstar_accum
60         integer nlbotmax, nstart, nlevs, nlphys, ndone         integer nlbotmax, nstart, nlevs, nlphys, ndone
61           integer nextra
62    c
63           nextra = ntry
64  c  c
65         do L = 1,Nr         do L = 1,Nr
66          do j = j1,j2          do j = j1,j2
# Line 86  c Line 89  c
89  c  c
90  c      do for each grid point:  c      do for each grid point:
91         do j = j1,j2         do j = j1,j2
92         do i = i1,i2+1         do i = i1,i2
93          Lbotij = Lbot(i,j,bi,bj)          Lbotij = Lbot(i,j,bi,bj)
94  c  c
95  c Find the maximum number of physics levels to fit in the bottom level  c Find the maximum number of physics levels to fit in the bottom level
# Line 166  c Make the one below thicker Line 169  c Make the one below thicker
169              endif              endif
170             endif             endif
171            enddo            enddo
172    C Need one more peice of logic - if we finished Lnew loop and
173    C  now we are done adding new physics layers, we need to be sure
174    C  that we are at the edge of a dynamics layer. if not, we need
175    C  to add one more layer.
176              if(nlphys.ge.ntry)then
177               if(abs(deltap-hfacC(i,j,L-1,bi,bj)*drF(L-1)).gt.0.001)then
178                nlphys = nlphys + 1
179                dpphys(i,j,nlphys,bi,bj) = hfacC(i,j,L-1,bi,bj)*drF(L-1)
180         .                                                          - deltap
181               endif
182              endif
183    
184           elseif(nlphys.eq.ntry)then           elseif(nlphys.eq.ntry)then
185  c Mostly done with new layers - make sure we end at dynamics edge,  c Mostly done with new layers - make sure we end at dynamics edge,
186  c      if not, make one more thinner (thinner than dyn grid) layer  c      if not, make one more thinner (thinner than dyn grid) layer
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22