/[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.5 by molod, Tue Jun 8 15:09:01 2004 UTC revision 1.10 by molod, Thu Jun 10 20:17:17 2004 UTC
# Line 31  c--------------------------------------- Line 31  c---------------------------------------
31  #include "fizhi_land_coms.h"  #include "fizhi_land_coms.h"
32  #include "chronos.h"  #include "chronos.h"
33  #include "gridalt_mapping.h"  #include "gridalt_mapping.h"
34    #include "GRID.h"
35    
36         integer myThid         integer myThid
37    
38         integer bi, bj         integer i,j,bi,bj
39         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
40         integer nymdb,nhmsb         integer nymdb,nhmsb
41         character*40 vegdata         character*40 vegdata
42         _RL pressure(Nrphys)         _RL pressure(Nrphys)
43           _RL lats(sNx,sNy,nSx,nSy), lons(sNx,sNy,nSx,nSy)
44    
45         im1 = 1-OLx         im1 = 1-OLx
46         im2 = sNx+OLx         im2 = sNx+OLx
# Line 52  c--------------------------------------- Line 54  c---------------------------------------
54         nhmsb = nhms0         nhmsb = nhms0
55    
56         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
57         call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,maxtyp,  
58       .            surftype,tilefrac,igrd,ityp,chfr )         do bj = myByLo(myThid), myByHi(myThid)
59           do bi = myBxLo(myThid), myBxHi(myThid)
60           do j = jm1,jm2
61           do i = im1,im2
62            lons(i,j,bi,bj) = xC(i,j,bi,bj)
63            lats(i,j,bi,bj) = yC(i,j,bi,bj)
64           enddo
65           enddo
66           enddo
67           enddo
68           call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,
69         .                 nSx*nPx,nSy*nPy,maxtyp,nchp,lons,lats,
70         .            surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
71    
72  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values (bottom-up)
73    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22