#include "CPP_OPTIONS.h" subroutine smooth_correl2Dw ( U fld_in,mask_in,xx_gen_file,mythid) IMPLICIT NONE #include "SIZE.h" #include "EEPARAMS.h" #include "GRID.h" #include "PARAMS.h" c#include "tamc.h" #include "smooth.h" #include "ctrl.h" #include "ecco_cost.h" _RL mask_in(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy) _RL fld_in(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RL tmpW character*(MAX_LEN_FNAM) xx_gen_file integer i,j,bi,bj integer itlo,ithi integer jtlo,jthi integer myThid jtlo = mybylo(mythid) jthi = mybyhi(mythid) itlo = mybxlo(mythid) ithi = mybxhi(mythid) DO bj = jtlo,jthi DO bi = itlo,ithi DO j = 1,sNy DO i = 1,sNx if ( xx_gen_file .EQ. xx_hflux_file ) then tmpW=whflux(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_sflux_file ) then tmpW=wsflux(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_tauu_file ) then tmpW=wtauu(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_tauv_file ) then tmpW=wtauv(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_atemp_file ) then tmpW=watemp(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_aqh_file ) then tmpW=waqh(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_precip_file ) then tmpW=wprecip(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_snowprecip_file ) then tmpW=wsnowprecip(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_swflux_file ) then tmpW=wswflux(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_swdown_file ) then tmpW=wswdown(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_lwflux_file ) then tmpW=wlwflux(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_lwdown_file ) then tmpW=wlwdown(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_evap_file ) then tmpW=wevap(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_apressure_file ) then tmpW=wapressure(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_uwind_file ) then tmpW=wuwind(i,j,bi,bj) elseif ( xx_gen_file .EQ. xx_vwind_file ) then tmpW=wvwind(i,j,bi,bj) else tmpW=0. WRITE(errorMessageUnit,'(2A)' ) & 'no weights implemented here for ',xx_gen_file STOP 'ABNORMAL END: S/R smooth_correl2Dw' endif if ((mask_in(i,j,1,bi,bj).NE.0.).AND.(tmpW.NE.0.)) then fld_in(i,j,bi,bj)=fld_in(i,j,bi,bj)/sqrt(tmpW) else fld_in(i,j,bi,bj)=fld_in(i,j,bi,bj)*0. endif ENDDO ENDDO ENDDO ENDDO _EXCH_XY_RL ( fld_in , myThid ) end