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 14224 – NEMO

Changeset 14224


Ignore:
Timestamp:
2020-12-19T15:49:56+01:00 (3 years ago)
Author:
smasson
Message:

trunk: minor bugfix in CANAL and with ln_dyn_trd

Location:
NEMO/trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/cfgs/SHARED/field_def_nemo-oce.xml

    r14146 r14224  
    183183    <field id="heatc"        long_name="Heat content vertically integrated"                 standard_name="integral_of_sea_water_potential_temperature_wrt_depth_expressed_as_heat_content"   unit="J/m2"       /> 
    184184    <field id="saltc"        long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3*kg/m2" /> 
     185    <field id="saltc2"       long_name="Salt content vertically integrated"                                                                                                                   unit="1e-3*kg/m2" > saltc * saltc </field > 
    185186 
    186187    <!-- EOS --> 

  • NEMO/trunk/src/OCE/DYN/dynatf.F90

    r14143 r14224  
    121121      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    122122      REAL(wp) ::   zue3a, zue3n, zue3b, zcoef    ! local scalars 
    123       REAL(wp) ::   zve3a, zve3n, zve3b, z1_2dt   !   -      - 
     123      REAL(wp) ::   zve3a, zve3n, zve3b           !   -      - 
    124124      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zue, zve, zwfld 
    125125      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zutau, zvtau 
    … …  
    320320      ENDIF 
    321321      IF( l_trddyn ) THEN                ! 3D output: asselin filter trends on momentum 
    322          zua(:,:,:) = ( puu(:,:,:,Kmm) - zua(:,:,:) ) * z1_2dt 
    323          zva(:,:,:) = ( pvv(:,:,:,Kmm) - zva(:,:,:) ) * z1_2dt 
     322         zua(:,:,:) = ( puu(:,:,:,Kmm) - zua(:,:,:) ) * r1_Dt 
     323         zva(:,:,:) = ( pvv(:,:,:,Kmm) - zva(:,:,:) ) * r1_Dt 
    324324         CALL trd_dyn( zua, zva, jpdyn_atf, kt, Kmm ) 
    325325      ENDIF 

  • NEMO/trunk/src/OCE/DYN/dynatf_qco.F90

    r14143 r14224  
    9898      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    9999      REAL(wp) ::   zue3a, zue3n, zue3b, zcoef    ! local scalars 
    100       REAL(wp) ::   zve3a, zve3n, zve3b, z1_2dt   !   -      - 
     100      REAL(wp) ::   zve3a, zve3n, zve3b           !   -      - 
    101101      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zue, zve 
    102102      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zua, zva 
    … …  
    234234      ENDIF 
    235235      IF( l_trddyn ) THEN                ! 3D output: asselin filter trends on momentum 
    236          zua(:,:,:) = ( puu(:,:,:,Kmm) - zua(:,:,:) ) * z1_2dt 
    237          zva(:,:,:) = ( pvv(:,:,:,Kmm) - zva(:,:,:) ) * z1_2dt 
     236         zua(:,:,:) = ( puu(:,:,:,Kmm) - zua(:,:,:) ) * r1_Dt 
     237         zva(:,:,:) = ( pvv(:,:,:,Kmm) - zva(:,:,:) ) * r1_Dt 
    238238         CALL trd_dyn( zua, zva, jpdyn_atf, kt, Kmm ) 
    239239      ENDIF 

  • NEMO/trunk/tests/CANAL/EXPREF/file_def_nemo-oce.xml

    r13472 r14224  
    1616     <field field_ref="ssh"  /> 
    1717     <field field_ref="socegrad"  /> 
    18      <field field_ref="eken_int"  /> 
     18     <field field_ref="ke_int"  /> 
    1919     <field field_ref="relvor"  /> 
    2020     <field field_ref="potvor"  /> 

  • NEMO/trunk/tests/CANAL/EXPREF/namelist_cfg

    r14053 r14224  
    231231&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection) 
    232232!----------------------------------------------------------------------- 
    233    ln_dynvor_ene = .false.  !  energy conserving scheme 
    234    ln_dynvor_ens = .false.  !  enstrophy conserving scheme 
    235    ln_dynvor_mix = .false.  !  mixed scheme 
    236    ln_dynvor_een = .false.  !  energy & enstrophy scheme 
    237 r_enT = .false.  !  energy conserving scheme (T-point) 
    238    ln_dynvor_eeT = .true.   !  energy conserving scheme (een using e3t) 
    239       nn_e3f_typ = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     233   ln_dynvor_ene = .false. !  energy    conserving scheme 
     234   ln_dynvor_ens = .false. !  enstrophy conserving scheme 
     235   ln_dynvor_mix = .false. !  mixed scheme 
     236   ln_dynvor_enT = .false. !  energy conserving scheme (T-point) 
     237   ln_dynvor_eeT =  .true. !  energy conserving scheme (een using e3t) 
     238   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     239   ! 
     240   nn_e3f_typ = 0          !  type of e3f (EEN, ENE, ENS, MIX only)  =0  e3f = mi(mj(e3t))/4 
     241   !                       !                                         =1  e3f = mi(mj(e3t))/mi(mj( tmask)) 
    240242/ 
    241243!----------------------------------------------------------------------- 

  • NEMO/trunk/tests/CANAL/MY_SRC/usrdef_istate.F90

    r14133 r14224  
    7676      CASE(-1)    ! stratif at rest 
    7777 
    78          ! sea level: 
    79          pssh(:,:) = 0. 
    8078         ! temperature: 
    8179         pts(:,:,1,jp_tem) = 25. !!30._wp 
    … …  
    146144         ! salinity:   
    147145         DO jk=1, jpkm1 
    148             pts(:,:,jk,jp_sal) = pssh(:,:) 
     146            pts(:,:,jk,jp_sal) = gphit(:,:) 
    149147         END DO 
    150148         ! velocities: 
Note: See TracChangeset for help on using the changeset viewer.