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) |
_RL pressure0(Nrphys+1) |
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) |
93 |
enddo |
enddo |
94 |
enddo |
enddo |
95 |
|
|
96 |
C Compute pressure profile to get methane and n2o values (bottom-up) |
C Compute pressure profile to get methane and n2o values |
97 |
|
C - First bottom-up (in mb) |
98 |
|
|
99 |
pressure0(1)=1000. |
pressure0(1)=1000. |
100 |
do L = 2,Nrphys+1 |
do L = 2,Nrphys+1 |
101 |
pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100. |
pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100. |
102 |
enddo |
enddo |
103 |
C Now flip pressure to a top-down array to send to fizhi init chem |
C Now flip pressure to a top-down array and average to mid levels |
104 |
do L = 1,Nrphys+1 |
C to send to fizhi init chem |
105 |
pressure(L)=pressure0(Nrphys+2-L) |
do L = 1,Nrphys |
106 |
|
pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2. |
107 |
enddo |
enddo |
108 |
|
|
109 |
call fizhi_init_chem(mythid, |
call fizhi_init_chem(mythid, |