/[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.1 by molod, Fri Jun 4 16:23:37 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 17  c                   now and interpolate Line 18  c                   now and interpolate
18  c      the other fixed parameters are:  c      the other fixed parameters are:
19  c         N2O, Methane (vary with space)  c         N2O, Methane (vary with space)
20  c         CO2, CFC11, CFC12, CFC22 (set to a global value)  c         CO2, CFC11, CFC12, CFC22 (set to a global value)
21  c   2) This routine will also read the data.fizhi namelist interface  c   3) For now, the fizhi package contains the alarms and clocks
22  c      to obtain information about diagnostic selection.  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 "land_SIZE.h"  #include "fizhi_land_SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29  #include "chemistry_coms.h"  #include "fizhi_chemistry_coms.h"
30  #include "earth_coms.h"  #include "fizhi_earth_coms.h"
31    #include "fizhi_land_coms.h"
32    #include "fizhi_ocean_coms.h"
33    #include "chronos.h"
34    #include "gridalt_mapping.h"
35    #include "GRID.h"
36    #include "PARAMS.h"
37    
38         integer myThid         integer myThid
39    
40         integer bi, bj         integer i,j,L,bi,bj
41         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim2, jdim2
42           integer nymdb,nhmsb
43           character*40 vegdata
44           _RL pressure(Nrphys)
45           _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
46           _RL fracland(sNx,sNy,Nsx,Nsy)
47    
48           idim2 = sNx+OLx
49           jdim2 = sNy+OLy
50           im1 = 1
51           im2 = sNx
52           jm1 = 1
53           jm2 = sNy
54           nymdb = nymd0
55           nhmsb = nhms0
56    
57         im1 = 1-OLx         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
        im2 = sNx+OLx  
        jm1 = 1-OLy  
        jm2 = sNy+OLy  
        idim1 = 1  
        idim2 = sNx  
        jdim1 = 1  
        jdim2 = sNy  
58    
59         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
60         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
61  C Dummy routine - nothing to do for now         do j = jm1,jm2
62           do i = im1,im2
63            lons(i,j,bi,bj) = xC(i,j,bi,bj)
64            lats(i,j,bi,bj) = yC(i,j,bi,bj)
65           enddo
66           enddo
67           enddo
68           enddo
69           call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,
70         .         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,
71         .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
72    
73    C And now fill the earth export landtype
74           do bj = myByLo(myThid), myByHi(myThid)
75           do bi = myBxLo(myThid), myBxHi(myThid)
76            call get_landfrac(im2,jm2,Nsx,Nsy,bi,bj,maxtyp,
77         .        surftype,tilefrac,fracland)
78            do j=jm1,jm2
79            do i=im1,im2
80             landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)
81             if(fracland(i,j,bi,bj).ge.0.3.and.surftype(i,j,1,bi,bj).ge.100)
82         .            landtype(i,j,bi,bj) = surftype(i,j,2,bi,bj)
83             if(sice(i,j,bi,bj).ne.0.0)landtype(i,j,bi,bj) = 101
84            enddo
85            enddo
86           enddo
87           enddo
88    
89    C Compute pressure profile to get methane and n2o values (bottom-up)
90    
91           pressure(1)=1000.
92           do L = 2,Nrphys+1
93            pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100.
94           enddo
95    
96           call fizhi_init_chem(mythid,
97         .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
98         .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
99         .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
100    
101           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         enddo
109         enddo         enddo
110    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22