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 |
|
|
35 |
integer myThid |
integer myThid |
36 |
|
|
38 |
integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2 |
integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2 |
39 |
integer nymdb,nhmsb |
integer nymdb,nhmsb |
40 |
character*40 vegdata |
character*40 vegdata |
41 |
|
_RL pressure(Nrphys) |
42 |
|
|
43 |
im1 = 1-OLx |
im1 = 1-OLx |
44 |
im2 = sNx+OLx |
im2 = sNx+OLx |
52 |
nhmsb = nhms0 |
nhmsb = nhms0 |
53 |
|
|
54 |
call fizhi_alarms(nymdb,nhmsb,deltaTClock) |
call fizhi_alarms(nymdb,nhmsb,deltaTClock) |
55 |
call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,maxtyp, |
call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy, |
56 |
. surftype,tilefrac,igrd,ityp,chfr ) |
. nSx*nPx,nSy*nPy,maxtyp,nchp,surftype,tilefrac,igrd,ityp,chfr ) |
57 |
|
|
58 |
C Compute pressure profile to get methane and n2o values (bottom-up) |
C Compute pressure profile to get methane and n2o values (bottom-up) |
|
do bj = myByLo(myThid), myByHi(myThid) |
|
|
do bi = myBxLo(myThid), myBxHi(myThid) |
|
59 |
|
|
60 |
enddo |
pressure(1)=1000. |
61 |
|
do L = 2,Nrphys+1 |
62 |
|
pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1) |
63 |
enddo |
enddo |
64 |
|
|
65 |
call fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs, |
call fizhi_init_chem(mythid, |
66 |
. ozone,nwatlats,nwatlevs,watlats,watlevs,water, |
. nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone, |
67 |
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
. nlatsq,nlevsq,ntimesq,latsq,levsq,stratq, |
68 |
|
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
69 |
|
|
70 |
return |
return |
71 |
end |
end |