igsm/src/obssolar.F


#include "ctrparam.h"

!       ==========================================================
!
!       OBSSOLAR.F:  This file contains subroitine obssolar.f
!                which reads observed solar constant
!
!       ----------------------------------------------------------
! 
!       Revision History:
!       
!       When    Who             What
!       ----    ----------      ------- 
!       070302  
!       
!
!       ==========================================================

      subroutine obssolar(S0const,YEAR1)
      parameter (nyd=2000,nyd1=nyd+1)
C
C
C     For simulations with observed solar constant using equivalent
C     S0conts  from file S0C_data
C     DATA are annual means, but they are
C       used as a data for the middle of the coresponding year
C
C
      dimension  solar(nyd),solar1(nyd1),yearS0(nyd)
      common /solardata/S0C_data
      character * 120 S0C_data
      logical first
      data first /.true./
      if(first)then
         print *,'  Solar constant from observation '
         open (unit=561,file=S0C_data,
     *     status='OLD',form='formatted')
          do 50 i=1,nyd
c        read (561,501,end=500),ighgyr(i),tmp,xco2(i)
            read (561,*,end=500),yearS0(i),solar(i)
c           print *,yearS0(i),solar(i)
   50    continue
  500    continue
         close (561)
         ny=i-1
         print *,'Last year for solardata is ',yearS0(ny)
         ny1=ny+1
         ny2=ny+2
  501    format(2f13.2)
         first=.false.
      endif
      year=year1
c        print *,' year=',year
      if(year.lt.yearS0(1))year=yearS0(1)
      if(year.gt.yearS0(ny))year=yearS0(ny)
      i=year-yearS0(1)+1
      dl=year-yearS0(i)
c     print *,' i=',i,' dl=',dl
c     print *,yearS0(i),yearS0(i+1)
c     print *,solar(i),solar(i+1)
      S0const=solar(i)*(1.-dl)+solar(i+1)*dl
c          print *,' year=',year,' S0x=',S0const
      return
      end
1
2	#include "ctrparam.h"
3
4	! ==========================================================
5	!
6	! OBSSOLAR.F: This file contains subroitine obssolar.f
7	! which reads observed solar constant
8	!
9	! ----------------------------------------------------------
10	!
11	! Revision History:
12	!
13	! When Who What
14	! ---- ---------- -------
15	! 070302
16	!
17	!
18	! ==========================================================
19
20	subroutine obssolar(S0const,YEAR1)
21	parameter (nyd=2000,nyd1=nyd+1)
22	C
23	C
24	C For simulations with observed solar constant using equivalent
25	C S0conts from file S0C_data
26	C DATA are annual means, but they are
27	C used as a data for the middle of the coresponding year
28	C
29	C
30	dimension solar(nyd),solar1(nyd1),yearS0(nyd)
31	common /solardata/S0C_data
32	character * 120 S0C_data
33	logical first
34	data first /.true./
35	if(first)then
36	print *,' Solar constant from observation '
37	open (unit=561,file=S0C_data,
38	* status='OLD',form='formatted')
39	do 50 i=1,nyd
40	c read (561,501,end=500),ighgyr(i),tmp,xco2(i)
41	read (561,*,end=500),yearS0(i),solar(i)
42	c print *,yearS0(i),solar(i)
43	50 continue
44	500 continue
45	close (561)
46	ny=i-1
47	print *,'Last year for solardata is ',yearS0(ny)
48	ny1=ny+1
49	ny2=ny+2
50	501 format(2f13.2)
51	first=.false.
52	endif
53	year=year1
54	c print *,' year=',year
55	if(year.lt.yearS0(1))year=yearS0(1)
56	if(year.gt.yearS0(ny))year=yearS0(ny)
57	i=year-yearS0(1)+1
58	dl=year-yearS0(i)
59	c print *,' i=',i,' dl=',dl
60	c print *,yearS0(i),yearS0(i+1)
61	c print *,solar(i),solar(i+1)
62	S0const=solar(i)(1.-dl)+solar(i+1)dl
63	c print *,' year=',year,' S0x=',S0const
64	return
65	end