Changeset 10251 for branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90
- Timestamp:
- 2018-10-29T15:20:26+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_AMM15_package/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90
r10249 r10251 19 19 20 20 PUBLIC trc_wri_my_trc 21 #if defined key_tracer_budget22 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:), SAVE :: trb_temp ! slwa23 #endif24 25 21 26 22 # include "top_substitute.h90" 27 23 CONTAINS 28 24 29 #if defined key_tracer_budget30 SUBROUTINE trc_wri_my_trc (kt, fl) ! slwa31 #else32 25 SUBROUTINE trc_wri_my_trc 33 #endif34 26 !!--------------------------------------------------------------------- 35 27 !! *** ROUTINE trc_wri_trc *** … … 37 29 !! ** Purpose : output passive tracers fields 38 30 !!--------------------------------------------------------------------- 39 #if defined key_tracer_budget40 INTEGER, INTENT( in ), OPTIONAL :: fl41 INTEGER, INTENT( in ) :: kt42 REAL(wp), DIMENSION(jpi,jpj,jpk) :: trpool !tracer pool temporary output43 #else44 INTEGER, INTENT( in ) :: kt45 #endif46 31 CHARACTER (len=20) :: cltra 47 INTEGER :: jn ,jk ! JC TODO jk defined here but may not be used32 INTEGER :: jn 48 33 !!--------------------------------------------------------------------- 49 34 50 35 ! write the tracer concentrations in the file 51 36 ! --------------------------------------- 52 53 54 #if defined key_tracer_budget55 IF( PRESENT(fl)) THEN56 ! depth integrated57 ! for strict budgetting write this out at end of timestep as an average between 'now' and 'after' at kt58 DO jn = jp_myt0, jp_myt159 IF(ln_trdtrc (jn))THEN60 trpool(:,:,:) = 0.5 * ( trn(:,:,:,jn) * fse3t_a(:,:,:) + &61 trb_temp(:,:,:,jn) * fse3t(:,:,:) )62 cltra = TRIM( ctrcnm(jn) )//"e3t" ! depth integrated output63 IF( kt == nittrc000 ) write(6,*)'output pool ',cltra64 DO jk = 1, jpk65 trpool(:,:,jk) = trpool(:,:,jk)66 END DO67 CALL iom_put( cltra, trpool)68 69 END IF70 END DO71 72 ELSE73 74 IF( kt == nittrc000 ) THEN75 ALLOCATE(trb_temp(jpi,jpj,jpk,jp_my_trc)) ! slwa76 ENDIF77 trb_temp(:,:,:,:)=trn(:,:,:,:) ! slwa save for tracer budget (unfiltered trn)78 79 80 END IF81 #else82 37 DO jn = jp_myt0, jp_myt1 83 38 cltra = TRIM( ctrcnm(jn) ) ! short title for tracer 84 39 CALL iom_put( cltra, trn(:,:,:,jn) ) 85 40 END DO 86 #endif87 41 ! 88 42 END SUBROUTINE trc_wri_my_trc … … 94 48 PUBLIC trc_wri_my_trc 95 49 CONTAINS 96 #if defined key_tracer_budget 97 SUBROUTINE trc_wri_my_trc (kt, fl) ! slwa 98 INTEGER, INTENT( in ), OPTIONAL :: fl 99 INTEGER, INTENT( in ) :: kt 100 #else 101 ! JC TODO Subroutine arguments (kt) inconsistent with earlier definition 102 SUBROUTINE trc_wri_my_trc (kt) 103 INTEGER, INTENT( in ) :: kt 104 #endif 50 SUBROUTINE trc_wri_my_trc ! Empty routine 105 51 END SUBROUTINE trc_wri_my_trc 106 52 #endif
Note: See TracChangeset
for help on using the changeset viewer.