Ignore:
Timestamp:
2013-11-20T09:48:15+01:00 (8 years ago)
Author:
clem
Message:

update lim3 for bdy purpose

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_LOCEAN_CMCC_INGV_MERC_UKMO_2013/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r4205 r4269  
    426426               !-Calculate stress tensor components zs1 and zs2  
    427427               !-at centre of grid cells (see section 3.5 of CICE user's guide). 
    428                zs1(ji,jj) = ( zs1(ji,jj) - dtotel*( ( 1._wp - alphaevp) * zs1(ji,jj) +   & 
    429                   &          ( delta / deltat(ji,jj) - zdd(ji,jj) / deltat(ji,jj) ) * zpresh(ji,jj) ) )  &        
    430                   &          / ( 1._wp + alphaevp * dtotel ) 
    431  
    432                zs2(ji,jj) = ( zs2(ji,jj) - dtotel * ( ( 1._wp - alphaevp ) * ecc2 * zs2(ji,jj) -   & 
    433                              zdt(ji,jj) / deltat(ji,jj) * zpresh(ji,jj) ) )   & 
    434                   &          / ( 1._wp + alphaevp * ecc2 * dtotel ) 
     428               !zs1(ji,jj) = ( zs1(ji,jj) - dtotel*( ( 1._wp - alphaevp) * zs1(ji,jj) +   & 
     429               !   &          ( delta / deltat(ji,jj) - zdd(ji,jj) / deltat(ji,jj) ) * zpresh(ji,jj) ) )  &        
     430               !   &          / ( 1._wp + alphaevp * dtotel ) 
     431 
     432               !zs2(ji,jj) = ( zs2(ji,jj) - dtotel * ( ( 1._wp - alphaevp ) * ecc2 * zs2(ji,jj) -   & 
     433               !              zdt(ji,jj) / deltat(ji,jj) * zpresh(ji,jj) ) )   & 
     434               !   &          / ( 1._wp + alphaevp * ecc2 * dtotel ) 
    435435 
    436436               ! new formulation from S. Bouillon to help stabilizing the code (no need of alphaevp) 
    437                !zs1(ji,jj) = ( zs1(ji,jj) + dtotel * ( ( zdd(ji,jj) / deltat(ji,jj) - delta / deltat(ji,jj) )  & 
    438                !   &         * zpresh(ji,jj) ) ) / ( 1._wp + dtotel ) 
    439                !zs2(ji,jj) = ( zs2(ji,jj) + dtotel * ( ecci * zdt(ji,jj) / deltat(ji,jj) * zpresh(ji,jj) ) )  & 
    440                !   &         / ( 1._wp + dtotel ) 
     437               zs1(ji,jj) = ( zs1(ji,jj) + dtotel * ( ( zdd(ji,jj) / deltat(ji,jj) - delta / deltat(ji,jj) )  & 
     438                  &         * zpresh(ji,jj) ) ) / ( 1._wp + dtotel ) 
     439               zs2(ji,jj) = ( zs2(ji,jj) + dtotel * ( ecci * zdt(ji,jj) / deltat(ji,jj) * zpresh(ji,jj) ) )  & 
     440                  &         / ( 1._wp + dtotel ) 
    441441 
    442442            END DO 
     
    472472 
    473473               !-Calculate stress tensor component zs12 at corners (see section 3.5 of CICE user's guide). 
    474                zs12(ji,jj) = ( zs12(ji,jj) - dtotel * ( (1.0-alphaevp) * ecc2 * zs12(ji,jj) - zds(ji,jj) /  & 
    475                   &          ( 2._wp * deltac(ji,jj) ) * zpreshc(ji,jj) ) )  & 
    476                   &          / ( 1._wp + alphaevp * ecc2 * dtotel )  
     474               !zs12(ji,jj) = ( zs12(ji,jj) - dtotel * ( (1.0-alphaevp) * ecc2 * zs12(ji,jj) - zds(ji,jj) /  & 
     475               !   &          ( 2._wp * deltac(ji,jj) ) * zpreshc(ji,jj) ) )  & 
     476               !   &          / ( 1._wp + alphaevp * ecc2 * dtotel )  
    477477 
    478478               ! new formulation from S. Bouillon to help stabilizing the code (no need of alphaevp) 
    479                !zs12(ji,jj) = ( zs12(ji,jj) + dtotel *  & 
    480                !   &          ( ecci * zds(ji,jj) / ( 2._wp * deltac(ji,jj) ) * zpreshc(ji,jj) ) )  & 
    481                !   &          / ( 1.0 + dtotel )  
     479               zs12(ji,jj) = ( zs12(ji,jj) + dtotel *  & 
     480                  &          ( ecci * zds(ji,jj) / ( 2._wp * deltac(ji,jj) ) * zpreshc(ji,jj) ) )  & 
     481                  &          / ( 1.0 + dtotel )  
    482482 
    483483            END DO ! ji 
Note: See TracChangeset for help on using the changeset viewer.