| 1 |
jscott |
1.1 |
|
| 2 |
|
|
#include "ctrparam.h" |
| 3 |
|
|
|
| 4 |
|
|
! ========================================================== |
| 5 |
|
|
! |
| 6 |
|
|
! FORCEDOZONE.F: Routine to get ozone deviations from |
| 7 |
|
|
! specified file and interpolate |
| 8 |
|
|
! to the appropriate data |
| 9 |
|
|
! |
| 10 |
|
|
! ---------------------------------------------------------- |
| 11 |
|
|
! |
| 12 |
|
|
! Revision History: |
| 13 |
|
|
! |
| 14 |
|
|
! When Who What |
| 15 |
|
|
! ---- ---------- ------- |
| 16 |
|
|
! 073100 Chien Wang repack based on CliChem3 & M24x11, |
| 17 |
|
|
! and add cpp. |
| 18 |
|
|
! |
| 19 |
|
|
! 110100 Added data for Jan. 2000, |
| 20 |
|
|
! dimension of o3f was changed from nm to nm+1 |
| 21 |
|
|
! |
| 22 |
|
|
! 011503 Chris Forest Converted to use GISS 1850-2050 ozone data |
| 23 |
|
|
! uses nlev = 14 to adjust ozone profile in three |
| 24 |
|
|
! layers above 10 hPa top of dynamics levels |
| 25 |
|
|
! (layer edges at 5, 2, and 0.00001 hPa) |
| 26 |
|
|
! |
| 27 |
|
|
! ========================================================== |
| 28 |
|
|
|
| 29 |
|
|
subroutine getforcedozone( jyear, jday ) |
| 30 |
|
|
|
| 31 |
|
|
parameter (nm=1728) |
| 32 |
|
|
c parameter (nm =252) |
| 33 |
|
|
|
| 34 |
|
|
#include "chem_para" |
| 35 |
|
|
#include "chem_com" |
| 36 |
|
|
|
| 37 |
|
|
common /o3data/o3_data |
| 38 |
|
|
character * 120 o3_data |
| 39 |
|
|
C |
| 40 |
|
|
c Routine to get ozone deviations from specified file and interpolate |
| 41 |
|
|
c to the appropriate date. |
| 42 |
|
|
C |
| 43 |
|
|
C Author: Chris E Forest. |
| 44 |
|
|
C Added Mar 18, 1998 |
| 45 |
|
|
C |
| 46 |
|
|
C tropmass = 28.97296245*1.e-3*0.8/P0 |
| 47 |
|
|
C pxxx = dpl(l)*part |
| 48 |
|
|
C |
| 49 |
|
|
C ULGAS(L,3)=o3(ILON,jyyy,l)/48.0 |
| 50 |
|
|
C & *pxxx*tropmass |
| 51 |
|
|
C |
| 52 |
|
|
C o3 = ppb(m) |
| 53 |
|
|
C |
| 54 |
|
|
C 48 = mol weight of o3 |
| 55 |
|
|
C |
| 56 |
|
|
C ULGAS = cm^3 (STP)/cm^2 |
| 57 |
|
|
|
| 58 |
|
|
logical first |
| 59 |
|
|
integer year0, nlevp3 ,yrmax |
| 60 |
|
|
! real date, o3f(nlat,nlev+3,nm), o3dummy ! 03/20/2006 |
| 61 |
|
|
real date, o3f(nlat,nlev+3,nm+12), o3dummy |
| 62 |
|
|
c real date, o3f(24,11,nm+1), o3dummy |
| 63 |
|
|
character*43 filen |
| 64 |
|
|
DIMENSION JDOFM(13) |
| 65 |
|
|
data first / .true. / |
| 66 |
|
|
DATA JDOFM /0,31,59,90,120,151,181,212,243,273,304,334,365/ |
| 67 |
|
|
data year0 / 1859 / |
| 68 |
|
|
|
| 69 |
|
|
c year0 = 1974 for o3hadley.dat |
| 70 |
|
|
c year0 = 1979 for o3giss.MR.dat |
| 71 |
|
|
c year0 = 1850 for o3.1850_2050.giss.46x11.dat |
| 72 |
|
|
c year0 = 1859 for o3.1859_2002.giss.46x11.dat |
| 73 |
|
|
c |
| 74 |
|
|
C filen = '/data11/ceforest/forcing/dummyozone.dat' |
| 75 |
|
|
c filen = '/data11/ceforest/forcing/o3giss.MR.dat' |
| 76 |
|
|
C filen = '/data11/ceforest/forcing/o3hadley.dat' |
| 77 |
|
|
|
| 78 |
|
|
if (first) then |
| 79 |
|
|
print *,'GetForcedOzone: reading data for 14 layers' |
| 80 |
|
|
print *,'nlat=',nlat,' nlev=',nlev |
| 81 |
|
|
open(unit=500, file=o3_data,status='old') |
| 82 |
|
|
c |
| 83 |
|
|
c GISS: nm = 252 , HadCM2: nm=265, 1850-2050: nm=2412 |
| 84 |
|
|
c |
| 85 |
|
|
c nlev = nlev + 3 to include 3 layers above 10 hPa. |
| 86 |
|
|
c |
| 87 |
|
|
|
| 88 |
|
|
nlevp3 = nlev + 3 |
| 89 |
|
|
|
| 90 |
|
|
do i=1,nm |
| 91 |
|
|
do k=1, nlevp3 |
| 92 |
|
|
do j = 1,nlat |
| 93 |
|
|
read(500,'(E12.4)') o3dummy |
| 94 |
|
|
o3f(j,k,i) = o3dummy*1.e9 ! convert to ppb(m) |
| 95 |
|
|
enddo |
| 96 |
|
|
enddo |
| 97 |
|
|
enddo |
| 98 |
|
|
close(500) |
| 99 |
|
|
C Added data for Jan. 2000 to allow simulation until Jan 1 2000 |
| 100 |
|
|
c NOT necessary with 1850-2050 GISS ozone data |
| 101 |
|
|
c do k=1, nlev |
| 102 |
|
|
c do j = 1,nlat |
| 103 |
|
|
c o3f(j,k,nm+1) = o3f(j,k,nm) |
| 104 |
|
|
c enddo |
| 105 |
|
|
c enddo |
| 106 |
|
|
! 03/20/2006 Data for anaditional year are added |
| 107 |
|
|
! really only data for January are used |
| 108 |
|
|
do m=1, 12 |
| 109 |
|
|
do k=1, nlev |
| 110 |
|
|
do j = 1,nlat |
| 111 |
|
|
o3f(j,k,nm+m) = o3f(j,k,nm-12+m) |
| 112 |
|
|
enddo |
| 113 |
|
|
enddo |
| 114 |
|
|
enddo |
| 115 |
|
|
first=.false. |
| 116 |
|
|
nyears=nm/12 |
| 117 |
|
|
yrmax=year0+nyears-1 |
| 118 |
|
|
print *,'Ozone data for ',nyears,' years ' |
| 119 |
|
|
print *,'from ',year0,' to ',yrmax |
| 120 |
|
|
endif |
| 121 |
|
|
|
| 122 |
|
|
|
| 123 |
|
|
! iy = jyear - year0 ! 03/17/06 |
| 124 |
|
|
jyear1=min(jyear,yrmax) |
| 125 |
|
|
iy = jyear1 - year0 |
| 126 |
|
|
! print *,'Ozone for year ',jyear1,iy |
| 127 |
|
|
if (iy.lt.0) then |
| 128 |
|
|
do i=1,nlat |
| 129 |
|
|
do j=1,nlevp3 |
| 130 |
|
|
o3dev(1,i,j) = 0. |
| 131 |
|
|
enddo |
| 132 |
|
|
enddo |
| 133 |
|
|
else |
| 134 |
|
|
C--- o3interp(o3f, o3dev, date, year0) |
| 135 |
|
|
|
| 136 |
|
|
jdd = jday-15 |
| 137 |
|
|
|
| 138 |
|
|
if (iy .eq. 0 .and. jdd.le.0) then |
| 139 |
|
|
do i=1,nlat |
| 140 |
|
|
do j=1,nlevp3 |
| 141 |
|
|
dl = ( 31. + float(jdd) )/31. |
| 142 |
|
|
o3dev(1,i,j) = o3f(i,j,1)*dl |
| 143 |
|
|
enddo |
| 144 |
|
|
enddo |
| 145 |
|
|
else |
| 146 |
|
|
im =12 |
| 147 |
|
|
do i=1,12 |
| 148 |
|
|
if (jdd.gt.jdofm(i) .and. jdd.le.jdofm(i+1)) then |
| 149 |
|
|
im = i |
| 150 |
|
|
endif |
| 151 |
|
|
enddo |
| 152 |
|
|
if (iy .ge. 1 .and. im.eq.12 .and. jdd.le.0 ) iy = iy-1 |
| 153 |
|
|
C if (jdd.le.0) iy = iy-1 |
| 154 |
|
|
C endif |
| 155 |
|
|
|
| 156 |
|
|
dl = ( float(jdd) - jdofm(im))/(jdofm(im+1) - jdofm(im)) |
| 157 |
|
|
if (im.eq.12.and.jdd.le.0) dl = ( 31. + float(jdd) )/31. |
| 158 |
|
|
|
| 159 |
|
|
imm = im + 12*iy |
| 160 |
|
|
! print *,jyear1,jday,iy,im,imm |
| 161 |
|
|
|
| 162 |
|
|
if(imm.lt.1.or.imm.ge.nm+1)then |
| 163 |
|
|
print *,' error in o3interp im=',imm |
| 164 |
|
|
stop 25 |
| 165 |
|
|
endif |
| 166 |
|
|
|
| 167 |
|
|
do j=1,nlat |
| 168 |
|
|
do k=1,nlevp3 |
| 169 |
|
|
o3dev(1,j,k) = o3f(j,k,imm)*(1.-dl)+o3f(j,k,imm+1)*dl |
| 170 |
|
|
enddo |
| 171 |
|
|
enddo |
| 172 |
|
|
C print *,jyear,jday,iy,im,dl,o3dev(0,0,9) |
| 173 |
|
|
endif |
| 174 |
|
|
endif |
| 175 |
|
|
|
| 176 |
|
|
c print *,'GetForcedOzone: JYEAR, JDAY, IY=', jyear,jday, iy |
| 177 |
|
|
c if (iy.ge.0) print *,'GetForcedOzone: imm=', imm |
| 178 |
|
|
|
| 179 |
|
|
return |
| 180 |
|
|
end |
| 181 |
|
|
|
| 182 |
|
|
|
Parent Directory
| 
Revision Log
| 
Revision Graph