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

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

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

revision 1.32 by molod, Fri Jul 30 20:37:50 2004 UTC revision 1.33 by molod, Wed Aug 4 18:23:44 2004 UTC
# Line 249  C compute pressures - all pressure are c Line 249  C compute pressures - all pressure are c
249        call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj),        call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj),
250       .                                                 pkl(1,1,1,bi,bj))       .                                                 pkl(1,1,1,bi,bj))
251    
252        if(bi.eq.1 .and. 2.eq.1 )then        if(1.eq.1 )then
253        print *,' In do fizhi, before fizhi driver - bi = ',bi        print *,' In do fizhi, before fizhi driver - bi = ',bi
254        do L = 1,Nrphys+1        do L = 1,Nrphys+1
255         do j = jm1,jm2         do j = jm1,jm2
# Line 257  C compute pressures - all pressure are c Line 257  C compute pressures - all pressure are c
257          tempij(i,j) = pephy(i,j,L,bi,bj)          tempij(i,j) = pephy(i,j,L,bi,bj)
258         enddo         enddo
259         enddo         enddo
260         print *,' pephy at level ',l,' ',tempij  c      print *,' pephy at level ',l,' ',tempij
261        enddo        enddo
262        do L = 1,Nrphys        do L = 1,Nrphys
263         do j = jm1,jm2         do j = jm1,jm2
# Line 265  C compute pressures - all pressure are c Line 265  C compute pressures - all pressure are c
265          tempij(i,j) = u(i,j,L)          tempij(i,j) = u(i,j,L)
266         enddo         enddo
267         enddo         enddo
268         print *,' u at level ',l,' ',tempij  c      print *,' u at level ',l,' ',tempij
269        enddo        enddo
270        do L = 1,Nrphys        do L = 1,Nrphys
271         do j = jm1,jm2         do j = jm1,jm2
# Line 273  C compute pressures - all pressure are c Line 273  C compute pressures - all pressure are c
273          tempij(i,j) = v(i,j,L)          tempij(i,j) = v(i,j,L)
274         enddo         enddo
275         enddo         enddo
276         print *,' v at level ',l,' ',tempij  c      print *,' v at level ',l,' ',tempij
277        enddo        enddo
278        do L = 1,Nrphys        do L = 1,Nrphys
279         do j = jm1,jm2         do j = jm1,jm2
# Line 289  C compute pressures - all pressure are c Line 289  C compute pressures - all pressure are c
289          tempij(i,j) = q(i,j,L,1)          tempij(i,j) = q(i,j,L,1)
290         enddo         enddo
291         enddo         enddo
292         print *,' q at level ',l,' ',tempij  c      print *,' q at level ',l,' ',tempij
293        enddo        enddo
294        do L = 1,Nrphys        do L = 1,Nrphys
295         do j = jm1,jm2         do j = jm1,jm2
# Line 297  C compute pressures - all pressure are c Line 297  C compute pressures - all pressure are c
297          tempij(i,j) = pl(i,j,L,bi,bj)          tempij(i,j) = pl(i,j,L,bi,bj)
298         enddo         enddo
299         enddo         enddo
300         print *,' pl at level ',l,' ',tempij  c      print *,' pl at level ',l,' ',tempij
301        enddo        enddo
302        do L = 1,Nrphys+1        do L = 1,Nrphys+1
303         do j = jm1,jm2         do j = jm1,jm2
# Line 305  C compute pressures - all pressure are c Line 305  C compute pressures - all pressure are c
305          tempij(i,j) = ple(i,j,L,bi,bj)          tempij(i,j) = ple(i,j,L,bi,bj)
306         enddo         enddo
307         enddo         enddo
308         print *,' ple at level ',l,' ',tempij  c      print *,' ple at level ',l,' ',tempij
309        enddo        enddo
310        do L = 1,Nrphys        do L = 1,Nrphys
311         do j = jm1,jm2         do j = jm1,jm2
# Line 313  C compute pressures - all pressure are c Line 313  C compute pressures - all pressure are c
313          tempij(i,j) = pkl(i,j,L,bi,bj)          tempij(i,j) = pkl(i,j,L,bi,bj)
314         enddo         enddo
315         enddo         enddo
316         print *,' pkl at level ',l,' ',tempij  c      print *,' pkl at level ',l,' ',tempij
317        enddo        enddo
318        do L = 1,Nrphys+1        do L = 1,Nrphys+1
319         do j = jm1,jm2         do j = jm1,jm2
# Line 321  C compute pressures - all pressure are c Line 321  C compute pressures - all pressure are c
321          tempij(i,j) = pkle(i,j,L,bi,bj)          tempij(i,j) = pkle(i,j,L,bi,bj)
322         enddo         enddo
323         enddo         enddo
324         print *,' pkle at level ',l,' ',tempij  c      print *,' pkle at level ',l,' ',tempij
325        enddo        enddo
326        do L = 1,Nrphys        do L = 1,Nrphys
327         do j = jm1,jm2         do j = jm1,jm2
# Line 329  C compute pressures - all pressure are c Line 329  C compute pressures - all pressure are c
329          tempij(i,j) = dpres(i,j,L,bi,bj)          tempij(i,j) = dpres(i,j,L,bi,bj)
330         enddo         enddo
331         enddo         enddo
332         print *,' dpres at level ',l,' ',tempij  c      print *,' dpres at level ',l,' ',tempij
333        enddo        enddo
334        do j = jm1,jm2        do j = jm1,jm2
335        do i = im1,im2        do i = im1,im2
336         tempij(i,j) = tgz(i,j,bi,bj)         tempij(i,j) = tgz(i,j,bi,bj)
337        enddo        enddo
338        enddo        enddo
339        print *,' tgz ',tempij  c     print *,' tgz ',tempij
340    
341        print *,' Just before fizhi driver call '        print *,' Just before fizhi driver call '
342    
# Line 377  C compute pressures - all pressure are c Line 377  C compute pressures - all pressure are c
377       . qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj),       . qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj),
378       . fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj))       . fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj))
379    
380        if(bi.eq.1 .and. (2.eq.1))then        if(1.eq.1)then
381        print *,' In do fizhi, after fizhi driver - bi = ',bi        print *,' In do fizhi, after fizhi driver - bi = ',bi
382        do L = 1,Nrphys        do L = 1,Nrphys
383         do j = jm1,jm2         do j = jm1,jm2
# Line 385  C compute pressures - all pressure are c Line 385  C compute pressures - all pressure are c
385          tempij(i,j) = turbu(i,j,L,bi,bj)          tempij(i,j) = turbu(i,j,L,bi,bj)
386         enddo         enddo
387         enddo         enddo
388         print *,' turbu at level ',l,' ',tempij  c      print *,' turbu at level ',l,' ',tempij
389        enddo        enddo
390        do L = 1,Nrphys        do L = 1,Nrphys
391         do j = jm1,jm2         do j = jm1,jm2
# Line 393  C compute pressures - all pressure are c Line 393  C compute pressures - all pressure are c
393          tempij(i,j) = turbv(i,j,L,bi,bj)          tempij(i,j) = turbv(i,j,L,bi,bj)
394         enddo         enddo
395         enddo         enddo
396         print *,' turbv at level ',l,' ',tempij  c      print *,' turbv at level ',l,' ',tempij
397        enddo        enddo
398        do L = 1,Nrphys        do L = 1,Nrphys
399         do j = jm1,jm2         do j = jm1,jm2
# Line 409  C compute pressures - all pressure are c Line 409  C compute pressures - all pressure are c
409          tempij(i,j) = turbq(i,j,L,1,bi,bj)/p(i,j,bi,bj)          tempij(i,j) = turbq(i,j,L,1,bi,bj)/p(i,j,bi,bj)
410         enddo         enddo
411         enddo         enddo
412         print *,' turbq at level ',l,' ',tempij  c      print *,' turbq at level ',l,' ',tempij
413        enddo        enddo
414        do L = 1,Nrphys        do L = 1,Nrphys
415         do j = jm1,jm2         do j = jm1,jm2
# Line 417  C compute pressures - all pressure are c Line 417  C compute pressures - all pressure are c
417          tempij(i,j) = moistu(i,j,L,bi,bj)          tempij(i,j) = moistu(i,j,L,bi,bj)
418         enddo         enddo
419         enddo         enddo
420         print *,' moistu at level ',l,' ',tempij  c      print *,' moistu at level ',l,' ',tempij
421        enddo        enddo
422        do L = 1,Nrphys        do L = 1,Nrphys
423         do j = jm1,jm2         do j = jm1,jm2
# Line 425  C compute pressures - all pressure are c Line 425  C compute pressures - all pressure are c
425          tempij(i,j) = moistv(i,j,L,bi,bj)          tempij(i,j) = moistv(i,j,L,bi,bj)
426         enddo         enddo
427         enddo         enddo
428         print *,' moistv at level ',l,' ',tempij  c      print *,' moistv at level ',l,' ',tempij
429        enddo        enddo
430        do L = 1,Nrphys        do L = 1,Nrphys
431         do j = jm1,jm2         do j = jm1,jm2
# Line 441  C compute pressures - all pressure are c Line 441  C compute pressures - all pressure are c
441          tempij(i,j) = moistq(i,j,L,1,bi,bj)/p(i,j,bi,bj)          tempij(i,j) = moistq(i,j,L,1,bi,bj)/p(i,j,bi,bj)
442         enddo         enddo
443         enddo         enddo
444         print *,' moistq at level ',l,' ',tempij  c      print *,' moistq at level ',l,' ',tempij
445        enddo        enddo
446        do L = 1,Nrphys        do L = 1,Nrphys
447         do j = jm1,jm2         do j = jm1,jm2
# Line 487  C compute pressures - all pressure are c Line 487  C compute pressures - all pressure are c
487          tempij(i,j) = duphy(i,j,L,bi,bj)          tempij(i,j) = duphy(i,j,L,bi,bj)
488         enddo         enddo
489         enddo         enddo
490         print *,' duphy at level ',l,' ',tempij  c      print *,' duphy at level ',l,' ',tempij
491        enddo        enddo
492        do L = 1,Nrphys        do L = 1,Nrphys
493         do j = jm1,jm2         do j = jm1,jm2
# Line 495  C compute pressures - all pressure are c Line 495  C compute pressures - all pressure are c
495          tempij(i,j) = dvphy(i,j,L,bi,bj)          tempij(i,j) = dvphy(i,j,L,bi,bj)
496         enddo         enddo
497         enddo         enddo
498         print *,' dvphy at level ',l,' ',tempij  c      print *,' dvphy at level ',l,' ',tempij
499        enddo        enddo
500        do L = 1,Nrphys        do L = 1,Nrphys
501         do j = jm1,jm2         do j = jm1,jm2
# Line 511  C compute pressures - all pressure are c Line 511  C compute pressures - all pressure are c
511          tempij(i,j) = dsphy(i,j,L,bi,bj)          tempij(i,j) = dsphy(i,j,L,bi,bj)
512         enddo         enddo
513         enddo         enddo
514         print *,' dsphy at level ',l,' ',tempij  c      print *,' dsphy at level ',l,' ',tempij
515        enddo        enddo
516        endif        endif
517    
Legend:
Removed from v.1.32  
changed lines
  Added in v.1.33
  ViewVC Help
Powered by ViewVC 1.1.22