- Timestamp:
- 2017-12-22T13:26:37+01:00 (6 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_update.F90
r8521 r9163 75 75 USE lib_mpp, ONLY: ctl_stop 76 76 USE par_kind, ONLY: wp 77 USE par_medusa, ONLY: jp_medusa, 77 USE par_medusa, ONLY: jp_medusa, jp_msa0, jp_msa1, & 78 78 jpalk, jpchd, jpchn, jpdet, jpdic, & 79 79 jpdin, jpdtc, jpfer, jpoxy, jppds, & … … 83 83 jpoxy_lc, jppds_lc, jpphd_lc, jpphn_lc, & 84 84 jpsil_lc, jpzme_lc, jpzmi_lc 85 USE par_oce, ONLY: jpi, jpim1, jpj, jpjm1 85 USE par_oce, ONLY: jpi, jpim1, jpj, jpjm1, jpk 86 86 USE par_trc, ONLY: jptra 87 87 USE sms_medusa, ONLY: friver_dep, & -
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsms_medusa.F90
r8521 r9163 23 23 USE trcsed_medusa 24 24 USE trcavg_medusa 25 !! for SMS trends 26 USE par_medusa, ONLY: jp_msa0, jp_msa1, jp_medusa 27 USE par_oce, ONLY: jpi, jpj, jpk 28 USE trd_oce, ONLY: jptra_sms, l_trdtrc 29 USE trdtrc 25 30 26 31 … … 47 52 !!---------------------------------------------------------------------- 48 53 INTEGER, INTENT(in) :: kt ! ocean time-step index 54 !! Loop variables 55 INTEGER :: jn 56 !! trend temporary array: 57 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrmed 58 49 59 50 60 # if defined key_debug_medusa … … 58 68 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 59 69 ENDIF 70 71 !! MEDUSA SMS trends: 72 IF( l_trdtrc ) THEN 73 CALL wrk_alloc( jpi, jpj, jpk, jp_medusa, ztrmed ) 74 ztrmed(:,:,:,:)=0.0 75 DO jn = 1, jp_medusa 76 ztrmed(:,:,:,jn) = tra(:,:,:,jp_msa0 + jn - 1) 77 END DO 78 END IF 60 79 61 80 CALL trc_avg_medusa( kt ) ! rolling average module … … 97 116 # endif 98 117 118 !! MEDUSA SMS trends: 119 IF( l_trdtrc ) THEN 120 DO jn = 1, jp_medusa 121 ztrmed(:,:,:,jn) = tra(:,:,:,jp_msa0 + jn - 1)-ztrmed(:,:,:,jn) 122 CALL trd_trc( ztrmed(:,:,:,jn), jn, jptra_sms, kt ) ! save trends 123 END DO 124 CALL wrk_dealloc( jpi, jpj, jpk, jp_medusa, ztrmed ) 125 END IF 126 127 99 128 END SUBROUTINE trc_sms_medusa 100 129
Note: See TracChangeset
for help on using the changeset viewer.