/[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.6 by molod, Tue Jun 8 22:26:08 2004 UTC revision 1.9 by molod, Wed Jun 9 20:33:37 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           real 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    
58           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,         call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,
69       .  nSx*nPx,nSy*nPy,maxtyp,nchp,surftype,tilefrac,igrd,ityp,chfr )       .                 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.6  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22