126 |
|
|
127 |
character*8 cname |
character*8 cname |
128 |
character*80 cdscrip |
character*80 cdscrip |
129 |
|
character*40 sicedata |
130 |
_RL fac1, fac2, lat0, lon0, timebc1, timebc2, timemod, undef |
_RL fac1, fac2, lat0, lon0, timebc1, timebc2, timemod, undef |
131 |
logical first, found, error |
logical first, found, error |
132 |
integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc,npxbc,npybc |
integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc,npxbc,npybc |
140 |
|
|
141 |
data first /.true./ |
data first /.true./ |
142 |
data error /.false./ |
data error /.false./ |
143 |
|
sicedata='sice19232.weekly.clim' |
144 |
|
|
145 |
c save header info |
c save header info |
146 |
save imbc,jmbc,npxbc,npybc,lat0,lon0,ndatebc,undef,nymdbc,nhmsbc |
save imbc,jmbc,npxbc,npybc,lat0,lon0,ndatebc,undef,nymdbc,nhmsbc |
178 |
if( iyear.ne.iyearbc ) then |
if( iyear.ne.iyearbc ) then |
179 |
|
|
180 |
close(iunit) |
close(iunit) |
181 |
open (iunit,form='unformatted',access='direct', |
open (iunit,file=sicedata,form='unformatted',access='direct', |
182 |
. recl=im2*jm2*nPgx*nPgy*4) |
. recl=im2*jm2*nPgx*nPgy*4) |
183 |
nrec = 1 |
nrec = 1 |
184 |
call bcheader (iunit, ndmax, nrec, |
call bcheader (iunit, ndmax, nrec, |
374 |
|
|
375 |
character*8 cname |
character*8 cname |
376 |
character*80 cdscrip |
character*80 cdscrip |
377 |
|
character*40 sstdata |
378 |
_RL fac1, fac2, lat0, lon0, timebc1, timebc2, timemod, undef |
_RL fac1, fac2, lat0, lon0, timebc1, timebc2, timemod, undef |
379 |
logical first, found, error |
logical first, found, error |
380 |
integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc,npxbc,npybc |
integer i,j,n,nn,iyear,iyearbc,nd,ndby3,imbc,jmbc,npxbc,npybc |
388 |
|
|
389 |
data first /.true./ |
data first /.true./ |
390 |
data error /.false./ |
data error /.false./ |
391 |
|
sstdata='sst19232.weekly.clim' |
392 |
|
|
393 |
c save header info |
c save header info |
394 |
save imbc,jmbc,npxbc,npybc,lat0,lon0,ndatebc,undef,nymdbc,nhmsbc |
save imbc,jmbc,npxbc,npybc,lat0,lon0,ndatebc,undef,nymdbc,nhmsbc |
425 |
if( iyear.ne.iyearbc ) then |
if( iyear.ne.iyearbc ) then |
426 |
|
|
427 |
close(iunit) |
close(iunit) |
428 |
open (iunit,form='unformatted',access='direct', |
open (iunit,file=sstdata,form='unformatted',access='direct', |
429 |
. recl=im2*jm2*nPgx*nPgy*4) |
. recl=im2*jm2*nPgx*nPgy*4) |
430 |
nrec = 1 |
nrec = 1 |
431 |
call bcheader (iunit, ndmax, nrec, |
call bcheader (iunit, ndmax, nrec, |