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 |
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 |
|
|