/[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.11 by molod, Tue Jun 15 21:18:18 2004 UTC revision 1.12 by molod, Tue Jul 13 21:11:08 2004 UTC
# Line 33  c--------------------------------------- Line 33  c---------------------------------------
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         _RL 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)         _RL fracland(sNx,sNy,Nsx,Nsy)
47    
        idim1 = 1-OLx  
48         idim2 = sNx+OLx         idim2 = sNx+OLx
        jdim1 = 1-OLy  
49         jdim2 = sNy+OLy         jdim2 = sNy+OLy
50         im1 = 1         im1 = 1
51         im2 = sNx         im2 = sNx
# Line 68  c--------------------------------------- Line 67  c---------------------------------------
67         enddo         enddo
68         enddo         enddo
69         call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,         call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,
70       .                 nSx*nPx,nSy*nPy,maxtyp,nchp,lons,lats,       .            nSx*nPx,nSy*nPy,maxtyp,nchp,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
# Line 79  C And now fill the earth export landtype Line 78  C And now fill the earth export landtype
78          do j=jm1,jm2          do j=jm1,jm2
79          do i=im1,im2          do i=im1,im2
80           landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)           landtype(i,j,bi,bj) = surftype(i,j,1,bi,bj)
81           if(fracland(i,j).ge.0.3.and.surftype(i,j,1,bi,bj).ge.100)           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)       .            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           if(sice(i,j,bi,bj).ne.0.0)landtype(i,j,bi,bj) = 101
84          enddo          enddo
# Line 97  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.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22