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

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

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

revision 1.4 by molod, Tue Jul 13 23:44:43 2004 UTC revision 1.5 by molod, Wed Jul 14 14:50:04 2004 UTC
# Line 6  C $Name$ Line 6  C $Name$
6        subroutine fizhi_step_diag(myThid,p,uphy,vphy,thphy,sphy,qq,pk,dp,        subroutine fizhi_step_diag(myThid,p,uphy,vphy,thphy,sphy,qq,pk,dp,
7       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,
8       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,
9       .  lwdt,swdt,lwdtclr,swdtclr,im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,bi,bj)       .  lwdt,swdt,lwdtclr,swdtclr,dlwdtg,
10         .  im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer)
11  C***********************************************************************  C***********************************************************************
12        implicit none        implicit none
13    
# Line 16  C*************************************** Line 17  C***************************************
17  #include "diagnostics.h"  #include "diagnostics.h"
18  #endif  #endif
19    
20        integer myThid,im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,bi,bj        integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer
21        real p(im2,jm2,Nsx,Nsy)        real p(im2,jm2,Nbi,Nbj)
22        real uphy(im2,jm2,Nrphys,Nsx,Nsy),vphy(im2,jm2,Nrphys,Nsx,Nsy)        real uphy(im2,jm2,Nrphys,Nbi,Nbj),vphy(im2,jm2,Nrphys,Nbi,Nbj)
23        real thphy(im2,jm2,Nrphys,Nsx,Nsy),sphy(im2,jm2,Nrphys,Nsx,Nsy)        real thphy(im2,jm2,Nrphys,Nbi,Nbj),sphy(im2,jm2,Nrphys,Nbi,Nbj)
24        real qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nsx,Nsy)        real qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj)
25        real dp(im2,jm2,Nrphys,Nsx,Nsy)        real dp(im2,jm2,Nrphys,Nbi,Nbj)
26        real radswt(im2,jm2,Nsx,Nsy),radswg(im2,jm2,Nsx,Nsy)        real radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)
27        real swgclr(im2,jm2,Nsx,Nsy),osr(im2,jm2,Nsx,Nsy)        real swgclr(im2,jm2,Nbi,Nbj),osr(im2,jm2,Nbi,Nbj)
28        real osrclr(im2,jm2,Nsx,Nsy),st4(im2,jm2,Nsx,Nsy)        real osrclr(im2,jm2,Nbi,Nbj),st4(im2,jm2,Nbi,Nbj)
29        real dst4(im2,jm2,Nsx,Nsy),tgz(im2,jm2,Nsx,Nsy)        real dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj)
30        real tg0(im2,jm2,Nsx,Nsy),radlwg(im2,jm2,Nsx,Nsy)        real tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj)
31        real lwgclr(im2,jm2,Nsx,Nsy)        real lwgclr(im2,jm2,Nbi,Nbj)
32        real turbu(im2,jm2,Nrphys,Nsx,Nsy),turbv(im2,jm2,Nrphys,Nsx,Nsy)        real turbu(im2,jm2,Nrphys,Nbi,Nbj),turbv(im2,jm2,Nrphys,Nbi,Nbj)
33        real turbt(im2,jm2,Nrphys,Nsx,Nsy),turbq(im2,jm2,Nrphys,Nsx,Nsy)        real turbt(im2,jm2,Nrphys,Nbi,Nbj)
34        real moistu(im2,jm2,Nrphys,Nsx,Nsy),moistv(im2,jm2,Nrphys,Nsx,Nsy)        real turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
35        real moistt(im2,jm2,Nrphys,Nsx,Nsy),moistq(im2,jm2,Nrphys,Nsx,Nsy)        real moistu(im2,jm2,Nrphys,Nbi,Nbj),moistv(im2,jm2,Nrphys,Nbi,Nbj)
36        real lwdt(im2,jm2,Nrphys,Nsx,Nsy),swdt(im2,jm2,Nrphys,Nsx,Nsy)        real moistt(im2,jm2,Nrphys,Nbi,Nbj)
37        real lwdtclr(im2,jm2,Nrphys,Nsx,Nsy)        real moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
38        real swdtclr(im2,jm2,Nrphys,Nsx,Nsy)        real lwdt(im2,jm2,Nrphys,Nbi,Nbj),swdt(im2,jm2,Nrphys,Nbi,Nbj)
39          real lwdtclr(im2,jm2,Nrphys,Nbi,Nbj)
40          real swdtclr(im2,jm2,Nrphys,Nbi,Nbj)
41          real dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)
42    
43        integer  i,j,L        integer  i,j,L
44        real pinv(im2,jm2), qbar(im2,jm2)        real pinv(im2,jm2), qbar(im2,jm2)
# Line 243  c -------------------------------------- Line 247  c --------------------------------------
247        do i=im1,im2        do i=im1,im2
248        qdiag(i,j,iswclr+l-1,bi,bj) = qdiag(i,j,iswclr+l-1,bi,bj) +        qdiag(i,j,iswclr+l-1,bi,bj) = qdiag(i,j,iswclr+l-1,bi,bj) +
249       .           swdtclr(i,j,l,bi,bj)*radswt(i,j,bi,bj)*       .           swdtclr(i,j,l,bi,bj)*radswt(i,j,bi,bj)*
250       .                           pk(i,j,l,bi,bj)*pinv(i,j,bi,bj)*86400       .                           pk(i,j,l,bi,bj)*pinv(i,j)*86400
251        enddo        enddo
252        enddo        enddo
253        endif        endif
# Line 339  c -------------------------------------- Line 343  c --------------------------------------
343        do j=jm1,jm2        do j=jm1,jm2
344        do i=im1,im2        do i=im1,im2
345        qdiag(i,j,ivdtmoist,bi,bj) = qdiag(i,j,ivdtmoist,bi,bj) +        qdiag(i,j,ivdtmoist,bi,bj) = qdiag(i,j,ivdtmoist,bi,bj) +
346       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
347        enddo        enddo
348        enddo        enddo
349        endif        endif
# Line 363  c -------------------------------------- Line 367  c --------------------------------------
367        do j=jm1,jm2        do j=jm1,jm2
368        do i=im1,im2        do i=im1,im2
369        qdiag(i,j,ivdtturb,bi,bj) = qdiag(i,j,ivdtturb,bi,bj) +        qdiag(i,j,ivdtturb,bi,bj) = qdiag(i,j,ivdtturb,bi,bj) +
370       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
371        enddo        enddo
372        enddo        enddo
373        endif        endif
# Line 388  c -------------------------------------- Line 392  c --------------------------------------
392        do j=jm1,jm2        do j=jm1,jm2
393        do i=im1,im2        do i=im1,im2
394        qdiag(i,j,ivdtradlw,bi,bj) = qdiag(i,j,ivdtradlw,bi,bj) +        qdiag(i,j,ivdtradlw,bi,bj) = qdiag(i,j,ivdtradlw,bi,bj) +
395       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
396        enddo        enddo
397        enddo        enddo
398        endif        endif
# Line 412  c -------------------------------------- Line 416  c --------------------------------------
416        do j=jm1,jm2        do j=jm1,jm2
417        do i=im1,im2        do i=im1,im2
418        qdiag(i,j,ivdtradsw,bi,bj) = qdiag(i,j,ivdtradsw,bi,bj) +        qdiag(i,j,ivdtradsw,bi,bj) = qdiag(i,j,ivdtradsw,bi,bj) +
419       . qbar(i,j)*radswt(i,j,bi,bj)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       . qbar(i,j)*radswt(i,j,bi,bj)*pinv(i,j)*pinv(i,j)*86400
420        enddo        enddo
421        enddo        enddo
422        endif        endif
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22