/[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.3 by edhill, Wed May 5 00:39:21 2004 UTC revision 1.4 by molod, Mon Nov 1 21:35:24 2004 UTC
# Line 53  c Line 53  c
53         integer nlperdyn(im1:im2,jm1:jm2,Nr,Nsx,Nsy)         integer nlperdyn(im1:im2,jm1:jm2,Nr,Nsx,Nsy)
54  c  c
55         integer i,j,L,Lbotij,Lnew         integer i,j,L,Lbotij,Lnew
56  c Require 12 bottom levels (300 mb worth) for the physics, the dp's are:  c Require 12 (or 15) levels near the surface (300 mb worth) for fizhi.
57         integer ntry  c   the dp's are in the dptry arrays:
58         data ntry /12/         integer ntry,ntry10,ntry40
59         _RL dptry(12), dptry_accum(12)         data ntry10 /12/
60         data dptry /300.000, 600.000,1000.000,1400.000,1700.000,2500.000,         data ntry40 /15/
61           _RL dptry(15),dptry10(12),dptry40(15)
62           _RL dptry_accum(15)
63           data dptry10/300.000,600.000,1000.000,1400.000,1700.000,2500.000,
64       .            2500.000,2500.000,2500.000,5000.000,5000.000,5000.000/       .            2500.000,2500.000,2500.000,5000.000,5000.000,5000.000/
65           data dptry40/300.000,600.000,1000.000,1400.000,1700.000,
66         .            2500.000,2500.000,2500.000,2500.000,2500.000,2500.000,
67         .            2500.000,2500.000,2500.000,2500.000/
68         _RL deltap, dpstar_accum         _RL deltap, dpstar_accum
69         integer nlbotmax, nstart, nlevs, nlphys, ndone         integer nlbotmax, nstart, nlevs, nlphys, ndone
70         integer nextra         integer nextra
71  c  c
72           if( (Nr.eq.10) .or. (Nr.eq.20) ) then
73            ntry = ntry10
74            do L = 1,ntry
75             dptry(L) = dptry10(L)
76            enddo
77           elseif(Nr.eq.40) then
78            ntry = ntry40
79            do L = 1,ntry
80             dptry(L) = dptry40(L)
81            enddo
82           else
83            print *,' Dont know how to make fizhi grid '
84            stop
85           endif
86    c
87         nextra = ntry         nextra = ntry
88  c  c
89         do L = 1,Nr         do L = 1,Nr
# Line 217  c     them as thin (1 mb) layers near th Line 238  c     them as thin (1 mb) layers near th
238          nlperdyn(i,j,Nr,bi,bj) = numlevphys          nlperdyn(i,j,Nr,bi,bj) = numlevphys
239         endif         endif
240  c END OF LOOP OVER GRID POINTS  c END OF LOOP OVER GRID POINTS
241    
242         enddo         enddo
243         enddo         enddo
244    
245           stop
246    
247         return         return
248         end         end
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22