/[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.9 by molod, Wed Jun 9 20:33:37 2004 UTC revision 1.14 by molod, Fri Jul 23 22:32:28 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"  #include "gridalt_mapping.h"
35  #include "GRID.h"  #include "GRID.h"
36    #include "PARAMS.h"
37    
38         integer myThid         integer myThid
39    
40         integer i,j,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)         _RL pressure(Nrphys)
45         real lats(sNx,sNy,nSx,nSy), lons(sNx,sNy,nSx,nSy)         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
46           _RL fracland(sNx,sNy,Nsx,Nsy)
47    
48         im1 = 1-OLx         idim2 = sNx+OLx
49         im2 = sNx+OLx         jdim2 = sNy+OLy
50         jm1 = 1-OLy         im1 = 1
51         jm2 = sNy+OLy         im2 = sNx
52         idim1 = 1         jm1 = 1
53         idim2 = sNx         jm2 = sNy
        jdim1 = 1  
        jdim2 = sNy  
54         nymdb = nymd0         nymdb = nymd0
55         nhmsb = nhms0         nhmsb = nhms0
56    
# Line 65  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,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
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)  C Compute pressure profile to get methane and n2o values (bottom-up)
90    
# Line 79  C Compute pressure profile to get methan Line 96  C Compute pressure profile to get methan
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         return         return
102         end         end
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22