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 |
56 |
. surftype,tilefrac,igrd,ityp,chfr ) |
. 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, |
call fizhi_init_chem(mythid, |