30 |
#include "fizhi_earth_coms.h" |
#include "fizhi_earth_coms.h" |
31 |
#include "fizhi_land_coms.h" |
#include "fizhi_land_coms.h" |
32 |
#include "chronos.h" |
#include "chronos.h" |
33 |
|
#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) |
43 |
|
_RL 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) |
|
call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,maxtyp, |
|
|
. surftype,tilefrac,igrd,ityp,chfr ) |
|
57 |
|
|
|
C Compute pressure profile to get methane and n2o values (bottom-up) |
|
58 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
59 |
do bi = myBxLo(myThid), myBxHi(myThid) |
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 |
enddo |
67 |
enddo |
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) |
73 |
|
|
74 |
|
pressure(1)=1000. |
75 |
|
do L = 2,Nrphys+1 |
76 |
|
pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1) |
77 |
|
enddo |
78 |
|
|
79 |
call fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs, |
call fizhi_init_chem(mythid, |
80 |
. ozone,nwatlats,nwatlevs,watlats,watlevs,water, |
. nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone, |
81 |
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
. nlatsq,nlevsq,ntimesq,latsq,levsq,stratq, |
82 |
|
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
83 |
|
|
84 |
return |
return |
85 |
end |
end |