;	Annual Primary Production with model_1:
;		--- Pb_opt is a function of T
;	the result is converted into logarithm scale:
;	0<count<246  =  (log10(Prim_Prod)+1.3367)/0.0357
;	for the case  0.005<Prim_Prod<2.83291e+07
;
@model_noT_cld_inh_ful.pro
img_name='Primary Production Map'

win_xsize=2048
win_ysize=1024
win_xsize1=win_xsize-1
win_ysize1=win_ysize-1
;;window,0,color=256,xsize=win_xsize/4,ysize=win_ysize/4,title='Global.PP'
;;restore,'annual.rgb'
;;tvlct,r,g,b

;	Input the Monthly Average Primary Production
;
Chl_file=strarr(14)
Irr_file=strarr(14)
Julian_date=intarr(14)
day_of_month=intarr(14)
Temperary_data=fltarr(win_xsize,win_ysize)
;;Annual_PP=fltarr(win_xsize,win_ysize)

;Chl_file(0)='/user5/czcs/full_fil_12.chl'
Chl_file(0)='/user5/newchl/new3_mean12.flt'
Irr_file(0)='/user5/seawifs/ful_8312.irr'
day_of_month(0)=31
Julian_date(0)=349

;Chl_file(1)='/user5/czcs/full_fil_01.chl'
Chl_file(1)='/user5/newchl/new3_mean01.flt'
Irr_file(1)='/user5/seawifs/ful_8401.irr'
day_of_month(1)=31
Julian_date(1)=15

;Chl_file(2)='/user5/czcs/full_fil_02.chl'
Chl_file(2)='/user5/newchl/new3_mean02.flt'
Irr_file(2)='/user5/seawifs/ful_8402.irr'
day_of_month(2)=28
Julian_date(2)=Julian_date(1)+day_of_month(1)

;Chl_file(3)='/user5/czcs/full_fil_03.chl'
Chl_file(3)='/user5/newchl/new3_mean03.flt'
Irr_file(3)='/user5/seawifs/ful_8403.irr'
day_of_month(3)=31
Julian_date(3)=Julian_date(2)+day_of_month(2)

;Chl_file(4)='/user5/czcs/full_fil_04.chl'
Chl_file(4)='/user5/newchl/new3_mean04.flt'
Irr_file(4)='/user5/seawifs/ful_8404.irr'
day_of_month(4)=30
Julian_date(4)=Julian_date(3)+day_of_month(3)

;Chl_file(5)='/user5/czcs/full_fil_05.chl'
Chl_file(5)='/user5/newchl/new3_mean05.flt'
Irr_file(5)='/user5/seawifs/ful_8405.irr'
day_of_month(5)=31
Julian_date(5)=Julian_date(4)+day_of_month(4)

;Chl_file(6)='/user5/czcs/full_fil_06.chl'
Chl_file(6)='/user5/newchl/new3_mean06.flt'
Irr_file(6)='/user5/seawifs/ful_8406.irr'
day_of_month(6)=30
Julian_date(6)=Julian_date(5)+day_of_month(5)

;Chl_file(7)='/user5/czcs/full_fil_07.chl'
Chl_file(7)='/user5/newchl/new3_mean07.flt'
Irr_file(7)='/user5/seawifs/ful_8307.irr'
day_of_month(7)=31
Julian_date(7)=Julian_date(6)+day_of_month(6)

;Chl_file(8)='/user5/czcs/full_fil_08.chl'
Chl_file(8)='/user5/newchl/new3_mean08.flt'
Irr_file(8)='/user5/seawifs/ful_8308.irr'
day_of_month(8)=31
Julian_date(8)=Julian_date(7)+day_of_month(7)

;Chl_file(9)='/user5/czcs/full_fil_09.chl'
Chl_file(9)='/user5/newchl/new3_mean09.flt'
Irr_file(9)='/user5/seawifs/ful_8309.irr'
day_of_month(9)=30
Julian_date(9)=Julian_date(8)+day_of_month(8)

;Chl_file(10)='/user5/czcs/full_fil_10.chl'
Chl_file(10)='/user5/newchl/new3_mean10.flt'
Irr_file(10)='/user5/seawifs/ful_8310.irr'
day_of_month(10)=31
Julian_date(10)=Julian_date(9)+day_of_month(9)

;Chl_file(11)='/user5/czcs/full_fil_11.chl'
Chl_file(11)='/user5/newchl/new3_mean11.flt'
Irr_file(11)='/user5/seawifs/ful_8311.irr'
day_of_month(11)=30
Julian_date(11)=Julian_date(10)+day_of_month(10)

;Chl_file(12)='/user5/czcs/full_fil_12.chl'
Chl_file(12)='/user5/newchl/new3_mean12.flt'
Irr_file(12)='/user5/seawifs/ful_8312.irr'
day_of_month(12)=31
Julian_date(12)=Julian_date(11)+day_of_month(11)

;Chl_file(13)='/user5/czcs/full_fil_01.chl'
Chl_file(13)='/user5/newchl/new3_mean01.flt'
Irr_file(13)='/user5/seawifs/ful_8401.irr'
day_of_month(13)=31
Julian_date(13)=15

;	Calculate the Annual Primary Productivity
;		from the average monthly

cell_area=19.*19.
pi=3.141592654
r_constant=pi/180.
;cell_area=cell_area*(2048./360.)*(2048./360.)

;	Annual_PP
Annual_PP=0.
for k=1,12 do begin
    model,Chl_file(k),Irr_file(k),Julian_date(k),Temperary_data
    Annual_PP=TEMPORARY(Annual_PP)+(Temperary_data > 0.)* $
	float(day_of_month(k))+(Temperary_data < 0.)
;;    tvscl,rebin(Temperary_data,win_xsize/4,win_ysize/4,/sample)
endfor
Annual_PP=(Annual_PP > 0.)-float(Annual_PP lt 0.)

;	Display the Annual_PP
window,0,color=256,xsize=win_xsize/4,ysize=win_ysize/4,title='Global.PP'
restore,'annual.rgb'
tvlct,r,g,b
tvscl,rebin(Annual_PP,win_xsize/4,win_ysize/4,/sample)

;	Save the Annual Primary Production
;
openw,2,'/user3/PP/Annu_Glb_noT_Cld_Inh_ful.dat'
writeu,2,Annual_PP
close,2

;	Calculate the Integrated Global Productivity
PP_global=0.
PP_lat=0
PP_lat=TOTAL((Annual_PP > 0.),1)
for j=0,win_ysize1 do begin
    r_lat=(float(j)/float(win_ysize1)-0.5)*pi
    PP_lat(j)=PP_lat(j)*COS(r_lat)
endfor
PP_global=TOTAL(PP_lat)*cell_area
PP_lat=0
Annual_PP=0
Temperary_data=0

print, 'Global_Annual_PP=',PP_global

;	Save the data to file 'annuPP_noT_cld_inh_ful.prn'
openw, 2, '/user3/PP/annuPP_noT_cld_inh_ful.prn'
printf,2, 'Annual_PP with no T, with clouds, with inhibition:'
printf,2, 'Global_Annual_PP=',PP_global
close,2

end