- Timestamp:
- 2020-04-03T18:54:55+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DIA/diamlr.F90
r12482 r12680 4 4 !! Management of the IOM context for multiple-linear-regression analysis 5 5 !!====================================================================== 6 !! History : ! 2019 (S. Mueller)6 !! History : 4.0 ! 2019 (S. Mueller) Original code 7 7 !!---------------------------------------------------------------------- 8 8 9 9 USE par_oce , ONLY : wp, jpi, jpj 10 10 USE phycst , ONLY : rpi 11 USE dom_oce , ONLY : adatrj 12 USE tide_mod 13 ! 11 14 USE in_out_manager , ONLY : lwp, numout, ln_timing 12 15 USE iom , ONLY : iom_put, iom_use, iom_update_file_name 13 USE dom_oce , ONLY : adatrj14 16 USE timing , ONLY : timing_start, timing_stop 15 17 #if defined key_iomput 16 18 USE xios 17 19 #endif 18 USE tide_mod19 20 20 21 IMPLICIT NONE 21 22 PRIVATE 22 23 23 LOGICAL, PUBLIC :: lk_diamlr = .FALSE. 24 LOGICAL, PUBLIC :: lk_diamlr = .FALSE. !: ===>>> NOT a DOCTOR norm name : use l_diamlr 25 ! lk_ is used only for logical controlled by a CPP key 24 26 25 27 PUBLIC :: dia_mlr_init, dia_mlr_iom_init, dia_mlr … … 42 44 !! 43 45 !!---------------------------------------------------------------------- 44 46 ! 45 47 lk_diamlr = .TRUE. 46 48 ! 47 49 IF(lwp) THEN 48 50 WRITE(numout, *) … … 50 52 WRITE(numout, *) '~~~~~~~~~~~~ multiple-linear-regression analysis' 51 53 END IF 52 54 ! 53 55 END SUBROUTINE dia_mlr_init 56 54 57 55 58 SUBROUTINE dia_mlr_iom_init … … 396 399 END SUBROUTINE dia_mlr_iom_init 397 400 401 398 402 SUBROUTINE dia_mlr 399 403 !!---------------------------------------------------------------------- … … 403 407 !! 404 408 !!---------------------------------------------------------------------- 405 406 409 REAL(wp), DIMENSION(jpi,jpj) :: zadatrj2d 410 !!---------------------------------------------------------------------- 407 411 408 412 IF( ln_timing ) CALL timing_start('dia_mlr') … … 411 415 ! (value of adatrj converted to time in units of seconds) 412 416 ! 413 ! A 2-dimensional field of constant value is sent, and subsequently used 414 ! directly or transformed to a scalar or a constant 3-dimensional field as 415 ! required. 417 ! A 2-dimensional field of constant value is sent, and subsequently used directly 418 ! or transformed to a scalar or a constant 3-dimensional field as required. 416 419 zadatrj2d(:,:) = adatrj*86400.0_wp 417 420 IF ( iom_use('diamlr_time') ) CALL iom_put('diamlr_time', zadatrj2d) 418 421 ! 419 422 IF( ln_timing ) CALL timing_stop('dia_mlr') 420 423 ! 421 424 END SUBROUTINE dia_mlr 422 425 426 !!====================================================================== 423 427 END MODULE diamlr
Note: See TracChangeset
for help on using the changeset viewer.