- Timestamp:
- 2017-04-29T17:24:54+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/dia25h.F90
r7953 r7990 9 9 USE dom_oce ! ocean space and time domain 10 10 USE zdf_oce ! ocean vertical physics 11 USE zdfgls , ONLY: mxln11 USE zdfgls , ONLY : hmxn 12 12 USE in_out_manager ! I/O units 13 13 USE iom ! I/0 library … … 107 107 IF( ln_zdfgls ) THEN 108 108 en_25h(:,:,:) = en(:,:,:) 109 rmxln_25h(:,:,:) = mxln(:,:,:)109 rmxln_25h(:,:,:) = hmxn(:,:,:) 110 110 ENDIF 111 111 #if defined key_lim3 || defined key_lim2 … … 169 169 IF( ln_zdfgls ) THEN 170 170 en_25h(:,:,:) = en_25h(:,:,:) + en(:,:,:) 171 rmxln_25h(:,:,:) = rmxln_25h(:,:,:) + mxln(:,:,:)171 rmxln_25h(:,:,:) = rmxln_25h(:,:,:) + hmxn(:,:,:) 172 172 ENDIF 173 173 cnt_25h = cnt_25h + 1 … … 217 217 zw3d(:,:,:) = vn_25h(:,:,:)*vmask(:,:,:) + zmdi*(1.0-vmask(:,:,:)) 218 218 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(:,:,:)) 220 220 CALL iom_put("vomecrtz25h", zw3d ) ! k-current 221 221 ! Write vertical physics 222 zw3d(:,:,:) = avt_25h(:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))222 zw3d(:,:,:) = avt_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 223 223 CALL iom_put("avt25h", zw3d ) ! diffusivity 224 zw3d(:,:,:) = avm_25h(:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))224 zw3d(:,:,:) = avm_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 225 225 CALL iom_put("avm25h", zw3d) ! viscosity 226 226 IF( ln_zdftke ) THEN 227 zw3d(:,:,:) = en_25h(:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))227 zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 228 228 CALL iom_put("tke25h", zw3d) ! tke 229 229 ENDIF 230 230 IF( ln_zdfgls ) THEN 231 zw3d(:,:,:) = en_25h(:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))231 zw3d(:,:,:) = en_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 232 232 CALL iom_put("tke25h", zw3d) ! tke 233 zw3d(:,:,:) = rmxln_25h(:,:,:)* tmask(:,:,:) + zmdi*(1.0-tmask(:,:,:))233 zw3d(:,:,:) = rmxln_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 234 234 CALL iom_put( "mxln25h",zw3d) 235 235 ENDIF … … 249 249 IF( ln_zdfgls ) THEN 250 250 en_25h(:,:,:) = en(:,:,:) 251 rmxln_25h(:,:,:) = mxln(:,:,:)251 rmxln_25h(:,:,:) = hmxn(:,:,:) 252 252 ENDIF 253 253 cnt_25h = 1
Note: See TracChangeset
for help on using the changeset viewer.