228 |
enddo |
enddo |
229 |
c |
c |
230 |
c All done adding layers - if we need more to make numlevphys, put |
c All done adding layers - if we need more to make numlevphys, put |
231 |
c them as thin (0.1 mb) layers near the top |
c them as thin (1 mb) layers near the top |
232 |
if(nlphys.lt.numlevphys)then |
if(nlphys.lt.numlevphys)then |
233 |
nlevs = numlevphys-nlphys |
nlevs = numlevphys-nlphys |
234 |
dpphys(i,j,nlphys,bi,bj) = dpphys(i,j,nlphys,bi,bj)-10. * nlevs |
dpphys(i,j,nlphys,bi,bj) = dpphys(i,j,nlphys,bi,bj)-100. * nlevs |
235 |
do Lnew = nlphys+1,numlevphys |
do Lnew = nlphys+1,numlevphys |
236 |
dpphys(i,j,Lnew,bi,bj) = 10. |
dpphys(i,j,Lnew,bi,bj) = 100. |
237 |
enddo |
enddo |
238 |
nlperdyn(i,j,Nr,bi,bj) = numlevphys |
nlperdyn(i,j,Nr,bi,bj) = numlevphys |
239 |
endif |
endif |