/[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.4 by molod, Mon Jun 7 21:36:45 2004 UTC revision 1.13 by molod, Thu Jul 22 22:33:55 2004 UTC
# Line 29  c--------------------------------------- Line 29  c---------------------------------------
29  #include "fizhi_chemistry_coms.h"  #include "fizhi_chemistry_coms.h"
30  #include "fizhi_earth_coms.h"  #include "fizhi_earth_coms.h"
31  #include "fizhi_land_coms.h"  #include "fizhi_land_coms.h"
32    #include "fizhi_ocean_coms.h"
33  #include "chronos.h"  #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         integer nymdb,nhmsb
43         character*40 vegdata         character*40 vegdata
44           _RL pressure(Nrphys)
45         im1 = 1-OLx         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
46         im2 = sNx+OLx         _RL fracland(sNx,sNy,Nsx,Nsy)
47         jm1 = 1-OLy  
48         jm2 = sNy+OLy         idim2 = sNx+OLx
49         idim1 = 1         jdim2 = sNy+OLy
50         idim2 = sNx         im1 = 1
51         jdim1 = 1         im2 = sNx
52         jdim2 = sNy         jm1 = 1
53           jm2 = sNy
54         nymdb = nymd0         nymdb = nymd0
55         nhmsb = nhms0         nhmsb = nhms0
56    
57         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
        call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,maxtyp,  
      .            surftype,tilefrac,igrd,ityp,chfr )  
58    
 C Compute pressure profile to get methane and n2o values (bottom-up)  
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           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,idim2,jdim2,Nsx,Nsy,
70         .            nSx*nPx,nSy*nPy,maxtyp,nchp,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         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)
94         enddo         enddo
95    
96         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
97       .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,       .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
98       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
99       .      Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22)       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
100    
101           print *,' Back from init chem '
102    
103         return         return
104         end         end
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22