Changeset 873


Ignore:
Timestamp:
05/17/19 16:31:11 (5 years ago)
Author:
ymipsl
Message:

Fix for tracer conservation when adding smooth physics tendency.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • codes/icosagcm/trunk/src/physics/physics_external.F90

    r872 r873  
    133133    
    134134          DO iq=1, nqtot 
    135             dq(:,:,iq)=((q(:,:,iq)-q0(:,:,iq))/itau_physics)/rhodz(:,:) 
     135            dq(:,:,iq)=((q(:,:,iq)-q0(:,:,iq))/itau_physics)*rhodz(:,:) 
    136136          ENDDO 
    137137          
     
    160160       CALL compute_rhodz(.TRUE., ps, rhodz) 
    161161       DO iq=1, nqtot 
    162          q(:,:,iq)=q(:,:,iq)+dq(:,:,iq)*rhodz(:,:) 
     162         q(:,:,iq)=q(:,:,iq)+dq(:,:,iq)/rhodz(:,:) 
    163163       ENDDO 
    164164     ENDDO 
Note: See TracChangeset for help on using the changeset viewer.