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 109 for trunk/NEMO/OPA_SRC/DYN – NEMO

Changeset 109 for trunk/NEMO/OPA_SRC/DYN


Ignore:
Timestamp:
2004-06-28T12:21:11+02:00 (20 years ago)
Author:
opalod
Message:

CT : UPDATE069 : Vorticity diagnostics has been added

Location:
trunk/NEMO/OPA_SRC/DYN
Files:
15 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/DYN/dynhpg.F90

    r106 r109  
    125125            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) + zuap 
    126126            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) + zvap 
    127 # if defined key_trddyn 
     127# if defined key_trddyn || defined key_trd_vor 
    128128            ! save the trend for diagnostics 
    129129            utrd(ji,jj,1,1) = zhpi(ji,jj,1) + zuap 
    … …  
    153153               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) + zuap 
    154154               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) + zvap 
    155 # if defined key_trddyn 
     155# if defined key_trddyn || defined key_trd_vor 
    156156               ! save the momentum trends for diagnostics 
    157157               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) + zuap 
    … …  
    238238            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) 
    239239            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) 
    240 # if defined key_trddyn 
     240# if defined key_trddyn || defined key_trd_vor 
    241241            ! save the momentum trends for diagnostics 
    242242            utrd(ji,jj,1,1) = zhpi(ji,jj,1) 
    … …  
    263263               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) 
    264264               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) 
    265 # if defined key_trddyn 
     265# if defined key_trddyn || defined key_trd_vor 
    266266               ! save the momentum trends for diagnostics 
    267267               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) 
    … …  
    293293               ! add the new one to the general momentum trend 
    294294               ua(ji,jj,iku) = ua(ji,jj,iku) + zhpi(ji,jj,iku) 
    295 # if defined key_trddyn 
     295# if defined key_trddyn || defined key_trd_vor 
    296296               ! save the momentum trends for diagnostics 
    297297               utrd(ji,jj,iku,1) = zhpi(ji,jj,iku) 
    … …  
    307307               ! add the new one to the general momentum trend 
    308308               va(ji,jj,ikv) = va(ji,jj,ikv) + zhpj(ji,jj,ikv) 
    309 # if defined key_trddyn 
     309# if defined key_trddyn || defined key_trd_vor 
    310310               ! save the momentum trends for diagnostics 
    311311               vtrd(ji,jj,ikv,1) = zhpj(ji,jj,ikv) 
    … …  
    394394            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) 
    395395            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) 
    396 # if defined key_trddyn 
     396# if defined key_trddyn || defined key_trd_vor 
    397397            ! save the momentum trends for diagnostics 
    398398            utrd(ji,jj,1,1) = zhpi(ji,jj,1) 
    … …  
    419419               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) 
    420420               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) 
    421 # if defined key_trddyn 
     421# if defined key_trddyn || defined key_trd_vor 
    422422               ! save the momentum trends for diagnostics 
    423423               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) 

  • trunk/NEMO/OPA_SRC/DYN/dynhpg_atsk.F90

    r106 r109  
    113113            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) + zuap 
    114114            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) + zvap 
    115 #   if defined key_trddyn 
     115#   if defined key_trddyn || defined key_trd_vor 
    116116            ! save the trend for diagnostics 
    117117            utrd(ji,jj,1,1) = zhpi(ji,jj,1) + zuap 
    … …  
    139139               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) + zuap 
    140140               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) + zvap 
    141 #   if defined key_trddyn 
     141#   if defined key_trddyn || defined key_trd_vor 
    142142               ! save the momentum trends for diagnostics 
    143143               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) + zuap 
    … …  
    226226            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) 
    227227            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) 
    228 #   if defined key_trddyn 
     228#   if defined key_trddyn || defined key_trd_vor 
    229229            ! save the momentum trends for diagnostics 
    230230            utrd(ji,jj,1,1) = zhpi(ji,jj,1) 
    … …  
    249249               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) 
    250250               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) 
    251 #   if defined key_trddyn 
     251#   if defined key_trddyn || defined key_trd_vor 
    252252               ! save the momentum trends for diagnostics 
    253253               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) 
    … …  
    272272               ! add the new one to the general momentum trend 
    273273               ua(ji,jj,iku) = ua(ji,jj,iku) + zhpi(ji,jj,iku) 
    274 #   if defined key_trddyn 
     274#   if defined key_trddyn || defined key_trd_vor 
    275275               ! save the momentum trends for diagnostics 
    276276               utrd(ji,jj,iku,1) = zhpi(ji,jj,iku) 
    … …  
    286286               ! add the new one to the general momentum trend 
    287287               va(ji,jj,ikv) = va(ji,jj,ikv) + zhpj(ji,jj,ikv) 
    288 #   if  defined key_trddyn 
     288#   if  defined key_trddyn || defined key_trd_vor 
    289289               ! save the momentum trends for diagnostics 
    290290               vtrd(ji,jj,ikv,1) = zhpj(ji,jj,ikv) 
    … …  
    376376            ua(ji,jj,1) = ua(ji,jj,1) + zhpi(ji,jj,1) 
    377377            va(ji,jj,1) = va(ji,jj,1) + zhpj(ji,jj,1) 
    378 #   if defined key_trddyn 
     378#   if defined key_trddyn || defined key_trd_vor 
    379379            ! save the momentum trends for diagnostics 
    380380            utrd(ji,jj,1,1) = zhpi(ji,jj,1) 
    … …  
    399399               ua(ji,jj,jk) = ua(ji,jj,jk) + zhpi(ji,jj,jk) 
    400400               va(ji,jj,jk) = va(ji,jj,jk) + zhpj(ji,jj,jk) 
    401 #   if defined key_trddyn 
     401#   if defined key_trddyn || defined key_trd_vor 
    402402               ! save the momentum trends for diagnostics 
    403403               utrd(ji,jj,jk,1) = zhpi(ji,jj,jk) 
    … …  
    411411 
    412412      IF(l_ctl) THEN         ! print sum trends (used for debugging) 
    413          WRITE(numout,*) ' hpg  - Ua: ', SUM(ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)) 
    414          WRITE(numout,*) '      - Va: ', SUM(va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
     413         WRITE(numout,*) ' hpg  - Ua: ', SUM(ua(2:nictl,2:njctl,1:jpkm1)*umask(2:nictl,2:njctl,1:jpkm1)), & 
     414         &                      ' Va: ', SUM(va(2:nictl,2:njctl,1:jpkm1)*vmask(2:nictl,2:njctl,1:jpkm1)) 
    415415         zuap = SUM( ua(2:nictl,2:njctl,1:jpkm1) * umask(2:nictl,2:njctl,1:jpkm1) ) 
    416416         zvap = SUM( va(2:nictl,2:njctl,1:jpkm1) * vmask(2:nictl,2:njctl,1:jpkm1) ) 

  • trunk/NEMO/OPA_SRC/DYN/dynkeg.F90

    r106 r109  
    100100               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
    101101               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    102 #if defined key_trddyn 
     102#if defined key_trddyn || defined key_trd_vor 
    103103               ! add to the general momentum trends 
    104104               utrd(ji,jj,jk,2) = zua 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r106 r109  
    208208               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
    209209               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    210 #if defined key_trddyn 
     210#if defined key_trddyn || defined key_trd_vor 
    211211               ! save the horizontal biharmonic diffusive trends 
    212212               utrd(ji,jj,jk,5) = zua 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

    r106 r109  
    117117               ua(ji,jj,jk) = ua(ji,jj,jk) + wk3(ji,jj,jk) 
    118118               va(ji,jj,jk) = va(ji,jj,jk) + wk4(ji,jj,jk) 
    119 #if defined key_trddyn  
     119#if defined key_trddyn  || defined key_trd_vor 
    120120               ! save the horizontal diffusive trends 
    121121               utrd(ji,jj,jk,3) = wk3(ji,jj,jk) 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_iso.F90

    r32 r109  
    233233               ua (ji,jj,jk) = ua (ji,jj,jk) + zuah 
    234234               va (ji,jj,jk) = va (ji,jj,jk) + zvah 
    235 #if defined key_trddyn 
     235#if defined key_trddyn || defined key_trd_vor 
    236236               ! save the trends for diagnostics 
    237237               utrd(ji,jj,jk,5) = zuah 

  • trunk/NEMO/OPA_SRC/DYN/dynldf_lap.F90

    r106 r109  
    107107               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
    108108               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    109 #if defined key_trddyn 
     109#if defined key_trddyn || defined key_trd_vor 
    110110               ! save the horizontal diffusive trends 
    111111               utrd(ji,jj,jk,5) = zua 

  • trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90

    r106 r109  
    157157            spgu(ji,jj) = zspgu + zegu 
    158158            spgv(ji,jj) = zspgv + zegv 
    159 #if defined key_trddyn 
     159#if defined key_trddyn || defined key_trd_vor 
    160160            ! save the surface pressure gradient trends  
    161161            utrd2(ji,jj,1) = zspgu 
    … …  
    175175               ua(ji,jj,jk) = ua(ji,jj,jk) + spgu(ji,jj) 
    176176               va(ji,jj,jk) = va(ji,jj,jk) + spgv(ji,jj) 
    177 #if defined key_trddyn 
     177#if defined key_trddyn || defined key_trd_vor 
    178178               ! save the surface pressure gradient trend for diagnostics 
    179179               utrd(ji,jj,jk,8) = spgu(ji,jj) 
    … …  
    315315            spgv(ji,jj) = z2dt * ztdgv 
    316316#endif 
    317 #if defined key_trddyn 
     317#if defined key_trddyn || defined key_trd_vor 
    318318            ! save the transport divergence gradient trends 
    319319            utrd2(ji,jj,2) = utrd2(ji,jj,2) + ztdgu 
    … …  
    333333               ua(ji,jj,jk) = (ua(ji,jj,jk) + spgu(ji,jj)) * umask(ji,jj,jk) 
    334334               va(ji,jj,jk) = (va(ji,jj,jk) + spgv(ji,jj)) * vmask(ji,jj,jk) 
    335 #if defined key_trddyn 
     335#if defined key_trddyn || defined key_trd_vor 
    336336               ! save the surface pressure gradient trend for diagnostics 
    337337               utrd(ji,jj,jk,8) = utrd(ji,jj,jk,8) + spgu(ji,jj)/z2dt 

  • trunk/NEMO/OPA_SRC/DYN/dynspg_fsc_atsk.F90

    r106 r109  
    160160            spgu(ji,jj) = zspgu + zegu 
    161161            spgv(ji,jj) = zspgv + zegv 
    162 #if defined key_trddyn 
     162#if defined key_trddyn || defined key_trd_vor 
    163163            ! save the surface pressure gradient trends  
    164164            utrd2(ji,jj,1) = zspgu 
    … …  
    176176               ua(ji,jj,jk) = ua(ji,jj,jk) + spgu(ji,jj) 
    177177               va(ji,jj,jk) = va(ji,jj,jk) + spgv(ji,jj) 
    178 #if defined key_trddyn 
     178#if defined key_trddyn || defined key_trd_vor 
    179179               ! save the surface pressure gradient trend for diagnostics 
    180180               utrd(ji,jj,jk,8) = spgu(ji,jj) 
    … …  
    338338            spgv(ji,jj) = z2dt * ztdgv 
    339339#endif 
    340 #if defined key_trddyn 
     340#if defined key_trddyn || defined key_trd_vor 
    341341            ! save the transport divergence gradient trends 
    342342            utrd2(ji,jj,2) = utrd2(ji,jj,2) + ztdgu 
    … …  
    354354               ua(ji,jj,jk) = (ua(ji,jj,jk) + spgu(ji,jj)) * umask(ji,jj,jk) 
    355355               va(ji,jj,jk) = (va(ji,jj,jk) + spgv(ji,jj)) * vmask(ji,jj,jk) 
    356 #if defined key_trddyn 
     356#if defined key_trddyn || defined key_trd_vor 
    357357               ! save the surface pressure gradient trend for diagnostics 
    358358               utrd(ji,jj,jk,8) = utrd(ji,jj,jk,8) + spgu(ji,jj)/z2dt 

  • trunk/NEMO/OPA_SRC/DYN/dynspg_rl.F90

    r78 r109  
    520520               ua(ji,jj,jk) = ua(ji,jj,jk) - spgu(ji,jj) 
    521521               va(ji,jj,jk) = va(ji,jj,jk) - spgv(ji,jj) 
    522 # if defined key_trddyn 
     522# if defined key_trddyn || defined key_trd_vor 
    523523               ! save the surface pressure gradient trend for diagnostics 
    524524               utrd(ji,jj,jk,8) = -spgu(ji,jj) 

  • trunk/NEMO/OPA_SRC/DYN/dynzad.F90

    r106 r109  
    115115               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
    116116               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    117 #   if defined key_trddyn 
     117#   if defined key_trddyn || defined key_trd_vor 
    118118               ! save the trends for diagnostics 
    119119               utrd(ji,jj,jk,6) = zua 
    … …  
    218218               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
    219219               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    220 #   if defined key_trddyn 
     220#   if defined key_trddyn || defined key_trd_vor 
    221221               ! save the trends for diagnostics 
    222222               utrd(ji,jj,jk,6) = zua 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_exp.F90

    r3 r109  
    6464      REAL(wp), DIMENSION(jpi,jpk) ::   & 
    6565         zwx, zwy, zwz, zww               ! temporary workspace arrays 
    66 #if defined key_trddyn 
     66#if defined key_trddyn || defined key_trd_vor 
    6767      INTEGER ::   & 
    6868         ikbu, ikbum1 , ikbv, ikbvm1      ! temporary integers 
    … …  
    9898               zwx(ji,jk) = ub(ji,jj,jk) 
    9999               zwz(ji,jk) = vb(ji,jj,jk) 
    100 #if defined key_trddyn 
     100#if defined key_trddyn || defined key_trd_vor 
    101101               utrd(ji,jj,jk,7) = ua(ji,jj,jk) 
    102102               vtrd(ji,jj,jk,7) = va(ji,jj,jk) 
    … …  
    131131         END DO   
    132132 
    133 #if defined key_trddyn 
     133#if defined key_trddyn || defined key_trd_vor 
    134134         ! diagnose the vertical diffusive momentum trends 
    135135 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

    r106 r109  
    7575      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    7676         zwi                             ! temporary workspace arrays 
    77 #if defined key_trddyn 
     77#if defined key_trddyn || defined key_trd_vor 
    7878      INTEGER ::   & 
    79          ikbu, ikbum1, ikbv, ikbvm1   &  ! temporary integers 
     79         ikbu, ikbum1, ikbv, ikbvm1     ! temporary integers 
    8080#endif 
    8181      !!---------------------------------------------------------------------- 
    … …  
    9494      IF( neuler == 0 .AND. kt == nit000 )   z2dt = rdt 
    9595      ! Normalization to obtain the general momentum trend ua 
    96 #if defined key_trddyn 
     96#if defined key_trddyn || defined key_trd_vor 
    9797      ! Save the previously computed trend 
    9898      DO jk = 1, jpkm1 
    … …  
    190190      END DO 
    191191 
    192 #if defined key_trddyn 
     192#if defined key_trddyn || defined key_trd_vor 
    193193      ! diagnose surface and bottom momentum fluxes 
    194194      DO jj = 2, jpjm1    
    … …  
    203203            ! subtract surface forcing and bottom friction trend from vertical 
    204204            ! diffusive momentum trend 
    205             utrd(ji,jj,1     ,7) = utrd(ji,jj,1     ,7) - tautrd(ji,jj,1) 
    206             utrd(ji,jj,ikbum1,7) = utrd(ji,jj,ikbum1,7) - tautrd(ji,jj,3) 
     205            utrd(ji,jj,1     ,7) = utrd(ji,jj,1     ,7) + tautrd(ji,jj,1) 
     206            utrd(ji,jj,ikbum1,7) = utrd(ji,jj,ikbum1,7) + tautrd(ji,jj,3) 
    207207         END DO 
    208208      END DO 
    … …  
    214214            DO ji = fs_2, fs_jpim1   ! vector opt. 
    215215               zua = ( ua(ji,jj,jk) - ub(ji,jj,jk) ) / z2dt 
    216 #if defined key_trddyn 
     216#if defined key_trddyn || defined key_trd_vor 
    217217               ! save the vertical diffusive momentum trend (general trend - previous one) 
    218218               utrd(ji,jj,jk,7) = zua - utrd(ji,jj,jk,7) 
    … …  
    309309      END DO 
    310310 
    311 #if defined key_trddyn 
     311#if defined key_trddyn || defined key_trd_vor 
    312312      ! diagnose surface and bottom momentum fluxes 
    313313      DO jj = 2, jpjm1    
    … …  
    322322            ! subtract surface forcing and bottom friction trend from vertical 
    323323            ! diffusive momentum trend 
    324             vtrd(ji,jj,1     ,7) = vtrd(ji,jj,1     ,7) - tautrd(ji,jj,2) 
    325             vtrd(ji,jj,ikbvm1,7) = vtrd(ji,jj,ikbvm1,7) - tautrd(ji,jj,4) 
     324            vtrd(ji,jj,1     ,7) = vtrd(ji,jj,1     ,7) + tautrd(ji,jj,2) 
     325            vtrd(ji,jj,ikbvm1,7) = vtrd(ji,jj,ikbvm1,7) + tautrd(ji,jj,4) 
    326326         END DO 
    327327      END DO 
    … …  
    333333            DO ji = fs_2, fs_jpim1   ! vector opt. 
    334334               zva = ( va(ji,jj,jk) - vb(ji,jj,jk) ) / z2dt 
    335 #if defined key_trddyn 
     335#if defined key_trddyn || defined key_trd_vor 
    336336               ! save the vertical diffusive momentum fluxes 
    337337               vtrd(ji,jj,jk,7) = zva - vtrd(ji,jj,jk,7) 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_imp_atsk.F90

    r106 r109  
    6969         zwx, zwy, zwz,               &  ! workspace 
    7070         zwd, zws, zwi, zwt 
    71 #if defined key_trddyn 
     71#if defined key_trddyn || defined key_trd_vor 
    7272      INTEGER ::   & 
    73          ikbu, ikbum1, ikbv, ikbvm1   &  ! temporary integers 
     73         ikbu, ikbum1, ikbv, ikbvm1     ! temporary integers 
    7474#endif 
    7575      !!---------------------------------------------------------------------- 
    … …  
    176176            DO ji = 2, jpim1 
    177177               zua = ( zwx(ji,jk) - ub(ji,jj,jk) ) / z2dt 
    178 #if defined key_trddyn 
     178#if defined key_trddyn || defined key_trd_vor 
    179179               ! save the vertical diffusive momentum trend 
    180180               utrd(ji,jj,jk,7) = zua - ua(ji,jj,jk) 
    … …  
    184184         END DO 
    185185 
    186 #if defined key_trddyn 
     186#if defined key_trddyn || defined key_trd_vor 
    187187         ! diagnose surface and bottom momentum fluxes 
    188188         DO ji = 2, jpim1 
    … …  
    285285            DO ji = 2, jpim1 
    286286               zva = ( zwx(ji,jk) - vb(ji,jj,jk) ) / z2dt 
    287 #if defined key_trddyn 
     287#if defined key_trddyn || defined key_trd_vor 
    288288               ! save the vertical diffusive momentum fluxes 
    289289               vtrd(ji,jj,jk,7) = zva - va(ji,jj,jk) 
    … …  
    293293         END DO 
    294294 
    295 #if defined key_trddyn 
     295#if defined key_trddyn || defined key_trd_vor 
    296296         ! diagnose surface and bottom momentum fluxes 
    297297         DO ji = 2, jpim1 

  • trunk/NEMO/OPA_SRC/DYN/dynzdf_iso.F90

    r32 r109  
    114114         zfuw, zdiu, zdju, zdj1u,      &  ! workspace 
    115115         zfvw, zdiv, zdjv, zdj1v 
    116 #if defined key_trddyn 
     116#if defined key_trddyn || defined key_trd_vor 
    117117      INTEGER ::   & 
    118118         ikbu, ikbum1 , ikbv, ikbvm1      ! temporary integers 
    … …  
    247247               ua(ji,jj,jk) = ua(ji,jj,jk) + zuav 
    248248               va(ji,jj,jk) = va(ji,jj,jk) + zvav 
    249 #if defined key_trddyn 
     249#if defined key_trddyn || defined key_trd_vor 
    250250               ! save the trends for diagnostics 
    251251               utrd(ji,jj,jk,5) = utrd(ji,jj,jk,5) + zuav 
    … …  
    290290            DO ji = 2, jpim1 
    291291               zua = ( zwx(ji,jk) - ub(ji,jj,jk) ) / z2dt 
    292 #if defined key_trddyn 
     292#if defined key_trddyn || defined key_trd_vor 
    293293               ! save the vertical diffusive momentum trend 
    294294               utrd(ji,jj,jk,7) = zua - ua(ji,jj,jk) 
    … …  
    298298         END DO   
    299299 
    300 #if defined key_trddyn 
     300#if defined key_trddyn || defined key_trd_vor 
    301301 
    302302         ! 1.4 diagnose surface and bottom momentum fluxes 
    … …  
    351351            DO ji = 2, jpim1 
    352352               zva = ( zwx(ji,jk) - vb(ji,jj,jk) ) / z2dt 
    353 #if defined key_trddyn 
     353#if defined key_trddyn || defined key_trd_vor 
    354354               ! save the vertical diffusive momentum fluxes 
    355355               vtrd(ji,jj,jk,7) = zva - va(ji,jj,jk) 
    … …  
    359359         END DO    
    360360 
    361 #if defined key_trddyn 
     361#if defined key_trddyn || defined key_trd_vor 
    362362         ! 2.4 diagnose surface and bottom momentum fluxes 
    363363         DO ji = 2, jpim1 
Note: See TracChangeset for help on using the changeset viewer.