NASA Ocean Color

ocssw V2022

         subroutine contst                                               
 c
 c***********************************************************************
 c       subroutine contstx checks for the convergance of upwelling and
 c       downwelling radiance in the all direction
 c
 c***********************************************************************
 c.......include common blocks
         implicit real*8 (a-h,o-z)
         include 'common_all.cmn'
 c***********************************************************************
         if(jpass.le.3)return
         do i=1,jjjj
              do j=1,jpart
                 d1px(1,i,j)=d1px(2,i,j)
                 d2px(1,i,j)=d2px(2,i,j)
              enddo
         enddo
 c
         do i=1,jjjj
            do j=1,jpart
               r1=geoiup(jpass-3,i,j)
               r2=geoiup(jpass-2,i,j)
               r3=geoiup(jpass-1,i,j)
               rat=(r3-r2)/(r2-r1)
               r4=r1+(r2-r1)/(1.0-rat)
               d1px(2,i,j)=r4
 c             if(jpass.eq.5)then
 c               write(6,100)i,j,r1,r2,r3,r4,d1px(2,i,j)
 100             format('i,j,r1,r2,r3,r4,d1px',2i4,1p5e11.3)
 c             endif
            enddo
         enddo
         do i=1,jjjj
            do j=1,jpart
               u1=geoidn(jpass-3,i,j)
               u2=geoidn(jpass-2,i,j)
               u3=geoidn(jpass-1,i,j)
               uat=(u3-u2)/(u2-u1)
               u4=u1+(u2-u1)/(1.0-uat)
               d2px(2,i,j)=u4
 c             if(jpass.eq.5)then
 c               write(6,200)i,j,u1,u2,u3,u4,d2px(2,i,j)
 200             format('i,j,u1,u2,u3,u4,d2px',2i4,1p5e11.3)
 c             endif
            enddo
         enddo
 c       compute the difference
         do i=1,jjjj
              do j=1,jpart
                 df1(i,j)=100.0d0*(d1px(2,i,j)-d1px(1,i,j))/d1px(2,i,j)
                 df2(i,j)=100.0d0*(d2px(2,i,j)-d2px(1,i,j))/d2px(2,i,j)
              enddo
         enddo
 c       compute the max value
         xmax1=0.0d0
         xmax2=0.0d0
         do i=1,jjjj
              do j=1,jpart
                 if(df1(i,j).gt.xmax1)xmax1=df1(i,j)
                 if(df2(i,j).gt.xmax2)xmax2=df2(i,j)
              enddo
         enddo
         d3=max(xmax1,xmax2)
 c
 c       write(6,2233)jpass,xmax1,xmax2
 c       write(6,2235)d3
 c
 c*******format statements***********************************************
 2233    format('jpass, upward & downward dd)',i3,1p2e12.3)
 2235    format('max. diff',1pe12.3)
 2237    format('factr,flxgsup,flxgsdn,flxcongs'/1p3e13.4,0pf8.5)
 c***********************************************************************
         return                                                          
         end                                                             
 c***********************************************************************

contst

subroutine contst

Definition: contst.f:2

max

#define max(A, B)

Definition: main_biosmap.c:61

and

===========================================================================V5.0.48(Terra) 03/20/2015 Changes shown below are differences from MOD_PR02 V5.0.46(Terra)============================================================================Changes noted for V6.1.20(Terra) below were also instituted for this version.============================================================================V6.1.20(Terra) 03/12/2015 Changes shown below are differences from MOD_PR02 V6.1.18(Terra)============================================================================Changes from v6.1.18 which may affect scientific output:A situation can occur in which a scan which contains sector rotated data has a telemetry value indicating the completeness of the sector rotation. This issue is caused by the timing of the instrument command to perform the sector rotation and the recording of the telemetry point that reports the status of sector rotation. In this case a scan is considered valid by L1B and pass through the calibration - reporting extremely high radiances. Operationally the TEB calibration uses a 40 scan average coefficient, so the 20 scans(one mirror side) after the sector rotation are contaminated with anomalously high radiance values. A similar timing issue appeared before the sector rotation was fixed in V6.1.2. Our analysis indicates the ‘SET_FR_ENC_DELTA’ telemetry correlates well with the sector rotation encoder position. The use of this telemetry point to determine scans that are sector rotated should fix the anomaly occured before and after the sector rotation(usually due to the lunar roll maneuver). The fix related to the sector rotation in V6.1.2 is removed in this version.============================================================================V6.1.18(Terra) 10/01/2014 Changes shown below are differences from MOD_PR02 V6.1.16(Terra)============================================================================Added doi attributes to NRT(Near-Real-Time) product.============================================================================V6.1.16(Terra) 01/27/2014 Changes shown below are differences from MOD_PR02 V6.1.14(Terra)============================================================================Migrate to SDP Toolkit 5.2.17============================================================================V6.1.14(Terra) 06/26/2012 Changes shown below are differences from MOD_PR02 V6.1.12(Terra)============================================================================Added the doi metadata to L1B product============================================================================V6.1.12(Terra) 04/25/2011 Changes shown below are differences from MOD_PR02 V6.1.8(Terra)============================================================================1. The algorithm to calculate uncertainties for reflective solar bands(RSB) is updated. The current uncertainty in L1B code includes 9 terms from prelaunch analysis. The new algorithm regroups them with the new added contributions into 5 terms:u1:the common term(AOI and time independent) and

Definition: HISTORY.txt:126

for

for(i=0;i< NROOTS;i++) s[i]

Definition: decode_rs.h:85