Changeset 6349


Ignore:
Timestamp:
2016-02-24T14:18:49+01:00 (5 years ago)
Author:
davestorkey
Message:

Temporary fix for bug in NEMO whereby grid cell thicknesses were written
to diagnostics at wrong time level. NB. This has now been fixed at the
latest version of 3.6_stable (ticket #1682). Only putting the fix in this
branch in order to rerun some old experiments.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/product_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r5962 r6349  
    594594         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   - interpolate scale factors and compute depths for next time step' 
    595595      ENDIF 
     596 
     597      ! Write outputs 
     598      ! ============= 
     599      CALL iom_put(     "e3t" , fse3t_n  (:,:,:) ) 
     600      CALL iom_put(     "e3u" , fse3u_n  (:,:,:) ) 
     601      CALL iom_put(     "e3v" , fse3v_n  (:,:,:) ) 
     602      CALL iom_put(     "e3w" , fse3w_n  (:,:,:) ) 
     603      CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) ) 
     604      IF( iom_use("e3tdef") )   & 
     605         CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 ) 
     606 
    596607      ! 
    597608      ! Time filter and swap of scale factors 
     
    666677      END DO 
    667678 
    668       ! Write outputs 
    669       ! ============= 
    670       CALL iom_put(     "e3t" , fse3t_n  (:,:,:) ) 
    671       CALL iom_put(     "e3u" , fse3u_n  (:,:,:) ) 
    672       CALL iom_put(     "e3v" , fse3v_n  (:,:,:) ) 
    673       CALL iom_put(     "e3w" , fse3w_n  (:,:,:) ) 
    674       CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) ) 
    675       IF( iom_use("e3tdef") )   & 
    676          CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 ) 
    677679 
    678680      ! write restart file 
Note: See TracChangeset for help on using the changeset viewer.