/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/wavebands_init_fixed.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/wavebands_init_fixed.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.2 by jahn, Thu Aug 9 21:18:53 2012 UTC
# Line 63  c local indeces Line 63  c local indeces
63  C band widths used to convert OASIM data to irradiation per nm  C band widths used to convert OASIM data to irradiation per nm
64  C put boundaries half-way between central values  C put boundaries half-way between central values
65  C but first and last boundary are at first and last "central" value  C but first and last boundary are at first and last "central" value
66        wb_width(1) = .5*(pwaves(2)-pwaves(1))        if (darwin_wavebands(1).EQ.0) then
67        do i=2,tlam-1          wb_width(1) = .5*(pwaves(2)-pwaves(1))
68          wb_width(i) = .5*(pwaves(i+1)-pwaves(i-1))          do i=2,tlam-1
69        enddo            wb_width(i) = .5*(pwaves(i+1)-pwaves(i-1))
70        wb_width(tlam) = .5*(pwaves(tlam)-pwaves(tlam-1))          enddo
71            wb_width(tlam) = .5*(pwaves(tlam)-pwaves(tlam-1))
72          else
73            do i=1,tlam
74              wb_width(i) = darwin_wavebands(i+1) - darwin_wavebands(i)
75            enddo
76          endif
77        wb_totalWidth = 0.0        wb_totalWidth = 0.0
78        do i=1,tlam        do i=1,tlam
79          wb_totalWidth = wb_totalWidth + wb_width(i)          wb_totalWidth = wb_totalWidth + wb_width(i)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22