Changeset 10640


Ignore:
Timestamp:
2019-02-06T13:33:10+01:00 (14 months ago)
Author:
jamesharle
Message:

Correction for un-initialised variable wn_25h in 25hour meaning
diagmostics (see ticket #2232).

Location:
NEMO/releases/release-4.0
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/cfgs/AMM12/EXPREF/file_def_nemo-oce.xml

    r9572 r10640  
    6767 
    6868     <file id="file15" name_suffix="_25hourm_grid_W" description="ocean W grid variables, 25h meaned" enabled=".TRUE." > 
    69      <field field_ref="vomecrtz25h"          name="vomerctz"  operation="instant" enabled=".TRUE."/> 
    70      <field field_ref="woce"                 name="vovecrtz"  operation="instant" enabled=".TRUE."/> 
     69     <field field_ref="vovecrtz25h"          name="vovecrtz"  operation="instant" enabled=".TRUE."/> 
    7170     <field field_ref="avt25h"               name="votkeavt"  operation="instant" enabled=".TRUE."/> 
    7271     <field field_ref="avm25h"               name="votkeavm"  operation="instant" enabled=".TRUE."/> 
  • NEMO/releases/release-4.0/src/OCE/DIA/dia25h.F90

    r10499 r10640  
    100100      un_25h  (:,:,:) = ub  (:,:,:) 
    101101      vn_25h  (:,:,:) = vb  (:,:,:) 
    102       wn_25h  (:,:,:) = wn  (:,:,:) 
    103102      avt_25h (:,:,:) = avt (:,:,:) 
    104103      avm_25h (:,:,:) = avm (:,:,:) 
     
    148147      ! local variable for debugging 
    149148      ll_print = ll_print .AND. lwp 
     149 
     150      ! wn_25h could not be initialised in dia_25h_init, so we do it here instead 
     151      IF( kt == nn_it000 ) THEN 
     152         wn_25h(:,:,:) = wn(:,:,:) 
     153      ENDIF 
    150154 
    151155      ! Sum of 25 hourly instantaneous values to give a 25h mean from 24hours every day 
     
    223227         CALL iom_put("vomecrty25h", zw3d  )   ! j-current 
    224228         zw3d(:,:,:) = wn_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
    225          CALL iom_put("vomecrtz25h", zw3d )   ! k-current 
     229         CALL iom_put("vovecrtz25h", zw3d )   ! k-current 
    226230         ! Write vertical physics 
    227231         zw3d(:,:,:) = avt_25h(:,:,:)*wmask(:,:,:) + zmdi*(1.0-tmask(:,:,:)) 
Note: See TracChangeset for help on using the changeset viewer.