/[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.10 by molod, Thu Jun 10 20:17:17 2004 UTC revision 1.11 by molod, Tue Jun 15 21:18:18 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"
# Line 40  c--------------------------------------- Line 41  c---------------------------------------
41         integer nymdb,nhmsb         integer nymdb,nhmsb
42         character*40 vegdata         character*40 vegdata
43         _RL pressure(Nrphys)         _RL pressure(Nrphys)
44         _RL lats(sNx,sNy,nSx,nSy), lons(sNx,sNy,nSx,nSy)         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
45           _RL fracland(sNx,sNy,Nsx,Nsy)
46    
47         im1 = 1-OLx         idim1 = 1-OLx
48         im2 = sNx+OLx         idim2 = sNx+OLx
49         jm1 = 1-OLy         jdim1 = 1-OLy
50         jm2 = sNy+OLy         jdim2 = sNy+OLy
51         idim1 = 1         im1 = 1
52         idim2 = sNx         im2 = sNx
53         jdim1 = 1         jm1 = 1
54         jdim2 = sNy         jm2 = sNy
55         nymdb = nymd0         nymdb = nymd0
56         nhmsb = nhms0         nhmsb = nhms0
57    
# Line 69  c--------------------------------------- Line 71  c---------------------------------------
71       .                 nSx*nPx,nSy*nPy,maxtyp,nchp,lons,lats,       .                 nSx*nPx,nSy*nPy,maxtyp,nchp,lons,lats,
72       .            surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)       .            surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
73    
74    C And now fill the earth export landtype
75           do bj = myByLo(myThid), myByHi(myThid)
76           do bi = myBxLo(myThid), myBxHi(myThid)
77            call get_landfrac(im2,jm2,Nsx,Nsy,bi,bj,maxtyp,
78         .        surftype,tilefrac,fracland)
79            do j=jm1,jm2
80            do i=im1,im2
81             landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)
82             if(fracland(i,j).ge.0.3.and.surftype(i,j,1,bi,bj).ge.100)
83         .            landtype(i,j,bi,bj) = surftype(i,j,2,bi,bj)
84             if(sice(i,j,bi,bj).ne.0.0)landtype(i,j,bi,bj) = 101
85            enddo
86            enddo
87           enddo
88           enddo
89    
90  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values (bottom-up)
91    
92         pressure(1)=1000.         pressure(1)=1000.
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22