Ignore:
Timestamp:
2017-01-11T17:15:08+01:00 (5 years ago)
Author:
davestorkey
Message:

Correct calculation of total tracer trend. New calculation uses ta-tn rather than
ta-tb. This is cleaner since ta and tn are both non-time-filtered whereas tb is
time-filtered. The new variable passes the sanity check whereby the time mean of
the total trend multiplied by the time-interval gives the change between the start
and end fields.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r7061 r7549  
    136136            CALL trd_tra( kt, 'TRA', jp_sal, jptra_zdfp, ztrds ) 
    137137         ENDIF 
    138          ! total trend (note slightly inaccurate because tsb is time-filtered and tsa isn't) 
     138         ! total trend for the non-time-filtered variables.  
    139139         ztrdt(:,:,:) = 0._wp 
    140140         ztrds(:,:,:) = 0._wp 
    141141         DO jk = 1, jpkm1 
    142             ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / r2dtra(jk)  
    143             ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / r2dtra(jk)  
     142            zfact = 1.0 / rdttra(jk) 
     143            ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsn(:,:,jk,jp_tem) ) * zfact  
     144            ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsn(:,:,jk,jp_sal) ) * zfact  
    144145         END DO 
    145146         CALL trd_tra( kt, 'TRA', jp_tem, jptra_tot, ztrdt ) 
Note: See TracChangeset for help on using the changeset viewer.