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

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

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

revision 1.13 by molod, Thu Jul 22 22:33:55 2004 UTC revision 1.16 by molod, Wed Jul 28 01:25:07 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "FIZHI_OPTIONS.h"
5         subroutine fizhi_init_fixed (myThid)         subroutine fizhi_init_fixed (myThid)
6  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
7  c  Routine to initialise the fizhi package.  c  Routine to initialise the fizhi package.
# Line 21  c   3) For now, the fizhi package contai Line 22  c   3) For now, the fizhi package contai
22  c      routines, so this routine will also initialize the alarms.  c      routines, so this routine will also initialize the alarms.
23  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
24         implicit none         implicit none
 #include "CPP_OPTIONS.h"  
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "fizhi_SIZE.h"  #include "fizhi_SIZE.h"
27  #include "fizhi_land_SIZE.h"  #include "fizhi_land_SIZE.h"
# Line 66  c--------------------------------------- Line 66  c---------------------------------------
66         enddo         enddo
67         enddo         enddo
68         enddo         enddo
69         call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,         call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,
70       .            nSx*nPx,nSy*nPy,maxtyp,nchp,nchpland,lons,lats,       .         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,
71       .            surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)       .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
72    
73  C And now fill the earth export landtype  C And now fill the earth export landtype
74         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 90  C Compute pressure profile to get methan Line 90  C Compute pressure profile to get methan
90    
91         pressure(1)=1000.         pressure(1)=1000.
92         do L = 2,Nrphys+1         do L = 2,Nrphys+1
93          pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)          pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100.
94         enddo         enddo
95    
96         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
# Line 98  C Compute pressure profile to get methan Line 98  C Compute pressure profile to get methan
98       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
99       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
100    
101         print *,' Back from init chem '         do bj = myByLo(myThid), myByHi(myThid)
102           do bi = myBxLo(myThid), myBxHi(myThid)
103            do j=jm1,jm2
104            do i=im1,im2
105             phis_var(i,j,bi,bj) = 0.
106            enddo
107            enddo
108           enddo
109           enddo
110    
111         return         return
112         end         end
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22