@@ -114,9 +114,6 @@ if ( .not.use_daily ) then
114114 temp1 (:,:) = temp2(:,:) ! SJL BUG fix: June 24, 2011
115115 else
116116 call READ_RECORD_ (' sst' , Date1, Udate1, temp1)
117- if ( use_ncep_sst .and . (.not . no_anom_sst) ) then
118- temp1 = temp1 + SST_ANOM_
119- endif
120117 call horiz_interp ( Interp%Hintrp, temp1, Interp%DATA1_)
121118 call CLIP_DATA_ (' sst' , Interp%DATA1_)
122119 Interp%Date1 = Date1
@@ -127,9 +124,6 @@ if ( .not.use_daily ) then
127124
128125 if (Date2 /= Interp%Date2) then
129126 call READ_RECORD_ (' sst' , Date2, Udate2, temp2)
130- if ( use_ncep_sst .and . (.not . no_anom_sst) ) then
131- temp2 = temp2 + SST_ANOM_
132- endif
133127 call horiz_interp ( Interp%Hintrp, temp2, Interp%DATA2_)
134128 call CLIP_DATA_ (' sst' , Interp%DATA2_)
135129 Interp%Date2 = Date2
@@ -145,11 +139,6 @@ if ( .not.use_daily ) then
145139! Nudging runs: (Note: NCEP SST updated only every 6 -hr)
146140! Compute SST anomaly from global SST datasets for subsequent forecast runs
147141!-------------------------------------------------------------------------------
148- if ( use_ncep_sst .and . no_anom_sst ) then
149- sst_anom = SST_NCEP_ - (temp1 + fmonth*(temp2 - temp1))
150- call horiz_interp (Interp%Hintrp, SST_NCEP_, sst)
151- call CLIP_DATA_ (' sst' , sst)
152- endif
153142
154143!! DEBUG CODE
155144 if (DEBUG) then
@@ -380,15 +369,7 @@ else
380369!-- SJL -------------------------------------------------------------
381370! Can NOT use ncep_sst to determine sea_ice For seasonal forecast
382371! Use climo sea ice for seasonal runs
383- if ( use_ncep_sst .and . use_ncep_ice ) then
384- where ( SST_NCEP_ <= (real(TFREEZE, FMS_AMIP_INTERP_KIND_)+real(tice_crit, FMS_AMIP_INTERP_KIND_)) )
385- sice = 1._lkind
386- elsewhere
387- sice = 0._lkind
388- endwhere
389- else
390- call READ_RECORD_ (' ice' , Date1, Udate1, sice)
391- endif
372+ call READ_RECORD_ (' ice' , Date1, Udate1, sice)
392373!--------------------------------------------------------------------
393374 call horiz_interp ( Interp%Hintrp, sice, Interp%DATA1_)
394375 call CLIP_DATA_ (' ice' , Interp%DATA1_)
@@ -401,15 +382,7 @@ else
401382 if (Date2 /= Interp%Date2) then
402383
403384!-- SJL -------------------------------------------------------------
404- if ( use_ncep_sst .and . use_ncep_ice ) then
405- where ( SST_NCEP_ <= (real(TFREEZE, FMS_AMIP_INTERP_KIND_)+real(tice_crit, FMS_AMIP_INTERP_KIND_)) )
406- sice = 1._lkind
407- elsewhere
408- sice = 0._lkind
409- endwhere
410- else
411- call READ_RECORD_ (' ice' , Date2, Udate2, sice)
412- endif
385+ call READ_RECORD_ (' ice' , Date2, Udate2, sice)
413386!--------------------------------------------------------------------
414387 call horiz_interp ( Interp%Hintrp, sice, Interp%DATA2_)
415388 call CLIP_DATA_ (' ice' , Interp%DATA2_)
0 commit comments