/[MITgcm]/MITgcm/pkg/fizhi/fizhi_init_veg.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_init_veg.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 18:32:21 2004 UTC revision 1.5 by molod, Mon Jun 7 20:26:43 2004 UTC
# Line 47  C                 10:  DARK DESERT Line 47  C                 10:  DARK DESERT
47  C                100:  OCEAN  C                100:  OCEAN
48  C***********************************************************************  C***********************************************************************
49        implicit none        implicit none
50    #include "EEPARAMS.h"
51    
52        integer mythid,im,jm,maxtyp,nchpmax,Nsx,Nsy        integer mythid,im,jm,maxtyp,nchpmax,Nsx,Nsy
53        integer surftype(im,jm,maxtyp,Nsx,Nsy)        integer surftype(im,jm,maxtyp,Nsx,Nsy)
# Line 55  C*************************************** Line 56  C***************************************
56        real chfr(nchpmax,Nsx,Nsy)        real chfr(nchpmax,Nsx,Nsy)
57        character*40 vegdata        character*40 vegdata
58        integer imdata,jmdata,Nsxdata,Nsydata        integer imdata,jmdata,Nsxdata,Nsydata
59        integer nchp,nchpland        integer nchp,nchpland,biglobal,bjglobal
60    
61        integer*4 im_32, jm_32, Nsx_32, Nsy_32        integer*4 im_32, jm_32, Nsx_32, Nsy_32
62        integer*4 iveg_32(im,jm,maxtyp,Nsx,Nsy)        integer*4 iveg_32(im,jm,maxtyp,Nsx,Nsy)
# Line 92  C*************************************** Line 93  C***************************************
93         call exit(102)         call exit(102)
94        ENDIF        ENDIF
95                                                                                                                                                                    
96        DO BJ = 1,Nsx        DO BJ = myByLo(myThid), myByHi(myThid)
97        DO BI = 1,Nsy        DO BI = myBxLo(myThid), myBxHi(myThid)
98    
99          biglobal=bi+(myXGlobalLo-1)/im
100          bjglobal=bj+(myYGlobalLo-1)/jm
101                                                                                                                                                                    
102        do k = 1,maxtyp        do k = 1,maxtyp
103        do j = 1,jm        do j = 1,jm
104        do i = 1,im        do i = 1,im
105         surftype(i,j,k,bi,bj) = iveg_32(i,j,k,bi,bj)         surftype(i,j,k,bi,bj) = iveg_32(i,j,k,biglobal,bjglobal)
106         tilefrac(i,j,k,bi,bj) = veg_32(i,j,k,bi,bj)         tilefrac(i,j,k,bi,bj) = veg_32(i,j,k,biglobal,bjglobal)
107        enddo        enddo
108        enddo        enddo
109        enddo        enddo
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22