New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2190 – NEMO

Changeset 2190


Ignore:
Timestamp:
2010-10-08T17:56:19+02:00 (14 years ago)
Author:
acc
Message:

NEMO branch DEV_r1924_nocs_latphys. Code changes to fix problems identified by NVTK. NVTK testing successful. Namelists updated and default setting for ln_traldf_grif now set to .false. (i.e. new code is inactive unless explicitly activated)

Location:
branches/DEV_r1924_nocs_latphys
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • branches/DEV_r1924_nocs_latphys/CONFIG/GYRE/EXP00/namelist

    r1741 r2190  
    401401   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    402402   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     403   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    403404                           !  Coefficient 
    404405   rn_aht_0         =  1000.    !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r1741 r2190  
    401401   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    402402   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     403   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    403404                           !  Coefficient 
    404405   rn_aht_0         =  1000.    !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/CONFIG/ORCA2_LIM/EXP00/namelist

    r1759 r2190  
    428428   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    429429   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     430   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    430431                           !  Coefficient 
    431432   rn_aht_0         =  2000.    !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist

    r1759 r2190  
    428428   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    429429   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     430   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    430431                           !  Coefficient 
    431432   rn_aht_0         =  2000.    !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r1798 r2190  
    9191   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    9292   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     93   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    9394!                               !  Coefficient 
    9495   aht0        =  2000.         !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/CONFIG/POMME/EXP00/namelist

    r1875 r2190  
    425425   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    426426   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
     427   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp") 
    427428                           !  Coefficient 
    428429   rn_aht_0         =   300.    !     horizontal eddy diffusivity for tracers [m2/s] 

  • branches/DEV_r1924_nocs_latphys/NEMO/OPA_SRC/LDF/ldfslp.F90

    r2163 r2190  
    558558        ! Depth of the w-point below ML base 
    559559        ! --------------------------------------- 
    560         DO jj = 1, jpjm1 
    561            DO ji = 1, fs_jpim1 
     560        DO jj = 1, jpj 
     561           DO ji = 1, jpi 
    562562              jk = nmln(ji,jj) 
    563563              zr_ml_basew(ji,jj)=1.0/fsdepw(ji,jj,jk+1) 
    … …  
    617617              DO jj = 1, jpjm1 
    618618                 DO ji = 1, fs_jpim1 
    619                     jk = nmln(ji,jj) 
    620                     zsx_ml_base(ji,jj,ip,kp)=zsx(ji,jj,jk+1-kp,ip,kp) 
     619                    jk = nmln(ji+ip,jj) 
     620                    zsx_ml_base(ji+ip,jj,1-ip,kp)=zsx(ji+ip,jj,jk+1-kp,1-ip,kp) 
    621621                 END DO 
    622622              END DO 
    … …  
    634634 
    635635                    DO ji = 1, fs_jpim1   ! vector opt. 
    636                        zfact = 1 - 1/(1 + (jk+1-kp)/nmln(ji,jj)) 
    637                        zsx(ji,jj,jk,ip,kp) = zfact*zsx(ji,jj,jk,ip,kp) + & 
    638                             & (1.0_wp-zfact)*(fsdepw(ji,jj,jk+1-kp)*zr_ml_basew(ji,jj))*zsx_ml_base(ji,jj,ip,kp)  
     636                       ! k index of uppermost point(s) of triad is jk+kp-1 
     637                       ! must be .ge. nmln(ji,jj) for zfact=1. 
     638                       !                   otherwise  zfact=0. 
     639                       zfact = 1 - 1/(1 + (jk+kp-1)/nmln(ji+ip,jj)) 
     640                       zsx(ji+ip,jj,jk,1-ip,kp) = zfact*zsx(ji+ip,jj,jk,1-ip,kp) + & 
     641                            & (1.0_wp-zfact)*(fsdepw(ji+ip,jj,jk+kp)*zr_ml_basew(ji+ip,jj))*zsx_ml_base(ji+ip,jj,1-ip,kp)  
    639642                    END DO 
    640643 
    … …  
    674677                       wslp2(ji+ip,jj,jk+kp)=wslp2(ji+ip,jj,jk+kp)+ & 
    675678                            & ((zvolf*ze3wr)/(e1t(ji+ip,jj)*e2t(ji+ip,jj)))*(zslope)**2 
    676                        psix_eiv(ji+ip,jj,jk+kp) = psix_eiv(ji+ip,jj,jk+kp) + 0.25_wp*zslope 
     679                       psix_eiv(ji,jj,jk+kp) = psix_eiv(ji,jj,jk+kp) + 0.25_wp*zslope 
    677680 
    678681                    END DO 
    … …  
    720723              DO jj = 1, jpjm1 
    721724                 DO ji = 1, fs_jpim1 
    722                     jk = nmln(ji,jj) 
    723                     zsy_ml_base(ji,jj,jp,kp)=zsy(ji,jj,jk+1-kp,jp,kp) 
     725                    jk = nmln(ji,jj+jp) 
     726                    zsy_ml_base(ji,jj+jp,1-jp,kp)=zsy(ji,jj+jp,jk+1-kp,1-jp,kp) 
    724727                 END DO 
    725728              END DO 
    … …  
    737740 
    738741                    DO ji = 1, fs_jpim1   ! vector opt. 
    739                        zfact = 1 - 1/(1 + (jk+1-kp)/nmln(ji,jj)) 
    740                        zsy(ji,jj,jk,jp,kp) = zfact*zsy(ji,jj,jk,jp,kp) +& 
    741                             & (1.0_wp-zfact)*(fsdepw(ji,jj,jk+1-kp)*zr_ml_basew(ji,jj))*zsy_ml_base(ji,jj,jp,kp)  
     742                       ! k index of uppermost point(s) of triad is jk+kp-1 
     743                       ! must be .ge. nmln(ji,jj) for zfact=1. 
     744                       !                   otherwise  zfact=0. 
     745                       zfact = 1 - 1/(1 + (jk+kp-1)/nmln(ji,jj+jp)) 
     746                       zsy(ji,jj+jp,jk,1-jp,kp) = zfact*zsy(ji,jj+jp,jk,1-jp,kp) + & 
     747                            & (1.0_wp-zfact)*(fsdepw(ji,jj+jp,jk+kp)*zr_ml_basew(ji,jj+jp))*zsy_ml_base(ji,jj+jp,1-jp,kp)  
    742748                    END DO 
    743749 
    … …  
    777783                       wslp2(ji,jj+jp,jk+kp)=wslp2(ji,jj+jp,jk+kp)+ & 
    778784                            & ((zvolf*ze3wr)/(e1t(ji,jj+jp)*e2t(ji,jj+jp)))*(zslope)**2 
    779                        psiy_eiv(ji,jj+jp,jk+kp) = psiy_eiv(ji,jj+jp,jk+kp) + 0.25_wp*zslope 
     785                       psiy_eiv(ji,jj,jk+kp) = psiy_eiv(ji,jj,jk+kp) + 0.25_wp*zslope 
    780786 
    781787                    END DO 

  • branches/DEV_r1924_nocs_latphys/NEMO/OPA_SRC/LDF/ldftra_oce.F90

    r2177 r2190  
    2020   LOGICAL , PUBLIC ::   ln_traldf_hor   = .FALSE.   !: horizontal (geopotential) direction 
    2121   LOGICAL , PUBLIC ::   ln_traldf_iso   = .TRUE.    !: iso-neutral direction 
    22    LOGICAL , PUBLIC ::   ln_traldf_grif  = .TRUE.    !: griffies skew flux 
     22   LOGICAL , PUBLIC ::   ln_traldf_grif  = .FALSE.   !: griffies skew flux 
    2323   REAL(wp), PUBLIC ::   rn_aht_0        = 2000._wp  !: lateral eddy diffusivity (m2/s) 
    2424   REAL(wp), PUBLIC ::   rn_ahtb_0       =    0._wp  !: lateral background eddy diffusivity (m2/s) 
    … …  
    2727 
    2828   REAL(wp), PUBLIC ::   aht0, ahtb0, aeiv0         !!: OLD namelist names 
     29   REAL(wp), PUBLIC ::   slpmax                     !: slope limit  
    2930 
    3031#if defined key_traldf_c3d 
    … …  
    4445   !!---------------------------------------------------------------------- 
    4546   LOGICAL, PUBLIC, PARAMETER               ::   lk_traldf_eiv   = .TRUE.   !: eddy induced velocity flag 
    46    REAL(wp), PUBLIC                         ::   slpmax                     !: slope limit  
    4747       
    4848# if defined key_traldf_c3d 

  • branches/DEV_r1924_nocs_latphys/NEMO/OPA_SRC/TRA/traldf.F90

    r2177 r2190  
    153153         WRITE(numout,*) '      horizontal (geopotential)   ln_traldf_hor   = ', ln_traldf_hor 
    154154         WRITE(numout,*) '      iso-neutral                 ln_traldf_iso   = ', ln_traldf_iso 
    155          WRITE(numout,*) '          iso-neutral (Griffies)      ln_traldf_grif  = ', ln_traldf_grif 
     155         WRITE(numout,*) '      iso-neutral (Griffies)      ln_traldf_grif  = ', ln_traldf_grif 
    156156      ENDIF 
    157157 

  • branches/DEV_r1924_nocs_latphys/NEMO/OPA_SRC/TRA/traldf_iso_grif.F90

    r2163 r2190  
    178178      END DO 
    179179      ! 
    180       psix_eiv(:,:,:) = psix_eiv(:,:,:)*fsaeiu(:,:,:) 
    181       psiy_eiv(:,:,:) = psiy_eiv(:,:,:)*fsaeiv(:,:,:) 
     180      DO jk=1,jpk 
     181      DO jj=2,jpjm1 
     182      DO ji=fs_2,fs_jpim1 
     183       psix_eiv(ji,jj,jk) = psix_eiv(ji,jj,jk)*fsaeiu(ji,jj,jk) 
     184       psiy_eiv(ji,jj,jk) = psiy_eiv(ji,jj,jk)*fsaeiv(ji,jj,jk) 
     185      END DO 
     186      END DO 
     187      END DO 
    182188 
    183189   END SUBROUTINE tra_ldf_iso_grif 
Note: See TracChangeset for help on using the changeset viewer.