/[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.2 by edhill, Wed May 5 00:39:21 2004 UTC revision 1.3 by molod, Tue Jun 8 15:25:27 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3           subroutine do_fizhi(uphy,vphy,thphy,sphy,pephy,
        subroutine do_fizhi(uphy,vphy,thphy,sphy,  
4       .                              ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       .                              ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
5       .im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,idim1,idim2,jdim1,jdim2,bi,bj,nchp,       .im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,idim1,idim2,jdim1,jdim2,bi,bj,nchp,
6       .                                         duphy,dvphy,dthphy,dsphy)       .                                         duphy,dvphy,dthphy,dsphy)
7  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
8  c Dummy routine to calculate physics increments - here set them to whatever  c Interface routine to calculate physics increments - call geos-3 driver
9    c Purpose of this routine is to pull things from model common blocks
10    c to send to physdrv. Also to take output from there and put into model
11    c common blocks.
12  c  c
13  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
14         implicit none         implicit none
# Line 17  c--------------------------------------- Line 19  c---------------------------------------
19         _RL uphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)         _RL uphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)
20         _RL vphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)         _RL vphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)
21         _RL thphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)         _RL thphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)
22           _RL pephy(im1:im2,jm1:jm2,Nrphys+1,Nsx,Nsy)
23         _RL sphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)         _RL sphy(im1:im2,jm1:jm2,Nrphys,Nsx,Nsy)
24         _RL ctmt(nchp,Nsx,Nsy),xxmt(nchp,Nsx,Nsy),yymt(nchp,Nsx,Nsy)         _RL ctmt(nchp,Nsx,Nsy),xxmt(nchp,Nsx,Nsy),yymt(nchp,Nsx,Nsy)
25         _RL zetamt(nchp,Nsx,Nsy)         _RL zetamt(nchp,Nsx,Nsy)
# Line 29  c--------------------------------------- Line 32  c---------------------------------------
32  c  c
33         integer i,j,L         integer i,j,L
34    
35         do L = 1,Nrphys         call physdrv(         )
         do j = jdim1,jdim2  
         do i = idim1,idim2  
          duphy(i,j,L,bi,bj) = 2.e-4*L  
          dvphy(i,j,L,bi,bj) = 0.5e-4*L  
          dthphy(i,j,L,bi,bj) = 4.e-4*L  
          dsphy(i,j,L,bi,bj) = 2.e-7*L  
         enddo  
         enddo  
        enddo  
36    
37         return         return
38         end         end
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22