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

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

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

revision 1.2 by molod, Tue Jun 15 16:06:03 2004 UTC revision 1.3 by molod, Thu Jun 24 19:57:02 2004 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5        subroutine lwrio (nymd,nhms,istrip,npcs,        subroutine lwrio (nymd,nhms,bi,bj,istrip,npcs,
6       .                  pz,tz,qz,plz,plze,pkz,pkht,oz,co2,       .                  pz,tz,qz,plz,plze,pkz,pkht,oz,co2,
7       .                  cfc11,cfc12,cfc22,       .                  cfc11,cfc12,cfc22,
8       .                  methane,n2o,emissivity,       .                  methane,n2o,emissivity,
# Line 19  C $Name$ Line 19  C $Name$
19    
20  c Input Variables  c Input Variables
21  c ---------------  c ---------------
22        integer nymd,nhms,istrip,npcs        integer nymd,nhms,istrip,npcs,bi,bj
23    
24        integer im,jm,lm                integer im,jm,lm        
25        real  ptop                      real  ptop              
# Line 279  C ************************************** Line 279  C **************************************
279  C ****            PASTE AND BUMP SOME DIAGNOSTICS                   ****  C ****            PASTE AND BUMP SOME DIAGNOSTICS                   ****
280  C **********************************************************************  C **********************************************************************
281    
282        IF(IOLR.GT.0)CALL PSTBMP(flx(1,1),QDIAG(1,1,IOLR),ISTRIP,        IF(IOLR.GT.0)CALL PSTBMP(flx(1,1),QDIAG(1,1,IOLR,bi,bj),ISTRIP,
283       .                                                      im*jm, 1,NN)       .                                                      im*jm, 1,NN)
284        IF(IOLRCLR.GT.0)CALL PSTBMP(flxclr(1,1),QDIAG(1,1,IOLRCLR),ISTRIP,        IF(IOLRCLR.GT.0)CALL PSTBMP(flxclr(1,1),QDIAG(1,1,IOLRCLR,bi,bj),
285       .                                                       im*jm,1,NN)       .                                                ISTRIP,im*jm,1,NN)
286        IF(IOZLW.GT.0)CALL PSTBMP(OZL(1,1),QDIAG(1,1,IOZLW),ISTRIP,        IF(IOZLW.GT.0)CALL PSTBMP(OZL(1,1),QDIAG(1,1,IOZLW,bi,bj),ISTRIP,
287       .                                                      im*jm,lm,NN)       .                                                      im*jm,lm,NN)
288    
289  C **********************************************************************  C **********************************************************************
# Line 310  C ************************************** Line 310  C **************************************
310        if(itgrlw.ne.0) then        if(itgrlw.ne.0) then
311        do j = 1,jm        do j = 1,jm
312        do i = 1,im        do i = 1,im
313        qdiag(i,j,itgrlw) = qdiag(i,j,itgrlw) + tgz(i,j)        qdiag(i,j,itgrlw,bi,bj) = qdiag(i,j,itgrlw,bi,bj) + tgz(i,j)
314        enddo        enddo
315        enddo        enddo
316        endif        endif
# Line 319  C ************************************** Line 319  C **************************************
319        do L = 1,lm        do L = 1,lm
320        do j = 1,jm        do j = 1,jm
321        do i = 1,im        do i = 1,im
322        qdiag(i,j,itlw+L-1) = qdiag(i,j,itlw+L-1) + tz(i,j,L)*pkz(i,j,L)        qdiag(i,j,itlw+L-1,bi,bj) = qdiag(i,j,itlw+L-1,bi,bj) +
323         .                                             tz(i,j,L)*pkz(i,j,L)
324        enddo        enddo
325        enddo        enddo
326        enddo        enddo
# Line 329  C ************************************** Line 330  C **************************************
330        do L = 1,lm        do L = 1,lm
331        do j = 1,jm        do j = 1,jm
332        do i = 1,im        do i = 1,im
333        qdiag(i,j,ishrad+L-1) = qdiag(i,j,ishrad+L-1) + qz(i,j,L)*1000        qdiag(i,j,ishrad+L-1,bi,bj) = qdiag(i,j,ishrad+L-1,bi,bj) +
334         .                                             qz(i,j,L)*1000
335        enddo        enddo
336        enddo        enddo
337        enddo        enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22