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 7990 for branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90 – NEMO

Ignore:
Timestamp:
2017-04-29T17:24:54+02:00 (7 years ago)
Author:
gm
Message:

#1880 (HPC-09): OPA remove avmu, avmv from zdf modules + move CALL tke(gls)_rst & gls_rst in zdftke(gls) + rename zdftmx and zdfqiao

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90

    r7953 r7990  
    99   USE dom_oce         ! ocean space and time domain 
    1010   USE zdf_oce         ! ocean vertical physics     
    11    USE zdfgls, ONLY: mxln 
     11   USE zdfgls   , ONLY : hmxn 
    1212   USE in_out_manager  ! I/O units 
    1313   USE iom             ! I/0 library 
     
    107107      IF( ln_zdfgls ) THEN 
    108108         en_25h(:,:,:) = en(:,:,:) 
    109          rmxln_25h(:,:,:) = mxln(:,:,:) 
     109         rmxln_25h(:,:,:) = hmxn(:,:,:) 
    110110      ENDIF 
    111111#if defined key_lim3 || defined key_lim2 
     
    169169         IF( ln_zdfgls ) THEN 
    170170            en_25h(:,:,:)        = en_25h(:,:,:) + en(:,:,:) 
    171             rmxln_25h(:,:,:)      = rmxln_25h(:,:,:) + mxln(:,:,:) 
     171            rmxln_25h(:,:,:)      = rmxln_25h(:,:,:) + hmxn(:,:,:) 
    172172         ENDIF 
    173173         cnt_25h = cnt_25h + 1 
     
    217217         zw3d(:,:,:) = vn_25h(:,:,:)*vmask(:,:,:) + zmdi*(1.0-vmask(:,:,:)) 
    218218         CALL iom_put("vomecrty25h", zw3d  )   ! j-current 
    219          zw3d(:,:,:) = wn_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     219         zw3d(:,:,:) = wn_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    220220         CALL iom_put("vomecrtz25h", zw3d )   ! k-current 
    221221         ! Write vertical physics 
    222          zw3d(:,:,:) = avt_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     222         zw3d(:,:,:) = avt_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    223223         CALL iom_put("avt25h", zw3d )   ! diffusivity 
    224          zw3d(:,:,:) = avm_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     224         zw3d(:,:,:) = avm_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    225225         CALL iom_put("avm25h", zw3d)   ! viscosity 
    226226         IF( ln_zdftke ) THEN 
    227             zw3d(:,:,:) = en_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     227            zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    228228            CALL iom_put("tke25h", zw3d)   ! tke 
    229229         ENDIF 
    230230         IF( ln_zdfgls ) THEN 
    231             zw3d(:,:,:) = en_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     231            zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    232232            CALL iom_put("tke25h", zw3d)   ! tke 
    233             zw3d(:,:,:) = rmxln_25h(:,:,:)*tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
     233            zw3d(:,:,:) = rmxln_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    234234            CALL iom_put( "mxln25h",zw3d) 
    235235         ENDIF 
     
    249249         IF( ln_zdfgls ) THEN 
    250250            en_25h(:,:,:) = en(:,:,:) 
    251             rmxln_25h(:,:,:) = mxln(:,:,:) 
     251            rmxln_25h(:,:,:) = hmxn(:,:,:) 
    252252         ENDIF 
    253253         cnt_25h = 1 
Note: See TracChangeset for help on using the changeset viewer.