/[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.2 by molod, Thu Jun 24 18:56:57 2004 UTC revision 1.7 by molod, Mon Jul 26 19:51:08 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3                                                                                          
4  #include "CPP_OPTIONS.h"  #include "FIZHI_OPTIONS.h"
5        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,
6       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,
7       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,
8       .  lwdt,swdt,lwdtclr,swdtclr,im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,bi,bj)       .  lwdt,swdt,lwdtclr,swdtclr,dlwdtg,
9         .  im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer)
10  C***********************************************************************  C***********************************************************************
11        implicit none        implicit none
12    
13    #ifdef ALLOW_DIAGNOSTICS
14    #include "SIZE.h"
15    #include "diagnostics_SIZE.h"
16  #include "diagnostics.h"  #include "diagnostics.h"
17    #endif
18    
19        integer myThid,im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,bi,bj        integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer
20        real p(im2,jm2,Nsx,Nsy)        _RL p(im2,jm2,Nbi,Nbj)
21        real uphy(im2,jm2,Nrphys,Nsx,Nsy),vphy(im2,jm2,Nrphys,Nsx,Nsy)        _RL uphy(im2,jm2,Nrphys,Nbi,Nbj)
22        real thphy(im2,jm2,Nrphys,Nsx,Nsy),sphy(im2,jm2,Nrphys,Nsx,Nsy)        _RL vphy(im2,jm2,Nrphys,Nbi,Nbj)
23        real qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nsx,Nsy)        _RL thphy(im2,jm2,Nrphys,Nbi,Nbj)
24        real dp(im2,jm2,Nrphys,Nsx,Nsy)        _RL sphy(im2,jm2,Nrphys,Nbi,Nbj)
25        real radswt(im2,jm2,Nsx,Nsy),radswg(im2,jm2,Nsx,Nsy)        _RL qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj)
26        real swgclr(im2,jm2,Nsx,Nsy),osr(im2,jm2,Nsx,Nsy)        _RL dp(im2,jm2,Nrphys,Nbi,Nbj)
27        real osrclr(im2,jm2,Nsx,Nsy),st4(im2,jm2,Nsx,Nsy)        _RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)
28        real dst4(im2,jm2,Nsx,Nsy),tgz(im2,jm2,Nsx,Nsy)        _RL swgclr(im2,jm2,Nbi,Nbj),osr(im2,jm2,Nbi,Nbj)
29        real tg0(im2,jm2,Nsx,Nsy),radlwg(im2,jm2,Nsx,Nsy)        _RL osrclr(im2,jm2,Nbi,Nbj),st4(im2,jm2,Nbi,Nbj)
30        real lwgclr(im2,jm2,Nsx,Nsy)        _RL dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj)
31        real turbu(im2,jm2,Nrphys,Nsx,Nsy),turbv(im2,jm2,Nrphys,Nsx,Nsy)        _RL tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj)
32        real turbt(im2,jm2,Nrphys,Nsx,Nsy),turbq(im2,jm2,Nrphys,Nsx,Nsy)        _RL lwgclr(im2,jm2,Nbi,Nbj)
33        real moistu(im2,jm2,Nrphys,Nsx,Nsy),moistv(im2,jm2,Nrphys,Nsx,Nsy)        _RL turbu(im2,jm2,Nrphys,Nbi,Nbj)
34        real moistt(im2,jm2,Nrphys,Nsx,Nsy),moistq(im2,jm2,Nrphys,Nsx,Nsy)        _RL turbv(im2,jm2,Nrphys,Nbi,Nbj)
35        real lwdt(im2,jm2,Nrphys,Nsx,Nsy),swdt(im2,jm2,Nrphys,Nsx,Nsy)        _RL turbt(im2,jm2,Nrphys,Nbi,Nbj)
36        real lwdtclr(im2,jm2,Nrphys,Nsx,Nsy)        _RL turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
37        real swdtclr(im2,jm2,Nrphys,Nsx,Nsy)        _RL moistu(im2,jm2,Nrphys,Nbi,Nbj)
38          _RL moistv(im2,jm2,Nrphys,Nbi,Nbj)
39          _RL moistt(im2,jm2,Nrphys,Nbi,Nbj)
40          _RL moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
41          _RL lwdt(im2,jm2,Nrphys,Nbi,Nbj)
42          _RL swdt(im2,jm2,Nrphys,Nbi,Nbj)
43          _RL lwdtclr(im2,jm2,Nrphys,Nbi,Nbj)
44          _RL swdtclr(im2,jm2,Nrphys,Nbi,Nbj)
45          _RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)
46    
47        integer  i,j,L        integer  i,j,L
48        real pinv(im2,jm2), qbar(im2,jm2)        _RL pinv(im2,jm2), qbar(im2,jm2)
49    
50  C **********************************************************************          C **********************************************************************        
51    
# Line 238  c -------------------------------------- Line 251  c --------------------------------------
251        do i=im1,im2        do i=im1,im2
252        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) +
253       .           swdtclr(i,j,l,bi,bj)*radswt(i,j,bi,bj)*       .           swdtclr(i,j,l,bi,bj)*radswt(i,j,bi,bj)*
254       .                           pk(i,j,l,bi,bj)*pinv(i,j,bi,bj)*86400       .                           pk(i,j,l,bi,bj)*pinv(i,j)*86400
255        enddo        enddo
256        enddo        enddo
257        endif        endif
# Line 334  c -------------------------------------- Line 347  c --------------------------------------
347        do j=jm1,jm2        do j=jm1,jm2
348        do i=im1,im2        do i=im1,im2
349        qdiag(i,j,ivdtmoist,bi,bj) = qdiag(i,j,ivdtmoist,bi,bj) +        qdiag(i,j,ivdtmoist,bi,bj) = qdiag(i,j,ivdtmoist,bi,bj) +
350       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
351        enddo        enddo
352        enddo        enddo
353        endif        endif
# Line 358  c -------------------------------------- Line 371  c --------------------------------------
371        do j=jm1,jm2        do j=jm1,jm2
372        do i=im1,im2        do i=im1,im2
373        qdiag(i,j,ivdtturb,bi,bj) = qdiag(i,j,ivdtturb,bi,bj) +        qdiag(i,j,ivdtturb,bi,bj) = qdiag(i,j,ivdtturb,bi,bj) +
374       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
375        enddo        enddo
376        enddo        enddo
377        endif        endif
# Line 383  c -------------------------------------- Line 396  c --------------------------------------
396        do j=jm1,jm2        do j=jm1,jm2
397        do i=im1,im2        do i=im1,im2
398        qdiag(i,j,ivdtradlw,bi,bj) = qdiag(i,j,ivdtradlw,bi,bj) +        qdiag(i,j,ivdtradlw,bi,bj) = qdiag(i,j,ivdtradlw,bi,bj) +
399       .      qbar(i,j)*pinv(i,j,bi,bj)*pinv(i,j,bi,bj)*86400       .      qbar(i,j)*pinv(i,j)*pinv(i,j)*86400
400        enddo        enddo
401        enddo        enddo
402        endif        endif
# Line 407  c -------------------------------------- Line 420  c --------------------------------------
420        do j=jm1,jm2        do j=jm1,jm2
421        do i=im1,im2        do i=im1,im2
422        qdiag(i,j,ivdtradsw,bi,bj) = qdiag(i,j,ivdtradsw,bi,bj) +        qdiag(i,j,ivdtradsw,bi,bj) = qdiag(i,j,ivdtradsw,bi,bj) +
423       . 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
424        enddo        enddo
425        enddo        enddo
426        endif        endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22