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 11922 for NEMO/branches/2019/dev_r11879_ENHANCE-05_SimonM-Harmonic_Analysis/src – NEMO

Ignore:
Timestamp:
2019-11-15T21:12:45+01:00 (4 years ago)
Author:
smueller
Message:

Addition of output of the elapsed model time for multiple-linear-regression analysis (ticket #2175)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11879_ENHANCE-05_SimonM-Harmonic_Analysis/src/OCE/DIA/diamlr.F90

    r11921 r11922  
    77   !!---------------------------------------------------------------------- 
    88 
     9   USE par_oce        , ONLY :   wp, jpi, jpj 
    910   USE in_out_manager , ONLY :   lwp, numout, ln_timing 
     11   USE iom            , ONLY :   iom_put, iom_use 
     12   USE dom_oce        , ONLY :   adatrj 
    1013   USE timing         , ONLY :   timing_start, timing_stop 
    1114   USE xios 
     
    6366      ! Get handle to multiple-linear-regression analysis configuration; if no 
    6467      ! configuration is found, disable diamlr 
    65       IF ( lk_diamlr .AND. xios_is_valid_fieldgroup( "diamlr_fields" ) ) THEN 
     68      IF ( lk_diamlr .AND. xios_is_valid_fieldgroup( "diamlr_fields" ) .AND. xios_is_valid_field( "diamlr_time" ) ) THEN 
    6669         CALL xios_get_handle("diamlr_fields",  slxhdl_fldgrp) 
    6770      ELSE 
     
    8386      !!---------------------------------------------------------------------- 
    8487 
     88      REAL, DIMENSION(jpi,jpj) ::   zadatrj2d 
     89 
    8590      IF( ln_timing )   CALL timing_start('dia_mlr') 
    8691 
    87       ! Update time 
     92      ! Update time to the continuous time since the start of the model run 
     93      ! (value of adatrj converted to time in units of seconds) 
     94      ! 
     95      ! A 2-dimensional field of constant value is sent, and subsequently used 
     96      ! directly or transformed to a scalar or a constant 3-dimensional field as 
     97      ! required. 
     98      zadatrj2d(:,:) = adatrj*86400.0_wp 
     99      IF ( iom_use('diamlr_time') ) CALL iom_put('diamlr_time', zadatrj2d) 
    88100       
    89101      IF( ln_timing )   CALL timing_stop('dia_mlr') 
Note: See TracChangeset for help on using the changeset viewer.