Ignore:
Timestamp:
2017-12-22T13:26:37+01:00 (3 years ago)
Author:
frrh
Message:

Add code from Julien Palmieri's Met Office GMED ticket 338.
This incorporates code from branches/NERC/dev_r5518_GO6_package_trdtrc
revisions 8454:9020 inclusive.

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  
    7575      USE lib_mpp,           ONLY: ctl_stop 
    7676      USE par_kind,          ONLY: wp 
    77       USE par_medusa,        ONLY: jp_medusa,                                & 
     77      USE par_medusa,        ONLY: jp_medusa, jp_msa0, jp_msa1,              & 
    7878                                   jpalk, jpchd, jpchn, jpdet, jpdic,        & 
    7979                                   jpdin, jpdtc, jpfer, jpoxy, jppds,        & 
     
    8383                                   jpoxy_lc, jppds_lc, jpphd_lc, jpphn_lc,   & 
    8484                                   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 
    8686      USE par_trc,           ONLY: jptra 
    8787      USE sms_medusa,        ONLY: friver_dep,                               & 
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcsms_medusa.F90

    r8521 r9163  
    2323   USE trcsed_medusa 
    2424   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 
    2530 
    2631 
     
    4752      !!---------------------------------------------------------------------- 
    4853      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 
    4959 
    5060# if defined key_debug_medusa 
     
    5868       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
    5969      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 
    6079 
    6180      CALL trc_avg_medusa( kt ) ! rolling average module 
     
    97116# endif 
    98117 
     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 
    99128   END SUBROUTINE trc_sms_medusa 
    100129    
Note: See TracChangeset for help on using the changeset viewer.