/[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.8 by molod, Mon Nov 27 16:09:55 2006 UTC revision 1.9 by molod, Sun Jul 8 20:24:51 2007 UTC
# Line 69  c   the dp's are in the dptry arrays: Line 69  c   the dp's are in the dptry arrays:
69         data bot_thick40/20000.000/         data bot_thick40/20000.000/
70         _RL deltap, dpstar_accum         _RL deltap, dpstar_accum
71         integer nlbotmax, nstart, nlevs, nlphys, ndone         integer nlbotmax, nstart, nlevs, nlphys, ndone
72           _RL thindp
73  c  c
74         if( (Nr.eq.10) .or. (Nr.eq.20) ) then         if( (Nr.eq.10) .or. (Nr.eq.20) ) then
75          ntry = ntry10          ntry = ntry10
# Line 76  c Line 77  c
77          do L = 1,ntry          do L = 1,ntry
78           dptry(L) = dptry10(L)           dptry(L) = dptry10(L)
79          enddo          enddo
80         elseif((Nr.eq.40).or.(Nr.eq.46)) then         elseif((Nr.eq.40).or.(Nr.eq.46).or.(Nr.eq.70)) then
81          ntry = ntry40          ntry = ntry40
82          bot_thick = bot_thick40          bot_thick = bot_thick40
83          do L = 1,ntry          do L = 1,ntry
# Line 86  c Line 87  c
87          print *,' Dont know how to make fizhi grid '          print *,' Dont know how to make fizhi grid '
88          stop          stop
89         endif         endif
90    
91           thindp=100.
92           if(Nr.eq.70)thindp=0.02
93  c  c
94         do L = 1,Nr         do L = 1,Nr
95          do j = j1,j2          do j = j1,j2
# Line 232  c  All done adding layers - if we need m Line 236  c  All done adding layers - if we need m
236  c     them as thin (1 mb) layers near the top  c     them as thin (1 mb) layers near the top
237         if(nlphys.lt.numlevphys)then         if(nlphys.lt.numlevphys)then
238          nlevs = numlevphys-nlphys          nlevs = numlevphys-nlphys
239          dpphys(i,j,nlphys,bi,bj) = dpphys(i,j,nlphys,bi,bj)-100. * nlevs          dpphys(i,j,nlphys,bi,bj)=dpphys(i,j,nlphys,bi,bj)-thindp*nlevs
240          do Lnew = nlphys+1,numlevphys          do Lnew = nlphys+1,numlevphys
241           dpphys(i,j,Lnew,bi,bj) = 100.           dpphys(i,j,Lnew,bi,bj) = thindp
242          enddo          enddo
243          nlperdyn(i,j,Nr,bi,bj) = numlevphys          nlperdyn(i,j,Nr,bi,bj) = numlevphys
244         endif         endif
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22