Changeset 9968


Ignore:
Timestamp:
2018-07-19T09:39:02+02:00 (2 years ago)
Author:
timgraham
Message:

Correction to fix in cpl_oasis3.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_merge_2017_GC_couple_pkg/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90

    r9924 r9968  
    8080   INTEGER, PUBLIC, PARAMETER ::   nmaxcat=5    ! Maximum number of coupling fields 
    8181   INTEGER, PUBLIC, PARAMETER ::   nmaxcpl=5    ! Maximum number of coupling fields 
     82 
     83   INTEGER          ::   nldi_save, nlei_save 
     84   INTEGER          ::   nldj_save, nlej_save 
    8285    
    8386   TYPE, PUBLIC ::   FLD_CPL               !: Type for coupling field information 
     
    335338      INTEGER                   , INTENT(  out) ::   kinfo     ! OASIS3 info argument 
    336339      INTEGER                   , INTENT(in   ) ::   kstep     ! ocean time-step in seconds 
    337       REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   pdata 
     340      REAL(wp), DIMENSION(:,:,:), INTENT(in   ) ::   pdata 
    338341      !! 
    339342      INTEGER                                   ::   jc,jm     ! local loop index 
     
    344347      ! snd data to OASIS3 
    345348       
    346       CALL lbc_lnk(pdata(:,:,:), ssnd(kid)%clgrid, ssnd(kid)%nsgn) 
    347349 
    348350      nldi_save = nldi   ;   nlei_save = nlei 
     
    408410      nldi_save = nldi   ;   nlei_save = nlei 
    409411      nldj_save = nldj   ;   nlej_save = nlej 
    410       IF( nimpp           ==      1 ) nldi = 1 
    411       IF( nimpp + jpi - 1 == jpiglo ) nlei = jpi 
    412       IF( njmpp           ==      1 ) nldj = 1 
    413       IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 
     412       
    414413      DO jc = 1, srcv(kid)%nct 
     414 
     415         IF( nimpp           ==      1 ) nldi = 1 
     416         IF( nimpp + jpi - 1 == jpiglo ) nlei = jpi 
     417         IF( njmpp           ==      1 ) nldj = 1 
     418         IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 
    415419         llfisrt = .TRUE. 
    416420 
     
    446450                     WRITE(numout,*) '     -     Sum value is ', SUM(pdata(:,:,jc)) 
    447451                     WRITE(numout,*) '****************' 
     452                     CALL FLUSH(numout) 
    448453                  ENDIF 
    449454                   
     
    459464         !--- check periodicity conditions (all cases) 
    460465         IF( .not. llfisrt )   CALL lbc_lnk( pdata(:,:,jc), srcv(kid)%clgrid, srcv(kid)%nsgn )    
    461  
    462          IF( nimpp           ==      1 ) nldi = 1 
    463          IF( nimpp + jpi - 1 == jpiglo ) nlei = jpi 
    464          IF( njmpp           ==      1 ) nldj = 1 
    465          IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 
    466466  
    467467      ENDDO 
Note: See TracChangeset for help on using the changeset viewer.