41 |
integer im1, im2, jm1, jm2, idim2, jdim2 |
integer im1, im2, jm1, jm2, idim2, jdim2 |
42 |
integer nymdb,nhmsb |
integer nymdb,nhmsb |
43 |
character*40 vegdata |
character*40 vegdata |
44 |
|
_RL pressure0(Nrphys) |
45 |
_RL pressure(Nrphys) |
_RL pressure(Nrphys) |
46 |
_RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy) |
_RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy) |
47 |
_RL fracland(sNx,sNy,Nsx,Nsy) |
_RL fracland(sNx,sNy,Nsx,Nsy) |
89 |
|
|
90 |
C Compute pressure profile to get methane and n2o values (bottom-up) |
C Compute pressure profile to get methane and n2o values (bottom-up) |
91 |
|
|
92 |
pressure(1)=1000. |
pressure0(1)=1000. |
93 |
do L = 2,Nrphys+1 |
do L = 2,Nrphys+1 |
94 |
pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100. |
pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100. |
95 |
|
enddo |
96 |
|
C Now flip pressure to a top-down array to send to fizhi |
97 |
|
do L = 1,Nrphys+1 |
98 |
|
pressure(L)=pressure0(Nrphys+2-L) |
99 |
enddo |
enddo |
100 |
|
|
101 |
call fizhi_init_chem(mythid, |
call fizhi_init_chem(mythid, |