2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
subroutine fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs, |
subroutine fizhi_init_chem(mythid,nozlats,nozlevs,ozlats,ozlevs, |
5 |
. ozone,nwatlats,nwatlevs,watlats,watlevs,water, |
. o3,nwatlats,nwatlevs,watlats,watlevs,water, |
6 |
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
. Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22) |
7 |
C*********************************************************************** |
C*********************************************************************** |
8 |
C Subroutine fizhi_init_chem - routine to read in the ozone and upper |
C Subroutine fizhi_init_chem - routine to read in the ozone and upper |
21 |
implicit none |
implicit none |
22 |
#include "CPP_EEOPTIONS.h" |
#include "CPP_EEOPTIONS.h" |
23 |
integer mythid,nozlevs,nozlats,nwatlevs,nwatlats,Nrphys |
integer mythid,nozlevs,nozlats,nwatlevs,nwatlats,Nrphys |
24 |
_RL oz(nozlats,nozlevs), water(nwatlats,nwatlevs) |
_RL o3(nozlats,nozlevs,12), water(nwatlats,nwatlevs,12) |
25 |
_RL ozlats(nozlats), ozlevs(nozlevs) |
_RL ozlats(nozlats), ozlevs(nozlevs) |
26 |
_RL watlats(nwatlats), watlevs(nwatlevs) |
_RL watlats(nwatlats), watlevs(nwatlevs) |
27 |
_RL pressure(Nrphys),methane(Nrphys),n2o(Nrphys) |
_RL pressure(Nrphys),methane(Nrphys),n2o(Nrphys) |
169 |
|
|
170 |
do time=1,ntime |
do time=1,ntime |
171 |
do lev=1,nlev |
do lev=1,nlev |
172 |
nrec = lev+(time-1)*nlev*2 ! Note: 2 quantities in Ozone Dataset |
C Note: 2 quantities in Ozone Dataset |
173 |
|
nrec = lev+(time-1)*nlev*2 |
174 |
read(ku,rec=nrec) o3 |
read(ku,rec=nrec) o3 |
175 |
do lat=1,nlat |
do lat=1,nlat |
176 |
oz(lat,nlev-lev+1,time) = o3(lat) |
oz(lat,nlev-lev+1,time) = o3(lat) |