Ignore:
Timestamp:
2019-05-10T18:43:09+02:00 (19 months ago)
Author:
acc
Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : Convert TOP routines in TOP/TRP directory and all knock on effects of these conversions. SETTE tested (GYRE_PISCES only)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/TRP/trcbbl.F90

    r10954 r10966  
    6161      IF( l_trdtrc )  THEN 
    6262         ALLOCATE( ztrtrd(jpi,jpj,jpk,jptra) ) ! temporary save of trends 
    63          ztrtrd(:,:,:,:)  = tra(:,:,:,:) 
     63         ztrtrd(:,:,:,:)  = tr(:,:,:,:,Krhs) 
    6464      ENDIF 
    6565 
     
    6767      IF( nn_bbl_ldf == 1 ) THEN 
    6868         ! 
    69          CALL tra_bbl_dif( trb, tra, jptra, Kmm )   
     69         CALL tra_bbl_dif( tr(:,:,:,:,Kbb), tr(:,:,:,:,Krhs), jptra, Kmm )   
    7070         IF( ln_ctl )   THEN 
    7171            WRITE(charout, FMT="(' bbl_dif')")  ;  CALL prt_ctl_trc_info(charout) 
    72             CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
     72            CALL prt_ctl_trc( tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    7373         ENDIF 
    7474         ! 
     
    7878      IF( nn_bbl_adv /= 0 ) THEN 
    7979         ! 
    80          CALL tra_bbl_adv( trb, tra, jptra, Kmm )   
     80         CALL tra_bbl_adv( tr(:,:,:,:,Kbb), tr(:,:,:,:,Krhs), jptra, Kmm )   
    8181         IF( ln_ctl )   THEN 
    8282            WRITE(charout, FMT="(' bbl_adv')")  ;  CALL prt_ctl_trc_info(charout) 
    83             CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
     83            CALL prt_ctl_trc( tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    8484         ENDIF 
    8585         ! 
     
    8888      IF( l_trdtrc )   THEN                      ! save the horizontal diffusive trends for further diagnostics 
    8989        DO jn = 1, jptra 
    90            ztrtrd(:,:,:,jn) = tra(:,:,:,jn) - ztrtrd(:,:,:,jn) 
     90           ztrtrd(:,:,:,jn) = tr(:,:,:,jn,Krhs) - ztrtrd(:,:,:,jn) 
    9191           CALL trd_tra( kt, Kmm, Krhs, 'TRC', jn, jptra_bbl, ztrtrd(:,:,:,jn) ) 
    9292        END DO 
Note: See TracChangeset for help on using the changeset viewer.