90 |
|
|
91 |
pressure(1)=1000. |
pressure(1)=1000. |
92 |
do L = 2,Nrphys+1 |
do L = 2,Nrphys+1 |
93 |
pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1) |
pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100. |
94 |
enddo |
enddo |
95 |
|
|
96 |
call fizhi_init_chem(mythid, |
call fizhi_init_chem(mythid, |
98 |
. nlatsq,nlevsq,ntimesq,latsq,levsq,stratq, |
. nlatsq,nlevsq,ntimesq,latsq,levsq,stratq, |
99 |
. Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22) |
. Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22) |
100 |
|
|
101 |
|
do bj = myByLo(myThid), myByHi(myThid) |
102 |
|
do bi = myBxLo(myThid), myBxHi(myThid) |
103 |
|
do j=jm1,jm2 |
104 |
|
do i=im1,im2 |
105 |
|
phis_var(i,j,bi,bj) = 0. |
106 |
|
enddo |
107 |
|
enddo |
108 |
|
enddo |
109 |
|
enddo |
110 |
|
|
111 |
return |
return |
112 |
end |
end |